Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,548,9250AC32.1% 73.2 / 11.1 57G382G (GGT→GGGmalQ4‑alpha‑glucanotransferase
Reads supporting (aligned to +/- strand):  ref base A (23/15);  new base C (8/10);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGTCTCGCCATACGGTATCCACCACAAACGCAGCATCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCGTCCTCTTTCACTTGCTGGGCATGTA  >  NC_000913/3548677‑3549169
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
gCCGTCTCGCCATACGGTATCCACCACAAACGCAGCATCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCCACCGGGCGGGGATAGCCGCCGGCGCGCCCACCCCTGCTTTCAGGAGAAATAGTTCACGGTTCCACCAGGTTTCCCCCCCAcc                                                                                                                                                                                                                                                     >  2:106929/1‑251 (MQ=255)
       cGCCATACGGTATCCACCACAAACGCAGCATCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAAACCCCAGTTCTGCCCCAACCGGCCGCGGATATCCGGCGGCGCGCCAACCGCTGCTTTCCGGCAATATAGTTCACGGTCACATCAGGGTTCCGCCCCACCTTCCGCt                                                                                                                                                                                                                                             >  1:235684/1‑251 (MQ=255)
                     caccacAAACGCAGCACCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAACCGGTTCATAGGCCCGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAg                                                                                                                                                                                                                               <  1:209986/251‑1 (MQ=255)
                          caAACGCAGCATCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCAc                                                                                                                                                                                                                          <  1:410052/251‑1 (MQ=255)
                           aaactcagcaccgtaaacacaTGGTACATACGTAATGCGCCGAAGTTTTGCATATTGGCACGCACCAGCGCGATAAACGGTTCATATGCACGCGCGGGGATGATATGCGGGTCCATTGGCGTTAACCCCCAGTTCTGCCCCACCGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGAAATATAGTCCACGGTCACACCAGGTTGCCGCCCCACCTCCCGCTACGCCAACCGCCAGATCACg                                                                                                                                                                                                                         <  2:92474/234‑1 (MQ=255)
                                 agcaTCGACATCACATTGTCAATTCTAAAGGCCCCGCAGTTTTGCATAGTGGCACGCAACAGCTCGAAAACCGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATGGGCGGTAATCCCCAGTTCTGCCCCACCGGGCCGAGGATACCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGTTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACaa                                                                                                                                                                                                                   <  2:552780/251‑1 (MQ=255)
                                 agcaTCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCCGTTCTGCCCCAACCGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCCGGCAATATAGTTCACGGTCCCACCAGGTTTCCGCCCCCCCCTCCGCTCCGCCAACCGccccgtcccgataca                                                                                                                                                                                                                    >  2:109179/1‑237 (MQ=255)
                                 agcaTCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGGTATCCGGCGGCGCGCCAACCGATGATTTCCGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCAACTTCCGCTACGCCAACCGCCAGAACACGATACa                                                                                                                                                                                                                    >  1:448824/1‑250 (MQ=255)
                                 agcaTCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACaa                                                                                                                                                                                                                   <  1:109179/251‑1 (MQ=255)
                                  gcaTCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCCCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAg                                                                                                                                                                                                                  <  1:523722/251‑1 (MQ=255)
                                       gacacacaaTGGTAAATTGGTAATGCGCGCAGTTTTTCATAATGGGCACGAAACAGCTAGATAAACGTTTCATGGTCAGGCGCGGTGATGATTGGCGGTTCCATTGGCGGAAACCCCCAGTTCGCCCCCACGGGGCGGAGGATATCCGGGGGCGCGCCACCCGATGTTTTCAGCAAATATAGTTCACGTTCACCCCAGTTTTCCCCCCCCCCTCCCCCTCCCCCAACCGCCAGATCACGATACAAGCCAAt                                                                                                                                                                                                             <  1:279192/243‑1 (MQ=255)
                                       gACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCCCCTTCCCCTACGCCAACCCCCAGATCACGATACAAACCaa                                                                                                                                                                                                              >  1:89089/1‑250 (MQ=255)
                                           tCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGGTATCCGGCGGCGCGCCCACCGATGCTTTCAGGCAATATAGTTTACGGTCACCCCAGGTTTCCGCCCCACCTTCCGCTACCCCCACCGCCAAGTCACGATACAAGCCAATCGGc                                                                                                                                                                                                         >  1:534330/1‑251 (MQ=255)
                                            cacaTGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCCGTTCTGCCCCCACGGGCCGAGGGTATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTTCCGGTCACACCCGGGTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGAACACGATACAAACCAATCGGc                                                                                                                                                                                                         >  1:504535/1‑250 (MQ=255)
                                             acaTGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGGAATCCCCCGTTCTGCCCCAACGGGCCGGGGGTATCCGGCGGCGCGCCCACCGATGCTTTCAGGCCATATAAGTCACAGTCACACCCGGGTTCCGCCCCCCCCTCCGCTCAGCCAACCCCCAGATCACACTACAAGCCCAACGGCAt                                                                                                                                                                                                       >  2:171691/1‑251 (MQ=255)
                                                 ggTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCTATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGGAATATAATTTACGGTCACACCCGGTTTCCGCCCCACCTTCCGCGTCACCAAACGCCCGATCACGATACAAGCCAAACGGCACTTc                                                                                                                                                                                                    >  1:61028/1‑250 (MQ=255)
                                                    cAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCAAAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATa                                                                                                                                                                                                 <  1:317179/250‑1 (MQ=255)
                                                           tAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTCCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCGTAGCCCTgg                                                                                                                                                                                          <  1:171691/250‑1 (MQ=255)
                                                             aTGCCCCGCAGTTTGGCATATTGGCACGCAACAGCTCGATAAACTGGTCATAGCCACGCGCGGTGATGATAGGTGGGGCCATTGGCGGTAATCCCCAGTCCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCACCCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCTCCTTCCGCTACGCCATCCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCtt                                                                                                                                                                                       <  1:201527/251‑1 (MQ=255)
                                                               gcgcCGCATTTTTGCATATGGGCACGCAACAGCTCGATAAACGGTCCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCGGAGGATATCCGGCGGCGCGCCACCCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTAt                                                                                                                                                                                     <  2:465683/251‑1 (MQ=255)
                                                                   cccAGTTTTGCATCTTGGCACGCAACAGCTCGATAATCGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTcc                                                                                                                                                                                 <  1:500300/249‑1 (MQ=255)
                                                                     gcattttggaatattgtcccgaaaCAGCGCGATAAACGGTTAATAGGCACGCGCGGTGAGGATATGCGGGTCCATGGGCGGTATCCCCCAGTTCGGCCCCACGGGCCCGAGGATATCCGGCGGCGCGCCAACCGCTGCTTTCAGGCAATATATTTCACGGTCACACCAGGTTCCCGCCCCCCCTTCCGCTCCGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTccc                                                                                                                                                                                <  2:407295/229‑1 (MQ=255)
                                                                     cAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCCATATAGTTCACGGTCACACCAGGTTTCCTCCCCACCTTCCGCTACGCCAACCGCCAGATCACGAAACAAGCCAATAGGCATTTCCAAGCCCTGGCTTATCTccc                                                                                                                                                                                >  1:412286/1‑250 (MQ=255)
                                                                      aGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGAAATCCCCAGTTCTGCCCCAACGGGCCGGGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCCCCGTGCCAACAGGTTTCCGCCCCCCCTTCCGCTACGCCAACCCCCAGATCACGATACAAACCAATCAGCATTTCCTAGCCCTGGCTTATCTCCcac                                                                                                                                                                              >  2:73613/1‑250 (MQ=255)
                                                                           tGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCCACCGATGCTTTCCGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCCTAGCCCTGGCTTATCTCCCCGCAggc                                                                                                                                                                         >  1:21479/1‑251 (MQ=255)
                                                                            gCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCCACGGGCCGGGGGTATCCGGCGGCGCGCCCACCGCTGGTTTCAGGCCATATATTTCACGGTCACACCCGGTTTTCCCCCCCCCCTCCGCCACCCCCACCGCCACATCCCGATACAATCCCATCCGCCTTTCATATCCCTGGCTTTTCCCCcagcaaggc                                                                                                                                                                        >  2:179521/1‑247 (MQ=255)
                                                                             gcatttgggCACGAAACATCCCGATAAACGGTTCAGAGGCACCCGCGGTGAGGATATACGGGTCAATGGGCGGTAATCCCCAGTTCGGCCCCAACGGACCGAGGATATCCGGCGGCGCGCCACCCGATGCTTTCAGGAAATATAGTCCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAggcg                                                                                                                                                                        <  1:154625/244‑1 (MQ=255)
                                                                                  ttGGCACGCAACAGCTCGATAAACGGTTCATAGGCCCGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCaa                                                                                                                                                                    <  1:290333/250‑1 (MQ=255)
                                                                                     gcccgcacccgCTCGACAACCGGTTCATAGGCACGGCGGGTGAGGATATGGGGGTCCATGGGCGGTAATCCCCATTACTGCCCCAACGGGCCGGGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGAAATATAGTCCACGGTCACACCAGGTTCCCCCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTg                                                                                                                                                                <  2:496167/241‑1 (MQ=255)
                                                                                                gtCGATAAACGGTGCATAGGCACGCGCGGTGAGGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAagcca                                                                                                                                                     <  1:264643/250‑1 (MQ=255)
                                                                                                    aTAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAAACCCCAGGTCTGCCCCAAACGGGCGAGGATATCCGGCGGCGCGCCCACCGATGCTTTCAGGCAATCTAAGTCACGGGCCCACCCGGTTTTCGCCCCCCCTTCCCCTACGCCCCCCGCCACATCCCGATACAAACCCATCTGCATTTTATAACCCTGGGTTATCTCCCCGCACGCGGCCCACTAGCTGTAagaccgccc                                                                                                                                                 >  2:350091/1‑244 (MQ=255)
                                                                                                    aTAAACGGTTCATAGGCACGCGCGGTGAGGATATGCGGGTCCATTGGCGGTAACCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAgccagcca                                                                                                                                                 <  1:344032/251‑1 (MQ=255)
                                                                                                     tAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCCGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCCCACCCGGTTTCCGCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAACCCATCGGCCTTTCCTAGCCCTGGCTTATCTCCCAGCAGGCGGCCAACCGGCTGTAagccagcccc                                                                                                                                                >  2:251073/1‑249 (MQ=255)
                                                                                                        aCGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGc                                                                                                                                             >  1:296963/1‑251 (MQ=255)
                                                                                                                tAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCCACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCCGATCACGATACAAACCAATCGGCATTTCATAGCCCTGGCTTTTCTCCCAGCAGGCGGCAAACTGGCTGTAAACCAGCCCCTGCaaacca                                                                                                                                       >  1:62578/1‑245 (MQ=255)
                                                                                                                      cgcgcgGTGATGATAGGCGGGTCAATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATaaaa                                                                                                                               <  1:60310/251‑1 (MQ=255)
                                                                                                                        cgcgGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCCCCTTCCGCTCCGCCCACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCCAACTGGCTGTAAGCCAGCCACTGCCACCCGAGGTaaaaa                                                                                                                              >  1:509752/1‑250 (MQ=255)
                                                                                                                          gggTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCGGCCCCAACGGGCCGAGGATTCCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACCCCAGGTTTCCGCCCCCCCTCCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCg                                                                                                                           <  1:574732/250‑1 (MQ=255)
                                                                                                                                               aTTGGCGTTAATCCCCATTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTCCACGGTCACACCAGTTTTCCGCCCCACCTGCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGttcttc                                                                                                       <  1:452601/251‑1 (MQ=255)
                                                                                                                                                           ccccAGTTCTGCCCCAACGGGCCGAGGGTGTCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAATTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACg                                                                                           >  1:418414/1‑251 (MQ=255)
                                                                                                                                                             ccAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGCACttaccc                                                                                          >  2:351746/1‑245 (MQ=255)
                                                                                                                                                                   cggCCCAAACGGGCCGCTGATACCCGGCGGCGCGCAAACCGTGCTTTTCAGGAAATTAATTTCCCGGTCCCCCCGGTTTCCCCCCCCCCCTTCCGCAACGCAAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTg                                                                                   <  2:367627/249‑1 (MQ=255)
                                                                                                                                                                         cAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAAt                                                                             >  1:88220/1‑251 (MQ=255)
                                                                                                                                                                                ggccgggggagacccGGCGGCGCCCCAACCTATTCTTTCAGCCAATAGTGTCCTCGTCCCCACCGGTTTCCCGCCCCCCCTCCCCCGACCCCACCCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCAc                                                                         <  2:135668/238‑1 (MQ=255)
                                                                                                                                                                                      aGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTg                                                                   >  1:159819/1‑251 (MQ=255)
                                                                                                                                                                                         atatCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGAt                                                                 >  1:259773/1‑250 (MQ=255)
                                                                                                                                                                                                         cAACCGATGCTTTCAGGCAATATAGTTCACGGTCACCCCAGGTTTCCGCCCCACCTCCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCa                                                <  1:289649/251‑1 (MQ=255)
                                                                                                                                                                                                             cGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGcc                                            >  2:559809/1‑251 (MQ=255)
                                                                                                                                                                                                               aTGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAgaga                                                                                                                                    <  1:210271/161‑1 (MQ=255)
                                                                                                                                                                                                               aTGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCACCCAgaga                                                                                                                                    >  2:210271/1‑161 (MQ=255)
                                                                                                                                                                                                                  cTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGccc                                       >  1:464261/1‑251 (MQ=255)
                                                                                                                                                                                                                              tataGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCAATGCAGGCCAGCCCCAACGCCTTTc                            >  2:130560/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCGTCctc                      >  2:199702/1‑251 (MQ=255)
                                                                                                                                                                                                                                             cccctgtttcccGCCCCCCCTTCCGCATCGCCAACCTCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCGTCCTCTTTCACTTg             <  1:106929/241‑1 (MQ=255)
                                                                                                                                                                                                                                             cacCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCGTCCTCTTTCACTTg             >  1:397135/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ccAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCGTCCTCTTTCACTTGCt           >  1:379127/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCAAGCCCCAGCGCATTTCGTCCTCTTTCCCTTGCTgg         >  2:95570/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         ccccccccCTCCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCGTCCTCTTTCACTTGCTGGGCATGTa  <  1:341589/249‑1 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
GCCGTCTCGCCATACGGTATCCACCACAAACGCAGCATCGACATCACATGGTCAATTCGTAATGCGCCGCAGTTTTGCATATTGGCACGCAACAGCTCGATAAACGGTTCATAGGCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGATATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGACGTCATCACGATGTTCTTCGCAGAACTGACGCACTTCTGGTGAATCCACGTTCTGATACATCTCTGGCCATGCAGGCCAGCCCCAGCGCATTTCGTCCTCTTTCACTTGCTGGGCATGTA  >  NC_000913/3548677‑3549169

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