Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913377,4020GA40.6% 48.7 / 74.2 65intergenic (+91/+133)yaiL/frmBDUF2058 domain‑containing protein YaiL/S‑formylglutathione hydrolase FrmB
Reads supporting (aligned to +/- strand):  ref base G (19/19);  new base A (14/12);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-02
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATAT  >  NC_000913/377153‑377647
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                            
aGCGAAGGTGAATCCGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAgg                                                                                                                                                                                                                                                              <  1:383228/250‑1 (MQ=255)
 gCGAGGGTGAAGACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACTGGATAAGGCg                                                                                                                                                                                                                                                            <  2:433583/251‑1 (MQ=255)
      gtgtaaaaCGCAATTATTACCGCCAGCGTTGCCGATAAAATTGCTCATCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTTCCTTATGATTTGATGTGGTAAACATAGATAAGTCTTGCGCAGATGCCTGATGCGCCGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCAc                                                                                                                                                                                                                                                      <  1:260896/244‑1 (MQ=255)
      ggTGAATTCGCTATTATTCCCGACAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGATCAGGATGAAGACGACCCGTATGCCGACTTAAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCCGATGCCTGATGCGACTCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGtttac                                                                                                                                                                                                                                                       <  2:330507/251‑4 (MQ=255)
            tACGCAATTATTCCCGCCAGCGTTGCCGATACAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGCTTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGTCGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGCAGTGTAGACCGGTTACGGCGTTCAcgccgc                                                                                                                                                                                                                                                 >  2:60930/1‑250 (MQ=255)
                 cATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                                                                                                                           >  1:303599/1‑251 (MQ=255)
                    tattCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCa                                                                                                                                                                                                                                         <  1:188288/250‑1 (MQ=255)
                                  ttGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGcc                                                                                                                                                                                                                        >  1:414685/1‑251 (MQ=255)
                                     ccGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTg                                                                                                                                                                                                                      >  2:400258/1‑250 (MQ=255)
                                      cGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCGTCCGGCATTTGTgcgctgatgcctg                                                                                                                                                                                                                       >  1:217272/1‑241 (MQ=255)
                                            aaTTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCACAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACg                                                                                                                                                                                                              <  1:319169/251‑1 (MQ=255)
                                                  tCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTgc                                                                                                                                                                                                        >  1:231546/1‑251 (MQ=255)
                                                   cAGCGCGATGCCAGCAGTATAGTGTTGCACACAGCGCTCAGCGCCGACGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGCTGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCtgacg                                                                                                                                                                                                        >  1:448783/1‑247 (MQ=38)
                                                    aGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCtgacgc                                                                                                                                                                                                       <  1:321044/251‑6 (MQ=38)
                                                        cgATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGGGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGACGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCtgacgcgtct                                                                                                                                                                                                   >  1:336010/1‑242 (MQ=38)
                                                          ataccAGCAGTATAGTGATGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCACAGTGCCTGATGATTTGATGTGGTAAACATAGATAAATCGTGCGCAGATGCCAGATGCGCCGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATCAGGCGTTCACGCCGCATCCGGCAGTCATGCGCCATACCTGATGCGACGCTTTGCGCGTCtt                                                                                                                                                                                                 >  2:315539/4‑250 (MQ=255)
                                                                gcagTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCatgc                                                                                                                                                                                          >  1:238767/1‑248 (MQ=37)
                                                                 cagTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGcc                                                                                                                                                                                         <  2:301043/251‑1 (MQ=37)
                                                                 cagTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGcc                                                                                                                                                                                         >  1:46228/1‑251 (MQ=37)
                                                                   gTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCt                                                                                                                                                                                        >  1:226031/1‑250 (MQ=37)
                                                                     ataGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTAc                                                                                                                                                                                      >  1:115675/1‑250 (MQ=37)
                                                                     ataGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCATATCATGCCTAc                                                                                                                                                                                      <  1:12194/250‑1 (MQ=37)
                                                                     ataGTGTTGAACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACg                                                                                                                                                                                     >  2:290791/1‑251 (MQ=255)
                                                                        gtgtTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGAGGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTGATGGGCTAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATTCGACGCGGACGCGTCTTATCATGCCTAcaaat                                                                                                                                                                                  <  2:448783/251‑5 (MQ=25)
                                                                         ggtTGCACAGCGCGCTCAGCGCCGAAGAGCAGGAGGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATGTGATGTGGTAAACATAGATAAGTCGTTCTCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCtt                                                                                                                                                                              <  2:231546/250‑1 (MQ=255)
                                                                          gtTGCACAGCGCGCGCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTGCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCtt                                                                                                                                                                              <  1:98586/250‑1 (MQ=255)
                                                                            tGCACAGCGCGCTCAGCGGCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCGGAGGATTTGATGTGGTAAACAGAGAGAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTAcaaatctg                                                                                                                                                                               <  2:86046/250‑8 (MQ=21)
                                                                                      gcTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTAcaaatctgtacgcgaaccg                                                                                                                                                                    >  1:78891/1‑233 (MQ=21)
                                                                                      gcTCAGCGACAAAGAGTAGGATGAATGCGGCCCGTATGCCTAGTTCAACGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCTCATATGCCGGATGCTACGCTATGCGCGTCTTCTCATGATTACGGGTTATTGTCGAAGTGTAGACCGGATAAGACGTTTACGCCGCTTCCGGCTTTTGTGCGCTGATGCCTGAGGCGACGCTGACGCGTCTTATCATGCCTActtatctgtactcgaatct                                                                                                                                                                    <  2:85960/251‑19 (MQ=255)
                                                                                       cTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGAACGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTAcaaatctgtacgcgaaccg                                                                                                                                                                    >  2:364456/1‑232 (MQ=21)
                                                                                         cAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGGGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTTCGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTTATGCCTGATGCGACCCTGACGCGTCTTATCATGCCTAcaaatctgtacgcgaaccgtag                                                                                                                                                                 <  1:328793/251‑22 (MQ=21)
                                                                                         cAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTACAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTAcaaatctgtacgcgaaccgt                                                                                                                                                                   >  1:390149/1‑230 (MQ=21)
                                                                                          aGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTAcaaatctgtacgcgaaccgtagg                                                                                                                                                                <  1:146486/251‑23 (MQ=21)
                                                                                           gcgcCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGGGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTAcaaatagtccgaaccgtaggccgg                                                                                                                                                               <  2:131173/251‑24 (MQ=21)
                                                                                             gcCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCgg                                                                                                                                                          <  1:168489/251‑1 (MQ=255)
                                                                                                gAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGAt                                                                                                                                                        >  1:173032/1‑250 (MQ=255)
                                                                                                             aagccgacCCGTGTGCCGACTTCAGAGTGCCTGAGGATTTGAGGTGGTAAACATAGACAAGTCGTTCTCAGATGCCTGATGCTACGCTATGCGCGTCTTATCGGGCCTACGGTTGATGGGCGAATTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTTCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                          <  2:142456/247‑1 (MQ=37)
                                                                                                                      cGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACG‑‑CGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCg                                                                                                                                   >  2:222855/1‑250 (MQ=14)
                                                                                                                       gTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                 >  2:339496/1‑250 (MQ=255)
                                                                                                                         aTGCCGACTTCAAAGTGCCTGAGGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCGGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACG‑‑CGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCgacct                                                                                                                               <  2:146633/251‑5 (MQ=2)
                                                                                                                           gCCGACTTCAAAGTGCATGATGATTTGATGTGGTAAACATAGATAAGTAGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATAGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGAGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCg                                                                                                                            >  1:262537/1‑251 (MQ=255)
                                                                                                                              gACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACG‑‑CGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCgacctgaaa                                                                                                                           >  1:403594/1‑242 (MQ=2)
                                                                                                                                   cAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACG‑‑CGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCgacctgaaaattctt                                                                                                                     >  2:98161/1‑237 (MQ=2)
                                                                                                                                    aaaGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACG‑‑CGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCgacctgaaaattctt                                                                                                                     >  1:190325/1‑236 (MQ=2)
                                                                                                                                          ccTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGc                                                                                                             >  2:36818/1‑251 (MQ=255)
                                                                                                                                             gatgatTTGATGTGGTAAACATAGATAAGTCGTGCGCAGTTGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAGGTTTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCat                                                                                                           <  1:401561/250‑1 (MQ=255)
                                                                                                                                                        gtgGTATACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTAtt                                                                                               >  1:118187/1‑251 (MQ=255)
                                                                                                                                                          ggTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCGTCGTGCTATTATCAACGCATATTCAGTTTATTgg                                                                                             >  2:90030/1‑251 (MQ=255)
                                                                                                                                                                 aTAGATTAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGt                                                                                       >  1:130216/1‑250 (MQ=255)
                                                                                                                                                                         gTCGTGGGCAGATGCCTGATGCGACGCTATGCGCGGTTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAatat                                                                              <  2:222019/251‑1 (MQ=255)
                                                                                                                                                                         gTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAatat                                                                              <  1:290791/251‑1 (MQ=255)
                                                                                                                                                                           cGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTGCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGc                                                                            <  1:364235/251‑1 (MQ=255)
                                                                                                                                                                                              cGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTgg                                                           >  1:211307/1‑248 (MQ=255)
                                                                                                                                                                                                     aTGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAaata                                                 <  1:280538/251‑1 (MQ=255)
                                                                                                                                                                                                        cgcgTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCATTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAataata                                              <  2:173032/251‑1 (MQ=255)
                                                                                                                                                                                                           gTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATGGCt                                           >  1:158597/1‑251 (MQ=255)
                                                                                                                                                                                                              ttATCAGGCCTACGGTTTATGGGCGAAGTGTAAACAGCATAAGCCGTTAACCCCGCCTCCGGCAGTCGAGCGCCATGCCTGATGCGACGATTTGCGCGTCTTATCAGGCCTACG‑‑GCTTAAGGGCGAAGTGTACGCCCGATAAGGCGTTAACGCCGCATCCGGCAGTCGCGCTATTATCAACGCATATTCAGTTTATTGGCGTGGAAGCCAATATGCACGCCAATAAAACTGGAGACAAAATAATAGCtggg                                        >  1:454310/1‑249 (MQ=255)
                                                                                                                                                                                                                         aCGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCt                             <  1:36818/251‑1 (MQ=255)
                                                                                                                                                                                                                         aCGGTTTATGGGCGAAGTGAAAACCGGATAAGGCTTCTACGCCGCATCCGGCAGTCTTGCGCCATGCTTGATGCGACGCGTTGCGCGGCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCTCCAATAAAACTGGAGACATAATAATAGCTGTGATCATAAccc                              <  2:421716/250‑1 (MQ=38)
                                                                                                                                                                                                                              ttATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGAGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTg                         >  2:119042/1‑250 (MQ=255)
                                                                                                                                                                                                                              ttATGGGAGAAGTGTAGACCGGAAAAGGCGTTCACGCCGCATCCGGCAGGCGGGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTAGCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTg                         <  2:17117/250‑1 (MQ=255)
                                                                                                                                                                                                                                  gggagAAGTGAAGCCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTTTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAAc                    <  2:153386/247‑1 (MQ=255)
                                                                                                                                                                                                                                        aGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGAtt               >  1:228067/1‑250 (MQ=255)
                                                                                                                                                                                                                                             aGACCGGATAGGGCGTTCACGCCGCATCCGGCAGGCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGt          <  2:130216/250‑1 (MQ=255)
                                                                                                                                                                                                                                              gACCGGAGAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCGTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAAGAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCt        <  1:60930/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    aTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGACAAGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTCCGGGCGAAATGTAGGCCGGATAAGGCGTTCACGCCGAATCCGGCAGACGAGCTATTATCAACACCTATTCAGTTTTTTGGCGCAGTAGGAGACATGCTCGGCAATAACACATGGGTACCAATACAAGCTTTGATTCTATCCCTCGTGATAACGGATTATGGTCTTGatat  >  2:79319/1‑251 (MQ=16)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                            
AGCGAAGGTGAATACGCCATTATTCCCGCCAGCGTTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATAT  >  NC_000913/377153‑377647

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