Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,977,227 T→A 15.1% intergenic (‑88/+39) motA ← / ← flhC proton conductor component of flagella motor/flagellar class II regulon transcriptional activator, with FlhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,977,2270TA15.1% 248.1 / 35.4 100intergenic (‑88/+39)motA/flhCproton conductor component of flagella motor/flagellar class II regulon transcriptional activator, with FlhD
Reads supporting (aligned to +/- strand):  ref base T (41/43);  new base A (8/7);  total (50/50)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.79e-01

Only 100 of 106 total aligned reads displayed.

AAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGG  >  NC_000913/1977129‑1977326
                                                                                                    |                                                                                                   
aagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑cgtc                                                                                                    >  1:2400809/1‑100 (MQ=255)
 agataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑cgtcg                                                                                                   >  1:1126535/1‑100 (MQ=255)
    taagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑GGTCGTTa                                                                                                <  2:1354330/100‑1 (MQ=255)
     aagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTAc                                                                                               >  1:958694/1‑100 (MQ=255)
       gCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCg                                                                                             >  2:658000/1‑100 (MQ=255)
         aagACATCACACTTCAACTGAGGACCATGAAAGGATATGCAATCGTAAGGAGTTAACGCGCGATTAGAGCAAAAAAAAGCAGCGGT‑A‑cgtcgt                                                                                                  <  1:767425/91‑1 (MQ=255)
            aCATCATCCTTCCACTGTTGACCATGACATGATGTTCAGGCGTCAGGCGTTAACGCGCGAGTGGGGCAAAAAAAAGCAGCGGTAA‑CGACG‑TACCgctgct                                                                                        <  1:1558827/100‑1 (MQ=37)
            aCATCATCCTTCCACTGTCTACCATGAAAGGATGTTTAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCgctg                                                                                          <  1:2371316/98‑1 (MQ=255)
             catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAA‑CGACG‑TACCgctgctt                                                                                       <  2:2347856/100‑2 (MQ=38)
             catcatCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGctgctg                                                                                       >  1:2402711/1‑100 (MQ=255)
              atcatcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAA‑CGACG‑TACCgctgct                                                                                        <  1:801137/98‑1 (MQ=38)
                catcCTTCTACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAA‑CGACG‑TACCgctgcttttt                                                                                    <  2:1100017/100‑5 (MQ=37)
                catcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGaa                                                                                    >  2:95721/1‑100 (MQ=255)
                catcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGaa                                                                                    >  2:55105/1‑100 (MQ=255)
                 atcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAAt                                                                                   <  1:1345791/100‑1 (MQ=255)
                  tcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAA‑CGACG‑TACCgctgcttttttt                                                                                  <  2:227329/100‑7 (MQ=37)
                   ccTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAA‑CGACG‑TACCgctgctttttttt                                                                                 <  1:1893808/100‑8 (MQ=37)
                   ccTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGt                                                                                 >  2:986762/1‑100 (MQ=255)
                    cTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGtt                                                                                >  2:2260300/1‑100 (MQ=255)
                     ttCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTCCCGCTGCTGGAATGTTg                                                                               <  2:1598123/100‑1 (MQ=255)
                     ttCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTg                                                                               <  1:1939739/100‑1 (MQ=255)
                      tCCACTGTTGAACATGACAGGATGTTCAGTCGTTAGGCCTTAACGCGCGATTGGGGGCAAAAAAAGCAGCGGT‑A‑CGGCGTTACCGGTGCTGGAATGTTg                                                                               >  1:834592/1‑99 (MQ=255)
                       ccACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTgcg                                                                             >  1:1838895/1‑100 (MQ=255)
                       ccACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTgcg                                                                             >  2:865566/1‑100 (MQ=255)
                        cACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTgcgc                                                                            >  1:1490016/1‑100 (MQ=255)
                         aCTGTTGACCATTACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATCGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                           <  1:2298871/100‑1 (MQ=255)
                         aCTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGTAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                           <  1:265415/100‑1 (MQ=255)
                         aCTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                           <  2:1357552/100‑1 (MQ=255)
                          cTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                           <  1:1117643/99‑1 (MQ=255)
                                aCCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGt                                                                    <  1:233171/100‑1 (MQ=255)
                                 ccATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTa                                                                   <  2:1510102/100‑1 (MQ=255)
                                   aTGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATc                                                                 <  2:1666216/100‑1 (MQ=255)
                                   aTGAAAGGATGTTCAGTCGTTAGGTGTTAACGCGAGATTGAGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg                                                                     <  1:327227/96‑1 (MQ=255)
                                     gACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACGGTATc                                                                 <  1:433703/98‑1 (MQ=255)
                                     gACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAg                                                               <  2:457456/100‑1 (MQ=255)
                                      aCAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGt                                                              >  2:2110357/1‑100 (MQ=255)
                                       cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                             >  1:2150561/1‑100 (MQ=255)
                                       cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                             <  1:2110343/100‑1 (MQ=255)
                                         ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                           >  2:2087284/1‑100 (MQ=255)
                                         ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                           >  2:551350/1‑100 (MQ=255)
                                         ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATAAGTTa                                                            >  1:1541843/1‑99 (MQ=255)
                                           aTGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAAc                                                         <  2:1524857/100‑1 (MQ=255)
                                              ttCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTACTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGc                                                      >  2:1151531/1‑100 (MQ=255)
                                               tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                     <  2:1428361/100‑1 (MQ=255)
                                               tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                     <  1:1349618/100‑1 (MQ=255)
                                               tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                     <  2:905224/100‑1 (MQ=255)
                                                cAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCt                                                    <  2:2080711/100‑1 (MQ=255)
                                                cAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCt                                                    >  1:811069/1‑100 (MQ=255)
                                                 aGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCATTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTg                                                   >  1:900636/1‑100 (MQ=255)
                                                     gtcTGGAGTTCATGCGCGATCGGATCATAAAAAAGCAGCGGT‑A‑CGTAGTTACCGCTGTTGGAATGTTGCGCATAACCGTAACAGTTAAACAGACTGTAc                                                >  1:1533610/1‑99 (MQ=255)
                                                     gtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTAct                                               >  1:1992079/1‑100 (MQ=255)
                                                      tcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActc                                              <  2:1126539/100‑1 (MQ=255)
                                                      tcAGGCGTTAACGCGCGATGGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActc                                              <  2:2191122/100‑1 (MQ=255)
                                                       cAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGAAATGTTGCGCCTCACCGTATCAGTTACACAGCCTGTAct                                               >  1:2222907/1‑98 (MQ=255)
                                                        aGGCGTTAACGCGCGATAGGGACAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctc                                            <  1:1739353/100‑1 (MQ=255)
                                                         ggCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctct                                           <  2:2311288/100‑1 (MQ=255)
                                                          gCGTTCAGGCGCGATTGTTTAAAAAAAAACTAGCGGT‑A‑CGTCGCTAGCGCTGTTGGAATGTGGCGCCTCACCGTATCAGTTAAACAGCCTGTActct                                             <  1:485805/97‑1 (MQ=255)
                                                           cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGt                                         <  1:511556/100‑1 (MQ=255)
                                                           cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGt                                         >  2:1539027/1‑100 (MQ=255)
                                                           cGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGt                                         >  1:1084360/1‑100 (MQ=255)
                                                            gTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                        >  1:1489138/1‑100 (MQ=255)
                                                            gTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                        >  1:177566/1‑100 (MQ=255)
                                                            gTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                        <  1:657997/100‑1 (MQ=255)
                                                            gTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                        <  2:2331558/100‑1 (MQ=255)
                                                             ttAACGCGCGATTGGGGCAAAAAAAAGCGGCGGT‑A‑CGTCGTTACCGGTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                        <  1:985537/99‑1 (MQ=255)
                                                             ttAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTc                                       <  2:793453/100‑1 (MQ=255)
                                                               aaCTCGCGATTTTGGCAAAAAAAAGCAGCGGT‑A‑CGTCTTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                        <  1:527782/97‑1 (MQ=255)
                                                                 cgcgcgATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATcc                                   <  1:1152720/100‑1 (MQ=255)
                                                                  gcgcgATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCca                                  >  2:2404285/1‑100 (MQ=255)
                                                                  gcgcgATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCca                                  >  1:601276/1‑100 (MQ=255)
                                                                  gcgcgATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCca                                  <  2:177567/100‑1 (MQ=255)
                                                                  gcgcgATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCca                                  <  2:393521/100‑1 (MQ=255)
                                                                    gcgATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagc                                <  1:1259530/100‑1 (MQ=255)
                                                                        ttGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGtt                            <  2:792148/100‑1 (MQ=255)
                                                                         tGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtg                           >  2:1590335/1‑100 (MQ=255)
                                                                           gggCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtgtg                         >  2:1369089/1‑100 (MQ=255)
                                                                            ggCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTgg                        >  2:2144550/1‑100 (MQ=255)
                                                                             gCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTggg                       >  1:4111/1‑100 (MQ=255)
                                                                              cAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGa                      <  1:526362/100‑1 (MQ=255)
                                                                              cAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGa                      <  1:1451006/100‑1 (MQ=255)
                                                                                aaaaaaaGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACGGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGa                      >  1:150318/1‑98 (MQ=255)
                                                                                aaaaaaaGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGata                    >  1:1092742/1‑100 (MQ=255)
                                                                                aaaaaaaGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGata                    >  2:672112/1‑100 (MQ=255)
                                                                                 aaaaaaGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGataa                   >  2:1048150/1‑100 (MQ=255)
                                                                                     cagcagcGGTAA‑CGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc               >  2:1210875/2‑100 (MQ=38)
                                                                                     cagcagcGGTAA‑CGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc               >  2:132226/2‑100 (MQ=38)
                                                                                     aaGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc               >  2:2387582/1‑100 (MQ=255)
                                                                                      agcagcGGTAA‑CGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg              >  1:2155845/1‑100 (MQ=38)
                                                                                      agcagcGGTAA‑CGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg              >  1:732022/1‑100 (MQ=38)
                                                                                       gcagcGGTAA‑CGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg             >  1:1026352/1‑100 (MQ=255)
                                                                                       gcagcGGTAA‑CGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg             >  2:219053/1‑100 (MQ=255)
                                                                                       gcagcGGTAA‑CGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg             >  2:8811/1‑100 (MQ=255)
                                                                                       gcagcGGT‑ATCGACG‑TACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg             >  1:1528557/1‑100 (MQ=255)
                                                                                       gcagcGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg             >  2:413642/1‑100 (MQ=255)
                                                                                           cGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAgg         <  2:815399/100‑1 (MQ=255)
                                                                                           cGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAgg         <  1:2560906/100‑1 (MQ=255)
                                                                                            ggT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGa        >  1:760389/1‑100 (MQ=255)
                                                                                              t‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGAtt      >  2:1391081/1‑100 (MQ=255)
                                                                                                    tcgtTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCCGCAGTTGTGGGATAATATCGGCAGGATTCTgg  >  2:389971/1‑100 (MQ=255)
                                                                                                    tcgtTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTgg  <  1:2387567/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
AAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGT‑A‑CGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGATTCTGG  >  NC_000913/1977129‑1977326

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