Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007792 1,817 C→G S213S (TCC→TCG)  SAUSA300_RS14745 → integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077921,8170CG81.5% 325.6 / 77.3 157S213S (TCC→TCG) SAUSA300_RS14745integrase
Reads supporting (aligned to +/- strand):  ref base C (14/15);  new base G (15/113);  total (29/128)
Fisher's exact test for biased strand distribution p-value = 3.42e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-04

Only 100 of 169 total aligned reads displayed.

AAAGTATCAATACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCTA‑CGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAACCAAATT  >  NC_007792/1720‑1916
                                                                                                  |                                                                                                    
aaaGTATCAATACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCTA‑CGGATTTTcccc                                                                                                    <  1:273699/100‑1 (MQ=40)
      tCAATACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCTA‑CGGATTTTCCCCagcact                                                                                              <  1:3201222/100‑6 (MQ=255)
         aTACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCTA‑CGGATTTTCCCCATGCCACgg                                                                                           <  1:3798824/100‑2 (MQ=255)
            cGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCTA‑CGGATTTTCCCCATGCCACGAAAt                                                                                        >  1:1021402/1‑100 (MQ=40)
                             aaaGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCTA‑CGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCaa                                                                       >  1:550683/1‑100 (MQ=35)
                                                       atataAAAAGAACCGTAGGTCTCTTCAGATCTA‑CGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGatatat                                             >  1:1093337/1‑100 (MQ=255)
                                                       atataAAAAGAACCGTAGGTCTCTTCAGATCTA‑CGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGatatat                                             >  1:5136614/1‑100 (MQ=255)
                                                              caGAACCGTAGGTCTCTTCAGATCTA‑CGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGAt                                      >  1:5994565/2‑100 (MQ=255)
                                                                    agcaGGTCTCTTCAGATCTA‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCta                                >  1:1261681/4‑100 (MQ=16)
                                                                      gtgtgctctTCCGATCTACCGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCtat                               <  1:1337568/95‑1 (MQ=11)
                                                                      caGGTCTCTTCAGATCTA‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCtatg                              <  1:5881196/99‑2 (MQ=255)
                                                                                  gATCTA‑C‑GATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGc                <  1:705394/100‑1 (MQ=255)
                                                                                   ctttA‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGc                <  1:2097247/97‑1 (MQ=255)
                                                                                   aTCTC‑CGGCTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGc                <  1:4132768/100‑1 (MQ=255)
                                                                                   aTCTA‑TGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGc                <  1:159495/100‑1 (MQ=255)
                                                                                   aTCTA‑CTGTTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGc                <  1:2984031/100‑1 (MQ=255)
                                                                                   aTCTA‑CCGCTTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGc                <  1:4687877/100‑1 (MQ=255)
                                                                                    tctcgtggATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTCACTTTGc                <  1:2441239/94‑1 (MQ=14)
                                                                                    tctactggtttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGc                <  1:5646159/91‑1 (MQ=14)
                                                                                    tCTA‑TGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCa               <  1:4241680/100‑1 (MQ=255)
                                                                                    tCTA‑CGTTTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCa               <  1:949018/100‑1 (MQ=255)
                                                                                    tCTA‑CGGTTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCa               <  1:3884300/100‑1 (MQ=255)
                                                                                    tCTA‑CGGGTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACCTTgta               <  1:575684/100‑3 (MQ=11)
                                                                                    tCTA‑CGGGCTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCa               <  1:616623/100‑1 (MQ=255)
                                                                                    tCTA‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGcc               <  1:4124101/100‑2 (MQ=25)
                                                                                    tCTA‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCa               <  1:4296533/100‑1 (MQ=255)
                                                                                    ctaa‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCa               <  1:1153451/97‑1 (MQ=255)
                                                                                    cgatctcgATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGc                <  1:2941940/93‑1 (MQ=255)
                                                                                     tctgtgggttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCa               <  1:25501/92‑1 (MQ=11)
                                                                                     cttgcGGGTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCa               <  1:5143850/96‑1 (MQ=255)
                                                                                     cTA‑CGTATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGccc              <  1:3362382/100‑3 (MQ=255)
                                                                                     cTA‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCaa              >  1:6026985/1‑100 (MQ=38)
                                                                                      tA‑CTGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:2782673/100‑1 (MQ=255)
                                                                                      tA‑CGGTTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:3479726/100‑1 (MQ=255)
                                                                                      tA‑CGGGTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:948808/100‑1 (MQ=255)
                                                                                      tA‑CGGGTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:977040/100‑1 (MQ=255)
                                                                                      tA‑CGGCTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:7364/100‑1 (MQ=255)
                                                                                      tA‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:1244189/100‑1 (MQ=38)
                                                                                      tA‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:4987125/100‑1 (MQ=255)
                                                                                      tA‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:1706861/100‑1 (MQ=25)
                                                                                      tA‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTCGCAAc             <  1:5711994/100‑1 (MQ=255)
                                                                                      tA‑CGGACTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:3786546/100‑1 (MQ=255)
                                                                                      tA‑CGAACTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTACAAc             <  1:5329893/100‑1 (MQ=14)
                                                                                      tA‑CAGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:701393/100‑1 (MQ=255)
                                                                                       tgcGGTTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCagg             <  1:1951322/98‑3 (MQ=255)
                                                                                       tcgtgATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:4378808/96‑1 (MQ=255)
                                                                                       tatgggttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:1489798/94‑1 (MQ=255)
                                                                                       tatggATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTCGCAAc             <  1:2622667/97‑1 (MQ=255)
                                                                                       tacgatttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:4815673/94‑1 (MQ=255)
                                                                                       tacattttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:4017700/94‑1 (MQ=255)
                                                                                       gctggATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:3037428/97‑1 (MQ=255)
                                                                                       gccGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCaat             <  1:2683004/98‑2 (MQ=255)
                                                                                       cttccgatctCCCCACGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:4385804/91‑1 (MQ=17)
                                                                                       ctacgtctttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:4237406/93‑1 (MQ=255)
                                                                                       actggtttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:1846702/94‑1 (MQ=255)
                                                                                       acgggATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACGTGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc             <  1:5158784/97‑1 (MQ=255)
                                                                                       a‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTTACTTTGCAAc             >  1:5905881/1‑100 (MQ=2)
                                                                                       a‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTacaaca            >  1:1357416/1‑100 (MQ=255)
                                                                                       a‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACa            >  1:1288781/1‑100 (MQ=38)
                                                                                       a‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACa            >  1:1644511/1‑100 (MQ=38)
                                                                                       a‑CGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATCAGCGATA‑TTTTTTACTTTGCAACa            >  1:4915104/1‑100 (MQ=255)
                                                                                       a‑CGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACa            >  1:754954/1‑100 (MQ=255)
                                                                                       a‑CGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACa            >  1:5850847/1‑100 (MQ=255)
                                                                                       a‑CGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTCGCAACa            >  1:1177886/1‑100 (MQ=255)
                                                                                       a‑CGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTCGCAACa            >  1:2468250/1‑100 (MQ=255)
                                                                                        tcgttttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACa            >  1:1142591/6‑100 (MQ=255)
                                                                                        gcGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACa            >  1:4239152/2‑100 (MQ=25)
                                                                                        cggagttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACa            <  1:4569992/95‑1 (MQ=255)
                                                                                        ccGGGTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACa            <  1:3066392/99‑1 (MQ=255)
                                                                                         tgggttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGccgtag           <  1:2263151/96‑6 (MQ=14)
                                                                                         ggctttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAg           <  1:748666/96‑1 (MQ=255)
                                                                                         gcgagtttCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAg           >  1:5546906/6‑100 (MQ=255)
                                                                                         cgtaTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAg           <  1:4949294/97‑1 (MQ=255)
                                                                                         cgtaTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAg           <  1:4593706/97‑1 (MQ=255)
                                                                                         cGGGTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAg           <  1:1353011/100‑1 (MQ=255)
                                                                                         cGGGTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAg           <  1:2412172/100‑1 (MQ=255)
                                                                                         cGGCTTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAg           <  1:3266262/100‑1 (MQ=255)
                                                                                         accagcctCCCCGCGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTTACTTTGCAACa            >  1:4489695/8‑100 (MQ=11)
                                                                                         aagATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAg           <  1:905988/98‑1 (MQ=255)
                                                                                          ggcttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATGTTAGTGATTAGCTATA‑TTTTTTACTTTGCAACatt          <  1:3854062/97‑3 (MQ=14)
                                                                                          ggcttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGa          <  1:4874640/97‑1 (MQ=255)
                                                                                             ttttCGGCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcc       <  1:1870260/100‑1 (MQ=255)
                                                                                             ttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcc       <  1:745663/100‑1 (MQ=38)
                                                                                             ttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcc       <  1:870997/100‑1 (MQ=255)
                                                                                             ttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcc       <  1:1759078/100‑1 (MQ=255)
                                                                                             ttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcc       <  1:1631920/100‑1 (MQ=255)
                                                                                             ttttCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcc       <  1:4112232/100‑1 (MQ=255)
                                                                                             tcttCCCCGTGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcc       <  1:1488134/98‑1 (MQ=255)
                                                                                             gtttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcc       <  1:5608449/99‑1 (MQ=25)
                                                                                             ctttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcc       <  1:921285/99‑1 (MQ=255)
                                                                                               ttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGCACCat     <  1:1464626/100‑2 (MQ=255)
                                                                                               ttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAACCat     <  1:204500/100‑2 (MQ=255)
                                                                                               ttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAACCat     <  1:986047/100‑2 (MQ=25)
                                                                                               ggccccATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAACCag     >  1:3120857/3‑99 (MQ=255)
                                                                                               ctCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAACCag     <  1:3957533/99‑2 (MQ=255)
                                                                                                 ccccATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAACCagcc   >  1:2632025/1‑97 (MQ=255)
                                                                                                 ccccATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAACCagcc   <  1:3085551/100‑4 (MQ=255)
                                                                                                 ccccATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAACCagca   <  1:3258291/100‑4 (MQ=255)
                                                                                                  cccGTGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAACCagcca  <  1:5418518/100‑5 (MQ=255)
                                                                                                  cccATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAACCagcca  <  1:5830369/100‑5 (MQ=255)
                                                                                                  |                                                                                                    
AAAGTATCAATACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCTA‑CGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAACCAAATT  >  NC_007792/1720‑1916

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑