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.

ATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950373‑3950562
                                                                                                |                                                                                              
attattATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttt                                                                                            <  1:534520/100‑1 (MQ=255)
attattATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttt                                                                                            <  1:1286279/100‑1 (MQ=255)
  tattatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttt                                                                                          >  2:2329400/1‑100 (MQ=255)
  tattatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑‑GGGtttttttt                                                                                        <  1:1334372/100‑1 (MQ=255)
   attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttttttt                                                                                         <  2:1325464/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGTTGCACTTGGACGAGGAAAGGCTTAGGGATCAAGCCTTAACTAACTAATACCCCCGCATTGAAATGTCCGTGGTGTTTTTTTg                                                                                       <  2:1154802/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTg                                                                                       <  1:2259039/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTg                                                                                       <  1:1264336/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTg                                                                                       <  1:1308205/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTg                                                                                       <  2:2289185/100‑1 (MQ=255)
     tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTg                                                                                       <  2:849301/100‑1 (MQ=255)
      atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGa                                                                                      >  1:375311/1‑100 (MQ=255)
      atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGa                                                                                      >  1:1060948/1‑100 (MQ=255)
        cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGAcc                                                                                    <  1:1696245/100‑1 (MQ=255)
        cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAcc                                                                                    <  1:203734/100‑1 (MQ=255)
        cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAc                                                                                     >  1:767340/1‑99 (MQ=255)
        cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTGACCt                                                                                   >  1:348861/1‑100 (MQ=255)
         ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCt                                                                                   <  1:1002564/100‑1 (MQ=255)
         ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCt                                                                                   >  1:384040/1‑100 (MQ=255)
         ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCt                                                                                   >  1:2377644/1‑100 (MQ=255)
          cACGGTGGGGGGTTTCAATTTGACGAGGAAGGGGTTAGACGTGAAGCCGTAAAGAACTTAGGCCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCtt                                                                                  <  2:943571/100‑1 (MQ=255)
          cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCCCACCGAAAGGTCCGGGTTTTTTTTTGACCTTa                                                                                 >  1:607113/1‑100 (MQ=255)
              gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaaa                                                                              <  1:1206867/100‑1 (MQ=255)
                 cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACa                                                                           >  1:1478020/1‑100 (MQ=255)
                 cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACa                                                                           <  1:910097/100‑1 (MQ=255)
                  ggggTTACACTTAGCCGAGGGAAGGGTTAGAGATCAAGCCTAAACAAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTTGACCTTAAAAACa                                                                           <  1:1762429/100‑1 (MQ=255)
                  ggggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACa                                                                           <  1:1785907/99‑1 (MQ=255)
                  ggggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACAt                                                                          >  1:2121898/1‑100 (MQ=255)
                  ggggATGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCGTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACa                                                                           <  1:494849/99‑1 (MQ=255)
                   gggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACa                                                                           <  1:2011724/98‑1 (MQ=255)
                   gggCTGCACTCGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATa                                                                         <  2:1574913/100‑1 (MQ=255)
                     gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTAACTTaaaaaaaaaaa                                                                       >  1:212847/1‑94 (MQ=255)
                     gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAc                                                                       >  2:692860/1‑100 (MQ=255)
                       tGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAcc                                                                      <  1:829579/99‑1 (MQ=255)
                       tGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCg                                                                     >  1:849818/1‑100 (MQ=255)
                        gCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCga                                                                    >  1:1196522/1‑100 (MQ=255)
                         cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  <  1:2327318/100‑1 (MQ=255)
                           cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCgagga                                                                 >  1:553976/1‑100 (MQ=255)
                             tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                               <  2:2385266/100‑1 (MQ=255)
                              ggACGAGGAAAGGCTCAGAGATCGAGCCTTAACGAACTAAGACACCCGCACCGAAAGGTCCGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                              <  1:1062418/99‑1 (MQ=255)
                               gACGAGGAAAGGTTTAAAGATCACGCGTTAGGGAAATAAAATGCCCGGACCGAAAGGTCAGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                              <  2:722779/99‑1 (MQ=255)
                               gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACGAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTGTGACCTTAAAAACATAACCGAGAAGCa                                                              >  2:344404/1‑99 (MQ=255)
                                   aGGGAAGGCCTAGTGATAAAGCCTGAGTGAAGCAAGCCACACGAAACGACAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGAc                                                           <  2:703334/98‑1 (MQ=255)
                                   aGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACaa                                                         >  1:1462545/1‑100 (MQ=255)
                                   aGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        >  2:150367/1‑100 (MQ=255)
                                   aGGAAAGGCTGAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                              <  1:241802/95‑1 (MQ=255)
                                    ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        >  1:1593422/1‑100 (MQ=255)
                                    ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        >  1:333534/1‑100 (MQ=255)
                                     gagaGGATTAGAGACCAAGCGTGAGGGGACGAAGGCGCCCGCACAGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  1:815292/97‑1 (MQ=255)
                                     gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCTAGGAGCAGACAATg                                                       >  2:1361065/1‑100 (MQ=255)
                                     gAAAGGCTTAGAGACCAAGCCGTAGCGACCGAAGACCCCCGCACAGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  1:2099518/100‑1 (MQ=255)
                                      aaaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGa                                                      >  1:2308521/1‑100 (MQ=255)
                                         gggttAGAGATCAAGCCTTAACGAACTAAGTCCCCCGCACCGAAAGGACCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  2:2372746/93‑1 (MQ=255)
                                         ggCTTAGAGATCAATCCTTACCGAACTGAGGCCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                    <  2:714037/99‑1 (MQ=255)
                                         ggCTTAAAGATCAAGCCTGAAGGAAGGAAGACCCCCGCACCGAAAGGGCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                   <  1:466358/100‑1 (MQ=255)
                                           tttAGAGACCAGGCCTGAACGAACTAGGACCCCGGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        <  1:1366162/92‑1 (MQ=255)
                                           cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 >  1:3535/1‑100 (MQ=255)
                                            ttAGAGATCAAGCCTTAACGAACTAAGACCACCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  1:219805/93‑1 (MQ=255)
                                             tAGAGATTAAGTTTTATCGAACTAAGTCCCCCGCACGGAAAGGTACGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                               <  2:2188938/100‑1 (MQ=255)
                                             tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTAAAATTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  2:357030/98‑1 (MQ=255)
                                             tAGAGATCAAGCCGTAACGGACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  1:816537/98‑1 (MQ=255)
                                               gagaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                             >  2:1624090/1‑100 (MQ=255)
                                                 gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                           >  1:1285731/1‑100 (MQ=255)
                                                 gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                           >  2:899103/1‑100 (MQ=255)
                                                  aTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaa                                          >  1:1393021/1‑100 (MQ=255)
                                                   tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTACTTTAAAAACTTAACCGAGGAGCATACAATGAATAACAGcaca                                           >  2:182201/1‑98 (MQ=255)
                                                   tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaa                                         >  2:1759965/1‑100 (MQ=255)
                                                   tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGGCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGGAAACAATGAATAACAGAACcaa                                         >  1:1917329/1‑100 (MQ=255)
                                                    cAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaaa                                        >  1:1715419/1‑100 (MQ=255)
                                                    cAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTCCTTTTTTGACCTCAAAACCATAACCGAGCAGCAGACCATGAATAACAGCACaaa                                         >  1:57292/1‑99 (MQ=255)
                                                       gCCTTAACGAATTTATACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                       <  2:1057451/98‑1 (MQ=255)
                                                         cTTAACGAACTAAGACCCCCGCGCCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   <  1:1659619/100‑1 (MQ=255)
                                                          ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATGACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                  >  1:1171729/1‑100 (MQ=255)
                                                          ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                  >  1:1245761/1‑100 (MQ=255)
                                                              cGAACTTAGCCCCTCGCACAGAAAGGTCCGGGGTATTTTTTTGACCTTAAAAACATAACCGAGGAGCAAACAATGAATAACAGCACAAAATTCTGTTtc                               <  1:2386109/99‑1 (MQ=255)
                                                               gaacatagtcccccGCACGGGAGGGGCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAACAACAGCACAAAATTCTGTTtc                               <  1:526491/89‑1 (MQ=255)
                                                                 aCTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                           <  1:1140405/100‑1 (MQ=255)
                                                                  cTAAGTCCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                           <  1:293853/99‑1 (MQ=255)
                                                                  cTAAGACCCCCGCACCGAAGGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                          <  1:2314322/100‑1 (MQ=255)
                                                                  cTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                          >  1:9776/1‑100 (MQ=255)
                                                                   ttggACGCCTGCAACGAAAGGGCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtctc                             <  2:820760/93‑1 (MQ=255)
                                                                        cccccGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAtt                       <  2:1715399/97‑1 (MQ=255)
                                                                         ttctgCACTGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  1:1540171/96‑1 (MQ=255)
                                                                           ccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGAc                 <  1:913899/100‑1 (MQ=255)
                                                                            cGCAGCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  2:2147825/97‑1 (MQ=255)
                                                                             gcgccgattggTCCGGGTTTTTTTGTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgg               <  1:420760/91‑1 (MQ=255)
                                                                             gCACCGAAAGGTCCGGGGTTTTTTTCTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAACTCTGTTTCTCAAGATTCAgg                   <  2:1468015/96‑1 (MQ=255)
                                                                               agcGAATGGTCCGGTTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  2:26527/98‑1 (MQ=255)
                                                                               aCCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  2:1591872/100‑1 (MQ=255)
                                                                                ccGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa            >  1:76411/1‑100 (MQ=255)
                                                                                  gacaGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggg              <  2:1285422/93‑1 (MQ=255)
                                                                                  gAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          >  1:1709937/1‑100 (MQ=255)
                                                                                  gAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          <  2:442845/100‑1 (MQ=255)
                                                                                  gAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          <  1:1598773/100‑1 (MQ=255)
                                                                                   aaaGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaact         <  1:1529695/100‑1 (MQ=255)
                                                                                     aGGTCCGGGGTTTTTTTTTTTTTTTTTAAACATAACCGAGGACCAGACAATGAATAACAGCACAAAATTCTTTTTCTCAAGATTCACGACGGGGaactaa       >  1:292940/1‑100 (MQ=255)
                                                                                      ggTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCCGTTTCTCAAGATTCAGGACGGGGaactaac      <  2:2380321/100‑1 (MQ=255)
                                                                                      ggTCCGGGGTTTTTTTTTGACCTGAAAAACATAACCGAGGAGCAGACAATGAATACCAGCACAAAACTCTGTTTCTCAAGATTCAGGACGGGGaacga        <  1:1620459/98‑3 (MQ=255)
                                                                                        tCCGGGTTTTTTTTTTGACCTTAACAACATAACCGAGGAGCAGGCAATGAATAACTGCACATAATTCTGTTTCTCAAGATTCAGGCCGGGGaactaac      <  2:866245/98‑1 (MQ=255)
                                                                                          cGGGGTTTTTTTTTGACCTGAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATg  <  1:1238302/100‑1 (MQ=255)
                                                                                                |                                                                                              
ATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  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: