Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,548,9250AC25.6% 166.2 / 12.0 83G382G (GGT→GGGmalQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base A (35/26);  new base C (4/17);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 2.62e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCACAAACGCAGCATCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCGTCC  >  NC_000913/3548699‑3549147
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               
accacAAACGCAGCATCGACATCCCATGGTCAATTCGTAACGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAAAGGTTCATCGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAAAGGGCCGGGGGAATCCCGCGGCGGGCCCACCGATGGTTTCCGGCAACATAGGTCACGGTCACTCCCAGTTTCCGCCCCCACCTTctcttc                                                                                                                                                                                                                       >  2:198310/1‑233 (MQ=255)
        cGCAGCATCGACATCACATGGTCAATTCGTAAGGCGCCGCAGTTGTGCATATTTGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCaa                                                                                                                                                                                                                 <  1:1286677/235‑1 (MQ=255)
                             tCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCCACCCCCAGATCACGATACAAcccca                                                                                                                                                                                         >  1:1073691/1‑236 (MQ=255)
                             tCAATTCGTAATGCGCCACAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCa                                                                                                                                                                                                      <  1:753152/225‑1 (MQ=255)
                              cAATTGGTAAGGCGCCGCGGTTTTGCATTTTGGCGCGCAACAGCTCGATAACCGGTTCATAGGCCCGCGCGGTGATGATAGGCGGGTCCATGGGCGTTAATCCCCAGTTCTGCCCAAACGGGCCGGGGATATCCGGCGGCGCGCCACCCGATGCTTTCAGGCAATATAGTTCCCGGTCACCCCGGGTTCCCGCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAAt                                                                                                                                                                                        <  2:900583/238‑1 (MQ=255)
                                   cGTAAGGCGCCGCATTTTGGCATATTGGCAGGAAACAGCCCGATAACCGGTTAATGGGCACGCGGGGGGATGATGTGCGGGTCCATGGGGGGTAACCCCCAGTTCTCCCCCACCGGGCCGGGGATACCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACCCCAGGTTTCCGCCCCCCCTTCCCCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGc                                                                                                                                                                                    <  2:1380392/237‑1 (MQ=255)
                                       aTGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGAAACAAGCCAATCGGCATTTCa                                                                                                                                                                             >  1:900583/1‑239 (MQ=255)
                                       aTGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCACTTCTGCCCCAACGGGCCGAGGATATCCCGCGGGGCGCCAACCGATGCTTTCAGGCCATATAGTTCACGGTCCAACCCGGGTTCCGCCCCACCCTCCGCCACGCCAACCGCCCGACCACGCTaacaaa                                                                                                                                                                                            >  2:753152/1‑220 (MQ=255)
                                                                                              ccgcgggGTGATGATATGGGGGTCCATGGGCGGTAACCCCCAGTTCTGCCCCACCGGGCCGAGGATATCCGGCGGCGCGCCACCCGATGCTTTCAGGAAATATAGTTCAGGTTCACCCCAGTTTTCCGCCCCCCCTTCCGCTACGCAACCCGCCAGATCACGATACAAGCCAATCGGCATTTCa                                                                                                                                                                             <  2:871519/183‑1 (MQ=255)
                                                                                              aCGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCCACGGGCCGAGGAGATCCGGCGGCGCGCCAACCGATGCTTTCCGGCAATATAGTTTACGGTCACACCAGGTTTCCCCCCCCCCCTCCCCTACGCCAACCGCCAGATCACGATACAAACCCACCCGCATTTCCAAGCCCTGGCTTTTCTCCCCCCACGCGGGCAACCTCCTGGAAGCCAGCCAcacac                                                                                                                      >  2:1295652/1‑235 (MQ=255)
                                                                                               cgcgcgGTGATGATATGCGGGTCCATTGGCGGTAATCCCCCGTTCTGCCCCAACGGGGCGCGGATATCCGGCGGCCCGCCAACCGATGCTTTCAGGCAATATAGTTCCCGGTCACACCAGGGTTCCGCCCCCCCTTCCGCTACGCCAACCGCCCGATCACGATACAAGACAATCGGCCTTTCATAGACCTGGCTTATCTCCCAGAAGGCGGCAAAATGGCTGTGCGCCAGCCACTGCa                                                                                                                      >  2:1122517/1‑238 (MQ=255)
                                                                                               cgcgcgGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGGTTTGCCCCAAAGGGCCGGGGGTATCCGGCGGGGCGCCAACCGATGCTTTCAGGCAATATAGGTTACGGTCACACCAGGTTTTCGCCCCCCCTTCCGCCACGCCAACCGCCAGATCACGATAAAAGCCAAACGGCATTTCCATACCCTGGCCTATTCCCCCGCCGGCCGCAACATGGCAGTAagccagcc                                                                                                                            >  2:1235227/1‑232 (MQ=255)
                                                                                                gcgcgGTGATGATGGGCGGGCCCATGGGCGGTAATCCCCAGTCCCGCCCCAACGGGCCGGGTATATCCGGCGGCGCGCCAACCGATGCTTCCAGGCAATATAGTTCAGGGTCACACCAGGTTCCCGCCCCCCCTTCCGCTACGCCACCCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAAc                                                                                                                    <  2:529049/239‑1 (MQ=255)
                                                                                                      tgatgaTATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCCCACCAGGTTTCCGCCCCACCTTCCGCTACGCCa                                                                                                                                                                                                                  <  1:43312/140‑1 (MQ=255)
                                                                                                      tgatgaTATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCCCCTTCCGCTACGCCAACCGCCa                                                                                                                                                                                                           <  1:721274/147‑1 (MQ=255)
                                                                                                        atgatATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGGTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCCGGc                                                                                                                                                                   >  1:871519/1‑184 (MQ=255)
                                                                                                            taggCGGGTCCATTGGGGGTAATCCCCAGTTCTGCCCCAACGGGCCGGGGATACCCGGCGGCGCGCCACCCGAGGCTTTCAGGCAATATAGTTCACGGTCACCCCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCa                                                                                                                                                                             <  1:1398300/167‑1 (MQ=255)
                                                                                                                cGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGGGGATATCCGGCGGCGCGCCCACCGATGCTTTCAGGCCATATAGTTCAAGGGCACAACAGGTTTTCGCCCCACCTTCCGCCACCCCAACCGCCAGATCACGATa                                                                                                                                                                                                 >  2:721274/1‑147 (MQ=255)
                                                                                                                 gggTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACCCCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGAc                                                                                                   <  1:1295652/239‑1 (MQ=255)
                                                                                                                      cATTGGCGGTAATCCCCCGTTCTGCCCCCACGGGGCGAGGGAAACCCGCGGCGCGCCCACCGATGCTTTTAGGCAATATAGTTCACGGTCACACCAGGTTTTCGCCCCACCTTCCGCTACCCCCACCGCCCGATCACGATACCAGCCAAACCGCATTTCATAGCCCCGGc                                                                                                                                                                   >  2:1398300/1‑170 (MQ=255)
                                                                                                                       aatgggcgttaTCCCCAGGTCCTCCCCAACCGGCCGGAGAATCTCCGGGGGGGCGCCAACCGTTGCTTTCAGGCAATATAGTTCACGGTCACCCCGGGTTCCCCCCCCACCTCCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGtca                                                                                               <  2:1140044/228‑1 (MQ=255)
                                                                                                                          ggcggaaccccccATTTCGCCCCAAACGGCCCGAGGATACCGGGGGGCGCCCCACCCGATTCTTACAGGAAATATAGTTCAGGGCAACACCAGGTTCCCCCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACg                                                                                          <  2:163543/230‑1 (MQ=255)
                                                                                                                                    cctattccgcccccAGCGGCCGGGGGCTATCGGGGGGCGCGCCAACCGAGCCTTCCAGGCAATAAAGTACCCGGTCACCCCGGGTTCCCGCCCCCCCTCCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCg                                                                                <  2:1073691/229‑1 (MQ=255)
                                                                                                                                      cAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACAAGAGATAAAAATCGACGTCATACCGATGgtcttaccc                                                                              >  2:693753/1‑235 (MQ=255)
                                                                                                                                         ttCTGCCCCAACGGGCCGAGGAAATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGAAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTAGCAGa                                                                            >  1:816630/1‑238 (MQ=255)
                                                                                                                                          tCTGCCCCAACGGCCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCCCGGTCACACCAGGTTTCCGCCCCACCTCCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAAc                                                                          <  1:1284504/239‑1 (MQ=255)
                                                                                                                                            tGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAAc                                                                                                                                                                                                                <  1:497108/104‑1 (MQ=255)
                                                                                                                                            tGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAAc                                                                                                                                                                                                                >  2:497108/1‑104 (MQ=255)
                                                                                                                                            tGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATc                                                                                                      >  2:297895/1‑209 (MQ=255)
                                                                                                                                            tGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATc                                                                                                      <  1:297895/209‑1 (MQ=255)
                                                                                                                                            tGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGa                                                                            >  2:79654/1‑235 (MQ=255)
                                                                                                                                              ccccAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCa                                                                                                                                                                                                           >  1:525530/1‑107 (MQ=255)
                                                                                                                                              ccccAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCa                                                                                                                                                                                                           <  2:525530/107‑1 (MQ=255)
                                                                                                                                              ccccAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGAc                                                                      >  2:1178635/1‑239 (MQ=255)
                                                                                                                                               cccAAGGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCCCCCCAGGTTTCCGCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACg                                                                     <  1:797907/239‑1 (MQ=255)
                                                                                                                                               ccaaccGGGCCGAGGTTCTCCGGCGGCGCGCCACCCGATGCTTTCAGGCAATATATTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACg                                                                     <  1:1094488/234‑1 (MQ=255)
                                                                                                                                                   aCGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACCCCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCa                                                                                                                                                                             <  1:473499/131‑1 (MQ=255)
                                                                                                                                                   aCGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCa                                                                                                                                                                             >  2:473499/1‑131 (MQ=255)
                                                                                                                                                   aCGGGCCGAGGATACCCGGCGGCGCGCCACCCGATGCTTTCAGGCAATATAGTCCACGGTCCCCCCAGGTTTCCGCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGc                                                                                                                                                                                    <  2:329080/125‑1 (MQ=255)
                                                                                                                                                        ccGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTAt                                                                                                                                                               >  2:1230596/1‑140 (MQ=255)
                                                                                                                                                        ccGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTAt                                                                                                                                                               <  1:1230596/140‑1 (MQ=255)
                                                                                                                                                              gATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGc                                                                                                                                                                         >  1:329080/1‑125 (MQ=255)
                                                                                                                                                                tatCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATcc                                                     <  1:79654/239‑1 (MQ=255)
                                                                                                                                                                 atCCGGCGGCGCGCCACCCGATGCTTTCAGGCAATATAGTTCACGGTCACCCCAGGTTTCCGCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCa                                                    <  1:582512/239‑1 (MQ=255)
                                                                                                                                                                  tCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGGGAATCCAc                                                   >  1:1279384/1‑239 (MQ=255)
                                                                                                                                                                  tCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGAAGCACTTCTGGGGAATCCAc                                                   >  1:1027207/1‑239 (MQ=255)
                                                                                                                                                                  cccGGCGCCGCCACACCCTAGGCTTTCAGGCAAAAAAGTCAACGGCCACACCAGCTTCCCGCCCCCCCTCCCGATCCGCCACCCGCCAGATCACGATACAAGCCAATCGGCa                                                                                                                                                                                   <  2:617804/111‑1 (MQ=255)
                                                                                                                                                                     ggcggcCCCCCACCCGATGCTTTCAGGCAATAGAGTTCCCGGTCACCCCAGGTTTCCGCCCCCCCTTCCTTTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGtt                                                <  1:1127356/239‑1 (MQ=255)
                                                                                                                                                                      gcggcgCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGt                                                 >  1:805939/1‑237 (MQ=255)
                                                                                                                                                                       cggcgCGCCACCCGATGCTTTCAGGCAATATAGTTCCCGGTCCCACCAGGTTTCCCCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTc                                               <  1:693753/238‑1 (MQ=255)
                                                                                                                                                                       cggccCCCCACCCGATGCTTCCAGGAAATAATGTTCACGGTCCCCCCAGGTTTCCCCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTc                                               <  1:1122517/238‑1 (MQ=255)
                                                                                                                                                                        ggcgCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGt                                                 >  1:475866/1‑235 (MQ=255)
                                                                                                                                                                           cggtaacccGTTGTTTTGAGACAATATAGTCAACGGCACCACCGGTTCTCCGCCCCTCCTCCCTCTCCGCCACCCGCCAGATCACGATGCAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCATCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATa                                          <  2:1026467/232‑1 (MQ=255)
                                                                                                                                                                            cgcCACCCGATCATTCCAGGCATTCTAGTCCGCGGTCACTCCAGGTTCCCGCCCCACCTCCCCCTCCGCCACCCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTg                                                                                                                        <  2:556774/160‑1 (MQ=255)
                                                                                                                                                                            cgcCAACCGATGCTTTCAGGCAATATAGTTCAAGGTCCCCCCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGcc                                                                                                                                                                        >  1:617804/1‑112 (MQ=255)
                                                                                                                                                                                caccggacccttCAGGCCATTATGGTCCCCGGTCACACCGGTTTCCCCCCCCCCCTCCCCCTCCCCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACAtc                                      <  2:805939/229‑1 (MQ=255)
                                                                                                                                                                                      tGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCCGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAgaga                                                                                                              >  1:556774/1‑160 (MQ=255)
                                                                                                                                                                                       gCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACAtct                                     >  1:53612/1‑232 (MQ=255)
                                                                                                                                                                                         tttCAGGCAATATAGTCCCCGGTCACCCCAGGTTTCCCCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATg                            <  1:1018771/239‑1 (MQ=255)
                                                                                                                                                                                         ttccAGGCAATATATTTCAGGGTCCCCCCAGTTTCCCCCCCCCCCTTCCCCTACGCCAACCGCCAGATCACGATACAAGCCAAt                                                                                                                                                                                        <  2:541321/81‑1 (MQ=255)
                                                                                                                                                                                          tccAGGAAATATATTCCACGGACACACCAGGTTCCCGCCCCACCTCCCGCTCCGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGc                           <  2:1279384/237‑1 (MQ=255)
                                                                                                                                                                                                 cAATATAGTTCACGTTCACCCCAGTTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCCCTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGcc                      <  1:198310/238‑1 (MQ=255)
                                                                                                                                                                                                    tatatgtCACGGTCACCCCAGTTTACCCCCCCACCTCCCGCACCGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCa                                                    <  2:1256753/199‑1 (MQ=255)
                                                                                                                                                                                                    tataGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCa                                                                                                                                                                             >  1:541321/1‑84 (MQ=255)
                                                                                                                                                                                                         ttCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTTCGCCAACCGCCa                                                                                                                                                                                                           <  1:613867/50‑1 (MQ=255)
                                                                                                                                                                                                         ttCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCa                                                                                                                                                                                                           >  2:613867/1‑50 (MQ=255)
                                                                                                                                                                                                         ttCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGccc                                                                                                                                                                       >  2:579215/1‑85 (MQ=255)
                                                                                                                                                                                                         ttCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGccc                                                                                                                                                                       <  1:579215/85‑1 (MQ=255)
                                                                                                                                                                                                          tCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCa                                                                                                                                                                                   >  2:979621/1‑73 (MQ=255)
                                                                                                                                                                                                          tCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCa                                                                                                                                                                                   <  1:979621/73‑1 (MQ=255)
                                                                                                                                                                                                           cACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAgcgc           >  2:994825/1‑239 (MQ=255)
                                                                                                                                                                                                              ggTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATa                                          >  1:1256753/1‑205 (MQ=255)
                                                                                                                                                                                                              ggTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCAtt        >  2:530459/1‑239 (MQ=255)
                                                                                                                                                                                                              ggTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCAtt        >  1:104674/1‑239 (MQ=255)
                                                                                                                                                                                                                tCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCa                                                                                                                                                                                   <  2:804435/67‑1 (MQ=255)
                                                                                                                                                                                                                tCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCa                                                                                                                                                                                   >  1:804435/1‑67 (MQ=255)
                                                                                                                                                                                                                 cccacaAGGTTCCCCCCCCACCTCCCGTAACCCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCg     <  1:1235227/237‑1 (MQ=255)
                                                                                                                                                                                                                 cacaCCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCCCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCg     <  1:1178635/239‑1 (MQ=255)
                                                                                                                                                                                                                          tttCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCa                                                                                                                                                                                                      >  1:1376296/1‑38 (MQ=255)
                                                                                                                                                                                                                          tttCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCa                                                                                                                                                                                                      <  2:1376296/38‑1 (MQ=255)
                                                                                                                                                                                                                              cGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCGTcc  >  1:1057823/1‑229 (MQ=255)
                                                                                                                                                                                                                                   ccACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCa                                                                                                                  >  2:573833/1‑113 (MQ=255)
                                                                                                                                                                                                                                   ccACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCa                                                                                                                  <  1:573833/113‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               
ACCACAAACGCAGCATCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCGTCC  >  NC_000913/3548699‑3549147

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