Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913377,4230GT38.7% 60.7 / 68.3 62intergenic (+112/+112)yaiL/frmBDUF2058 domain‑containing protein YaiL/S‑formylglutathione hydrolase FrmB
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base T (14/10);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTC  >  NC_000913/377187‑377669
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                            
ttGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGcc                                                                                                                                                                                                                                              >  1:414685/1‑251 (MQ=255)
   ccGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTg                                                                                                                                                                                                                                            >  2:400258/1‑250 (MQ=255)
    cGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCGTCCGGCATTTGTgcgctgatgcctg                                                                                                                                                                                                                                             >  1:217272/1‑241 (MQ=255)
          aaTTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCACAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACg                                                                                                                                                                                                                                    <  1:319169/251‑1 (MQ=255)
                tCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTgc                                                                                                                                                                                                                              >  1:231546/1‑251 (MQ=255)
                 cAGCGCGATGCCAGCAGTATAGTGTTGCACACAGCGCTCAGCGCCGACGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGCTGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCtgacg                                                                                                                                                                                                                              >  1:448783/1‑247 (MQ=38)
                  aGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCtgacgc                                                                                                                                                                                                                             <  1:321044/251‑6 (MQ=38)
                      cgATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGGGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGACGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCtgacgcgtct                                                                                                                                                                                                                         >  1:336010/1‑242 (MQ=38)
                        ataccAGCAGTATAGTGATGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCACAGTGCCTGATGATTTGATGTGGTAAACATAGATAAATCGTGCGCAGATGCCAGATGCGCCGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATCAGGCGTTC‑ACGCCGCATCCGGCAGTCATGCGC‑C‑ATACCTGATGCGACGCTTTGCGCGTCtt                                                                                                                                                                                                                       >  2:315539/4‑250 (MQ=255)
                              gcagTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCatgc                                                                                                                                                                                                                >  1:238767/1‑248 (MQ=37)
                               cagTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGcc                                                                                                                                                                                                               <  2:301043/251‑1 (MQ=37)
                               cagTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGcc                                                                                                                                                                                                               >  1:46228/1‑251 (MQ=37)
                                 gTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCt                                                                                                                                                                                                              >  1:226031/1‑250 (MQ=37)
                                   ataGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTAc                                                                                                                                                                                                            >  1:115675/1‑250 (MQ=37)
                                   ataGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCATATCATGCCTAc                                                                                                                                                                                                            <  1:12194/250‑1 (MQ=37)
                                   ataGTGTTGAACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACg                                                                                                                                                                                                           >  2:290791/1‑251 (MQ=255)
                                      gtgtTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGAGGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTGATGGGCTAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATTCGACGCGG‑ACGCGTCTTATCATGCCTAcaaat                                                                                                                                                                                                        <  2:448783/251‑5 (MQ=25)
                                       ggtTGCACAGCGCGCTCAGCGCCGAAGAGCAGGAGGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATGTGATGTGGTAAACATAGATAAGTCGTTCTCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCt‑t                                                                                                                                                                                                    <  2:231546/250‑1 (MQ=255)
                                        gtTGCACAGCGCGCGCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTGCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCt‑t                                                                                                                                                                                                    <  1:98586/250‑1 (MQ=255)
                                          tGCACAGCGCGCTCAGCGGCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCGGAGGATTTGATGTGGTAAACAGAGAGAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTAcaaatctg                                                                                                                                                                                                     <  2:86046/250‑8 (MQ=21)
                                                    gcTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTAcaaatctgtacgcgaaccg                                                                                                                                                                                          >  1:78891/1‑233 (MQ=21)
                                                    gcTCAGCGACAAAGAGTAGGATGAATGCGGCCCGTATGCCTAGTTCAACGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCTCATATGCCGGATGCTACGCTATGCGCGTCTTCTCATGATTACGGGTTATTGTCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCTTCCGGCTTTTGTGCGCTG‑ATGCCTGAGGCGACGCTG‑ACGCGTCTTATCATGCCTActtatctgtactcgaatct                                                                                                                                                                                          <  2:85960/251‑19 (MQ=255)
                                                     cTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGAACGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTAcaaatctgtacgcgaaccg                                                                                                                                                                                          >  2:364456/1‑232 (MQ=21)
                                                       cAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGGGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTTCGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGC‑TTATGCCTGATGCGAC‑CCTGACGCGTCTTATCATGCCTAcaaatctgtacgcgaaccgtag                                                                                                                                                                                       <  1:328793/251‑22 (MQ=21)
                                                       cAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTACAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTAcaaatctgtacgcgaaccgt                                                                                                                                                                                         >  1:390149/1‑230 (MQ=21)
                                                        aGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTAcaaatctgtacgcgaaccgtagg                                                                                                                                                                                      <  1:146486/251‑23 (MQ=21)
                                                         gcgcCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGGGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTAcaaatagtccgaaccgtaggccgg                                                                                                                                                                                     <  2:131173/251‑24 (MQ=21)
                                                           gcCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCgg                                                                                                                                                                                <  1:168489/251‑1 (MQ=255)
                                                              gAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGAt                                                                                                                                                                              >  1:173032/1‑250 (MQ=255)
                                                                           aagccgacCCGTGTGCCGACTTCAGAGTGCCTGAGGATTTGAGGTGGTAAACATAGACAAGTCGTTCTCAGATGCCTGATGCTACGCTATGCGCGTCTTATCGGGCCTACGGTTGATGGGCGAATTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTTCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                                                <  2:142456/247‑1 (MQ=37)
                                                                                    cGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTACAAATCTGTACG‑‑CGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCg                                                                                                                                                         >  2:222855/1‑250 (MQ=14)
                                                                                     gTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                                       >  2:339496/1‑250 (MQ=255)
                                                                                       aTGCCGACTTCAAAGTGCCTGAGGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCGGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTACAAATCTGTACG‑‑CGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCgacct                                                                                                                                                     <  2:146633/251‑5 (MQ=2)
                                                                                         gCCGACTTCAAAGTGCATGATGATTTGATGTGGTAAACATAGATAAGTAGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATAGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGAGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCg                                                                                                                                                  >  1:262537/1‑251 (MQ=255)
                                                                                            gACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTACAAATCTGTACG‑‑CGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCgacctgaaa                                                                                                                                                 >  1:403594/1‑242 (MQ=2)
                                                                                                 cAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTACAAATCTGTACG‑‑CGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCgacctgaaaattctt                                                                                                                                           >  2:98161/1‑237 (MQ=2)
                                                                                                  aaaGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGACGTTT‑ACGCCGCATCCGGCATTTGTGCGCTG‑ATGCCTGATGCGACGCTG‑ACGCGTCTTATCATGCCTACAAATCTGTACG‑‑CGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCgacctgaaaattctt                                                                                                                                           >  1:190325/1‑236 (MQ=2)
                                                                                                        ccTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGc                                                                                                                                   >  2:36818/1‑251 (MQ=255)
                                                                                                           gatgatTTGATGTGGTAAACATAGATAAGTCGTGCGCAGTTGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAGGTTTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCat                                                                                                                                 <  1:401561/250‑1 (MQ=255)
                                                                                                                      gtgGTATACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTAtt                                                                                                                     >  1:118187/1‑251 (MQ=255)
                                                                                                                        ggTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCGTCGTGCTATTATCAACGCATATTCAGTTTATTgg                                                                                                                   >  2:90030/1‑251 (MQ=255)
                                                                                                                               aTAGATTAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGt                                                                                                             >  1:130216/1‑250 (MQ=255)
                                                                                                                                       gTCGTGGGCAGATGCCTGATGCGACGCTATGCGCGGTTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAatat                                                                                                    <  2:222019/251‑1 (MQ=255)
                                                                                                                                       gTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAatat                                                                                                    <  1:290791/251‑1 (MQ=255)
                                                                                                                                         cGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTGC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGc                                                                                                  <  1:364235/251‑1 (MQ=255)
                                                                                                                                                            cGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑Tgg                                                                                 >  1:211307/1‑248 (MQ=255)
                                                                                                                                                                   aTGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAaata                                                                       <  1:280538/251‑1 (MQ=255)
                                                                                                                                                                      cgcgTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCATTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAataata                                                                    <  2:173032/251‑1 (MQ=255)
                                                                                                                                                                         gTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAATAATGGCt                                                                 >  1:158597/1‑251 (MQ=255)
                                                                                                                                                                            ttATCAGGCCTACGGTTTATGGGCGAAGTGTAAACAGCATAAGCCGTTA‑ACCCCGCCTCCGGCAGTCGAGCGC‑C‑ATGCCTGATGCGACGATTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TAAGGGCGAAGTGTACGCCCGATAAGGCGTTAACGCCGCATCCGGCAGTCGCGCTATTATCAACGCATATTCAGTTTATTGGCGTGGAAGCCAATATGCACGCCAATAAAAC‑TGGAGACAAAATAATAGCtggg                                                              >  1:454310/1‑249 (MQ=255)
                                                                                                                                                                                       aCGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAATAATAGCTGTGATCATAACCCt                                                   <  1:36818/251‑1 (MQ=255)
                                                                                                                                                                                       aCGGTTTATGGGCGAAGTGAAAACCGGATAAGGC‑TTCTACGCCGCATCCGGCAGTCTTGCGC‑C‑ATGCTTGATGCGACGCGTTGCGCGGCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCTCCAATAAAAC‑TGGAGACATAATAATAGCTGTGATCATAAccc                                                    <  2:421716/250‑1 (MQ=38)
                                                                                                                                                                                            ttATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGAGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAATAATAGCTGTGATCATAACCCTCGTg                                               >  2:119042/1‑250 (MQ=255)
                                                                                                                                                                                            ttATGGGAGAAGTGTAGACCGGAAAAGGCGTTC‑ACGCCGCATCCGGCAGGCGGGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTAGCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAATAATAGCTGTGATCATAACCCTCGTg                                               <  2:17117/250‑1 (MQ=255)
                                                                                                                                                                                                gggagAAGTGAAGCCCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTTTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAAc                                          <  2:153386/247‑1 (MQ=255)
                                                                                                                                                                                                      aGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGAtt                                     >  1:228067/1‑250 (MQ=255)
                                                                                                                                                                                                           aGACCGGATAGGGCGTTC‑ACGCCGCATCCGGCAGGCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGt                                <  2:130216/250‑1 (MQ=255)
                                                                                                                                                                                                            gACCGGAGAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCGTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAATAAGAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCt                              <  1:60930/251‑1 (MQ=255)
                                                                                                                                                                                                                  aTAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGA‑C‑AAGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TCCGGGCGAAATGTAGGCCGGATAAGGCGTTCACGCCGAATCCGGCAGACGAGCTATTATCAACACCTATTCAGTTTTTTGGCGCAGTAGGAGACATGCTCGGCAATAACACATGGGTACCAA‑TACAAGCTTTGATTCTATCCCTCGTGATAACGGATTATGGTCTTGatat                        >  2:79319/1‑251 (MQ=16)
                                                                                                                                                                                                                                     aTCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATCTTCAGTTTATTGGCGTGGTCGGCAAAATGCTCGCCAATAAAAC‑TGGAGCCAAAATTATAGCTGGGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTTGCAGATCt      >  2:179752/1‑251 (MQ=38)
                                                                                                                                                                                                                                         gtctgTCGTGGGC‑C‑ATGCCTGATTGGACGCTTTTCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTc  <  2:329857/247‑1 (MQ=255)
                                                                                                                                                                                                                                           cAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAg                                                                                                           <  1:102272/145‑1 (MQ=255)
                                                                                                                                                                                                                                           cAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAg                                                                                                           >  2:102272/1‑145 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                            
TTGCCGATAAAATTGCTCAGCGCGATGCCAGCAGTATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTC‑ACGCCGCATCCGGCAGTCGTGCGC‑C‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACG‑‑GCT‑TACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAAC‑TGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTC  >  NC_000913/377187‑377669

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