Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,098,1550AC25.0% 173.6 / 11.2 80L140F (TTA→TTCrhaSDNA‑binding transcriptional activator RhaS
Reads supporting (aligned to +/- strand):  ref base A (36/24);  new base C (0/20);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 6.07e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGAT  >  NC_000913/4098019‑4098293
                                                                                                                                        |                                                                                                                                          
tcGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTa                                                                                                                                            <  2:1350588/137‑1 (MQ=255)
tcGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTa                                                                                                                                            >  1:1350588/1‑137 (MQ=255)
  gCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTc                                                                                                                                       <  2:943570/140‑1 (MQ=255)
      ggCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGAcc                                                                                                                                   >  2:1136620/1‑140 (MQ=255)
       gCTGAATCAGTGGCTGCAACAAGAGCTGGATGGGCAGTATCCGTCTCCCTGGCCAGTTACCCACGGGGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCg                                                                                                                                  <  2:1236553/140‑1 (MQ=255)
        cTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGc                                                                                                                                 <  1:461664/140‑1 (MQ=255)
          gAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCa                                                                                                                               <  2:1316069/140‑1 (MQ=255)
          gAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCa                                                                                                                               <  2:747085/140‑1 (MQ=255)
             tCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTc                                                                                                                            <  1:1129787/140‑1 (MQ=255)
               atttccggcCACAAGATCTGGATGGGCAGTATCCGTCTCACTGGCCCGTTAACCACAGCGTATTGCAGCAGGTGCAACATCTGGTTGCACAGATGGAACAGCCGGAAGGGGAAAATGATTTCCCCTCAACCGCCAGTcgc                                                                                                                          <  2:556424/133‑1 (MQ=255)
                 ttGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCga                                                                                                                        >  1:1138186/1‑140 (MQ=255)
                 ttGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCga                                                                                                                        >  2:342245/1‑140 (MQ=255)
                    ctgcCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGAt                                                                                                                     <  1:253392/140‑1 (MQ=255)
                     tgcCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACTACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATc                                                                                                                    <  1:742108/140‑1 (MQ=255)
                     tgcCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTc                                                                                                                            <  1:766126/132‑1 (MQ=255)
                     tgcCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTc                                                                                                                            >  2:766126/1‑132 (MQ=255)
                      gcCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCt                                                                                                                   <  1:398950/140‑1 (MQ=255)
                         acaAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTcgcg                                                                                                                         >  1:155948/1‑131 (MQ=255)
                         acaAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTcgcg                                                                                                                         <  2:155948/131‑1 (MQ=255)
                         acaAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGt                                                                                                                >  2:316524/1‑140 (MQ=255)
                          caAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATATACCCTCGACCGCCAGTCGCGAGATCTTGtt                                                                                                               <  2:419331/140‑1 (MQ=255)
                           aaGAGCTGGAGGGGCAGTATCCGTCTCACTGGGGCGTTACCCACAGCGTATTGCGGAAGGTGGGACTGCTGGTTGCACAGAGGAAACAGCGGGAAGGGGAAAATGATTTTCCCGCGACCGCCAGTCGCGAGATCTTGttt                                                                                                              <  2:674943/140‑1 (MQ=255)
                             gagCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTAt                                                                                                            <  2:259875/140‑1 (MQ=255)
                                        gCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGc                                                                                                 >  2:1213877/1‑140 (MQ=255)
                                         cAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGAAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGct                                                                                                <  2:1257496/140‑1 (MQ=255)
                                            tATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCtt                                                                                             >  2:630685/1‑140 (MQ=255)
                                             aTCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtg                                                                                            >  1:912229/1‑140 (MQ=255)
                                             aTCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtg                                                                                            >  1:822129/1‑140 (MQ=255)
                                              tCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtgc                                                                                           <  2:625288/140‑1 (MQ=255)
                                                    tcACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTa                                                                                     >  1:1314086/1‑140 (MQ=255)
                                                      aCTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTaaa                                                                                   >  2:294519/1‑140 (MQ=255)
                                                      aCTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTaaa                                                                                   >  1:219271/1‑140 (MQ=255)
                                                      aCTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTaaa                                                                                   >  1:863551/1‑140 (MQ=255)
                                                        tGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAg                                                                                 <  2:641035/140‑1 (MQ=255)
                                                        gggCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGTTTCCACAGATGAAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAg                                                                                 <  2:1314086/139‑1 (MQ=255)
                                                          gcgcGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCa                                                                               <  1:1093543/140‑1 (MQ=255)
                                                          gcgcGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCa                                                                               <  2:863551/140‑1 (MQ=255)
                                                          gcgcGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCa                                                                               <  2:1138186/140‑1 (MQ=255)
                                                          gcgcGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCa                                                                               <  1:905435/140‑1 (MQ=255)
                                                          gcgcGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCa                                                                               <  1:1266919/140‑1 (MQ=255)
                                                           cgcgTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAg                                                                              >  1:1332610/1‑140 (MQ=255)
                                                             cgTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGtt                                                                            <  2:1075090/140‑1 (MQ=255)
                                                             cgTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGtt                                                                            <  2:467306/140‑1 (MQ=255)
                                                                  aCCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAg                                                                       <  2:1156771/140‑1 (MQ=255)
                                                                   ccACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAgg                                                                      <  2:768419/140‑1 (MQ=255)
                                                                    cacaGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGga                                                                     >  1:327333/1‑140 (MQ=255)
                                                                      caGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGgaga                                                                   <  2:219271/140‑1 (MQ=255)
                                                                      caGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGgaga                                                                   <  2:1163611/140‑1 (MQ=255)
                                                                      caGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGgaga                                                                   <  1:342245/140‑1 (MQ=255)
                                                                      caGCGTATGGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAGGGGGAAAATGATTTACCCTCGACCGCCAGTCGCTAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGgaga                                                                   <  2:821099/140‑1 (MQ=255)
                                                                      aaGCGTATTGCAGCAAGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGgaga                                                                   <  2:577221/139‑1 (MQ=255)
                                                                       aGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGaa                                                                  >  2:961221/1‑140 (MQ=255)
                                                                       aGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGaa                                                                  >  1:408691/1‑140 (MQ=255)
                                                                        gCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAAc                                                                 <  1:1136620/140‑1 (MQ=255)
                                                                         cGTATTGCAGCAGGTGCGACAGCTTGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAAcc                                                                >  2:1354800/1‑140 (MQ=255)
                                                                            aTTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTgg                                                             >  1:925727/1‑140 (MQ=255)
                                                                                   cagGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAgc                                                     >  1:423927/1‑140 (MQ=255)
                                                                                     ggTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAgcg                                                    >  1:438045/1‑140 (MQ=255)
                                                                                       tGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCa                                                  <  2:408691/140‑1 (MQ=255)
                                                                                        gCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCAt                                                 >  2:477907/1‑140 (MQ=255)
                                                                                        gCGACAGCTGGTGCCAAAGAGGGAACAGCAGGAAGGGAAAAATGTTTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCAt                                                 <  2:925727/140‑1 (MQ=255)
                                                                                         cGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATc                                                <  1:1213877/140‑1 (MQ=255)
                                                                                          gACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCa                                               <  1:961221/140‑1 (MQ=255)
                                                                                               ttGGTTGCACAGATGGAACAGCAGGAAGGGAAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGtct                                          <  2:327333/139‑1 (MQ=255)
                                                                                                tGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGtctc                                         >  1:767685/1‑140 (MQ=255)
                                                                                                  gctGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCaa                                       >  1:1237246/3‑140 (MQ=255)
                                                                                                    tGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACt                                     <  2:923625/140‑1 (MQ=255)
                                                                                                    tGCACAGATGGAACAGCAGGAAGGGAAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACt                                     <  2:335568/140‑1 (MQ=255)
                                                                                                              gAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctggc                           <  1:47250/140‑1 (MQ=255)
                                                                                                                    cagGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggct                     >  2:777945/1‑140 (MQ=255)
                                                                                                                     agGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggctg                    <  1:630685/140‑1 (MQ=255)
                                                                                                                     agGAAGGGAAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggctg                    <  1:1354800/140‑1 (MQ=255)
                                                                                                                       gAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTgga                  >  1:330549/1‑140 (MQ=255)
                                                                                                                         gggggAAAATGATTTCCCCCCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTggagg                <  2:482665/139‑1 (MQ=255)
                                                                                                                         aGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTggagg                >  2:801863/1‑140 (MQ=255)
                                                                                                                         aGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTggagg                >  1:1208253/1‑140 (MQ=255)
                                                                                                                         aGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTggagg                >  1:988519/1‑140 (MQ=255)
                                                                                                                            ggAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGAcc             >  1:1211582/1‑140 (MQ=255)
                                                                                                                              aaaaTGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCACCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCAt           >  1:425616/1‑140 (MQ=255)
                                                                                                                                aaTGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCAttt         <  1:316524/140‑1 (MQ=255)
                                                                                                                                  tGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTg       <  1:777945/140‑1 (MQ=255)
                                                                                                                                     tttACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCg    >  2:619021/1‑140 (MQ=255)
                                                                                                                                     tttACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCg    >  1:383807/1‑140 (MQ=255)
                                                                                                                                     tttACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCg    >  2:1205339/1‑140 (MQ=255)
                                                                                                                                     tttACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCg    >  1:832914/1‑140 (MQ=255)
                                                                                                                                       tACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGAt  <  2:988519/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGAT  >  NC_000913/4098019‑4098293

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: