Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,977,224:1 +A 15.1% intergenic (‑85/+42) 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,2231.A15.1% 538.2 / 103.2 166intergenic (‑84/+43)motA/flhCproton conductor component of flagella motor/flagellar class II regulon transcriptional activator, with FlhD
Reads supporting (aligned to +/- strand):  ref base . (73/68);  new base A (17/8);  total (90/76)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 170 total aligned reads displayed.

CTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGAT  >  NC_000913/1977124‑1977321
                                                                                                     |                                                                                                  
cTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGt                                                                                                     <  1:1927296/100‑1 (MQ=255)
 taataaGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTa                                                                                                   >  2:862414/1‑100 (MQ=255)
    taagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAcgt                                                                                                >  2:912159/1‑100 (MQ=255)
     aagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAcgtc                                                                                               <  2:3784941/100‑1 (MQ=255)
     aagataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAGCAGCGGTAcgtcg                                                                                              <  2:262738/100‑1 (MQ=255)
      agctaagCTCGACATAATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCTTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAGCAGCGGTAcg                                                                                                 >  1:4126260/4‑96 (MQ=255)
      agataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAcgtcg                                                                                              >  2:738568/1‑100 (MQ=255)
       gataagCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAATGCGCGATTGGGGCAAAAACAAGCAGCGGTAcgtcgt                                                                                             >  1:491027/1‑100 (MQ=255)
           agCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGCTTGGGGCAAAAAAAAGCAGCGGTACGTCGttcc                                                                                          >  1:1000988/1‑97 (MQ=255)
            gCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCg                                                                                        >  2:1286207/1‑100 (MQ=255)
            gCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCg                                                                                        <  2:1614235/100‑1 (MQ=255)
              aCGACATCACCCTTTCACAATTGACCAGGATATGATGAACAGTATTCAGGCGTTAACGCGCAAAGGGCGCAAAAAAAAGGAGCGATACGTCGATAca                                                                                         >  1:1389207/1‑96 (MQ=255)
                gACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCgctgc                                                                                    >  1:2647985/1‑100 (MQ=255)
                   atcatcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTgg                                                                                 >  2:2705593/1‑100 (MQ=255)
                    tcatcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGa                                                                                <  2:2509818/100‑1 (MQ=255)
                     catcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAGGTCGTTACCGCTGCTGGaa                                                                               >  2:4100677/1‑100 (MQ=255)
                      atcCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAAt                                                                              <  2:52052/100‑1 (MQ=255)
                        ccTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGt                                                                            <  2:1247236/100‑1 (MQ=255)
                         cTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACGTACCgctgcttttttttg                                                                           <  1:3006528/100‑9 (MQ=25)
                         cTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGtt                                                                           >  2:403811/1‑100 (MQ=255)
                         cTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGtt                                                                           >  1:512165/1‑100 (MQ=255)
                         cTTCCACTGTTGAACATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACGTACCgctgcttttttttg                                                                           <  1:809219/100‑9 (MQ=25)
                          ttCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTAACGACGTACCgctgcttttttttgc                                                                          <  1:2191417/100‑10 (MQ=25)
                           tCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTgc                                                                         <  1:2255645/100‑1 (MQ=255)
                            ccACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAGAGCAGTGGTACGTCGTTACCGCTGCTGGAATGTTgcg                                                                        >  1:2973724/1‑100 (MQ=255)
                             cACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTgcgc                                                                       <  1:3191574/100‑1 (MQ=255)
                             cACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTgc                                                                         >  1:2670014/1‑98 (MQ=255)
                              aCTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                      <  2:2742308/100‑1 (MQ=255)
                                tGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTc                                                                    <  2:1074449/100‑1 (MQ=255)
                                tGTTAACCATGACAGAATGTTCAGTCGTCAGGTGTTAACGCGAGAAAAAAAAAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGcc                                                                      <  1:3164821/98‑1 (MQ=255)
                                   tGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCAcc                                                                 <  1:3273999/100‑1 (MQ=255)
                                    gACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg                                                                >  2:2817769/1‑100 (MQ=255)
                                    gACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCg                                                                >  2:1434168/1‑100 (MQ=255)
                                     aCCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGt                                                               >  2:792578/1‑100 (MQ=255)
                                       cATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTAt                                                             <  2:2296615/100‑1 (MQ=255)
                                            naGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                        >  2:3230546/2‑100 (MQ=255)
                                            cAGGCTGTTCAGTCGGCAGGCGTTAACGCGCGATTGGTGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAg                                                          <  1:3185037/98‑1 (MQ=255)
                                            cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                        >  1:662935/1‑100 (MQ=255)
                                            cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGtt                                                        >  2:3342547/1‑100 (MQ=255)
                                            cAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGt                                                         >  2:293495/1‑99 (MQ=255)
                                              ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                      >  2:1061766/1‑100 (MQ=255)
                                              ggATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTaa                                                      >  1:3624353/1‑100 (MQ=255)
                                               gATGTTCAGCCGTCAGGCGTTCACGCGCGATTGGGGCAAAAAAAAGCACCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCATCACCGTACCAGTTaaa                                                     >  1:3500116/1‑100 (MQ=255)
                                                 tGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACa                                                   <  1:847664/100‑1 (MQ=255)
                                                  gTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAg                                                  >  1:494424/1‑100 (MQ=255)
                                                  gTTCAGTAGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAg                                                  >  1:1521910/1‑100 (MQ=255)
                                                    tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                <  1:4363352/100‑1 (MQ=255)
                                                    tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGcc                                                <  1:4100647/100‑1 (MQ=255)
                                                    tCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGACTCACCGTATCAGTTAAACAGcc                                                >  1:2179996/1‑100 (MQ=255)
                                                     cAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCt                                               <  2:4026147/100‑1 (MQ=255)
                                                     cAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCt                                               <  1:2167413/100‑1 (MQ=255)
                                                      aGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTCAACAGCCTg                                              >  1:1652773/1‑100 (MQ=255)
                                                      aGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTg                                              <  2:275132/100‑1 (MQ=255)
                                                         cgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTAc                                           >  2:2230658/1‑100 (MQ=255)
                                                         cgtcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTAc                                           <  2:3188870/100‑1 (MQ=255)
                                                          gtcAGGCGTTAACGCGCGAATGGGGCAAAAAAAAGCAGCGTTACGTCGTTACCGCTGCTGGAATGGAGCGCCTCACCGTATCAATTAAACAGCATATAc                                           >  1:886996/1‑99 (MQ=255)
                                                           tcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActc                                         >  1:24522/1‑100 (MQ=255)
                                                           tcAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActc                                         <  2:288056/100‑1 (MQ=255)
                                                             aGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctc                                       <  2:3619457/100‑1 (MQ=255)
                                                             aGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctc                                       >  2:2536291/1‑100 (MQ=255)
                                                              ggCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctct                                      >  1:4008438/1‑100 (MQ=255)
                                                              ggCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctct                                      >  2:3223993/1‑100 (MQ=255)
                                                               gcgcTAACGCGCGATTTGGGCAAATAAACGCAGCGGTACGTCGTTACCGCTGCTGGAATCTTGCGCCTCACCGTATCAGTTAAACAGCCTGTActctc                                       <  1:3707325/98‑1 (MQ=255)
                                                               gCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTg                                     >  2:2955454/1‑100 (MQ=255)
                                                                 gTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                   >  1:3606137/1‑100 (MQ=255)
                                                                 gTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                   <  1:3723327/100‑1 (MQ=255)
                                                                 gTTAACGCGCGATTGAGGCAAAAAAAAGCAGCGGTACGTCTTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGtt                                   <  1:429274/100‑1 (MQ=255)
                                                                 gTTAACGCGCGAGTGGGGCAAAAAAAAGCAGCGGTACGTAGTTACCGCTGCTGGAATGGTGCGACTAACCGTATCAGTTAAACAGCCTGTActctc                                       >  1:4095736/1‑96 (MQ=255)
                                                                   tAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCa                                 >  1:2601866/1‑100 (MQ=255)
                                                                    aaCGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCAt                                <  2:4121323/100‑1 (MQ=255)
                                                                     aCGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATc                               >  2:3598849/1‑100 (MQ=255)
                                                                     aCGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATc                               <  1:173392/100‑1 (MQ=255)
                                                                      cgcgcgATTGGGGCAAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATc                               >  1:1570238/1‑100 (MQ=255)
                                                                      cgcgcgATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATcc                              >  1:3424857/1‑100 (MQ=255)
                                                                       gcgcgATTGGGGCAAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATcc                              <  2:1963603/100‑1 (MQ=255)
                                                                         gcgATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCcagc                           <  2:2930573/100‑1 (MQ=255)
                                                                              tGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtg                      >  2:2805770/1‑100 (MQ=255)
                                                                               ggggCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTtgt                     >  2:3377789/1‑100 (MQ=255)
                                                                                 ggCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTgg                   >  2:2342378/1‑100 (MQ=255)
                                                                                 ggCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTgg                   >  2:1319466/1‑100 (MQ=255)
                                                                                  gCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTggg                  <  1:222784/100‑1 (MQ=255)
                                                                                  agaaaaaaaaaGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTgg                   <  1:2845057/97‑1 (MQ=255)
                                                                                     aaaaaaaaGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTggg                  <  1:4277193/98‑1 (MQ=255)
                                                                                     aaaaaaaaGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGat                >  2:3618121/1‑100 (MQ=255)
                                                                                        aaaaaGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTACGCCTAACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGata               >  1:1193137/1‑98 (MQ=255)
                                                                                         aaaaGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGataata            <  2:2468616/100‑1 (MQ=255)
                                                                                          aaaGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAatat           >  2:1019479/1‑100 (MQ=255)
                                                                                           cagcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc          >  1:2304111/2‑100 (MQ=38)
                                                                                           cagcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc          >  2:3431643/2‑100 (MQ=38)
                                                                                           cagcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc          >  2:3551173/2‑100 (MQ=38)
                                                                                           cagcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATc          >  2:2265792/2‑100 (MQ=38)
                                                                                            agcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg         >  1:2612233/1‑100 (MQ=38)
                                                                                            agcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg         >  2:1852355/1‑100 (MQ=38)
                                                                                            agcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg         >  2:2139004/1‑100 (MQ=38)
                                                                                            agcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg         >  1:328623/1‑100 (MQ=38)
                                                                                            agcagcGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCg         <  2:4171874/100‑1 (MQ=255)
                                                                                             gcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg        >  1:2092852/1‑100 (MQ=255)
                                                                                             gcagcGGTAACGACGTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg        >  2:1642105/1‑100 (MQ=255)
                                                                                             gcagcGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCgg        >  2:2043475/1‑100 (MQ=255)
                                                                                                   gTACGTCGTTACCGCTGCGGGAATGTTGCGCCGCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCa      <  1:3661453/96‑1 (MQ=255)
                                                                                                     |                                                                                                  
CTAATAAGATAAGCACGACATCATCCTTCCACTGTTGACCATGACAGGATGTTCAGTCGTCAGGCGTTAACGCGCGATTGGGGCAAAAAAAAGCAGCGGTACGTCGTTACCGCTGCTGGAATGTTGCGCCTCACCGTATCAGTTAAACAGCCTGTACTCTCTGTTCATCCAGCAGTTGTGGGATAATATCGGCAGGAT  >  NC_000913/1977124‑1977321

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: