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,4680GT98.0% 306.1 / ‑1.4 103intergenic (+48/‑39)ilvL/ilvXilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (50/50);  total (53/50)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.72e-01

Only 100 of 110 total aligned reads displayed.

TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950371‑3950562
                                                                                                  |                                                                                              
tgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑ttt                                                                                              <  2:941640/100‑1 (MQ=255)
 gATAATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTATGACTTCCGCACCGAAAGGTCCGGGG‑tttt                                                                                             <  2:1534046/100‑1 (MQ=255)
   ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttttt                                                                                           >  1:1434378/1‑100 (MQ=255)
     attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCACCGCACCGAAAGGTCCGGGG‑tttttttt                                                                                         >  1:1316438/1‑100 (MQ=255)
     agtatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttttttt                                                                                         <  2:1838662/98‑1 (MQ=255)
      ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGt‑ttttttttt                                                                                        <  1:2212999/100‑1 (MQ=255)
      ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑Ttttttttt                                                                                        <  2:470202/100‑1 (MQ=255)
       tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTg                                                                                       <  1:2810205/100‑1 (MQ=255)
         tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGAc                                                                                     >  1:1387966/1‑100 (MQ=255)
         tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGAc                                                                                     >  1:1748983/1‑100 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGAcat                                                                                   >  2:2000496/1‑98 (MQ=255)
           ccACCGTACGGGGCTAAACATGGACGAGGAAAGGCATAGAAATCAAGCCTTAACGAACTAAGACCCCCGCACCAAAAAATCCGGGG‑Tatttttt                                                                                         >  2:760979/1‑94 (MQ=255)
            cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCtt                                                                                  >  1:2396576/1‑100 (MQ=255)
             aCCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTa                                                                                 >  2:523827/1‑100 (MQ=255)
                gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaaa                                                                               <  1:2906711/99‑1 (MQ=255)
                    ggggCTGCACTTGGACGAGGAAAGGCTTAGAGTTCAAGCCTTAACGAATTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACa                                                                           <  1:1597694/99‑1 (MQ=255)
                       gattaaCTTGGACAAGAAAAGGCTTAGATATCAAGCCATAAAGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTCTTTTGACCTTAAAAACATAAc                                                                       <  1:437663/95‑1 (MQ=255)
                         tGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAAc                                                                       <  1:2857215/98‑1 (MQ=255)
                          gCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCga                                                                    >  1:502239/1‑100 (MQ=255)
                           cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgag                                                                   >  1:437081/1‑100 (MQ=255)
                           cACTTGGACGAGAAAAGGCTTAGAGATCAAGCCTTAACGAACTAATACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCga                                                                    <  2:1458571/99‑1 (MQ=255)
                            aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  >  1:1715306/1‑100 (MQ=255)
                            aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  <  2:1702684/100‑1 (MQ=255)
                             cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagga                                                                 >  1:350567/1‑100 (MQ=255)
                               tGGACGAGGAAATGCTTCGAGATCAAGCCTTAACGAACTTAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  <  2:427899/97‑1 (MQ=255)
                                ggACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGAAATTAAAAAAATAACCAAGGAGCa                                                              >  1:513763/1‑100 (MQ=255)
                                   cGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                             <  1:1604184/98‑1 (MQ=255)
                                   cGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGAc                                                           <  2:877518/100‑1 (MQ=255)
                                     aGGAAAGGCTTAGAGATCGAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACaa                                                         >  1:103912/1‑100 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCGTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        <  2:280230/100‑1 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACATTAAAAACATAACCGAGGAGCAGACAAt                                                        >  1:1475926/1‑100 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTATAAACTTAAAAAAATAAAAAAGGAacaaacaa                                                         >  1:2111479/1‑99 (MQ=255)
                                        aaaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGa                                                      <  2:2297737/100‑1 (MQ=255)
                                         aaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGaa                                                     <  2:1895795/100‑1 (MQ=255)
                                           ggCTTGGGGATCAAGCCTGAGCGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGa                                                      <  2:756893/97‑1 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                    >  2:111859/1‑100 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACTTTAAAAACATAACCGAGGAGCAGACAATGact                                                    >  2:940394/1‑97 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGAAATTAAAAACATAACAGAGGAGCAGAAAATGAATa                                                   >  1:2907480/1‑100 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGT‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 >  1:686608/1‑100 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  2:1298104/100‑1 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAAACACGGAGAAGACAATGAATAAc                                                 >  1:2991583/1‑100 (MQ=255)
                                              ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAATGGTCCGGGG‑TTTTTTTTTGACCATAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  1:940400/99‑1 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACTCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                               <  1:1454386/100‑1 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGTCCTTAAAAACATAACCGAGCAGCAAACAATGAATAAAAg                                               >  2:642162/1‑100 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                               >  2:84736/1‑100 (MQ=255)
                                                agatATCAAGTATTAACGTAGTAAGACCCCCGCACAGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGc                                              <  1:303529/100‑1 (MQ=255)
                                                 gataTCAACACTTAACGAAGTGAGGCCCCCGCTCTGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGc                                              <  1:1579445/96‑1 (MQ=255)
                                                   gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                           >  2:1519570/1‑100 (MQ=255)
                                                    aTCAAGCCTTAACGAACTAAGACCCCCGCCCCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAAACAAAAAACAAACAATGAATAACAGCACaa                                          >  1:471916/1‑100 (MQ=255)
                                                     tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaa                                         <  2:648212/100‑1 (MQ=255)
                                                        aGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                       >  2:937106/1‑99 (MQ=255)
                                                         gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAAAATAACCGAAGAGAAGACAATGAATAACAGCACAAAATTc                                     >  1:437504/1‑100 (MQ=255)
                                                         gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAACAGACAATGAATAAGAGAACAAAAtt                                      >  1:1492811/1‑98 (MQ=255)
                                                          ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAACACAAAAt                                       >  2:1815885/1‑97 (MQ=255)
                                                          ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCAAGGAGCAAACAATGAAAAACAGCACAAAATTc                                     >  1:2571404/1‑99 (MQ=255)
                                                           cTTAACGAACTAGGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   >  1:633458/1‑100 (MQ=255)
                                                            ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                                    <  1:2696815/98‑1 (MQ=255)
                                                            ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   <  1:3016912/99‑1 (MQ=255)
                                                             tAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGtt                                 >  2:2917711/1‑100 (MQ=255)
                                                              aaggAACTAAGACCCCCGCACCGAAAGGACCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                       <  2:239367/90‑1 (MQ=255)
                                                              aaCGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGttt                                >  1:501516/1‑100 (MQ=255)
                                                               aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                               >  1:2885742/1‑100 (MQ=255)
                                                                cGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                              >  1:720063/1‑100 (MQ=255)
                                                                cGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                              >  2:330826/1‑100 (MQ=255)
                                                                cGAACTAAGACCCCCGCAACGAAAGGTCCGGGG‑TTTTTTTTGGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                              <  2:2477748/100‑1 (MQ=255)
                                                                  aaCTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                            <  2:2148361/100‑1 (MQ=255)
                                                                    cTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACTTTAAAAAAATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTGTCAAg                          >  2:769247/1‑100 (MQ=255)
                                                                    cTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                           <  1:553715/99‑1 (MQ=255)
                                                                    cTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                          >  1:1125707/1‑100 (MQ=255)
                                                                     tAATACCCCCGGACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                            <  1:1357659/97‑1 (MQ=255)
                                                                      aaGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGAACAAAATTCTGTTTATCAAGAt                        >  1:647754/1‑100 (MQ=255)
                                                                      aaGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAt                        >  2:512998/1‑100 (MQ=255)
                                                                        gACCTCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAt                        <  1:1400499/98‑1 (MQ=255)
                                                                        gACCCCCGCACCGAAAGGTCCGGGT‑TTTTTTTTTGACCTTAAAAACATAACCAAGGAGCAAACAATGAATAAAAGCACAAAATTCTGTTTCTCAAGATTc                      >  1:2595408/1‑100 (MQ=255)
                                                                        gACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                      >  1:1818307/1‑100 (MQ=255)
                                                                         aCCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                     <  2:1178076/100‑1 (MQ=255)
                                                                          cccccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                    <  1:980875/100‑1 (MQ=255)
                                                                          cccccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                    <  2:1434364/100‑1 (MQ=255)
                                                                           ttccgatctgAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  2:1454853/91‑1 (MQ=255)
                                                                            cgcGCACTGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  1:585185/97‑1 (MQ=255)
                                                                             ccGCACTGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                     <  1:2941501/96‑1 (MQ=255)
                                                                             ccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGAc                 >  1:1783496/1‑100 (MQ=255)
                                                                              cGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACg                >  1:2330811/1‑100 (MQ=255)
                                                                                 gccGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  2:2495631/99‑1 (MQ=255)
                                                                                  ccGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  1:154544/99‑1 (MQ=255)
                                                                                  ccGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggg              >  2:1541293/1‑98 (MQ=255)
                                                                                  ccGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa            >  1:2502468/1‑100 (MQ=255)
                                                                                  ccGAAAGGTCCGGGG‑TTTTTTGTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  2:741238/99‑1 (MQ=255)
                                                                                   cGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaa           <  1:2653731/100‑1 (MQ=255)
                                                                                    gAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          <  1:1004387/100‑1 (MQ=255)
                                                                                    gAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          >  1:536411/1‑100 (MQ=255)
                                                                                       aGGGCCGGGG‑TTTTTTTTTGTCCTTAGAAACATAGCCGAGGAGCAGACAATGAATAACAGTACAATATTCTGTTTCTCAAGCTTCAGGACggg              <  2:1579419/93‑1 (MQ=255)
                                                                                        ggggcGGGG‑TTTTTTTTTGACCTTAAAAAGATAACCGAGGAGCAGACAATGAATAAGAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaac      <  2:1984949/96‑1 (MQ=255)
                                                                                        ggTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaac      <  2:73569/100‑1 (MQ=255)
                                                                                         gTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaact     <  2:2885710/100‑1 (MQ=255)
                                                                                         gTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaact     >  1:1949742/1‑100 (MQ=255)
                                                                                         gTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaac      <  1:2218854/99‑1 (MQ=255)
                                                                                         gTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACAGGGaactaa       >  1:754257/1‑98 (MQ=255)
                                                                                          tCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACTGAGGAGTAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAactaacta    <  1:109497/100‑1 (MQ=255)
                                                                                            ccggg‑gTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATg  >  2:1731172/3‑100 (MQ=255)
                                                                                                  |                                                                                              
TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950371‑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: ‑