Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,851,6540TG30.9% 111.1 / 21.2 69G218G (GGT→GGGhypBhydrogenase isoenzymes nickel incorporation protein HypB
Reads supporting (aligned to +/- strand):  ref base T (11/36);  new base G (21/0);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.87e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGAAGGCGACCAGCAAACCGTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTG  >  NC_000913/2851409‑2851900
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                
aTTGAAGGCGACCAGCAAACCGTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAg                                                                                                                                                                                                                                                             <  1:311896/250‑1 (MQ=255)
    aaGGCGACCAGCAAACCGTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACaa                                                                                                                                                                                                                                                         >  1:40519/1‑250 (MQ=255)
       gCGAACAGCAAACCGTGAACGATGCCGCACGCATTCGTGTTCCCGGCACACCAGCGATTCAGGTGAACACCGGTAAATGCTGCCATCTGGACGCACAGATGAGTGCCGACGCCGCCCCGCGTCTGCCACGGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCa                                                                                                                                                                                                                                                     <  2:485172/251‑1 (MQ=255)
         gACCAGCAAACCGTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGACTTCCCGTAACGAAGGTGAAGACAAACCAcc                                                                                                                                                                                                                                                   >  1:133787/1‑250 (MQ=255)
            cAGCAAACCGTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGaa                                                                                                                                                                                                                                                <  2:359750/251‑1 (MQ=255)
                       gAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCCCTTCCGAAGGGGAAAGCAAAACCCTGAAATATCCGCat                                                                                                                                                                                                                                      >  1:246377/1‑250 (MQ=255)
                          cGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTAGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCCGTACCGCAGGGGCAGACaaaaccctggaataaccgactattt                                                                                                                                                                                                                                  >  1:374499/1‑229 (MQ=255)
                                cgcacgcaTTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCCTTTCCCAAGGGGAAGAAAAACCACTGAAAAATCCGCataattttgttc                                                                                                                                                                                                                            >  1:508743/1‑242 (MQ=255)
                                     gcaTTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGTTTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTgccgccgcc                                                                                                                                                                                                                       <  1:201405/251‑1 (MQ=255)
                                           gcgcTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTTCGGTACCGAAGGGGAAAGAAAAACACTGAAAAATCCCCATATGTTTGCCGCCGCCCCGcgg                                                                                                                                                                                                                 >  1:500908/1‑249 (MQ=255)
                                                aCCGGCACACTAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGGCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCAGATGTTTGCCGCCGCCTCGCTGAtg                                                                                                                                                                                                             <  1:448447/249‑1 (MQ=255)
                                                           cgcgacttaggtgaaccacagTAAACGCTGCCATCTTGACGCACAATTGATTGGCCATGCCGCTCCCGGTCTTGCACTGGACGATTAATGTAATATGTTTTTTTAAAACGTTGGCAACCTCGTTTGCCCGGCCAACTTCGATCTCGGTGAAAAACACAAAGGGGCGGTTCTTTCCGTTTCCCAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAg                                                                                                                                                                                                <  1:200010/231‑1 (MQ=255)
                                                           aGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAg                                                                                                                                                                                                <  1:329238/251‑1 (MQ=255)
                                                             cGATTCAGGTGAACACCGGTAAAGGCGGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGtt                                                                                                                                                                                              <  1:286210/251‑1 (MQ=255)
                                                                     gTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCGCAACAAAGTTGACCTGtt                                                                                                                                                                                     <  1:68201/251‑1 (MQ=255)
                                                                        aaCACCGGTAAAGGCTGCCATCTTGACGCACAGAGGATGGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGcc                                                                                                                                                                                  <  2:154497/251‑1 (MQ=255)
                                                                                   aGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCCAAGGGGAAGACAAACCACTGAAATATCCCCATATGTTTGCCGCCGCCTCGCTGATGCCGCTCAACAAAATTGACCTGTTGCCGTATCTCCACt                                                                                                                                                                      >  1:212767/1‑251 (MQ=255)
                                                                                        gCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTTCGTTACCGAAGGGGAAGACAAAACAAAGAAATATCCCCATATGTTTTCCGCCCCCCCGCTGATGCTGCTCAAAAAAGATGACCTGTTGCCGTATTTCAACttttac                                                                                                                                                                 >  1:252917/1‑248 (MQ=255)
                                                                                          cactcagacgcgacgatgctTGCCGACGCCGCACCGGGTGTGCCACTGGTCGAGAACGGTATTCTGTTTATCGAAAGCGTTGGCACCCTCGTATGCCCGGCCAGCTTCGATCTGGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCTACTTTGACGt                                                                                                                                                              <  2:83130/237‑1 (MQ=255)
                                                                                             cTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTAACGAAAGTGAAGACAAACCACTGGGATAACCGCATATGTTTGCCGCCGCCCCGCTGATGCTGCTCAACAAAGATGAACAGTTGCCGCATCTCAAATTTTACGTTgt                                                                                                                                                          >  1:162954/1‑250 (MQ=255)
                                                                                               ggACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGtt                                                                                                                                                             <  2:247381/246‑1 (MQ=255)
                                                                                                    cacaGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCCTTCCGTTTCCGAAGGGGAAGACAAACCACTGAAAAATCCGCATATTTTTGCCGCCGCCCCGCTGATGCCGCTCAACAAAGTTAGACTGGTGCCGTATCTCACCATTGACGTTGGAAAAggc                                                                                                                                                   >  1:366735/1‑249 (MQ=255)
                                                                                                       gtggtattgaggtcgccgcACCGCGTCTGCCCCTAGACGAAAACGGTATTCCGATTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCACCTTCGATCGCGGGGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTGGGGCGCCGCCTCGCTGATGCTGCTCAAGAAAGTTGACCTGTTGTCGTAACTCAACGTTGACGGTGAAAAGTGCATc                                                                                                                                               <  2:173869/238‑1 (MQ=255)
                                                                                                        gaggatggCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCg                                                                                                                                              <  2:301071/244‑1 (MQ=255)
                                                                                                             ttgccgagtccccACCGCGTCTGCCACCGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGt                                                                                                                                                    <  2:90395/228‑1 (MQ=255)
                                                                                                             ttGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGTTTTCCGTTTCCGCAGGGGAAAACAAACCACTGGAAAATCCGCATATGTTTGCCGCCGCCCCGCTGATGCTTGTCAACAAAGTTGACCTGTTGCCGTATCTTAAATTTGTCGTTTAAAAgcagattgggca                                                                                                                                          >  1:90395/1‑239 (MQ=255)
                                                                                                                       cgcACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGGTACCCCAGGGGAAGAAAAAACACTGAAATAACCGCCTAAGTTTTGCGGCGCCTCGCTGATGCTTCTCAACAAAATTGACCTGTTTCCGTCTCTCCACTTTTAAGTTGGAAAGAGCactgcctccccccgcca                                                                                                                                 >  2:52511/1‑233 (MQ=255)
                                                                                                                          cccGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCa                                                                                                                            <  1:104929/250‑1 (MQ=255)
                                                                                                                              ccttcgcgtcTGGCCGATAACGGGATTCTAGTTAGCGAAAACGTCGGCAACCGCGTGTGCCGGGCCAGATTCGATCTCTGTAAAAAACACAAAGTGGCGGCGCTTTCCGGTACCGAAGGTGAATACAAACCTCTGAAATATCCGCATACTTTTGCCGCTGCCTCGCTGTTTCTGCTGAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAATTGCACCGCCTGCGCCCGCGATGTCTATcc                                                                                                                        <  2:520203/242‑1 (MQ=255)
                                                                                                                                   gCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTTCCGAAGGGGAAGACAAAACACCGAAATAACCGCATATGTTTGCCGCCGCCTCGCTGGTGGTGCTCAACCAAGATGACCTGTTGACGGATTTCAACTTTTACGGTTAAAAAGGCATCGCCTGCGCCCGCCGAGTCAAACCAGaa                                                                                                                    >  1:107553/1‑250 (MQ=255)
                                                                                                                                   gCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTATCGAAGGGGGAGACAAACAACAGAAAAATCCGCCTATGTTTGCCGCCCCCTCCCGGAAGCCGCCCAACAAAATTTACCCTTTTCCGTATCTCACCTTTGTCTTTTAAAAAAGCATCGCCCGCGCCCGCGAAGGCAAACCaaaa                                                                                                                    >  2:217347/1‑247 (MQ=255)
                                                                                                                                           aCGATAACGGTATTCTGTTTATCGAAAACGTTGGCAGCCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAatc                                                                                                          <  1:238326/251‑1 (MQ=255)
                                                                                                                                               tAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAg                                                                                                                      <  1:186594/236‑1 (MQ=255)
                                                                                                                                                        tCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGAGCTTTCCGTTACAGAAGGGGCAGACAAACCACTGTAATATTCCCATATTTTTGCCGCCCCCTCGCTGATGCTGCTCAAAAAAATTTGCCTGTTTGCGTATCTCAACTTTGTCGTTTAAAAATGCATCGCCCGCCCCCCCGAAGTCAATCCCGAGATTTAAAACATCCTTTTTTccc                                                                                            >  1:107882/1‑250 (MQ=255)
                                                                                                                                                          tGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTTCGTTACCGGAGGGGAAGAAAAAACACTGGAATATCCGCATATGTTTGGCGCCGCCTCGCTGATGCTGCTCCACCAAGTTGTCCTGTTGCCGTATCTCCACTTTTACGTTTAAAAGAGGATCGCCCTCGCCCCCGAAGTCAATACAGAAAATTAAATCATCCTTATTTCCGcc                                                                                          >  1:220975/1‑251 (MQ=255)
                                                                                                                                                           gTTTATCGAAAACTTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGcc                                                                                          <  2:246377/250‑1 (MQ=255)
                                                                                                                                                           gTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGCGCCTTTCGTTTCCGAAGGGGAAAAAAAACCAACGAAATATTCCCATATGTTTGTCGCCCCCCCGCTGGATCTGCTCCACAAAGTTGACCTTTTGCCGCATCTCAACCTTGACGTTTGAAAATGCATCGCCCGCCCCCCCGAAGACAAACCACAAATTTGAAatatcacttttttccccc                                                                                           >  2:207056/1‑235 (MQ=255)
                                                                                                                                                                 cGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAACCACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAgcg                                                                                   <  2:459857/251‑1 (MQ=255)
                                                                                                                                                                  gAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTTCGTTTCCGCAGGGGAAAAAAAACCACCGAAAAATCCGCCTATTTTTTCCCCCCCCCCCCTGGTGCTTCTCAACAAAGTTGAACTGTTTCCGGATCTCCACTTTGTCGTTGAAAAGTGCCTCGCCCGGCCCCCCGGAGTCAAACCCGAAATTTAAATCATCCTTTTTTCCGCCCCCAccgg                                                                                  >  2:136554/1‑251 (MQ=255)
                                                                                                                                                                          tGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGCGCCTTTCCTTTAAGCAAGGGAAAACCAAACACCTAGATATCCGCATATGTTTGCCGCCGCCCCGCTGGATGTTGTCAACAAAAGTGACCTGTTGCCGGATCTCCCACTTTACGTTGAAAAGTGCATCGCCCGTCCCCGCGACGTCAACCCCCAAAATGAAAACATCCTTATTTACGccccccgccacacaaggg                                                                          >  2:123050/1‑235 (MQ=255)
                                                                                                                                                                               aCCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGAGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAAGCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGAc                                                                    <  2:357579/251‑1 (MQ=255)
                                                                                                                                                                                   cGTATGCCCGGCCAGCTTCGATCTCGGTCAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGATGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATGGAAATCACCCTTATTTCCGCCACCAGCGGCGAAGGGATGAACCAGt                                                                <  2:108434/251‑1 (MQ=255)
                                                                                                                                                                                    gTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCCGTTCCCAAGGGGAAGACAAACCCCTGAAAAATCCGCATATGTTTTCCGCCCCCCCGCTGGAGCTGCTCAACAAAGTTGACCTGTTGCCGTATTTCAACATTGACGTTGAAAAATGCATCGCCTGCGCCCCCCAAGTCAATCCAGAAATTTAAATACTCCTTTTTTCCCCCCCCCGCGGCGAAGGGATgggacagtg                                                               >  2:305093/1‑244 (MQ=255)
                                                                                                                                                                                        gCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCt                                                            <  1:450432/250‑1 (MQ=255)
                                                                                                                                                                                        gCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTGCCGCCACCAGCGGCGAAGGGATGGACCAGTGGcg                                                            <  2:52168/250‑2 (MQ=255)
                                                                                                                                                                                           cGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGCAGGGGAAGAAAAACCAATTAAATATTCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAATTGACCTGTTGCCGGATCTCAACCTTGACGGTGTAAAGAGCCTCGCCCGCGCCCGCGGCGTCCAACCAGAAATTGAAAATATCCTTATTTTCCCCCCCAGCGGCGAAGGGGTGGGCCCATGGCTGaa                                                         >  2:383685/1‑250 (MQ=255)
                                                                                                                                                                                             gCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCAGCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAActg                                                      <  2:302019/251‑1 (MQ=255)
                                                                                                                                                                                              ccAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTCCTTTCCGTTTCCGCAGGGGAAAACAAACCACTGAAAAATAAGCATATGTTTTCCGCCGCCTCGCTGAAGCTTCACAAAAAAAGTGACCTGTTTCCCTATTCCAACATTGTCGTTGTAAAGAGCATCGCCCGCGCCCCCGCAGACAATACCGAAAATGAAATCACTCCTATTTCCCCCCCCAGCGGCGCAAGGGTGGAcccgaggcgtgactgg                                                     >  2:384946/1‑237 (MQ=255)
                                                                                                                                                                                                    tCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGGTACCGAAGAGGGAGACAAAACACCTAAATTTTCGGATATGATTGCCGCCGCCTCGCTGATTCCGCTCCAAAAAGATGACCCGTTTCCGTATCTCAACTTTGACGTTGAAAAGGGCATCGCCTGCGCCCCCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGGCAACCGCGGGGGAGGGAAGGGCAAGTGGCTGGActggccgggg                                               >  1:181135/1‑246 (MQ=255)
                                                                                                                                                                                                      gATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGCTAACGCAAGGTAAGACAAACAACTGAAAAAACCGCCTATGTTTGACGCCGCCCCGCTGGATCTGCTCAAAAAAATTTACCTGTTGCCGTATCTCAACTTTTAAGTTGAAAAGAGGAACGCCCGCGCCCGCGGAGGCAATCCAAAAATTTAAAATAACCTTATTTCCCCCACCCACGGCGAAGGGATGGGCCAGTGGCCGGACTTGCTGGGGac                                             >  2:389916/1‑251 (MQ=255)
                                                                                                                                                                                                           cGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGac                                             <  1:478550/246‑1 (MQ=255)
                                                                                                                                                                                                           cGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGAcaca                                          <  1:341527/249‑1 (MQ=255)
                                                                                                                                                                                                           cGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGGGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTTCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGc                                        >  1:23340/1‑251 (MQ=255)
                                                                                                                                                                                                                  aaaCACAAAGTGGCGGTGCTTTCCGTTACCGAAGGGGAAGACAAACCACTGAAAAATCCCCCTATGTTTTCCGCCCCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGCATCTCAACTTTTACGTTGAAAAATGCATCGCCCCCGCCCCCGCAGACAATCCCGAAATTGTAAACCTCCCTTTTTCCGCCCCCAGCGGCGAAGGGGAGGGCCAGTGGCTGAACTGGGTGGGGACACACCGAGGTGc                                 >  2:380519/1‑251 (MQ=255)
                                                                                                                                                                                                                      acaAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAAGCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAg                             <  2:350617/251‑1 (MQ=255)
                                                                                                                                                                                                                           gTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGtt                        <  2:252917/251‑1 (MQ=255)
                                                                                                                                                                                                                            gggCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAAGATCCGCATATGTGTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTc                       <  2:415693/250‑1 (MQ=255)
                                                                                                                                                                                                                             ggCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTcc                      >  1:324055/1‑251 (MQ=255)
                                                                                                                                                                                                                             ggCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTc                       <  1:122862/250‑1 (MQ=255)
                                                                                                                                                                                                                              gCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAActg                                                      <  1:53777/218‑1 (MQ=255)
                                                                                                                                                                                                                                ggTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCgg                   >  1:71366/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGCCGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGc                  >  2:418615/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCTGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGcc                 >  1:77331/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCt                                                            <  2:521055/201‑1 (MQ=255)
                                                                                                                                                                                                                                           cGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCt                                                            >  1:521055/1‑201 (MQ=255)
                                                                                                                                                                                                                                                 ccGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGccccc     >  2:142436/1‑247 (MQ=255)
                                                                                                                                                                                                                                                   gAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCAtt   <  1:508322/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTg  >  1:125040/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     aaGGTGAAGACAAACCACTGAAAAATCCGCTTATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTg  <  2:366735/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                
ATTGAAGGCGACCAGCAAACCGTGAACGATGCCGCACGCATTCGCGCTACCGGCACACCAGCGATTCAGGTGAACACCGGTAAAGGCTGCCATCTTGACGCACAGATGATTGCCGACGCCGCACCGCGTCTGCCACTGGACGATAACGGTATTCTGTTTATCGAAAACGTTGGCAACCTCGTATGCCCGGCCAGCTTCGATCTCGGTGAAAAACACAAAGTGGCGGTGCTTTCCGTTACCGAAGGTGAAGACAAACCACTGAAATATCCGCATATGTTTGCCGCCGCCTCGCTGATGCTGCTCAACAAAGTTGACCTGTTGCCGTATCTCAACTTTGACGTTGAAAAGTGCATCGCCTGCGCCCGCGAAGTCAATCCAGAAATTGAAATCATCCTTATTTCCGCCACCAGCGGCGAAGGGATGGACCAGTGGCTGAACTGGCTGGAGACACAGCGATGTGCATAGGCGTTCCCGGCCAGATCCGCACCATTG  >  NC_000913/2851409‑2851900

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