Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,097,4840CT21.9% 121.6 / 13.6 64S73F (TCC→TTC) hemWheme chaperone HemW
Reads supporting (aligned to +/- strand):  ref base C (12/38);  new base T (13/1);  total (25/39)
Fisher's exact test for biased strand distribution p-value = 4.33e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGCTGGACGCTTTACGTAATGGTTAAATTACCTCCGCTGAGTCTCTACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAAGCTGAAACGACTTGGGCGTATTCATGGCCCGCAAGAAGCGAAACGCGCGGCGAAGCTGGCGAGCG  >  NC_000913/3097247‑3097726
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
gctgctgGACGCTTTACGTAATGGTTAAATTACCTCCGCTGAGTCTCTACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATg                                                                                                                                                                                                                                        <  1:416219/251‑1 (MQ=255)
       gACGCTTTACGTAATGGTTAAATTACCTCCGCTGAGTCTCTACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTTCGGCCCGGCG‑CTGCAAACgc                                                                                                                                                                                                                                 >  1:12157/1‑251 (MQ=255)
           cTTTACGTAATGGTTAAATTACCTCCGCTGAGCCTCTACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGCAGTAAAGACAATCTTTATTGGCGGTTGGT‑CGCCGGGCCCGGTTTTCGGCCCGGCggttgacaactcgtt                                                                                                                                                                                                                              >  2:233945/1‑237 (MQ=255)
            tttACGTCATGGTTAAATTACCTCCGCTGAGTCTCTACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGGACGCCCGGCCCGCTTTTCGGCCCGGCG‑ATGCAAAcaccgcgg                                                                                                                                                                                                                            >  1:391550/1‑244 (MQ=255)
                 gTAATGGTTAAATTACCTCCGCTGAGTCTCTACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGA‑GGGGCGCCGAGCCTGCTTTCCGGCCCGGCG‑AGGCAAACGCTGCTGGGCgg                                                                                                                                                                                                                       >  1:528255/1‑251 (MQ=255)
                    aTGGTTAAATTACCTCCGCTGAGTCTCTACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCCAGGCTGCTTTCCCGCCCGGGC‑ATGCAAACGCCGCTggcgggcgg                                                                                                                                                                                                                    >  1:485898/1‑244 (MQ=255)
                               aCCTCCGCTGAGTCTCTACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGtt                                                                                                                                                                                                         <  1:261097/251‑1 (MQ=255)
                                        gAGTCTCTACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCt                                                                                                                                                                                                  <  1:263360/249‑1 (MQ=255)
                                        gAGTCTCTACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTgg                                                                                                                                                                                                <  1:352186/251‑1 (MQ=255)
                                            ttctACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATTTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGTAAACGCTGCGGGACGGCGTGCGTGCGCGTTTGCCGCTGGCat                                                                                                                                                                                             <  2:418447/249‑2 (MQ=255)
                                              ctaatgtCACAGCCCGGGGTGCGCGCAGAAATGCCCGTTATGCGAGTTCAACTCACACGCGTTGAAAGTAGAAGTGCCGCACGACGATTAGGATGAGCAACGGCTTAACGAGCGGGACAACGATGTGGCTTACGCGCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCT‑ATGCAAACGCTTCTTGACGGCGTGCGTGCGTGttttgcgcgggcag                                                                                                                                                                                             <  2:486489/242‑12 (MQ=255)
                                                     cacaTCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAg                                                                                                                                                                                    <  1:164976/250‑1 (MQ=255)
                                                     cacaGCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGGTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCAGCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGTCGGT‑GGTACTCCGAGCCTGCTGTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCTCGTTTGCCGCTGGCAGCGGATGCAg                                                                                                                                                                                    <  2:262033/250‑1 (MQ=255)
                                                         tCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGGACGCCGAGCCTGCTTTTCGGCCCCGCG‑AGGCAAAAGCTGCTGGGCGGCGGGCGTGGGCGGTTGCCGCCGGGAGCGGGTGGAGaaaat                                                                                                                                                                               >  1:195326/1‑249 (MQ=255)
                                                               ggtgCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGCGCCCGCTTTCCGGCCCCGCG‑ATGCAAACACTTCTGGACGGCGTGCGTGCGCGGTTGCCGCTGGCCGCGGGTGCCAAAATTACCATg                                                                                                                                                                         >  1:508367/1‑251 (MQ=255)
                                                                        aGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTCCGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AAc                                                                                                                                                               <  2:27006/251‑1 (MQ=255)
                                                                              gCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCAGCGGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGc                                                                                                                                                         <  2:528255/251‑1 (MQ=255)
                                                                                ccGTACTGCGATTTCAACTCTCACTCGGTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGTCAACGATGTGGCTTACGCTCAGGGCCGGGAAGTAAAGACAATCTGTATTGGCGGA‑GGTACGCCGAGCCTTCTTTCCTGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTTCGATGGAAGCG‑CACCCTTGCAc                                                                                                                                                       <  2:174093/251‑1 (MQ=255)
                                                                                              cAACTCTCACGCGTTGAAAGTAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGGTGTGGCTTACGCTCAGGGCCGCTAAGTAAATACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGTTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATc                                                                                                                                         <  2:426192/251‑1 (MQ=255)
                                                                                                     cccgcgTGGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCtt                                                                                                                                     <  2:486466/246‑1 (MQ=255)
                                                                                                     cACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGAAAACGATGTGGCTTCCGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCTCTTTGt                                                                                                                                  <  1:112664/251‑1 (MQ=255)
                                                                                                      aCGCGTTGAAAGGAGAAGTGCCGCACGATGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGAGGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGG‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑ACCCCTGGCACGGTAGAAGCCGATCGCTTTGTc                                                                                                                                 <  2:499892/251‑1 (MQ=255)
                                                                                                       cgcgTTTAAAGGAGAAGTGCCGCACGACGATTATTTCCCGCATCTGCTTAACGATCTGGTCAACGATGTGGCTTACTCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCg                                                                                                                                <  2:336259/251‑1 (MQ=255)
                                                                                                                aGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGc                                                                                                                       <  1:203958/251‑1 (MQ=255)
                                                                                                                 ggAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGGAAGCCGAGGCTGCTTTCCGGCCCCGGG‑AAGCAAACGCTGCGGGACGGCGGGCGGGGGCGGTTGCCGCTGGCAGCGGGTGCAGACAATAATATGGGTGCG‑AACACTGGCACGGGAGGAGACGAACGATTTGTCGCTTAtaaacg                                                                                                                      >  1:274990/1‑246 (MQ=255)
                                                                                                                  gagaAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTGTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGttt                                                                                                                                                                                                        <  1:282411/169‑1 (MQ=255)
                                                                                                                     aaGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCt                                                                                                                  <  2:12157/251‑1 (MQ=255)
                                                                                                                     aaGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCt                                                                                                                  <  2:391550/251‑1 (MQ=255)
                                                                                                                               cgacgaTTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAAcc                                                                                                        <  1:233945/251‑1 (MQ=255)
                                                                                                                                  cgaTTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTTCGGCCCCGCG‑ATGCCAACGCTGCTGGGCGGGGGGGGTGCGCGTTTGCCGCCGGCAGCGGATGCAGGAATTACTATGGGAACG‑GAACCTGGCACGGGAGAAGGCGGTCGCTTTGTTGGTTATTAGCGGGATGGTGTGGACCGcc                                                                                                     >  1:535085/1‑250 (MQ=255)
                                                                                                                                   gaTTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTTCGGCCCCGCG‑ATGCAAACGCTGCTGGGCGGCGTGCGTGCGCGTTTGCCGCCGGCAGCGGGTGCAGAAATTACTATGGGAGCG‑AACCCTGGCACGGGAGAAGCCGATTACTTTGTCGATTATCAGCGTGCTGGGGTGAACCCCAt                                                                                                    >  1:129070/1‑251 (MQ=255)
                                                                                                                                   gaTTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTCCTTTCCGGCCCCGCG‑CTGCAAACGCCGCTGGACGGCGTGCGTGCGCGGTTGCCGCCGGCAGCGGGTGCAGGAATTACTATGGAAGCG‑AACCCCTGCCCGGTAGAAGCCGGTCGCTTTGTCGATTTTCAGCGGGCTGGTGTGGACCGCAt                                                                                                    >  1:226683/1‑251 (MQ=255)
                                                                                                                                       aTGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTTCTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATctct                                                                                                <  2:386406/251‑1 (MQ=255)
                                                                                                                                       aTGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGG‑GGTACGCCCAGCCCGCTTTTCGGCCCGGCG‑ATGCCAACGCTGCTGGGCCGGGTGCGGGCGCGTTTGCCGCTGGCAGCGGGTGTGGAAATTACTCTGGGAGCG‑AACCCCGGCACGGGAGAAGCCGATCGCTTTTTCGATTTTCAGCGTGCTGGGGGGAACCGCacttct                                                                                                >  1:473674/1‑246 (MQ=255)
                                                                                                                                          ttCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGGAAGCCCAGACCGCTTTCCCGCCCGGCG‑ATGCAAACACTGCTGGACGGCGGGCGTGCGCGGTTGCCGCTGGCCGCGGAGGCAGGAATTTATATGGGAACG‑GACCCCGGCCCGGTGGCAGCCGAACGCGTTTTTGATTTTCAGCGTGCTGGGGTGAACCGccgctctctt                                                                                             >  2:283378/1‑242 (MQ=255)
                                                                                                                                                atcggCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGCCAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTAGGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGtgtg                                                                                       <  1:397919/247‑1 (MQ=255)
                                                                                                                                                aTCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTTCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGGGCGTGCGCGGTTGCCGCCGGCCGCGGATGCAGAAAATACTATGGAAGCG‑AACCCTGGCACGGGAGAAGCCGGTCGCTTTGTCGATTATCAGCGGGCTGGTGTGAAACGCctctctctttgtggg                                                                                       >  1:209550/1‑242 (MQ=255)
                                                                                                                                                    gCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGc                                                                                      <  2:160838/248‑1 (MQ=255)
                                                                                                                                                      ttAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGGACGCCGAGCCCGCTTTTCGGCCCGGCG‑GTGCAAACGCTGCTGGACGGCGGGCGTGCGCGGTTGGCGCTGGCAGCGGAGGCAGAAATTACTATGGGAGCG‑GACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATTAGCGTGCTGGTGTGGACCGCATCTCTATTTGTGTGGagagg                                                                                 >  1:354947/1‑250 (MQ=255)
                                                                                                                                                        aaCGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGG‑GGTACGCCGAGCCTGCTGTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AATCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGtt                                                                                <  1:231073/250‑1 (MQ=255)
                                                                                                                                                          cGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGGGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGGCGGCGGGCGTGCGCGGTTGTCGCTGGCCGCGGGTGCAGAAATTACTATGGAAGCG‑AACCCCGGCACGGTAGAAACCGATCGCTTTGTTGGTTATTAGCGGGCTGGGGGGAAACGCATCTCTATTtgtggagatcggaca                                                                              >  1:19657/1‑239 (MQ=255)
                                                                                                                                                          cGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTAGTGGAGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGgtgt                                                                                        <  2:19657/240‑1 (MQ=255)
                                                                                                                                                                ggACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCgaa                                                                        <  1:47472/250‑1 (MQ=255)
                                                                                                                                                                 gACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGGACGCCGAGCCTGCTTTCCGGCCCCGCG‑GTGCAAACGCTGCTGGACGGCGGGCGGGCGCGGTTGCCGCTGGCAGCGGGTGCAGGAAATACTATGGGAGCG‑AACCCCGGCACGGGAGAAGCCGATCGCCTTTTCGGTTATCAGCGTGCTGGTGTGGACCGCATCTCTCTTGGTGTGCAGAGTTTTTGCga                                                                         >  1:44183/1‑248 (MQ=255)
                                                                                                                                                                  aCAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGGACGCCCAGCCTGCTTTTCGGCCCCGCG‑ATGCAAACGCTGCTGGACGGCGGGCGGGCGCGTTTGCCGCCGGCAGCGGGTGCAGAAATTAATTTGGAAGCG‑GACCCCGGCACGGGAGAAGGCGATCGACTTGTCGGGTATCAGCGGGCTGGTGTGGAACGCATCTTTATTTGGGGGCCGAGGTTTAGCgaaga                                                                      >  1:415984/1‑250 (MQ=255)
                                                                                                                                                                   cAACGAGGGGGCTTTCGCTCGAGGCTGTGAAGTAAAGACAATCATTATTGGCGGT‑GGGACGCCGAGCCTGCTTTTCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGGCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTTGTTTGCAGAGTTTTAGCGAAGaaa                                                                    <  2:314196/251‑1 (MQ=255)
                                                                                                                                                                       gATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑AGTAAGCCGAGGCCGCTTTTCGGCCCCGGG‑GTGCAAACCCCGCTGGGCGGGGTGCGTGCGGGTTTTCCGCTGCCCGCGGGTGCCCAAATTACTATGGGAGCG‑AAACCTAGCACGGTAGAAGCCGATCGCTTTTTCGATTATTAACGCGCCGGTGTGAACCGCCTCTCTCTTGGGGTGGAGAGTTTTAGCGCAGaaaacct                                                                >  1:448602/1‑248 (MQ=255)
                                                                                                                                                                            ggCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAAGCTGAAAc                                                           <  1:223706/251‑1 (MQ=255)
                                                                                                                                                                                    cTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGC‑GGGGAGACGGACCTGCTTTCCCGGCCGGCGGTTGCAAA‑GCTGCCGGAAGGCGTGCGCGTGCGTTTTGCGCCTGCCGCGGATGCCGAAATTACTATGGAAGCG‑AACCCTGGCGCGGGAGAAGCCGGTCGCTTTTTCGAATATCAACGTGCCGGTGGGGAACGCATCTCTATTGGTGTGGAGGGATTTTGCGGAGGAAAGCTGaaaagaatttgg                                                   >  1:104100/1‑242 (MQ=255)
                                                                                                                                                                                       ggggCCGTGAAGTAAAGACAATCTTTAGTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTGTGTCGATTATCAGCGTGCTGGTGTGAACCGCAGCTCTATTGGTGTGCAGAGTCTTAGCGAAGAAAAGCTGAAACGACGTGGGCGt                                                <  2:485898/250‑1 (MQ=255)
                                                                                                                                                                                       aGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACACCGACCCCGCTTTCCCGCCCGGCG‑CAGCAAACGCTGCTGGGCGGCGGGCGTGCGCGTTTGCCGCTGGCCGCGGGTGCAGAAATTATTACAGGAGCGCAACCC‑GGGACGGTAGAAGCCGAACGCTCTGTCGGTTATCAGCGGGCTGGTGTGGACCGCACCTCTATTGGTGTGCAGAGTTTTTGCGACGAAAAACGGAAACGACTTgggggg                                                >  1:504122/1‑248 (MQ=255)
                                                                                                                                                                                                aaGTAAAGACATGCTTTATTGGCGGG‑GGTACGCCGAGCCTGCTTGCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAAGCTGAAACGACTTGGGCGTATTCATGGc                                       <  2:473674/251‑1 (MQ=255)
                                                                                                                                                                                                 aGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAAGCTGAAACGACTTGGGCGTATTCATGGcc                                      <  1:168711/251‑1 (MQ=255)
                                                                                                                                                                                                 aGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAAGCTGAAACGACTTGGGCGTATTCATGGcc                                      <  1:153254/251‑1 (MQ=255)
                                                                                                                                                                                                  gTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTTCCGCCCGGGG‑GTGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGGAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGGTTTAGCGAAGAAAAGCTGAAACg                                                          >  1:318508/1‑230 (MQ=255)
                                                                                                                                                                                                  gTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAAGCTGAAACg                                                          <  2:318508/230‑1 (MQ=255)
                                                                                                                                                                                                           aTCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTTCCGCCCCGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTTCCGCTGGCAGCGGGTGCAGAAATTAATTTGGAAGCGGAACCC‑GGCACGGGAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGGGTGAACCCCCTCTCTATTTGGGTGACGGGTTTTAGGGAAGAAAAGCTGCAACGACTTGGGCGTGTTCATGGCCCGCaagaagg                            >  1:332797/1‑250 (MQ=255)
                                                                                                                                                                                                               ttATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAATCTGAAACGACTTGGGCGTATTCATGGCCCGCAAGAAGc                            <  1:71592/247‑1 (MQ=255)
                                                                                                                                                                                                                tATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAAGCTGAAACGACTTGGGCGTATTCATGGCCCGCAAGAAGCGaaa                        <  1:62889/250‑1 (MQ=255)
                                                                                                                                                                                                                           ggtACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAAGCTGAAACGACTTGGGCGTATTCATGGCCCGCAAGAAGCGAAACGCGCGGCGaa             <  2:195326/251‑1 (MQ=255)
                                                                                                                                                                                                                                gCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCt                                                                                                                  <  1:478468/145‑1 (MQ=255)
                                                                                                                                                                                                                                gCCGAGCCTGCTTTCAGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGCTCGCTTTGTCGATTATCAGCGTGCt                                                                                                                  >  2:478468/1‑145 (MQ=255)
                                                                                                                                                                                                                                  cGAGCCTGATTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGAAGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAAGCTGAAACGACTTGGGCGTATTCATGGCCCGCAAGCAGCGAAACGCGCGGCGAAGCTGGCg      >  2:261394/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ccTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAAGCTGAAACGACTTGGGCGTATTCATGGCCCGCAAGAAGCGAAACGCGCGGCGAAGCTGGCGAGCg  <  1:283378/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
GCTGCTGGACGCTTTACGTAATGGTTAAATTACCTCCGCTGAGTCTCTACATTCACATCCCGTGGTGCGTGCAGAAATGCCCGTACTGCGATTTCAACTCTCACGCGTTGAAAGGAGAAGTGCCGCACGACGATTATGTTCAGCATCTGCTTAACGATCTGGACAACGATGTGGCTTACGCTCAGGGCCGTGAAGTAAAGACAATCTTTATTGGCGGT‑GGTACGCCGAGCCTGCTTTCCGGCCCGGCG‑ATGCAAACGCTGCTGGACGGCGTGCGTGCGCGTTTGCCGCTGGCAGCGGATGCAGAAATTACTATGGAAGCG‑AACCCTGGCACGGTAGAAGCCGATCGCTTTGTCGATTATCAGCGTGCTGGTGTGAACCGCATCTCTATTGGTGTGCAGAGTTTTAGCGAAGAAAAGCTGAAACGACTTGGGCGTATTCATGGCCCGCAAGAAGCGAAACGCGCGGCGAAGCTGGCGAGCG  >  NC_000913/3097247‑3097726

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