Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913377,4020GA21.2% 126.9 / 29.1 67intergenic (+91/+133)yaiL/frmBDUF2058 domain‑containing protein YaiL/S‑formylglutathione hydrolase FrmB
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base A (10/4);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTT  >  NC_000913/377153‑377642
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    
aGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGc                                                                                                                                                                                                                                                     >  1:299968/1‑251 (MQ=255)
  cGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGt                                                                                                                                                                                                                                                   >  1:314228/1‑251 (MQ=255)
   gAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGTTAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGtt                                                                                                                                                                                                                                                  >  2:69646/1‑251 (MQ=255)
       gTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCAcg                                                                                                                                                                                                                                              >  2:39720/1‑251 (MQ=255)
       gTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTAcg                                                                                                                                                                                                                                              >  2:107227/1‑251 (MQ=255)
       gTGAATACGCCATTATTCCCGACAGCGTTGCCGATAAAATTGATCAGCGCGATGCCAGCAGCATAGCGTTGCACAGCGCGCTCAGCGCAGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTAAAAGTGCCTGATGATTTGATGTGGGAAACATAGATAAGTCGTGCGCAGATGCATGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAGGTGTAGACCGGATAAGACGtttagg                                                                                                                                                                                                                                              >  2:107229/1‑247 (MQ=255)
          aaTACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCAcgccg                                                                                                                                                                                                                                           >  2:106014/1‑251 (MQ=255)
              cGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATc                                                                                                                                                                                                                                       <  1:296270/251‑1 (MQ=255)
               gCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGAGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGCTTTGATGTGGCAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGAAAAGACGTTTACGCCGCATcc                                                                                                                                                                                                                                      >  2:346978/1‑251 (MQ=255)
                ccATTATCCCCGCCAGCGTTGCCGATAAAATTTCTCAGCGCGATGCCAGCAGTATAGTTTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTAGGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                                                                                                                     <  1:72600/251‑1 (MQ=255)
                   ttattCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGc                                                                                                                                                                                                                                   >  1:7896/1‑250 (MQ=255)
                           gCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcg                                                                                                                                                                                                                          >  2:293009/1‑251 (MQ=255)
                              aGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTTTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCa                                                                                                                                                                                                                      <  2:63945/251‑1 (MQ=255)
                                cGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATg                                                                                                                                                                                                                    >  1:14411/1‑251 (MQ=255)
                                    gCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcgctgatgcct                                                                                                                                                                                                                 >  1:235619/1‑243 (MQ=255)
                                            aaTTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACg                                                                                                                                                                                                        <  1:201376/251‑1 (MQ=255)
                                             aTTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACg                                                                                                                                                                                                        <  1:167524/251‑1 (MQ=39)
                                                gCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCttt                                                                                                                                                                                                    >  1:146842/1‑251 (MQ=255)
                                                 cTCAGCGCGATGCCAGCAGTATAGTGTTGCCCAGCGCGCTCAGCGCCGAAGCGCAGGATGAAGACGACCCGTATGCCGCCTTCAAAGTGCCTGATGATTTGATGTGGTAAACCTAGATACGTCGCGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGTAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcgttgctgcctgatgcgtcgctg                                                                                                                                                                                                     >  2:310302/1‑229 (MQ=255)
                                                  tCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGACGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGc                                                                                                                                                                                                                   >  1:356447/1‑234 (MQ=255)
                                                  tCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTGATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCAGCCGGCAGTCGTGCGCCATGc                                                                                                                                                                                                                   <  2:356447/234‑1 (MQ=255)
                                                            gCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATc                                                                                                                                                                                        <  1:219985/251‑1 (MQ=255)
                                                            gCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTAATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATc                                                                                                                                                                                        <  2:244637/251‑1 (MQ=255)
                                                                   gTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTa                                                                                                                                                                                 <  2:321941/251‑1 (MQ=37)
                                                                         ggtggCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCtt                                                                                                                                                                        <  2:14411/247‑1 (MQ=255)
                                                                          gtTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTa                                                                                                                                                                       <  1:92903/251‑1 (MQ=255)
                                                                            tGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACg                                                                                                                                                                     >  1:1430/1‑251 (MQ=255)
                                                                                   cgcgcTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTAGGCCGACTTCAAAGTGCCTGATGATTTGATGGGGTAAACATAGATAAGTCGTGCGCAGATTCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTATACCGGAGAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTAcaaatctgtacgcga                                                                                                                                                                  >  1:2297/1‑236 (MQ=21)
                                                                                             gcCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGCCGCTTTGCGCGCCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCgg                                                                                                                                                    >  2:148301/1‑251 (MQ=255)
                                                                                             gcCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCgg                                                                                                                                                    <  1:38450/251‑1 (MQ=255)
                                                                                             gcCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTAcaaatctgtacgcgaaccgtaggccg                                                                                                                                                       <  2:162140/251‑26 (MQ=18)
                                                                                                                 cgacCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                <  1:313345/251‑1 (MQ=255)
                                                                                                                   ccccGTATGCCGGCTTCAAAGTGCCGGATGATTGGATGTGGTAAACATGTAAAAGTCGGGCGCAGACCCCTGACGAGACGCTTTGCGCTTCTTACCAGCCTGCCGTTTTATGGGTCAAGTGTGGACCGGATAAGCCTTTCACGCCCCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTGGCGCGTCTGATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                                              <  2:146063/250‑1 (MQ=255)
                                                                                                                         atgccgccttgaatgggCCTTATGATTGGATGGGGTAAACATAGAAAAGTCGTGCGCAGATGCGTGAGCGGACGCTATGCGATTCTGATCAGGCCTAGGGTGTATGGGCGAAGGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCGGTCGGGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGGGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAg                                                                                                                         <  2:299968/234‑1 (MQ=21)
                                                                                                                           gCCGACTTCAAAGTGCCTGATGATTTGAGGTGGAAAACATAGATAAGTCGTGCGCCGATGCCGGGTGGGACGCTATGCGCGTCTTATAGGGCCAACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTc                                                                                                                       <  2:146842/250‑1 (MQ=37)
                                                                                                                                 ttCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACG‑‑CGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCgacctgaaaattc                                                                                                                 >  1:302673/1‑239 (MQ=2)
                                                                                                                                  tcaaagggccgtagtatTTGATGTGGTAAAAATAGAGAAGTCGTGCTCAGATGCGTGATGGGGCGCTAGGCGTGTTTTAGCGGCGACCCGTTGTATGGGTGATGTGTAGGCGGGTTATGGCTTTCACGCCGCTTCCGGCAGGCGTGCGCCATGCCTGATGCGACTCTTTGCGCGTCTTATCAGGCCTACG‑‑GTTTGCGGGCGATGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCtat                                                                                                                <  2:318879/235‑1 (MQ=12)
                                                                                                                                                     gaggtgGTAAACATAGAAAAGTCGTGCGCAGAAGCCGGATGCGAGGCGATGCGCGTGTTATAGGGCCTACGGTTTCTGGGCGAAGGGTAGCCGGGAAAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGttt                                                                                            <  2:90703/248‑1 (MQ=25)
                                                                                                                                                         tgGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGAGACGCTTTTCGCGTCTTATCAGGCCTAAG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCCCGCCGCATACAGCAGGCGTGCTATTATCAACGCATATTCAGtttcttt                                                                                        >  2:126145/1‑247 (MQ=37)
                                                                                                                                                                  tAGATAAGTCGTGCGCAGATGCCTGATGCGACGCGATGCGCGTCTTATCAGGCTTACGGTTTAGTGGCGAGGTGTAGACCGGATAAGGCGGTCACGCCGCAGCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCGTCCGGCAGTCGTGCTATTATCAACGCat                                                                                                     <  1:313568/229‑1 (MQ=255)
                                                                                                                                                                  tAGATAAGTCGTGCGAAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGAGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCAGTCGTGCTATTATCAACGCat                                                                                                     >  2:313568/1‑229 (MQ=255)
                                                                                                                                                                      tAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGACGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCa                                                                            >  2:39348/1‑250 (MQ=255)
                                                                                                                                                                      tAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCa                                                                            >  1:360675/1‑250 (MQ=255)
                                                                                                                                                                         gTCGTGCGCAGATGCCTGATGCGGCGCTAGGAGCGTCTTATCAGGCCAACGGTTTATGGGCGAAGTGTAGACGGGATAAGGCGTTCACGCCTCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAatat                                                                        <  2:320303/251‑1 (MQ=255)
                                                                                                                                                                                gcAGATGCCTGATGAGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCAACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGCGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGccattttatcgcc                                                                 >  2:52447/1‑239 (MQ=255)
                                                                                                                                                                                   gATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAAt                                                              <  1:106014/251‑1 (MQ=255)
                                                                                                                                                                                        cTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGTCGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTAGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCGATaaaa                                                          >  2:74713/1‑250 (MQ=255)
                                                                                                                                                                                           aTGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTgg                                                      >  1:250905/1‑251 (MQ=255)
                                                                                                                                                                                            ggCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGga                                                     <  2:1430/250‑1 (MQ=255)
                                                                                                                                                                                               gACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGAc                                                  >  1:188847/1‑251 (MQ=255)
                                                                                                                                                                                                 aCGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACa                                                 <  1:74592/251‑1 (MQ=255)
                                                                                                                                                                                                  ggtgaggggcgtgtgagcgggCCTAGGGTTTAGGGGGGAAGGGTAGACGGGATAAGGGGTTCACGCCGCATCCGGCGGTGGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGTGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACac                                                <  2:106780/232‑2 (MQ=37)
                                                                                                                                                                                                    cTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACaaa                                               >  1:317852/1‑250 (MQ=255)
                                                                                                                                                                                                       tGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCACGCCAATAAAACTGGAGACAAaata                                            >  1:212346/1‑250 (MQ=255)
                                                                                                                                                                                                                    aGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCAt                              <  1:292889/251‑1 (MQ=255)
                                                                                                                                                                                                                        cTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAcc                          >  1:29583/1‑251 (MQ=255)
                                                                                                                                                                                                                           agggTTTATGGGCAAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCt                        <  2:263111/249‑1 (MQ=255)
                                                                                                                                                                                                                            cGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTc                       <  1:223884/251‑1 (MQ=255)
                                                                                                                                                                                                                              gTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTTCG‑‑GCTTACGGGCGAAGTGTCGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTGGTGCTATTATCAACGCATATTCAGCTTATTGGCGTGGGAGGAAAGATGCTCGCCAATAAAACTGGAGACAAAATAATAACTGTGATGATaaacatcg                      >  2:241789/1‑244 (MQ=38)
                                                                                                                                                                                                                              gTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGt                     >  1:50468/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGa                   >  1:2851/1‑250 (MQ=255)
                                                                                                                                                                                                                                  tATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGAt                  >  1:143866/1‑251 (MQ=255)
                                                                                                                                                                                                                                             tgtAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAAc       >  1:308340/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gtAGACCGGCTAAGGCGTTCACGCCGCATCCGGCGGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACg      <  1:10546/251‑1 (MQ=255)
                                                                                                                                                                                                                                                aGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATTGCTGTGATCATAACCCTCGTGATAACGGATTAACGt     <  2:314228/250‑1 (MQ=255)
                                                                                                                                                                                                                                                aGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTc    >  1:151692/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  aCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCtt  <  2:360675/251‑1 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    
AGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTT  >  NC_000913/377153‑377642

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