Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913417,5330TG44.0% 60.9 / 26.7 77V131G (GTG→GGG) phoBresponse regulator in two‑component regulatory system with PhoR
Reads supporting (aligned to +/- strand):  ref base T (23/19);  new base G (33/0);  total (58/19)
Fisher's exact test for biased strand distribution p-value = 1.71e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAG‑AG‑CCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGC  >  NC_000913/417309‑417748
                                                                                                                                                                                                                                |                                                                                                                                                                                                                          
aCCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAgg                                                                                                                                                                                                              >  1:2883171/1‑239 (MQ=255)
aCCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATAACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGAAGAGGGGATTGAGATGCAgg                                                                                                                                                                                                              >  1:1074711/1‑239 (MQ=255)
    ggcggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGTGGTTGAGATGCAGGGAtt                                                                                                                                                                                                          >  1:1241692/1‑239 (MQ=255)
    ggcggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGGtt                                                                                                                                                                                                          >  1:3410626/1‑239 (MQ=255)
    ggcggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGGGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGCATTTCGCCAATGGCGGTGGAAGAGGGGATTtaga                                                                                                                                                                                                                    >  1:1669340/1‑229 (MQ=255)
    cgcggcTCCGGCATCCCGTTCATCAAACACCTCAAGAGCGCGTCGATGTCCCGGGCTACTCCAGCGGCGATGTTGACCGCCGGAGGGGCAGAAGAAGATCCCGTGCGCGGCCTTGCAACCGGCGCGGATGCCTATATCACCAAGCCGTTTTCGCCTTCGGAGCTGGGGGCGCGAATCAAAGCGGAAACGCGCCGTATGTCGCCGACGGCGGGGGGAAGGGTGGTGTGGGTGCAGGGAtt                                                                                                                                                                                                          >  2:817691/2‑239 (MQ=255)
     gcggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCATTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGCGGAGAGGGGGATGGGATGCGGGGGTTa                                                                                                                                                                                                         >  2:2025113/1‑239 (MQ=255)
      cggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGGGAAGCAGGGGTTaa                                                                                                                                                                                                        >  1:1597452/1‑239 (MQ=255)
         cTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGCGGAGATTGAGATGCAGGGATTAAg                                                                                                                                                                                                       >  2:3533257/1‑237 (MQ=255)
         cTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGCGGAAGAGGGGATTGAGAGGCAGGGATTAAg                                                                                                                                                                                                       >  1:2920515/1‑237 (MQ=255)
          tCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGAGGGGGGGGGTGGGAAGCAgggggttaggc                                                                                                                                                                                                     >  2:1437293/1‑230 (MQ=255)
                   cAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGCGGTGATTGAGATGCAGGGGTTAAG‑TCTCGACCCGa                                                                                                                                                                                           >  1:3085294/1‑238 (MQ=255)
                   cAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGAc                                                                                                                                                                                          >  1:2624925/1‑239 (MQ=255)
                      ttCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGAc                                                                                                                                                                                               >  1:253864/1‑231 (MQ=255)
                      ttCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACa                                                                                                                                                                                         >  1:2743876/1‑237 (MQ=255)
                                acCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGa                                                                                                                                                                                                >  1:266876/1‑220 (MQ=255)
                                  cTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGAt                                                                                                                                                                           >  1:961868/1‑239 (MQ=255)
                                  cTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGAt                                                                                                                                                                           >  1:2118284/1‑239 (MQ=255)
                                     aaGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGGTGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATggc                                                                                                                                                                        >  1:1525460/1‑239 (MQ=255)
                                          cgAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGAGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGAGGCAGGGATTAAG‑GCCCGACCCGCCAATCCACCGAGTGATggcggggg                                                                                                                                                                   >  1:1825131/1‑237 (MQ=255)
                                            aGTCGATGACCCGGGCTATTTCAGTGGTGATGTTGACCGCCAGAGGGGAAGGAGAAGACCGCGTGCGTGTCCTTGTAACCGGCGCGGATGACTATCTCCCCACGCCCGTTTCGGCGAAGGAGCTGGTGGCGCGCATCAAAGCGACAATGCGCCGTATTTCGCCAATGGCGGTGGTAGAGGGTATTGAGATGCAGGGGATAAG‑GATCGACGCGAGATCCTACCGCGTGATggcgggcgcg                                                                                                                                                                 >  2:2993446/1‑237 (MQ=255)
                                            aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATggcgggggga                                                                                                                                                                 >  1:2355128/1‑235 (MQ=255)
                                            aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAAGGGGGATTGAGATGCAGGGATTAAGTTCCCGA‑CCGGCATCTCACCGAGGGagagggggggga                                                                                                                                                                 >  2:2679792/1‑228 (MQ=255)
                                            aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATTTCACCGAGCGAGGGCGGGCGaa                                                                                                                                                                 >  1:2694539/1‑239 (MQ=255)
                                            aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCTCCGAGGGATggcgggggga                                                                                                                                                                 >  2:2115860/1‑235 (MQ=255)
                                            aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGCGaa                                                                                                                                                                 >  1:3469067/1‑239 (MQ=255)
                                            aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑GCTCGACCCGACATCTCACCGAGGGATggcgggggga                                                                                                                                                                 >  1:1442925/1‑235 (MQ=255)
                                            aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGGAGAGGGGCTTGAGGGGCCGGGGTTAAG‑GCTCGACCCGACATCTCACCGAGGGGTtgcgggggga                                                                                                                                                                 >  1:1587187/1‑235 (MQ=255)
                                            aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGGAGAGGGGATTGGGATGCCGGGGTTAAG‑TCTCGAACCCACATCCCACCGGGGGATggcgggggga                                                                                                                                                                 >  2:324830/1‑235 (MQ=255)
                                                             aaaccAGTGGGGATGTGGACCGCAAGAGAGAAAGAAGAAGAGCGCGTGCGCGTCCTTGAAGCCGGGGCGAAGGACTATATAACCAAGACATTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAAGGCGAAGAATTTCGCCAATGGCGGTGGAAGAGGTGATCGAGATGCAGGCAGTAAA‑GCTCGACCCGACATCTCACCGAATGATGGCGGGCGAAG‑AG‑GCGCAGGAGATggg                                                                                                                                              <  2:1233017/236‑1 (MQ=255)
                                                                    tggtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑t                                                                                                                                                                                                     >  1:3200749/1‑179 (MQ=255)
                                                                    tggtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGGGATGGCGGGGGGAG‑AG‑CCCCCGGGGGTGGGGCCgccc                                                                                                                                       >  1:3489163/1‑236 (MQ=255)
                                                                    tggtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTaa                                                                                                                                                                                                        >  1:256338/1‑177 (MQ=255)
                                                                         aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAg                                                                                                                                                                                                               >  1:1146534/1‑165 (MQ=255)
                                                                                ccGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTaa                                                                                                                                                                                                        >  1:3345552/1‑165 (MQ=255)
                                                                                              gaagaagaTCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGACCCGAAATCTCACCGAGTGATGGCGGGGGGAG‑AGCCCCCGGGGGGTGGGGCCG‑CAGAAATTAAACTGCGGGCCtttttttt                                                                                                             >  1:2065739/1‑237 (MQ=255)
                                                                                                         cgTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCa                                                                                                                                                                                    >  1:1094632/1‑159 (MQ=255)
                                                                                                         cgTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGGGAAA‑AG‑CCCCCGGGGATGGGGGCGCCAGAAATTTAAATGCTGCCCTTTTTTTTGAAGCATc                                                                                                     >  1:3352664/1‑236 (MQ=255)
                                                                                                         cgTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGCGGAG‑AG‑CCGCCGGGGATGGGGCCCGCAGAATTTTAACTGCTGCACTTTTTTTTGACGCAt                                                                                                      >  1:829771/1‑235 (MQ=255)
                                                                                                          gTGCGCGGCATTGAAGCCGGGGCGGATGACTATGTCACCAAGCCGTTTGCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGACGGTGGAGGAGGTGATTGAGATGCAGGTATTAAG‑ACGCGACCCGACAGCTCGCCGAGTGATGGCGGGCGAAG‑AG‑ACGCTGGAGGTGGGGCCGACAGAATTTAATCTGCTGCGCTTTTTTATGACGCATccag                                                                                                  <  2:1825131/238‑3 (MQ=255)
                                                                                                                      gAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGGTTGAGATGCAGGGATTAAG‑TCTCGACCCGACCTCTCACCGAGTGATGGCGGGGGGAG‑GG‑CGGCGGGGGAGGGGGGCGGCCGAATTTATACAGTGCCGCTTTTTTTTTTGGCATCCTCGGCGCGGGTGCa                                                                                      >  2:29795/1‑238 (MQ=255)
                                                                                                                                  gaagcatttatCGCCATGGCGTTTTCGCCGGAGAGGCTGGGTGCGTAAATCGACGCAGTGATACGCCGTATTTCGCCCATGGCTGTGGAAAAGATGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAG‑AG‑CCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGCCGCTTCCTGAGCGCGTGTACAGCCGCGagca                                                                            <  2:643308/228‑1 (MQ=255)
                                                                                                                                          tatCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAA‑TCTCGACCCGa                                                                                                                                                                                           >  1:1852689/1‑119 (MQ=255)
                                                                                                                                                  aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑t                                                                                                                                                                                                     >  1:247200/1‑101 (MQ=255)
                                                                                                                                                  aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGCGTGATGGCGGGCGGAG‑GG‑CCGCCTGGGATGGGGCCGACAGAAATTAAAATTCTGCACTTTTTTATGACGCGTCCTGAgc                                                                                               >  2:2182704/1‑201 (MQ=255)
                                                                                                                                                  aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGGGGAG‑AG‑ACGCCGGGGAAGGGGCCGGCAGAATTTTAAATGATGCACTTTTTTTTGACGCATCCTGAGCGCGTGTACAGCCGCGGGCAGCTGTTAAAACACGCCTgg                                                         >  2:2752290/1‑239 (MQ=255)
                                                                                                                                                  aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGGGAAA‑AG‑CCGCt                                                                                                                                                       >  1:1167831/1‑145 (MQ=255)
                                                                                                                                                  aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGGGATGGCGGGCGCAG‑AG‑ACCCCGGGGGTGGGGCCGGCAGAAATTTAAATGCTGCACCTTTTTTTTTAGCCTCCTGCGGGCGCGTTCAACCCCGAGCAGATGTTTAATCACATCttt                                                         >  2:169421/1‑237 (MQ=255)
                                                                                                                                                  aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCa                                                                                                                                                                                    >  1:1306535/1‑118 (MQ=255)
                                                                                                                                                  aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGGGGAGAAG‑CGGCGGGGGGGGGGGCCGAAGGAATTTAAAATCTGGCCCTTTTTTTTGA‑GGATACTGGACGCGGGTGAAGCCGCCGGCAGCCGTTTAACCACgcccgt                                                         >  2:3424866/1‑234 (MQ=255)
                                                                                                                                                  aGCCGTTTTCGCCGAAGGAGCTGGAGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCAcc                                                                                                                                                                                  >  1:886933/1‑120 (MQ=255)
                                                                                                                                                    ccGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGCGCGAAG‑AG‑CCg                                                                                                                                                         >  2:1422617/1‑141 (MQ=255)
                                                                                                                                                    ccGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGAGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCa                                                                                                                                                                                    >  1:931618/1‑116 (MQ=255)
                                                                                                                                                           cgacGAAGGAGGTGGTGGCGTGAATCAAAGCGGTAAGGCGCCGTATTTCGCCGATGGCGGTGGAAGAGGTGATTGAGGTGCAGGGATTAAG‑GTTCGTCCCTACATCTACCCGAGTGATGGCGGGCGAAG‑AG‑TAGCTGGATATTGGGCCGCGTGAATTTAATCTGCGGCACTTTTTTAAGACGCATTCTGACCGCATGTTCAGCCTCAAACAGCTGTTAAACCACGTCTTGGGAACTatc                                                <  2:590921/236‑3 (MQ=255)
                                                                                                                                                              cGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGGTTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGGGGAG‑AG‑CCGCTTGGGAGGGGGGCGCCACAAATTTAACAGCTGCACTTTTTTTTTACGCaac                                                                                                     >  2:2886791/1‑181 (MQ=255)
                                                                                                                                                                                  aTCAAAGTGGTAAGGCGCCTTATTTCGCCAATGCCGGTGGAAGGGGTGATTGAGATTCAGGGATTAAG‑TCTCGACCCTACATCTCACCGAGTGATGGCGGGCGAAG‑AG‑CCTCTGGAGATGGGGCCGACGGAATCTAACCTGCTGCACTTTTTTTCGACGCTCGCTGAGCGCGGGTGCAGCCAAGAGCAGGTGTTAAACCACGTCTGGGGAACTAACGTGTATATGGaaagtcggtcggg                         <  2:1142973/239‑12 (MQ=255)
                                                                                                                                                                                       aGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGt                                                                                                                                                                                                                            <  1:1815154/42‑1 (MQ=255)
                                                                                                                                                                                            tAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTaa                                                                                                                                                                                                        >  1:608489/1‑57 (MQ=255)
                                                                                                                                                                                            tAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTaa                                                                                                                                                                                                        >  1:2387979/1‑57 (MQ=255)
                                                                                                                                                                                              aTGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGGTTGAGGTGCAggg                                                                                                                                                                                                             >  1:3395844/1‑50 (MQ=255)
                                                                                                                                                                                                    cGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGCGGAG‑AG‑GCGCAGGGGATgg                                                                                                                                               >  1:1756951/1‑103 (MQ=255)
                                                                                                                                                                                                     gTATTTCGCCAATGGGGGTGGAAGTGGTGATGTGGATGCAGGGATTAAG‑CCTCGACCCGACATCTCCCCGAGTGATGGCGGGCGAAG‑AG‑CCGCTGGAGATGGGGCCGACAGAATTTAACCTGCTGCACTTTTTTATGCCGCATCCTGGGCGCGGGTACAGCCGCAAGCAGCTGTTAAACCACGACTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGTTGTCCACATTCGTCGc      <  2:2743876/239‑1 (MQ=255)
                                                                                                                                                                                                      tATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGTCATCTCACCGAGTGATGGCGGGGGGAG‑AG‑CCGCCGGGGATgggg                                                                                                                                             >  1:1192191/1‑103 (MQ=255)
                                                                                                                                                                                                         ttCGCCAATGGCGGTGGGAGGGGGGATTGAGATGCAGGGGTTAAG‑TCTCGACCCGACATCTCACCGAGGGATGGCGGGCGGAG‑AG‑CCGCTGGAGCATGGGCCCACAGGAATTTAACCGCGGCACCTTTTTTTTACGCATCCCTGGCGCGGGTAACACCCCGCGCGCATGTTTAAACACGTGTGTGGGACTAACGTGTATGTGGGGGACCGCACGGGCGCAGTCCACAATCGTCGCCTGc  >  2:1988528/1‑239 (MQ=255)
                                                                                                                                                                                                             ccAATGGCGGTGGAAGAGGTGAt                                                                                                                                                                                                                         <  2:1146534/23‑1 (MQ=255)
                                                                                                                                                                                                               aaTGGCGGTGGAAGAGGTGATTg                                                                                                                                                                                                                       <  2:3395844/23‑1 (MQ=255)
                                                                                                                                                                                                                aTGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGGGGAG‑AG‑CCGCTGGAGATGGGGCCGGCAGAATTTTAAATGCTGCCCTTTTTTTTTAAGCCTCCCGAGCGCGGGTAAAggc                                                                                   >  2:3019262/1‑149 (MQ=255)
                                                                                                                                                                                                                    cGGTGGAAGAGGTGATTGAGATg                                                                                                                                                                                                                  <  2:256338/23‑1 (MQ=255)
                                                                                                                                                                                                                    cGGTGGAAGAGGTGATTGAGATg                                                                                                                                                                                                                  <  2:608489/23‑1 (MQ=255)
                                                                                                                                                                                                                    cGGTGGAAGAGGTGATTGAGATg                                                                                                                                                                                                                  <  2:3345552/23‑1 (MQ=255)
                                                                                                                                                                                                                    cGGTGGAAGAGGTGATTGAGATg                                                                                                                                                                                                                  <  2:2387979/23‑1 (MQ=255)
                                                                                                                                                                                                                      gTGGAAGAGGTGATTGAGATGCa                                                                                                                                                                                                                <  2:3200749/23‑1 (MQ=255)
                                                                                                                                                                                                                      gTGGAAGAGGTGATTGAGATGCa                                                                                                                                                                                                                <  2:1597452/23‑1 (MQ=255)
                                                                                                                                                                                                                      gTGGAAGAGGTGATTGAGATGCa                                                                                                                                                                                                                <  2:3444352/23‑1 (MQ=255)
                                                                                                                                                                                                                      gTGGAAGAGGTGATTGAGATGCa                                                                                                                                                                                                                <  2:247200/23‑1 (MQ=255)
                                                                                                                                                                                                                           agagGTGATTGAGATGCAGGGAt                                                                                                                                                                                                           <  2:266876/23‑1 (MQ=255)
                                                                                                                                                                                                                            gagGTGATTGAGATGCAGGGAtt                                                                                                                                                                                                          <  2:253864/23‑1 (MQ=255)
                                                                                                                                                                                                                               gTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAG‑AG‑CCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGGGCGCGTGTACAGCCGCGAGCAGCTGTTAAGCCACGTCTGGGGAACt                                                   <  1:1805413/168‑1 (MQ=255)
                                                                                                                                                                                                                                tGATTGAGATGCAGGGATTAAG‑t                                                                                                                                                                                                     <  2:1852689/23‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                          
ACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAG‑TCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAG‑AG‑CCGCTGGAGATGGGGCCGACAGAATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTACAGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTGGAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGC  >  NC_000913/417309‑417748

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑