Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913417,5330TG48.1% 35.0 / 22.8 54V131G (GTG→GGG) phoBresponse regulator in two‑component regulatory system with PhoR
Reads supporting (aligned to +/- strand):  ref base T (17/11);  new base G (26/0);  total (43/11)
Fisher's exact test for biased strand distribution p-value = 3.05e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGA‑GCCGCTGGAGATGGGGCCGACA‑G‑AATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTAC‑AGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTG‑GAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAA‑GCACTGGAGCCCGGCGG  >  NC_000913/417303‑417770
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
gATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGAAGAGGGGGTTggga                                                                                                                                                                                                                                             >  1:301573/1‑239 (MQ=255)
  tGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGCTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGAAGAGGGGAGTGAGGTg                                                                                                                                                                                                                                           >  2:1443382/1‑239 (MQ=255)
  tGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGAt                                                                                                                                                                                                                                            >  1:753368/1‑238 (MQ=255)
       ccTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGAAGAGGGGGTTGGGATGGAggg                                                                                                                                                                                                                                      >  1:592133/1‑239 (MQ=255)
          ggcggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGGtt                                                                                                                                                                                                                                   >  1:1048203/1‑239 (MQ=255)
           gcggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGGGATTGAAATGCAGGGATTa                                                                                                                                                                                                                                  >  2:1701846/1‑239 (MQ=255)
           gcggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGAAGAGGGGGTTGAGATGCAGGGATTa                                                                                                                                                                                                                                  >  1:2034692/1‑239 (MQ=255)
            cggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGGGGGGATTGGGAAGCAAGGAGTaa                                                                                                                                                                                                                                 >  2:1596357/1‑239 (MQ=255)
            cggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGAAGAGGGGATTGAGATGCAGGGATTaa                                                                                                                                                                                                                                 >  1:2339784/1‑239 (MQ=255)
                            ttCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGTAGGAGCTGGTGGCGCGAATCAACGCGGTAATGCGCCGCATTTCGCCAATGGCGGCGGCAGGGGGGATTGAAATGCCAGGATTAGTACTCGCACCGACAt                                                                                                                                                                                                                  >  2:337331/1‑238 (MQ=255)
                               atcaAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGACGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCACTGGCGGTGGAAGAGGTGATTTGGATGCCGGGGTTAAATTTCGACCCGCCAtcgc                                                                                                                                                                                                               >  2:518939/1‑236 (MQ=255)
                                      acCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGGTTAAGTCTCGACCCGACATCTCACCGAGTg                                                                                                                                                                                                       >  1:1884675/1‑239 (MQ=255)
                                                 gAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAACGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGGGGTTGAGAGGCAGGGTGTAAGTCTCGACCCGACGTCTCTCCGAGTgcaggtggggg                                                                                                                                                                                             >  2:862221/1‑228 (MQ=255)
                                                  aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTTTCGACCCGGCCTCTCACCGAGGGATGGCGGGGGac                                                                                                                                                                                           >  2:200544/1‑238 (MQ=255)
                                                  aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGaa                                                                                                                                                                                           >  1:2155827/1‑239 (MQ=255)
                                                                    ttCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGGAGAGGGGATTGAGATGCAGGGGTTAAGTCTCGACCCCACATCTCACCGAGGGATGGCGGGGGAAGA‑GGCCCCGGTGGTgggg                                                                                                                                                                        >  2:1070245/1‑239 (MQ=255)
                                                                         gtggtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGGGGAGA‑GGCGCCGGGGATGGGGGCggc                                                                                                                                                                   >  1:1152748/1‑237 (MQ=255)
                                                                            gtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCCGGGGTTAAGTCTCGACCCGACATCCCACCGAGTGATGGGGGGGGGAGAGGCCCCTGGGGAGGGGGcgcccaa                                                                                                                                                                 >  2:1080474/1‑233 (MQ=255)
                                                                            gtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGGTTAAGGCTCGACCCGACATCTCACCGAGGGATGGCGGGGGAAGA‑GGCCCTGGGGAGGGGGGCGGCA‑G‑a                                                                                                                                                              >  1:574239/1‑239 (MQ=255)
                                                                            gtgATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGAAATCTCACCGAGTGATGGCGGGGGGAAG‑GCCGCCGGAGATGGGGGCGACA‑G‑a                                                                                                                                                              >  1:1028571/1‑239 (MQ=255)
                                                                               aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGACGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTTGGATGCAGGGATTAAGt                                                                                                                                                                                                                               >  1:509155/1‑174 (MQ=255)
                                                                               aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGGGATGGCGGGGGAGGA‑GCCGCTGGAGAGGGGGGCGACA‑C‑aagt                                                                                                                                                           >  1:1221543/1‑237 (MQ=255)
                                                                               aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGAACCGACATCTCACCGAGTGATGGCGGGGGGAGA‑GCCGCCGGGGGTGGGGGCCACA‑G‑aaat                                                                                                                                                           >  1:2268133/1‑237 (MQ=255)
                                                                               aTGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGGGAAGA‑GCCGCCGGAGATGGGGCCCGCA‑G‑aaat                                                                                                                                                           >  1:1874050/1‑237 (MQ=255)
                                                                                                               cgTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGGGGTTGAGATGCAGGGATTAAGTCTCCAACCGAGATCTCCCCGAGTGATGGCGGGCGGAGA‑GACGCCTGTGAGGGGGCCGGCG‑G‑AGATTAAAATGCGTCACTTTTTTTTCATGCAt                                                                                                                               >  2:1640533/1‑235 (MQ=255)
                                                                                                                                  ggCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCa                                                                                                                                                                                                                                         >  1:179022/1‑113 (MQ=255)
                                                                                                                                  ggCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGa                                                                                                                                                                                                                     >  1:795030/1‑133 (MQ=255)
                                                                                                                                  ggCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGGGAAAA‑GCCCCTGGGGAt                                                                                                                                                                            >  1:1491363/1‑173 (MQ=255)
                                                                                                                                    cgcgGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGGAGA‑GCCGCCGggga                                                                                                                                                                             >  1:1712478/1‑170 (MQ=255)
                                                                                                                                                   cACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGGAGA‑GCCGCCGGAGAt                                                                                                                                                                            >  1:1886044/1‑156 (MQ=255)
                                                                                                                                                     ccAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCa                                                                                                                                                                                                                                         >  1:891390/1‑94 (MQ=255)
                                                                                                                                                       aaGCCGTTTTCGCCGAAGTAGCGGGTGGCGCGAATCAACGCGGTAATGCGCCTTATTTCGCCAATGGCGGTGGAAGAGGTGATGGAGATGCAGGGATGAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGA‑GCCGCTGGAGATGGGGCCGACA‑G‑AATTTAAACTGCTGCACTTTTTTATGGCGCATCCTGAGCGCGTGTAC‑AGCCGCGAGCAGCTGTTAAACCACGTCTg                                                                                  <  2:570461/239‑1 (MQ=255)
                                                                                                                                                        aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGGGATTGAGATGCAGGGGTTAAGTCTCGACCCGACATCTCACCGGGTGATGGCGGGGGGAGA‑GGCGCTGGGGATGGGGGCGGCAGG‑AATTTAAA‑TGGTGCACTTTTTTTTTATGCAACCCGAGGGGGTGTAC‑AGCCGCGGGCCGGTGTTAAACACCGTCttg                                                                                 >  2:582316/1‑237 (MQ=255)
                                                                                                                                                        aGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGAc                                                                                                                                                                                                                         >  1:644978/1‑107 (MQ=255)
                                                                                                                                                                                     cgAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGGGGAGA‑AGCGCCGGGGGTGGGGGCGAGA‑G‑AAATTTAACTTGTGCACTTTTTTTTTAAGCCTCCTGGGCGCGGGGACAACCCGCG‑GGAGGCGTTTAAACACGCCTCGGGGACTAACCGGTATTTGTGAAAACCGcc                                                    >  2:130246/1‑238 (MQ=255)
                                                                                                                                                                                     cgAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGGAGAGGGGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGGGGAGG‑AGCGCCGGGGGGGGGGGCGGCA‑GAAATTTAAATggctgc                                                                                                                                                >  2:132639/1‑149 (MQ=255)
                                                                                                                                                                                       aaTCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGGTTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGGAGA‑GGCGCTGGAGATGGGGCCCGCA‑G‑AATTTAAAATGCTGCCCTTTTTTTTGAAGCATCCt                                                                                                                            >  2:514036/1‑166 (MQ=255)
                                                                                                                                                                                                         ccGTATTTCGCCAATGGCGGTGGGAGAGGGGATTGAGATGCAGGGGTTAAg                                                                                                                                                                                                                                >  1:31216/1‑51 (MQ=255)
                                                                                                                                                                                                         ccGTATTTCGCCAATGGCGGTGGAAGAGGGGATTGAGATGCAGGGATTAAg                                                                                                                                                                                                                                >  1:679045/1‑51 (MQ=255)
                                                                                                                                                                                                            tATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCa                                                                                                                                                                                                              >  1:1564640/1‑66 (MQ=255)
                                                                                                                                                                                                            tATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGt                                                                                                                                                                                                        >  1:191123/1‑72 (MQ=255)
                                                                                                                                                                                                                  gCCAATGGCGGTGGAAGAGGTGa                                                                                                                                                                                                                                                   <  2:891390/23‑1 (MQ=255)
                                                                                                                                                                                                                  gCCAATGGCGGTGGAAGAGGTGa                                                                                                                                                                                                                                                   <  2:179022/23‑1 (MQ=255)
                                                                                                                                                                                                                          cGGTGGAAGAGGTGATTGAGATg                                                                                                                                                                                                                                           <  2:1462631/23‑1 (MQ=255)
                                                                                                                                                                                                                           ggTGGAAGAGGTGATTGAGATGc                                                                                                                                                                                                                                          <  2:31216/23‑1 (MQ=255)
                                                                                                                                                                                                                           ggTGGAAGAGGTGATTGAGATGc                                                                                                                                                                                                                                          <  2:679045/23‑1 (MQ=255)
                                                                                                                                                                                                                            gTGGAAGAGGTGATTGAGATGCa                                                                                                                                                                                                                                         <  2:592133/23‑1 (MQ=255)
                                                                                                                                                                                                                            gTGGAAGAGGTGATTGAGATGCa                                                                                                                                                                                                                                         <  2:753368/23‑1 (MQ=255)
                                                                                                                                                                                                                            gTGGAAGAGGTGATTGAGATGCa                                                                                                                                                                                                                                         <  2:509155/23‑1 (MQ=255)
                                                                                                                                                                                                                            gTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGGGAAGA‑GGCCCCGGGGATGGGGGCGGCA‑GAAATTTAAAC‑GTTGCCCCTTTTTTTGAAGCAACCt                                                                                                                            >  2:2298935/1‑129 (MQ=255)
                                                                                                                                                                                                                                  gagGTGATTGAGATGCAGGGAtt                                                                                                                                                                                                                                   <  2:644978/23‑1 (MQ=255)
                                                                                                                                                                                                                                     gTGATTGAGATGCAGGGATTAAGTCTCGTCCCGACTTCTCACCGAGTGATTGCGGGCGAATA‑GCCGCTGGAGATGGGGCCGAc                                                                                                                                                                   <  2:806241/83‑1 (MQ=255)
                                                                                                                                                                                                                                     gTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGA‑GCCGCTGGGGATGGGGCCGACA‑G‑AATTTAAACTGCTGCCCTTTTTTATTACGCATCCTGGGCGCGGGTAC‑AGCCGCGAGCAGCTGTTAAACCACGGCTGGGGGACTAACGTGTATTTG‑GGAGACCCCACGGTCGATGTCCACATTCGTCGCCCGCGTAAA‑GCCATGGAGCCcggccg  >  2:801240/1‑237 (MQ=255)
                                                                                                                                                                                                                                     gTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGA‑GACGCTGGAGAAGGGCGCGGCA‑G‑GATTTAAACTGGCGCCCTTTTTTATGACGGCTCCTGAGCGCGTGTAC‑AGCCGCGAGCAGCCGTTAAAACACGGCTGGGGGACTAACGTGTATGTG‑GAAGACCGCACGGTCGATGTCCACATTCGTACCCCGCGGAAATGCATTGGAcccgcggg   >  2:779455/1‑234 (MQ=255)
                                                                                                                                                                                                                                      tGATTGAGATGCAGGGATTAAGt                                                                                                                                                                                                                               <  2:1048203/23‑1 (MQ=255)
                                                                                                                                                                                                                                      tGATTGAGATGCAGGGATTAAGt                                                                                                                                                                                                                               <  2:795030/23‑1 (MQ=255)
                                                                                                                                                                                                                                      tGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGA‑GCCGCTGGAGATGGGGCCGACA‑G‑AATTTAAACTGCTGCACTTTTTTATGAc                                                                                                                                   <  1:999610/112‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
GATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGAAGAGGTGATTGAGATGCAGGGATTAAGTCTCGACCCGACATCTCACCGAGTGATGGCGGGCGAAGA‑GCCGCTGGAGATGGGGCCGACA‑G‑AATTTAAACTGCTGCACTTTTTTATGACGCATCCTGAGCGCGTGTAC‑AGCCGCGAGCAGCTGTTAAACCACGTCTGGGGAACTAACGTGTATGTG‑GAAGACCGCACGGTCGATGTCCACATTCGTCGCCTGCGTAAA‑GCACTGGAGCCCGGCGG  >  NC_000913/417303‑417770

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