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.

TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950371‑3950562
                                                                                                  |                                                                                              
tgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGttt                                                                                              <  2:941640/100‑1 (MQ=255)
 gATAATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTATGACTTCCGCACCGAAAGGTCCGGGGtttt                                                                                             <  2:1534046/100‑1 (MQ=255)
   ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttttt                                                                                           >  1:1434378/1‑100 (MQ=255)
     attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCACCGCACCGAAAGGTCCGGGGtttttttt                                                                                         >  1:1316438/1‑100 (MQ=255)
     agtatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttttttt                                                                                         <  2:1838662/98‑1 (MQ=255)
      ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGtttttttttt                                                                                        <  1:2212999/100‑1 (MQ=255)
      ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttttt                                                                                        <  2:470202/100‑1 (MQ=255)
       tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTg                                                                                       <  1:2810205/100‑1 (MQ=255)
         tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAc                                                                                     >  1:1387966/1‑100 (MQ=255)
         tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAc                                                                                     >  1:1748983/1‑100 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAcat                                                                                   >  2:2000496/1‑98 (MQ=255)
           ccACCGTACGGGGCTAAACATGGACGAGGAAAGGCATAGAAATCAAGCCTTAACGAACTAAGACCCCCGCACCAAAAAATCCGGGGTatttttt                                                                                         >  2:760979/1‑94 (MQ=255)
            cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCtt                                                                                  >  1:2396576/1‑100 (MQ=255)
             aCCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTa                                                                                 >  2:523827/1‑100 (MQ=255)
                gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaa                                                                               <  1:2906711/99‑1 (MQ=255)
                    ggggCTGCACTTGGACGAGGAAAGGCTTAGAGTTCAAGCCTTAACGAATTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACa                                                                           <  1:1597694/99‑1 (MQ=255)
                       gattaaCTTGGACAAGAAAAGGCTTAGATATCAAGCCATAAAGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTCTTTTGACCTTAAAAACATAAc                                                                       <  1:437663/95‑1 (MQ=255)
                         tGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAc                                                                       <  1:2857215/98‑1 (MQ=255)
                          gCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCga                                                                    >  1:502239/1‑100 (MQ=255)
                           cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgag                                                                   >  1:437081/1‑100 (MQ=255)
                           cACTTGGACGAGAAAAGGCTTAGAGATCAAGCCTTAACGAACTAATACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCga                                                                    <  2:1458571/99‑1 (MQ=255)
                            aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  >  1:1715306/1‑100 (MQ=255)
                            aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  <  2:1702684/100‑1 (MQ=255)
                             cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagga                                                                 >  1:350567/1‑100 (MQ=255)
                               tGGACGAGGAAATGCTTCGAGATCAAGCCTTAACGAACTTAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  <  2:427899/97‑1 (MQ=255)
                                ggACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAAATTAAAAAAATAACCAAGGAGCa                                                              >  1:513763/1‑100 (MQ=255)
                                   cGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                             <  1:1604184/98‑1 (MQ=255)
                                   cGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGAc                                                           <  2:877518/100‑1 (MQ=255)
                                     aGGAAAGGCTTAGAGATCGAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACaa                                                         >  1:103912/1‑100 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCGTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        <  2:280230/100‑1 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACATTAAAAACATAACCGAGGAGCAGACAAt                                                        >  1:1475926/1‑100 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTATAAACTTAAAAAAATAAAAAAGGAacaaacaa                                                         >  1:2111479/1‑99 (MQ=255)
                                        aaaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGa                                                      <  2:2297737/100‑1 (MQ=255)
                                         aaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGaa                                                     <  2:1895795/100‑1 (MQ=255)
                                           ggCTTGGGGATCAAGCCTGAGCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGa                                                      <  2:756893/97‑1 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                    >  2:111859/1‑100 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACTTTAAAAACATAACCGAGGAGCAGACAATGact                                                    >  2:940394/1‑97 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAAATTAAAAACATAACAGAGGAGCAGAAAATGAATa                                                   >  1:2907480/1‑100 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 >  1:686608/1‑100 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  2:1298104/100‑1 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAACACGGAGAAGACAATGAATAAc                                                 >  1:2991583/1‑100 (MQ=255)
                                              ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAATGGTCCGGGGTTTTTTTTTGACCATAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                 <  1:940400/99‑1 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACTCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                               <  1:1454386/100‑1 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGTCCTTAAAAACATAACCGAGCAGCAAACAATGAATAAAAg                                               >  2:642162/1‑100 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                               >  2:84736/1‑100 (MQ=255)
                                                agatATCAAGTATTAACGTAGTAAGACCCCCGCACAGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGc                                              <  1:303529/100‑1 (MQ=255)
                                                 gataTCAACACTTAACGAAGTGAGGCCCCCGCTCTGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGc                                              <  1:1579445/96‑1 (MQ=255)
                                                   gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                           >  2:1519570/1‑100 (MQ=255)
                                                    aTCAAGCCTTAACGAACTAAGACCCCCGCCCCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAACAAAAAACAAACAATGAATAACAGCACaa                                          >  1:471916/1‑100 (MQ=255)
                                                     tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaa                                         <  2:648212/100‑1 (MQ=255)
                                                        aGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                       >  2:937106/1‑99 (MQ=255)
                                                         gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAAAATAACCGAAGAGAAGACAATGAATAACAGCACAAAATTc                                     >  1:437504/1‑100 (MQ=255)
                                                         gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAACAGACAATGAATAAGAGAACAAAAtt                                      >  1:1492811/1‑98 (MQ=255)
                                                          ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAACACAAAAt                                       >  2:1815885/1‑97 (MQ=255)
                                                          ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCAAGGAGCAAACAATGAAAAACAGCACAAAATTc                                     >  1:2571404/1‑99 (MQ=255)
                                                           cTTAACGAACTAGGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   >  1:633458/1‑100 (MQ=255)
                                                            ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                                    <  1:2696815/98‑1 (MQ=255)
                                                            ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   <  1:3016912/99‑1 (MQ=255)
                                                             tAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGtt                                 >  2:2917711/1‑100 (MQ=255)
                                                              aaggAACTAAGACCCCCGCACCGAAAGGACCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                       <  2:239367/90‑1 (MQ=255)
                                                              aaCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGttt                                >  1:501516/1‑100 (MQ=255)
                                                               aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                               >  1:2885742/1‑100 (MQ=255)
                                                                cGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                              >  1:720063/1‑100 (MQ=255)
                                                                cGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                              >  2:330826/1‑100 (MQ=255)
                                                                cGAACTAAGACCCCCGCAACGAAAGGTCCGGGGTTTTTTTTGGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                              <  2:2477748/100‑1 (MQ=255)
                                                                  aaCTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                            <  2:2148361/100‑1 (MQ=255)
                                                                    cTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACTTTAAAAAAATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTGTCAAg                          >  2:769247/1‑100 (MQ=255)
                                                                    cTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                           <  1:553715/99‑1 (MQ=255)
                                                                    cTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                          >  1:1125707/1‑100 (MQ=255)
                                                                     tAATACCCCCGGACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                            <  1:1357659/97‑1 (MQ=255)
                                                                      aaGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGAACAAAATTCTGTTTATCAAGAt                        >  1:647754/1‑100 (MQ=255)
                                                                      aaGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAt                        >  2:512998/1‑100 (MQ=255)
                                                                        gACCTCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAt                        <  1:1400499/98‑1 (MQ=255)
                                                                        gACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCAAGGAGCAAACAATGAATAAAAGCACAAAATTCTGTTTCTCAAGATTc                      >  1:2595408/1‑100 (MQ=255)
                                                                        gACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                      >  1:1818307/1‑100 (MQ=255)
                                                                         aCCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                     <  2:1178076/100‑1 (MQ=255)
                                                                          cccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                    <  1:980875/100‑1 (MQ=255)
                                                                          cccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                    <  2:1434364/100‑1 (MQ=255)
                                                                           ttccgatctgAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  2:1454853/91‑1 (MQ=255)
                                                                            cgcGCACTGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  1:585185/97‑1 (MQ=255)
                                                                             ccGCACTGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                     <  1:2941501/96‑1 (MQ=255)
                                                                             ccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGAc                 >  1:1783496/1‑100 (MQ=255)
                                                                              cGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACg                >  1:2330811/1‑100 (MQ=255)
                                                                                 gccGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  2:2495631/99‑1 (MQ=255)
                                                                                  ccGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  1:154544/99‑1 (MQ=255)
                                                                                  ccGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggg              >  2:1541293/1‑98 (MQ=255)
                                                                                  ccGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa            >  1:2502468/1‑100 (MQ=255)
                                                                                  ccGAAAGGTCCGGGGTTTTTTGTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg             <  2:741238/99‑1 (MQ=255)
                                                                                   cGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaa           <  1:2653731/100‑1 (MQ=255)
                                                                                    gAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          <  1:1004387/100‑1 (MQ=255)
                                                                                    gAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac          >  1:536411/1‑100 (MQ=255)
                                                                                       aGGGCCGGGGTTTTTTTTTGTCCTTAGAAACATAGCCGAGGAGCAGACAATGAATAACAGTACAATATTCTGTTTCTCAAGCTTCAGGACggg              <  2:1579419/93‑1 (MQ=255)
                                                                                        ggggcGGGGTTTTTTTTTGACCTTAAAAAGATAACCGAGGAGCAGACAATGAATAAGAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaac      <  2:1984949/96‑1 (MQ=255)
                                                                                        ggTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaac      <  2:73569/100‑1 (MQ=255)
                                                                                         gTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaact     <  2:2885710/100‑1 (MQ=255)
                                                                                         gTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaact     >  1:1949742/1‑100 (MQ=255)
                                                                                         gTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaac      <  1:2218854/99‑1 (MQ=255)
                                                                                         gTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACAGGGaactaa       >  1:754257/1‑98 (MQ=255)
                                                                                          tCCGGGGTTTTTTTTTGACCTTAAAAACATAACTGAGGAGTAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAactaacta    <  1:109497/100‑1 (MQ=255)
                                                                                            ccggggTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATg  >  2:1731172/3‑100 (MQ=255)
                                                                                                  |                                                                                              
TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  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: