Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,098,1550AC24.7% 181.4 / 13.7 81L140F (TTA→TTC) rhaSDNA‑binding transcriptional activator RhaS
Reads supporting (aligned to +/- strand):  ref base A (41/20);  new base C (0/20);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 2.94e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGA  >  NC_000913/4098016‑4098292
                                                                                                                                           |                                                                                                                                         
ttCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTa                                                                                                                                           <  1:850998/140‑1 (MQ=255)
  ctcGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTAcc                                                                                                                                         <  1:778038/140‑1 (MQ=255)
   tcGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTAccc                                                                                                                                        >  2:79875/1‑140 (MQ=255)
   tcGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTAccc                                                                                                                                        >  2:577350/1‑140 (MQ=255)
    cGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCt                                                                                                                                       <  1:615331/140‑1 (MQ=255)
    cGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGACAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCt                                                                                                                                       >  2:1056064/1‑140 (MQ=255)
         ggCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGG‑AAATGATTTACCCTCGACCg                                                                                                                                 <  1:688068/140‑1 (MQ=255)
         ggCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGAAGGAAGGGG‑AAATGATTTACCCTCGACCg                                                                                                                                 <  1:1103898/140‑1 (MQ=255)
             gAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCa                                                                                                                              <  2:908980/140‑1 (MQ=255)
              aaTCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAg                                                                                                                             <  1:382440/140‑1 (MQ=255)
                 cAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGAGGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTcg                                                                                                                          <  2:527208/140‑1 (MQ=255)
                  ggTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGGGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTcgc                                                                                                                         <  2:360604/139‑1 (MQ=255)
                  aGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTcgc                                                                                                                         <  1:632809/140‑1 (MQ=255)
                  aGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTcgc                                                                                                                         <  2:326102/140‑1 (MQ=255)
                  aGTTGCTGCCACAAGAGCTGGATGGACAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTcgc                                                                                                                         <  1:1056064/140‑1 (MQ=255)
                    ttGCTGCCACAAGAGCCGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCga                                                                                                                       >  1:1012024/1‑140 (MQ=255)
                       ctgcCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACGGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGAt                                                                                                                    <  2:186557/140‑1 (MQ=255)
                         gcCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCt                                                                                                                  >  2:1064014/1‑140 (MQ=255)
                          ccACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTCACCCTCGACCGCCAGTCGCGAGATCtt                                                                                                                 >  1:1153345/1‑140 (MQ=255)
                               agagCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTa                                                                                                            <  2:141362/140‑1 (MQ=255)
                               agagCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTa                                                                                                            <  1:412248/140‑1 (MQ=255)
                                  gCTGGATGGGCAGTACCCGTCTCACGGGCGCTTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGAAACAGCAGGAAGGGAAAAATGATTCCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGc                                                                                                         <  2:1153345/140‑1 (MQ=255)
                                   cTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCa                                                                                                        <  1:415341/140‑1 (MQ=255)
                                      gATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAAtt                                                                                                     <  1:997915/140‑1 (MQ=255)
                                            cAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGct                                                                                               <  2:944782/140‑1 (MQ=255)
                                            cAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGct                                                                                               >  2:256175/1‑140 (MQ=255)
                                                aTCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtg                                                                                           <  2:351038/140‑1 (MQ=255)
                                                aTCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGCGATCTTGTTTATGCAATTACTGCTCTtg                                                                                           >  2:605621/1‑140 (MQ=255)
                                                aTCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtg                                                                                           >  2:461434/1‑140 (MQ=255)
                                                aTCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtg                                                                                           >  1:841353/1‑140 (MQ=255)
                                                 tCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtgc                                                                                          <  1:79875/140‑1 (MQ=255)
                                                   cGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtgctg                                                                                        <  2:178407/140‑1 (MQ=255)
                                                   cGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtgctg                                                                                        <  1:824273/140‑1 (MQ=255)
                                                    gTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCCACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtgctgc                                                                                       >  1:1049375/1‑140 (MQ=255)
                                                     tctcACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCg                                                                                      >  1:24109/1‑140 (MQ=255)
                                                         aCTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTaaa                                                                                  >  1:538084/1‑140 (MQ=255)
                                                           tGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAg                                                                                <  1:697562/140‑1 (MQ=255)
                                                           tGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAg                                                                                >  2:312760/1‑140 (MQ=255)
                                                             gcgcGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCa                                                                              <  1:585024/140‑1 (MQ=255)
                                                             gcgcGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCa                                                                              >  2:668111/1‑140 (MQ=255)
                                                              cgcgTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAg                                                                             >  1:470872/1‑140 (MQ=255)
                                                               gcgTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGt                                                                            >  2:312855/1‑140 (MQ=255)
                                                                cgTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGtt                                                                           >  2:723346/1‑140 (MQ=255)
                                                                      ccACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGAAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAgg                                                                     <  2:1012024/140‑1 (MQ=255)
                                                                       cgaaGCGTATTGCAGCAGGTGCACAAGCTGGTTGCACAGATGGATCAGCAGGAAGGGAAAAATGTTTTCCCCTGGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGga                                                                    <  2:102209/137‑1 (MQ=255)
                                                                            cGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGga                                                                    >  2:225284/1‑135 (MQ=255)
                                                                            cGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGga                                                                    <  1:225284/135‑1 (MQ=255)
                                                                               aTTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTgg                                                            >  1:789010/1‑140 (MQ=255)
                                                                               aTTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTgg                                                            >  1:733393/1‑140 (MQ=255)
                                                                               aTTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTgg                                                            >  1:990126/1‑140 (MQ=255)
                                                                                 tGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGga                                                                    <  2:779585/130‑1 (MQ=255)
                                                                                 tGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGga                                                                    >  1:779585/1‑130 (MQ=255)
                                                                                 tGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGaa                                                          >  1:102209/1‑140 (MQ=255)
                                                                                   cagcagGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGTACCTGGaaaa                                                        >  2:712795/1‑140 (MQ=255)
                                                                                    agcagGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAAc                                                       <  1:312855/140‑1 (MQ=255)
                                                                                      cagGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCGAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAg                                                     >  2:1018156/1‑140 (MQ=255)
                                                                                      cagGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAg                                                     >  1:439712/1‑140 (MQ=255)
                                                                                        ggTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGgaga                                                                  <  1:451192/125‑1 (MQ=255)
                                                                                        ggTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGgaga                                                                  >  2:451192/1‑125 (MQ=255)
                                                                                        ggTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAgcg                                                   >  1:1093402/1‑140 (MQ=255)
                                                                                           gggACAGTTGGTTGCACAGATGAAACACAAGGCGGGGCAAAATCATTTAACTGCGACCGACAGTCGCGAGATCTTGTTTATGCAATTACTGATCTTGCTGCGTAAAAGCAGTTTGAAGGAGCACATGGAACACAGCGCAt                                                <  2:1093402/138‑1 (MQ=255)
                                                                                           gCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCAt                                                <  1:790796/140‑1 (MQ=255)
                                                                                             gACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCa                                              <  1:723346/140‑1 (MQ=255)
                                                                                                  ttGGTTGCACAGATGGAACAGAAGGAAGGGAAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGtct                                         <  2:538084/139‑1 (MQ=255)
                                                                                                     gTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCaa                                      >  2:535624/1‑140 (MQ=255)
                                                                                                       tGCACAGATGGAACAGCAGGAAGGGGAAAATGTTTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACt                                    <  1:1153099/140‑1 (MQ=255)
                                                                                                       tGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACt                                    <  1:461434/140‑1 (MQ=255)
                                                                                                              aTGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTct                             >  1:124271/1‑140 (MQ=255)
                                                                                                                 gAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctggc                          <  2:7300/140‑1 (MQ=255)
                                                                                                                 gAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctggc                          <  2:841353/140‑1 (MQ=255)
                                                                                                                 aaaCAGCAGGAGGGGAAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctggc                          <  1:577350/139‑1 (MQ=255)
                                                                                                                      gcagGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctggcctggc                     <  2:24109/140‑1 (MQ=255)
                                                                                                                      gcagGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctggcctggc                     <  1:535624/140‑1 (MQ=255)
                                                                                                                      gcagGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctggcctggc                     >  2:707374/1‑140 (MQ=255)
                                                                                                                       cagGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggct                    >  2:49248/1‑140 (MQ=255)
                                                                                                                        agGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggctg                   >  2:353772/1‑140 (MQ=255)
                                                                                                                         ggAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggctgg                  >  2:367020/1‑140 (MQ=255)
                                                                                                                         ggAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggctgg                  >  2:886325/1‑140 (MQ=255)
                                                                                                                               ggAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGAcc            >  2:494598/1‑140 (MQ=255)
                                                                                                                                gAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCa           >  1:1110677/1‑140 (MQ=255)
                                                                                                                                 aaaaTGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCAt          >  1:772413/1‑140 (MQ=255)
                                                                                                                                    aTGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCa           >  2:702360/1‑136 (MQ=255)
                                                                                                                                    aTGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCa           <  1:702360/136‑1 (MQ=255)
                                                                                                                                         ttACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGa  >  1:820634/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                         
TTCTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGA  >  NC_000913/4098016‑4098292

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: ‑