Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,933,4090TC31.6% 182.4 / 10.0 79E399G (GAA→GGA) edd6‑phosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (46/8);  new base C (0/25);  total (46/33)
Fisher's exact test for biased strand distribution p-value = 5.56e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGTACAACATCAGAAAGGTCAGAGAAGTCA  >  NC_000913/1933189‑1933642
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
accaccacaACCGCGGCCCAGAAACCCCGCTTCAAGGCCCGGCATACGTTCATGCGGCTTTAAAAAACCACCCGCTGGCGCTCCAACCCCGTGTTTCACACCCGGCACGGCAGGGGTTTTAATACCCGCGCGCCCCAGGTTACCGCTTAACCCTTTTGCCCCACCATGATGAGAGAAAGGTTGTTCGAAGAAAGCGTTCACTTGGCTGTAGAGTGTTTTTTCCGCCCCTCCCCGc                                                                                                                                                                                                                             <  2:928786/230‑1 (MQ=255)
                      acaCCCGCTTCAAAGGCCGGCAAAACGTCATGCTGGTTTTAAAAAACCACCGCTGGCGCTCCAATCCCCTGGTTCCCACCCGGCACGGAAGAGTTTTTCATAACCGCACGGCCCAGGTTACCGTTTAACACTTTTCCCCCACCATGAGGGGAGAAAGGTTGTTCGAAGGAAGCGACCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCa                                                                                                                                                                                                   <  2:146793/237‑1 (MQ=255)
                          ccGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGt                                                                                                                                                                                               >  2:1374024/1‑239 (MQ=255)
                          ccGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGt                                                                                                                                                                                               >  1:931648/1‑239 (MQ=255)
                               tCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAg                                                                                                                                                                                          >  1:167932/1‑239 (MQ=255)
                                 aacggcagacATACCGCCATGTCGACTTACAAAAACACCCGCGGGCGATTCAATCCCTGTTTTCCCACCCGGCACGGCAGGGGTTTTCAAACCCGCCCGGCCCGGTTACCCGCTTACCTCTTTGTCCCCCCAATGAGGAAAGAAGGGTTGTTGGAAGGAAGCGATCACATTGGTGGCGATTGATTTCCCCGCCCGTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAggg                                                                                                                                                                                        <  2:74265/230‑1 (MQ=255)
                                       ctgCATACCGTCTTGTTGTTTTTAAAAACCACCCGCTGTCCCTTCAATCCCCTGTTTCGCACCCGGCCCGGCGGAGTTTTTCATAACCGCACGCCCCATTTTACCGCTAAGCCTTTTGGCCCCCCCATGATGAGAGAAAGGTTGTTCAAGGGAAGCGATCACATTGCTGTCGGGTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATaa                                                                                                                                                                                   <  2:1043871/236‑1 (MQ=255)
                                                           ttaaaaaaaCCCCCGCTGCCCTTCCAATGCCCTGTTTCACACCCGCCCCGGCGTGGTTTTCCAAACCCGCACGGCCCAGTTTCCCCCTTAACTCTTTTGCCCCACCATGATGAGAGAAGGTTTGTTCGAAGGAAGCGATCACTTTGCTTTCGAGTGTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGcc                                                                                                                                                              <  2:219351/235‑1 (MQ=255)
                                                            tcaaaaaaCACCCGCTGGCGTTACACTCCCCTGGTTCCCACCCGGCACGGCAGAGTTTTTCATACCCGCACGGCCAGGGTAACCGCTTAACCTTTTTGCCCCACCATGAGGAGAGAAGGGTTGTCCAAGGGAAGCGATCACATTGCTGTGGAGTGATTTTCCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCa                                                                                                                                                             <  2:1299622/236‑1 (MQ=255)
                                                              cAAAAACAACCGCTGGCGCTTCAATGACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCAcc                                                                                                                                                           >  1:219351/1‑239 (MQ=255)
                                                                               cttttCAACCCCCTGGTTCCCAACGGCCAGGGCGGAGTTTTCCAAAACGCCAGGCCCCGGGTTCCCGCTAAACCTTTTTCCCCCACCAGGATGAGAGAAAGGTTGTCCGAAGGAAGCGATCCCTTGGCTGTCGAGTGATTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCa                                                                                                                                            <  2:1325159/235‑1 (MQ=255)
                                                                                               gTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCCGGCCt                                                                                                                                >  1:409796/1‑235 (MQ=255)
                                                                                               gTTCTCAACCGGCACGGCAGAGGTTTTCATAAACGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGCAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGGATAACGAGACAGACCAAAACCTGCCACCGTATTGACATATTCATGCAGCAGGCCt                                                                                                                                >  2:1131732/1‑235 (MQ=255)
                                                                                                          gcacggcaGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTc                                                                                                                  >  2:342559/1‑238 (MQ=255)
                                                                                                                caGAGGTTTTCATACCCCCCCGGCCAGGGTTCCCGCTTAACCCTTTTGTCCCACCATGATGAAAGAAGGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGGGATTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACgaga                                                                                                                                                                              <  2:373669/172‑1 (MQ=255)
                                                                                                                        ctCAAAACCGCGGGGCCCTGTTCACCGCTAAACGTTTTTGCCCCCCGATTAGGAGAGAAGGGTTGTACGAAGAAAGCGACCACATTGTTGTGGGGTTTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCa                                                                                                                                                                                            <  2:1494913/149‑1 (MQ=255)
                                                                                                                         tccaaacccGCGCGGCCAAGTTTCCCGCTAAACGCTTTGTCCCCACGAGGAGAAAGAAAGTTGTTGTAGAGGGAGGCCATCTCTTTGCTGTCGGTTGTTTTCCCCCCCCCTTCCCGCCAGTCCAGTTCCCCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGt                                                                                                  <  2:1589011/232‑1 (MQ=255)
                                                                                                                           cATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCaaaa                                                                                                                                                                    >  1:373669/1‑172 (MQ=255)
                                                                                                                                 cGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCa                                                                                                             >  1:1178330/1‑221 (MQ=255)
                                                                                                                                   cACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAAc                                                                                                                                                                                  >  1:1494913/1‑150 (MQ=255)
                                                                                                                                   cACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCa                                                                                            >  1:1253797/1‑236 (MQ=255)
                                                                                                                                        cccGGGTTCCCTCTTACCATTTTTCCCCCACATTGATGAGAAAAAGTTTGTCCAAAGTAAGCCATCACTTTGCTGCCGAGTATTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCa                                                                                                                                                                                                   <  2:540961/128‑1 (MQ=255)
                                                                                                                                        cccAGTTTCCCGCTTAACCCTTTTGTCCCCCCATGATGAGAAAAAGTTTTTTCGAAGGAAGCGATCACTTTGCTGTGGATTATTTTTCCCCCCCCTTCCCGCCAGTCCAGtt                                                                                                                                                                                                                   <  2:1210916/112‑1 (MQ=255)
                                                                                                                                        cccAGGTTCCCGCTTAACATTTTTGCCCCCCCATGATGAGAAAAAGGTTGTTCAAAGGAAGCGATCACATTGCTGTCGAGTGATTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGtt                                                                                                                                                                                              <  1:849438/133‑1 (MQ=255)
                                                                                                                                        cccAGGTTCCCGATAAACTCTTTTGCCCCACAAGGATGAGAGAAAGTTTGTTCGAAGGAAGCGATCACATTGCTTTCGGGTGATTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGtt                                                                                                                                                                                              <  2:1547003/133‑1 (MQ=255)
                                                                                                                                                   cGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGc                                                                                                                                                                                                                                        <  2:241579/82‑1 (MQ=255)
                                                                                                                                                   cGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGc                                                                                                                                                                                                                                        >  1:241579/1‑82 (MQ=255)
                                                                                                                                                   cGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCattatt                                                                                                                                                                                                         >  1:1210916/1‑112 (MQ=255)
                                                                                                                                                   cGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCa                                                                                                                                                                                                   <  1:566886/118‑1 (MQ=255)
                                                                                                                                                   cGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCa                                                                                                                                                                                                   >  2:566886/1‑118 (MQ=255)
                                                                                                                                                   cGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAgg                                                                                                                                                                                         >  1:540961/1‑128 (MQ=255)
                                                                                                                                                   cGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGtata                                                                                                                                                                                    >  2:849438/1‑133 (MQ=255)
                                                                                                                                                   cGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGtata                                                                                                                                                                                    >  1:1547003/1‑133 (MQ=255)
                                                                                                                                                   cGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGAc                                                                                                                                                   >  1:163233/1‑166 (MQ=255)
                                                                                                                                                   cGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTgg                                                                          >  2:1483259/1‑239 (MQ=255)
                                                                                                                                                   cGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTgg                                                                          >  1:1255916/1‑239 (MQ=255)
                                                                                                                                                     cTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGtt                                                                        >  1:1208042/1‑239 (MQ=255)
                                                                                                                                                        aaCACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTATGCCACCTGCCGCCTGGAAGTGGTTAat                                                                     >  1:699018/1‑239 (MQ=255)
                                                                                                                                                          ctcTTTTGCCCCACCATGATGAGAGAAAGGTTTTTGGAAGGAAGCGATCACTTTGCTGTCGGGTAATTTTCCCCCCCCTTCCCGccagtc                                                                                                                                                                                                                        <  2:740798/88‑1 (MQ=255)
                                                                                                                                                             ttttGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCCCATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCACAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCg                                                                 >  2:1226673/1‑238 (MQ=255)
                                                                                                                                                                  tCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCgg                                                            >  1:485335/1‑238 (MQ=255)
                                                                                                                                                                    ccaccaTGATGAGAGAAAGTTTGTTCGAAGGAAGCGATCACTTTGCTGTGGAGTTTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTc                                                                                                                                                                                             <  2:684162/107‑1 (MQ=255)
                                                                                                                                                                    ccaccaTGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCa                                                                                                                                                                                                              >  1:740798/1‑90 (MQ=255)
                                                                                                                                                                    ccaccaTGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGAcc                                                         >  1:1120566/1‑239 (MQ=255)
                                                                                                                                                                       ccaTGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGt                                                       >  1:1259842/1‑238 (MQ=255)
                                                                                                                                                                       ccaTGAGCAGAGAAGGTTTGTTCGAAGGAAGCGACCACTTTGGTGTCGTGTTTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGtt                                                      <  2:1255916/239‑1 (MQ=255)
                                                                                                                                                                        ccttgatgaGAAAAAGTTTTTTCGAAGGAAGCGATCACTTTGCTTTCGAGTGTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCa                                                                                                                                                                       <  2:1242187/123‑1 (MQ=255)
                                                                                                                                                                              tgaGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACgaga                                                                                                                                                                              >  2:268623/1‑113 (MQ=255)
                                                                                                                                                                              tgaGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACgaga                                                                                                                                                                              <  1:268623/113‑1 (MQ=255)
                                                                                                                                                                               gagagaAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATaa                                                                                                                                                                                   >  1:684162/1‑107 (MQ=255)
                                                                                                                                                                                  agaAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCa                                                                                                                                                             >  1:1242187/1‑126 (MQ=255)
                                                                                                                                                                                             ttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCa                                  >  2:117169/1‑239 (MQ=255)
                                                                                                                                                                                                  aggaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGtt                                                                                                                                                                                                                   >  2:474174/1‑57 (MQ=255)
                                                                                                                                                                                                  aggaagCGATCACATTGCTGTCGAGTGATTTTCCCGCCCCTTCCCGCCAGTCCAGtt                                                                                                                                                                                                                   <  1:474174/57‑1 (MQ=255)
                                                                                                                                                                                                    gactcacgcacaTTGCTGTAGAGTGTTTGGTCCGCCCCTTCCCGCCAGTCCAGTCCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCCCCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGTa                            <  2:409796/230‑1 (MQ=255)
                                                                                                                                                                                                    gaagCGATCACATTGCTGTCGAGTGATTTTTCTGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGCTCAAGGGCATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGATCGGGCAGAGACGAGCCATAAGCGGGac                           >  2:119609/1‑239 (MQ=255)
                                                                                                                                                                                                    gaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCa                                                                                                                 <  1:89457/153‑1 (MQ=255)
                                                                                                                                                                                                    gaagCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCa                                                                                                                 >  2:89457/1‑153 (MQ=255)
                                                                                                                                                                                                      agCCATCCCAGTGCTGTGAGTTACTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGTac                           <  2:316104/237‑1 (MQ=255)
                                                                                                                                                                                                       gCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGCCGTGCCATCAGCGGTacaac                        >  1:885509/1‑239 (MQ=255)
                                                                                                                                                                                                       gCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGGGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGTacaac                        >  2:19076/1‑239 (MQ=255)
                                                                                                                                                                                                            cactTTGCTGTGGAGTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCCCCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGTACAACATCAg                   <  2:1208042/239‑1 (MQ=255)
                                                                                                                                                                                                               tttGCTGTCGTTTTTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCattatt                                                                                                                                                                                                         <  2:753049/53‑1 (MQ=255)
                                                                                                                                                                                                               tttGCAGTCGCGTGTTTTTCCATCCCCTACCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGtt                                                      <  1:999091/200‑1 (MQ=255)
                                                                                                                                                                                                               gttGCTGTAGGGTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAg                                 <  2:314382/221‑1 (MQ=255)
                                                                                                                                                                                                               aTTGCTGTCGGTTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCa                                                                                                                                                                                                   <  2:370535/60‑1 (MQ=255)
                                                                                                                                                                                                                  gCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGTACAACATCAGAAAgg              >  1:1534393/1‑238 (MQ=255)
                                                                                                                                                                                                                   cTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCa                                                                                                                                                                                                              >  2:1458632/1‑45 (MQ=255)
                                                                                                                                                                                                                    tGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAgg                                                                                                                                                                                         >  1:652957/1‑65 (MQ=255)
                                                                                                                                                                                                                      tCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCa                                                                                                                                                                                                   <  1:1153100/53‑1 (MQ=255)
                                                                                                                                                                                                                      tCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCa                                                                                                                                                                                                   >  2:1153100/1‑53 (MQ=255)
                                                                                                                                                                                                                      tCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGTACAAAATCAGAAAGGCCaga         >  1:244296/1‑239 (MQ=255)
                                                                                                                                                                                                                        gAGTGATTTGCCCGCCCCTCCCCGTCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGTACAACATCAGAAAGGTCAgaga       <  2:931648/239‑1 (MQ=255)
                                                                                                                                                                                                                         cgtatcttttttCCCCCCTTCACCGCCAGTCCAGTTCACCATTATTCGGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGTACAACATCAGAAAGGTCaga         <  2:485335/232‑1 (MQ=255)
                                                                                                                                                                                                                         aGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGt                                                                                                                                                                                               >  1:753049/1‑54 (MQ=255)
                                                                                                                                                                                                                         aGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAgg                                                                                                                                                                                         >  1:370535/1‑60 (MQ=255)
                                                                                                                                                                                                                         aGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTTAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCCGGAAGTGGTTAATATCGGCCTGACCGTCTGGGCAgaga                                            >  2:999091/1‑201 (MQ=255)
                                                                                                                                                                                                                         aGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGGacaaca                       >  1:314382/1‑222 (MQ=255)
                                                                                                                                                                                                                           tGATTTTCCCGCCCCTTCCCGACAGTCCAGTTCACCAGTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGATGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGTACAACATCAGAAAGGTCAGAGAAg     <  2:1255570/238‑1 (MQ=255)
                                                                                                                                                                                                                             aTTTTTCCGCCCCTTCCCGCAAGTCCCGTTCACCATTATTCAGCCATGGTTCAAGGTCATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCGTTGAGCAGTTCACGCACCAGAACCGGTACGCCCCCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGGAGAGACGTGCCATCAGCGGCACAACATCAGAAAGGCAAGAGAAGTCa  >  2:1223843/1‑239 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
AACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCACCGTATTGACATCTTCATGCAGCAGGCCTGCTTTGAGCAGTTCACGCACCAGAACCGGTACGCCACCTGCCGCCTGGAAGTGGTTAATATCGGCCGGACCGTTCGGGTAGAGACGTGCCATCAGCGGTACAACATCAGAAAGGTCAGAGAAGTCA  >  NC_000913/1933189‑1933642

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