Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,950,468 G→T 100% 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,4680GT94.6% 328.8 / 8.1 108intergenic (+48/‑39)ilvL/ilvXilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base T (42/60);  total (44/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 110 total aligned reads displayed.

ATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950373‑3950562
                                                                                                |                                                                                              
attattATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTtttt                                                                                            <  1:534520/100‑1 (MQ=255)
attattATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTtttt                                                                                            <  1:1286279/100‑1 (MQ=255)
  tattatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTtttttt                                                                                          >  2:2329400/1‑100 (MQ=255)
  tattatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑‑G‑GGtttttttt                                                                                        <  1:1334372/100‑1 (MQ=255)
   attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑Gtttttttt                                                                                         <  2:1325464/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGTTGCACTTGGACGAGGAAAGGCTTAGGGATCAAGCCTTAACTAACTAATACCCCCGCATTGAAATGTCCGTG‑GTGTTTTTTTg                                                                                       <  2:1154802/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTg                                                                                       <  1:2259039/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTg                                                                                       <  1:1264336/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTg                                                                                       <  1:1308205/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTg                                                                                       <  2:2289185/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTg                                                                                       <  2:849301/100‑1 (MQ=255)
      atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGa                                                                                      >  1:375311/1‑100 (MQ=255)
      atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGa                                                                                      >  1:1060948/1‑100 (MQ=255)
        cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGAcc                                                                                    <  1:1696245/100‑1 (MQ=255)
        cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGAcc                                                                                    <  1:203734/100‑1 (MQ=255)
        cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGAc                                                                                     >  1:767340/1‑99 (MQ=255)
        cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GG‑GTTTTTTTTTGACCt                                                                                   >  1:348861/1‑100 (MQ=255)
         ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCt                                                                                   <  1:1002564/100‑1 (MQ=255)
         ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCt                                                                                   >  1:384040/1‑100 (MQ=255)
         ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCt                                                                                   >  1:2377644/1‑100 (MQ=255)
          cACGGTGGGGGGTTTCAATTTGACGAGGAAGGGGTTAGACGTGAAGCCGTAAAGAACTTAGGCCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCtt                                                                                  <  2:943571/100‑1 (MQ=255)
          cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCCCACCGAAAGGTCC‑GG‑GTTTTTTTTTGACCTTa                                                                                 >  1:607113/1‑100 (MQ=255)
              gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTaaaa                                                                              <  1:1206867/100‑1 (MQ=255)
                 cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACa                                                                           >  1:1478020/1‑100 (MQ=255)
                 cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACa                                                                           <  1:910097/100‑1 (MQ=255)
                  ggggTTACACTTAGCCGAGGGAAGGGTTAGAGATCAAGCCTAAACAAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTTGACCTTAAAAACa                                                                           <  1:1762429/100‑1 (MQ=255)
                  ggggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACa                                                                           <  1:1785907/99‑1 (MQ=255)
                  ggggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACAt                                                                          >  1:2121898/1‑100 (MQ=255)
                  ggggATGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCGTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACa                                                                           <  1:494849/99‑1 (MQ=255)
                   gggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACa                                                                           <  1:2011724/98‑1 (MQ=255)
                   gggCTGCACTCGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATa                                                                         <  2:1574913/100‑1 (MQ=255)
                     gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTTAACTTaaaaaaaaaaa                                                                       >  1:212847/1‑94 (MQ=255)
                     gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAAc                                                                       >  2:692860/1‑100 (MQ=255)
                       tGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAAcc                                                                      <  1:829579/99‑1 (MQ=255)
                       tGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCg                                                                     >  1:849818/1‑100 (MQ=255)
                        gCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCga                                                                    >  1:1196522/1‑100 (MQ=255)
                         cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GG‑GGTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  <  1:2327318/100‑1 (MQ=255)
                           cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCgagga                                                                 >  1:553976/1‑100 (MQ=255)
                             tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                               <  2:2385266/100‑1 (MQ=255)
                              ggACGAGGAAAGGCTCAGAGATCGAGCCTTAACGAACTAAGACACCCGCACCGAAAGGTCC‑GG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                              <  1:1062418/99‑1 (MQ=255)
                               gACGAGGAAAGGTTTAAAGATCACGCGTTAGGGAAATAAAATGCCCGGACCGAAAGGTCAGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                              <  2:722779/99‑1 (MQ=255)
                               gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACGAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTGTGACCTTAAAAACATAACCGAGAAGCa                                                              >  2:344404/1‑99 (MQ=255)
                                   aGGGAAGGCCTAGTGATAAAGCCTGAGTGAAGCAAGCCACACGAAACGACAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGAc                                                           <  2:703334/98‑1 (MQ=255)
                                   aGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACaa                                                         >  1:1462545/1‑100 (MQ=255)
                                   aGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GG‑GGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        >  2:150367/1‑100 (MQ=255)
                                   aGGAAAGGCTGAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                              <  1:241802/95‑1 (MQ=255)
                                    ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        >  1:1593422/1‑100 (MQ=255)
                                    ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        >  1:333534/1‑100 (MQ=255)
                                     gagaGGATTAGAGACCAAGCGTGAGGGGACGAAGGCGCCCGCACAGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  1:815292/97‑1 (MQ=255)
                                     gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCTAGGAGCAGACAATg                                                       >  2:1361065/1‑100 (MQ=255)
                                     gAAAGGCTTAGAGACCAAGCCGTAGCGACCGAAGACCCCCGCACAGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  1:2099518/100‑1 (MQ=255)
                                      aaaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGa                                                      >  1:2308521/1‑100 (MQ=255)
                                         gggttAGAGATCAAGCCTTAACGAACTAAGTCCCCCGCACCGAAAGGACCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  2:2372746/93‑1 (MQ=255)
                                         ggCTTAGAGATCAATCCTTACCGAACTGAGGCCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                    <  2:714037/99‑1 (MQ=255)
                                         ggCTTAAAGATCAAGCCTGAAGGAAGGAAGACCCCCGCACCGAAAGGGCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                   <  1:466358/100‑1 (MQ=255)
                                           tttAGAGACCAGGCCTGAACGAACTAGGACCCCGGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        <  1:1366162/92‑1 (MQ=255)
                                           cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 >  1:3535/1‑100 (MQ=255)
                                            ttAGAGATCAAGCCTTAACGAACTAAGACCACCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  1:219805/93‑1 (MQ=255)
                                             tAGAGATTAAGTTTTATCGAACTAAGTCCCCCGCACGGAAAGGTACGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                               <  2:2188938/100‑1 (MQ=255)
                                             tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTAAAATTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  2:357030/98‑1 (MQ=255)
                                             tAGAGATCAAGCCGTAACGGACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  1:816537/98‑1 (MQ=255)
                                               gagaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                             >  2:1624090/1‑100 (MQ=255)
                                                 gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                           >  1:1285731/1‑100 (MQ=255)
                                                 gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                           >  2:899103/1‑100 (MQ=255)
                                                  aTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaa                                          >  1:1393021/1‑100 (MQ=255)
                                                   tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTTACTTTAAAAACTTAACCGAGGAGCATACAATGAATAACAGcaca                                           >  2:182201/1‑98 (MQ=255)
                                                   tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaa                                         >  2:1759965/1‑100 (MQ=255)
                                                   tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGGCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGGAAACAATGAATAACAGAACcaa                                         >  1:1917329/1‑100 (MQ=255)
                                                    cAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaaa                                        >  1:1715419/1‑100 (MQ=255)
                                                    cAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTCCTTTTTTGACCTCAAAACCATAACCGAGCAGCAGACCATGAATAACAGCACaaa                                         >  1:57292/1‑99 (MQ=255)
                                                       gCCTTAACGAATTTATACCCCCGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                       <  2:1057451/98‑1 (MQ=255)
                                                         cTTAACGAACTAAGACCCCCGCGCCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   <  1:1659619/100‑1 (MQ=255)
                                                          ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATGACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                  >  1:1171729/1‑100 (MQ=255)
                                                          ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                  >  1:1245761/1‑100 (MQ=255)
                                                              cGAACTTAGCCCCTCGCACAGAAAGGTCCGGG‑GTATTTTTTTGACCTTAAAAACATAACCGAGGAGCAAACAATGAATAACAGCACAAAATTCTGTTtc                               <  1:2386109/99‑1 (MQ=255)
                                                               gaacatagtcccccGCACGGGAGGGGCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAACAACAGCACAAAATTCTGTTtc                               <  1:526491/89‑1 (MQ=255)
                                                                 aCTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                           <  1:1140405/100‑1 (MQ=255)
                                                                  cTAAGTCCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                           <  1:293853/99‑1 (MQ=255)
                                                                  cTAAGACCCCCGCACCGAAGGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                          <  1:2314322/100‑1 (MQ=255)
                                                                  cTAAGACCCCCGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                          >  1:9776/1‑100 (MQ=255)
                                                                   ttggACGCCTGCAACGAAAGGGCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtctc                             <  2:820760/93‑1 (MQ=255)
                                                                        cccccGCACCGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAtt                       <  2:1715399/97‑1 (MQ=255)
                                                                         ttctgCACTGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  1:1540171/96‑1 (MQ=255)
                                                                           ccGCACCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGAc                 <  1:913899/100‑1 (MQ=255)
                                                                            cGCAGCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  2:2147825/97‑1 (MQ=255)
                                                                             gcgccgattggTCCGGG‑TTTTTTTGTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgg               <  1:420760/91‑1 (MQ=255)
                                                                             gCACCGAAAGGTCCGGG‑GTTTTTTTCTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAACTCTGTTTCTCAAGATTCAgg                   <  2:1468015/96‑1 (MQ=255)
                                                                               agcGAATGGTCCGGT‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  2:26527/98‑1 (MQ=255)
                                                                               aCCGAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  2:1591872/100‑1 (MQ=255)
                                                                                ccGAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa            >  1:76411/1‑100 (MQ=255)
                                                                                  gacaGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggg              <  2:1285422/93‑1 (MQ=255)
                                                                                  gAAAGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          >  1:1709937/1‑100 (MQ=255)
                                                                                  gAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          <  2:442845/100‑1 (MQ=255)
                                                                                  gAAAGGTCCGGG‑GTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          <  1:1598773/100‑1 (MQ=255)
                                                                                   aaaGGTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaact         <  1:1529695/100‑1 (MQ=255)
                                                                                     aGGTCCGGG‑GTTTTTTTTTTTTTTTTTAAACATAACCGAGGACCAGACAATGAATAACAGCACAAAATTCTTTTTCTCAAGATTCACGACGGGGaactaa       >  1:292940/1‑100 (MQ=255)
                                                                                      ggTCCGGG‑TTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCCGTTTCTCAAGATTCAGGACGGGGaactaac      <  2:2380321/100‑1 (MQ=255)
                                                                                      ggTCCGGG‑GTTTTTTTTTGACCTGAAAAACATAACCGAGGAGCAGACAATGAATACCAGCACAAAACTCTGTTTCTCAAGATTCAGGACGGGGaacga        <  1:1620459/98‑3 (MQ=255)
                                                                                        tCCGGG‑TTTTTTTTTTGACCTTAACAACATAACCGAGGAGCAGGCAATGAATAACTGCACATAATTCTGTTTCTCAAGATTCAGGCCGGGGaactaac      <  2:866245/98‑1 (MQ=255)
                                                                                          cGGG‑GTTTTTTTTTGACCTGAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATg  <  1:1238302/100‑1 (MQ=255)
                                                                                                |                                                                                              
ATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGG‑GGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950373‑3950562

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