Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,468 G→T intergenic (+48/‑39) ilvL → / → ilvX ilvG operon leader peptide/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,4680GT95.6% 295.1 / 5.4 90intergenic (+48/‑39)ilvL/ilvXilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (46/40);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950371‑3950562
                                                                                                  |                                                                                              
tgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑ttt                                                                                              <  2:1938846/100‑1 (MQ=255)
 gATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttt                                                                                             >  1:1319722/1‑100 (MQ=255)
     attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttttttt                                                                                         >  1:949300/1‑100 (MQ=255)
      ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑Ttttttttt                                                                                        <  2:1286488/100‑1 (MQ=255)
      ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑Ttttttttt                                                                                        >  2:1387221/1‑100 (MQ=255)
      ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑Ttttttttt                                                                                        >  2:1434123/1‑100 (MQ=255)
       tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑Tttttttttt                                                                                       <  2:1727886/100‑2 (MQ=255)
       tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTg                                                                                       <  2:1851638/100‑1 (MQ=255)
       tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTg                                                                                       <  2:1474184/100‑1 (MQ=255)
       tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTg                                                                                       <  2:315607/100‑1 (MQ=255)
        atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGa                                                                                      >  1:765723/1‑100 (MQ=255)
          ccccccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGGGATCAAGCCTTAACGAACTAAGCCCCCCGCACCTAAAGGTCCGGGG‑TTTTTTTTTg                                                                                       <  2:2118575/97‑1 (MQ=255)
          cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑GTTTTTTTTGAcc                                                                                    <  2:1999056/100‑1 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCt                                                                                   >  1:855080/1‑100 (MQ=255)
            cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCtt                                                                                  >  1:129853/1‑100 (MQ=255)
            cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGATTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCtt                                                                                  >  2:1542013/1‑100 (MQ=255)
              ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaa                                                                                >  1:685150/1‑100 (MQ=255)
              ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaa                                                                                >  2:1490146/1‑100 (MQ=255)
               cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaaa                                                                               <  2:706524/100‑1 (MQ=255)
                gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaaaa                                                                              <  2:2253539/100‑1 (MQ=255)
                   cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACa                                                                           >  2:877873/1‑100 (MQ=255)
                    ggggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACAt                                                                          <  2:1771706/100‑1 (MQ=255)
                     gggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATa                                                                         <  1:1006877/100‑1 (MQ=255)
                     gggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATa                                                                         <  2:144007/100‑1 (MQ=255)
                       gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAAc                                                                       <  2:368062/100‑1 (MQ=255)
                        cTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAAcc                                                                      >  2:388320/1‑100 (MQ=255)
                          gCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCga                                                                    >  1:1274319/1‑100 (MQ=255)
                          gCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCga                                                                    >  2:1557663/1‑100 (MQ=255)
                            aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  >  1:1747764/1‑100 (MQ=255)
                            aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  >  1:408887/1‑100 (MQ=255)
                            aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  <  1:1434127/100‑1 (MQ=255)
                            aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  >  2:368623/1‑100 (MQ=255)
                            aCTTGGACAAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GGG‑GTTTTTTTTGACCTTAAAAACATAACCgagga                                                                 >  1:313047/1‑100 (MQ=255)
                              ttGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgaggag                                                                <  1:2245948/100‑1 (MQ=255)
                               tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                               <  1:386796/100‑1 (MQ=255)
                               tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                               <  1:1595763/100‑1 (MQ=255)
                               tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                               >  2:2185229/1‑100 (MQ=255)
                               tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                               >  2:1489244/1‑100 (MQ=255)
                               tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                               >  1:1019564/1‑100 (MQ=255)
                               tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAAGTCCGGGG‑TTTTTTTTTGACAATCATAACGTAAACgaggcgg                                                               >  1:427922/1‑97 (MQ=255)
                                tgACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                              <  1:2309769/99‑1 (MQ=255)
                                 gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCaa                                                             >  2:2309760/1‑99 (MQ=255)
                                 gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                             >  1:679122/1‑100 (MQ=255)
                                 gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                             >  2:1481311/1‑100 (MQ=255)
                                  aCGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                            <  1:1490150/100‑1 (MQ=255)
                                   cGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGAc                                                           <  1:1481315/100‑1 (MQ=255)
                                     aGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACaa                                                         >  2:2320078/1‑100 (MQ=255)
                                     aGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACaa                                                         >  2:2320281/1‑100 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        >  2:628123/1‑100 (MQ=255)
                                       gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        >  1:2027914/1‑100 (MQ=255)
                                        aaaGGCTTAGAGAGCAATCCTTAACTAACGAATACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAACCATCCCCGCGGAGAAGAAACTGa                                                      <  2:2182671/100‑1 (MQ=255)
                                         aaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                    <  2:313045/100‑1 (MQ=255)
                                          aGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                    >  1:1932221/1‑100 (MQ=255)
                                          aGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                    <  1:1719213/100‑1 (MQ=255)
                                          aGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                    <  1:1557667/100‑1 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                   <  1:1489248/100‑1 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCACCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                   >  2:1530407/1‑100 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 >  2:1855192/1‑100 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 >  1:1487253/1‑100 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  1:2320290/100‑1 (MQ=255)
                                              ttAGAGATCAAGCCTTAACGAACTAAGCCCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  1:2320087/99‑1 (MQ=255)
                                              ttAGAGATCAAGCCTTAACGAACTAAGA‑CCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                               >  1:2076785/1‑100 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                               >  1:510611/1‑100 (MQ=255)
                                                agagATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGc                                              >  2:476194/1‑100 (MQ=255)
                                                   gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                           >  1:504088/1‑100 (MQ=255)
                                                   gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                           >  1:554139/1‑100 (MQ=255)
                                                         gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                     >  1:2335791/1‑100 (MQ=255)
                                                          ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                                    >  1:1425245/1‑100 (MQ=255)
                                                           cTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   <  2:876016/100‑1 (MQ=255)
                                                            ttAACGAACTAAGCCCCCCGCGCCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                  <  1:690864/100‑1 (MQ=255)
                                                             tAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGtt                                 >  2:1008774/1‑100 (MQ=255)
                                                             tAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGtt                                 <  2:303343/100‑1 (MQ=255)
                                                              aaCGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGttt                                >  1:722644/1‑100 (MQ=255)
                                                              aaCGAACTAACACCCCCGAACCGAAAGTTCCGTGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                  <  2:2102198/98‑1 (MQ=255)
                                                               aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                               >  2:2059888/1‑100 (MQ=255)
                                                               aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                               <  2:855078/100‑1 (MQ=255)
                                                                 gAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                              <  1:1008776/99‑1 (MQ=255)
                                                                  aaCTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                            <  2:1319719/100‑1 (MQ=255)
                                                                   aCTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                           <  1:1530411/100‑1 (MQ=255)
                                                                        gACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                      >  1:801183/1‑100 (MQ=255)
                                                                        gACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                      >  1:1759047/1‑100 (MQ=255)
                                                                         aCCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                     <  2:504086/100‑1 (MQ=255)
                                                                          cccccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                      <  2:685149/98‑1 (MQ=255)
                                                                           ctccGCTCCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  1:569667/98‑1 (MQ=255)
                                                                           ccccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  1:2185236/100‑1 (MQ=255)
                                                                            cccACACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTAAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                     <  2:765722/97‑1 (MQ=255)
                                                                             ccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGTCCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGag                 >  2:569666/1‑99 (MQ=255)
                                                                              cGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACg                >  1:549226/1‑100 (MQ=255)
                                                                                  ccGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa            >  1:2141522/1‑100 (MQ=255)
                                                                                    gAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          <  1:732405/100‑1 (MQ=255)
                                                                                     aaaGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaact         <  1:877875/100‑1 (MQ=255)
                                                                                        ggTCCGGGG‑GTTTTTTTTTTTTTTTTTAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaac      >  2:29325/1‑100 (MQ=255)
                                                                                            cGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATg  >  2:1601428/1‑100 (MQ=255)
                                                                                                  |                                                                                              
TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950371‑3950562

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