Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,632,9550AC27.1% 114.5 / 10.7 59G372G (GGT→GGGguaBinosine 5'‑monophosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (36/7);  new base C (0/16);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 2.23e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGCCGACTTTAACCGC  >  NC_000913/2632712‑2633191
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
cAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGCCCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCgg                                                                                                                                                                                                                                           >  2:247539/1‑249 (MQ=255)
     cgtcgATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttct                                                                                                                                                                                                                                  <  1:9329/251‑1 (MQ=255)
      gtcgATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCGGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttctt                                                                                                                                                                                                                                 <  2:496272/251‑1 (MQ=255)
      gtcgATAGTACCACAGCCGGTCAGACCCATACAGGAGAGCAGGCCGCCCATCCGCTGGTGAATGATATCTTTCAGGCGACCTTTCTAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCCGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttctt                                                                                                                                                                                                                                 >  1:264864/1‑251 (MQ=255)
           tAGTACCACAGCCGGTCAGACCCATACAGGCGCGCAGGACGCCCATATGCTGGTGAATGATCTCTTTCAGGCGACATTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGt                                                                                                                                                                                                                             >  1:64199/1‑250 (MQ=255)
                 cacaGCCGGTCAGACCCATGCAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCGTTCGGGCGGCCTGTTTAGGCTACGCGACCTTCGATGCCTTCCGGCACCAGTTTGTCGGCAGCTTTATCGCTCTGGAAATAACGGTCAGAGTAACCGTTGGACGTGGCGCCCAGGGAGCCGATACCACGGTAAGATTTGGAGGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGcc                                                                                                                                                                                                                      <  1:289219/251‑1 (MQ=255)
                  acaGCCGGTCCGACCCCTACAGTAGCGCAGGCCGCCCATCTGCTGGGGACTGATATCTTTCAGGCGACATTTATAGGCTACGCGACCATCGCTCCCTTACGGCTCCAGTTTGTCGGCAGCGTTATCGCTCTGGAAGTAACGGTCAGAGGAACCTTTGGACATCGCGCCCCGGGAACCCATACCACGGTGAGATTTGTAACAACGGCCCTGTGAGATTTCGCTTTCACCCGGTGATCCTTCCGTACCCGccc                                                                                                                                                                                                                     >  2:77649/1‑250 (MQ=255)
                     gCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGc                                                                                                                                                                                                                   >  1:219927/1‑250 (MQ=255)
                       gggggaGACCCATACAGGGGCGCAGCCCCCCCATCTGCGGGGGAATGATCTTTTTCGGGGGACTTTTATAGGCTACGCGCCCTGCGATACCTCCCGGCACAAGTTGGTCGGCAGCGTTAGCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGACCCCATACCACGGTAAGTTTTGTAAGAACGGCCCGGGTAGAGTCCGATTTCACCCGGAGTTTCTTCGGTACCCGCCAGCATg                                                                                                                                                                                                                <  2:121492/246‑1 (MQ=255)
                         gTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGACTTCACCCGGAGATTCTTCAGTACCCGCCAGCATgg                                                                                                                                                                                                               >  2:527725/1‑250 (MQ=255)
                           cAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAAc                                                                                                                                                                                                            >  1:208720/1‑251 (MQ=255)
                                                              ggTGAATGTTCTCTTTCAGCCCCCTTTTATAGGTACCCCGCCTTCCGATCCCTCCCGGCCCCAGTTTGTCGCCAGCGTTATCGCTCTGGAAATAACGTGCAGAGGAACCTTGGGACATGGCGCCCGGGGAACCCATCCCACGGTAAGTTTTGTAAGACCGGCCCTTGTAGAGTTCGTTTTCCCCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAg                                                                                                                                                                           <  2:511572/249‑1 (MQ=255)
                                                                  atttatttctTTCAGGCGACCTTTATAGGCTCCGCTTCCTTCGATCCCTCCCGGCACCAGTGTGCCGGCAGCGTTATCGCTTTGGAAATACCGGGCAGTGGATCCGTTGGACACCGCTCCAAGGGAACCCATACCACGGTAAGATTTGAAAGATGGCCCCTGGGAGTTTTTGATTTCCCCCGGAGATTTTTCAGTTCCCGCTAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTg                                                                                                                                                                      <  1:68689/243‑1 (MQ=255)
                                                                     gATCTCTTTCAGGCGACCTTTATAGGCTACGCGCCCTTCGATCCCTCCCGGCACCAGTTTGTCGCCAGCTTTATCGCTCTGAAAATACCGGCCAGAGGAACCTTTGGACACCGCGCCCAGGGAACCCATCCCACGGTAAGATTTGTAAGAACGGCCCGGGTATAGTGCGATTTCCCCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCg                                                                                                                                                                   <  1:287520/250‑1 (MQ=255)
                                                                          ctTTCAGGCGACTTTTAAAGGCTACGCGACCTTCGATCCCTTCCGGCACCAGTTTGTCGGCAGCGTTACCGCTCTGGAAATAACGGTCAGAGGAACCTGTGGACATCGCGCCCAGGAAACCCATCCCACGGTAATATTTGAAGGAAGGGCCCTGGTAGTGTTCGATTTCACCGGGAGAGTCTGCATTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATg                                                                                                                                                                <  1:521703/248‑1 (MQ=255)
                                                                          ctTTCAGGCGACCTTTATAGGCTCCGCGACCTTCGTTCCCTTCCGGCACCATGTTGGCGGCAGCGTTACCGCTCTGAAAATACCGGGCAGAGTACCCTTTGAACACCCCGCCCAGTAACCCAATCCCCCGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTCCGATTCCCCCCGGAAATTCTTCAGTCCCCGCCAGCATGGAACCTTCCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGAt                                                                                                                                                                 <  1:387876/247‑1 (MQ=255)
                                                                            ttcaggcgccttttgtagtcggcgcgACCTTCGATACCTTCGGGCACCAGTTCTCCGGCAGCGTTATCCCTCGGGAAAAAACGGTCAGAGAAACCTTTGGACACTGCGCCCAGGGAACCCAGCCCCAGGTAAGATTTGTTAGACCGCCCCTGGTAGAGTCTGTTTTCCCCCGGAGATTCTCCGATCCCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTc                                                                                                                                                              <  2:439359/226‑1 (MQ=255)
                                                                             tCAGGCGCCCTTTATAGGCTACGCCACCTTCGATCCCTCCCGGCCCCATTTGGTGGACAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCACACCACGGTAAGATTTGTAAGAACGCCCCTGTTAGAGTTCGGTTTCCCCCGGAGATTCTTCGGCCCCCGCCAGCATGGAACCTACCATCACCGCGCTTGCTCCAGCGGCGATAGCTGTGGCGATGTCGcc                                                                                                                                                           <  2:564482/250‑1 (MQ=255)
                                                                              cAGGCGACATTTATAGGCTACGCGACCTTCGATACCTACCCCCACCCGTTCGTCGGCAGCGTTATCGCTCAGGAAATAACGGTCAGAGGAACCATTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATATGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGACAGCGGCGATAGCTTTGGCGATGTCGCCgg                                                                                                                                                         >  1:570084/1‑251 (MQ=255)
                                                                               tggCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGACATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTCAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCgg                                                                                                                                                         >  2:473997/2‑250 (MQ=255)
                                                                                ggCGACTTTTATAGGCTACGCCACCTTCAATCCCTTCCGGCACCAGTTGTGCGGCAGCGTTAACCCTCGGAAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCGGGGAACCAATACCCCGGAAGGATTTGTAAGACCGCCCCTGGTAGAGTTCGATTTCCCCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGgag                                                                                                                                                       <  1:438901/251‑1 (MQ=255)
                                                                                     ccTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGAGTTATCGCTCTGGAAATAACGGTCAGAGGCACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGCTTTGTAAGAACGGCCCTGGCAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATCGCTTTGGCGATGTCGCCGGAGAAGCg                                                                                                                                                  >  1:70123/1‑251 (MQ=255)
                                                                                          ataGGCTACGCGACCTTCGATACCTTCCGGCACCAGTATGTCGGCAGCGCTATCGCTCTGGAAATAACGGTCAGAGGCACCTTCGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGCACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATAc                                                                                                                                             >  1:41559/1‑251 (MQ=255)
                                                                                           aagacTACGCGACCGTCGATACCGCCCGGCACCAGTTTGCCGGCAGCGTTTTCGCTTGGGAAATAAGGGTCAGAGAACCCTTTGGACACCGCCCCCGGGGACCCCATCCCACGGTAAGTTTGGTAAGAACGGCCCGGGTAGAGTGCGATTTCCCCGGGAGATGCTTCAGTCCCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGAGAGCGTTGGCGATGTCGCCGGAGAAGCGAATAcc                                                                                                                                            <  2:312024/247‑1 (MQ=255)
                                                                                           aaGGCTACGCGCCCTTCGATCCCTTCCGGCACCAGTTTGCCGGCAGCGTTATCGCTCTGAAAATAACGGTCAGAGGACCCTTGGGACATCGCGCCCGGGGAACCCATCCCCCGGTAAGATTTGTAAGAACGGCCCTGTTAGAGTTGGTTTTCACCCGGAGATTCTTCAGTCCCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATAcc                                                                                                                                            <  1:579546/250‑1 (MQ=255)
                                                                                               ttACGCGACCTTCGATAACTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGACCCTTTGGACATCGCGCCCAGGAAACCCATACAACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGACGCGAATACCGcc                                                                                                                                         >  2:89859/2‑250 (MQ=255)
                                                                                               cTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATg                                                                                                                                                                >  1:440566/1‑227 (MQ=255)
                                                                                                 aCGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATc                                                                                                                                      >  2:128903/1‑251 (MQ=255)
                                                                                                    cgACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGCAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGc                                                                                                                                   >  2:406503/1‑251 (MQ=255)
                                                                                                     gACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCg                                                                                                                                  >  2:221011/1‑251 (MQ=255)
                                                                                                           cGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAAcc                                                                                                                            >  1:349454/1‑251 (MQ=255)
                                                                                                             gATACCTTGCGGCACAAGTATGTCGGCACCGATATCGCTCTGGAAATACCGGCCAGAGGAACCTTTGGACATCGCGCCCCGGGAACCCATACCACGGTAAGATTTGTAAGAACGGACCTGGTAGAGTTCGATTTCACCCGGCGCTTCTTCAGTACCCGCCAGCATGGAACCTACCATCCCCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGCTGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAAc                                                                                                                             >  2:575758/1‑249 (MQ=255)
                                                                                                                 ccTTCCGGCACCAGTTTGCCGGGAGGTTTATCGCTCTGGAAATAACGTACATAGAAACTTTTGTACTGCGCCCCCAGGGACCCCATCCCGTGGAAAGATTTGAAAAAACGGCCCGGGTAGAGTTCGATTTCCCCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGGGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAAt                                                                                                                       <  1:43403/251‑1 (MQ=255)
                                                                                                                      cGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGCAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGCGAAGCGAATACCGCCATCCGCGATAACCGGAATACCgg                                                                                                                  >  1:429954/1‑251 (MQ=255)
                                                                                                                       ggCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCttt                                                                                                                                                                       >  2:442529/1‑197 (MQ=255)
                                                                                                                       ggCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACACCGCGCCCAGGGAACCCATCCCCCGGTAAGATTTGTAAGAACGGCCCGGGTAGAGTTCGATTTCCCCCGGAGATTCTTAAGTCCCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCttt                                                                                                                                                                       <  1:442529/197‑1 (MQ=255)
                                                                                                                                       gacgctttaccGCTCTGAAAAAACCGGTCAGAGGAACCTTGGGACATCGCGCCCAGGAAACCCATCCCACGGTAAGATTTGTAAGACCGGCCCGGGTAGATTTCGATTTCCCCCGGAGATTCTTCATTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTc                                                                                                  <  1:52531/241‑1 (MQ=255)
                                                                                                                                        cAGCGTTATCTCTCGGGAAAAAGCGGTCAGAGAACCCTTTGGACATCCCCCCCAGGGAACCAATCCCACGAAAATATTTAAAAGAACGCCCCGGGTGGAGTCCGATTCCCCCGGAGGCTTCTCCAGTCCCCGCCAGCATGGAACCTACCATCCCCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCt                                                                                                 <  2:455076/251‑1 (MQ=255)
                                                                                                                                            gTTATCGCACTGGAAATACCGGTCAGAGGCACCTTTGCACATTGCGCCTAGGGAACCCCTACCACGGTAAGATTCGTAAGAACGGACCTGGTCGCGTTCGATTTCACCCGGAGATTCCTCAGTACCCGCCCGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTg                                                                                             >  2:547771/1‑251 (MQ=255)
                                                                                                                                               aTCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGCAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGt                                                                                          >  1:220225/1‑251 (MQ=255)
                                                                                                                                                         aaataagggtcgggggttccTTGGGACACCGCGCCGGGGGAACCCATACCGCGGAAAGATTTGAAAGAACGGCCTGGTGAGAGTCCGTTTCCACCGGGGGATTCTTCAGTCCCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATCCCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGc                                                                                <  2:8036/233‑1 (MQ=255)
                                                                                                                                                             aaCGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGt                                                                             >  1:258903/1‑250 (MQ=255)
                                                                                                                                                                        gAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGGAATCTGCGGAac                                                                 >  2:204275/1‑251 (MQ=255)
                                                                                                                                                                         acccttgggccaccgcgCCCGGGGAACCCATCCCCCGGTAAGTTTTGAAAGAACGGCCCGGGTAGAGTTCGATTTCCCCCGGAATTTCTTCAGTCCCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAaca                                                                <  1:128903/238‑1 (MQ=255)
                                                                                                                                                                                 gACACCGCGCCCAGGGACCCCACCCCAGGGTAAGTGTTTAAAGACCGGCCCTGTAAGATTTCAATTTCCCCCGGAATTTCTGCAGTCCCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGcc                                                        <  2:349454/251‑1 (MQ=255)
                                                                                                                                                                                        cgcCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCATGGTCGAGTTCGATTTCACACGGAGATCCTTCAGTACCCGCCAGCATGGAAACTACCATCACCGCGCTTGCGACAGCGGCGCTAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCCGTCACg                                                 >  1:289468/1‑251 (MQ=255)
                                                                                                                                                                                                        acccACGGTAAGTTTTGAAAGACCGCCCCGGGTAGAGTCCGTTTCCCCCGGGAGATTCTCCAGTCCCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAg                                  <  1:221011/249‑1 (MQ=255)
                                                                                                                                                                                                              ggTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATag                              >  1:474741/1‑249 (MQ=255)
                                                                                                                                                                                                                 aaGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGACTACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACCGCGGTAATCTGCGGCACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGca                          >  2:116126/1‑249 (MQ=255)
                                                                                                                                                                                                                   gATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAgg                       >  1:107991/1‑251 (MQ=255)
                                                                                                                                                                                                                          aaGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACAGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGTTAAAGCCAGGGCCAATg                >  1:317916/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaaccGGCCCGGGTAGGTTTCGTTTTCCCCGGGAGATTCTCAAGTCCCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGc               <  1:440339/247‑1 (MQ=255)
                                                                                                                                                                                                                           aGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGc               >  1:192098/1‑251 (MQ=255)
                                                                                                                                                                                                                             aaCGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGCCg             >  1:257792/1‑251 (MQ=255)
                                                                                                                                                                                                                              aCGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGCCgc            >  1:71011/1‑250 (MQ=255)
                                                                                                                                                                                                                                ggCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGCCGACt          >  1:444205/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGCCGACTTTAAcc    >  2:90173/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGGGTTGTACAGATaga                             >  2:206457/1‑225 (MQ=255)
                                                                                                                                                                                                                                        tAGAGTTCGATTTCACCCGGAGCTTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGACATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGCCGACTTTAACCGc  >  1:413144/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
CAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAGAAGCGAATACCGCCATCAGCGATAACCGGAATACCGGTGCCTTCCAGGGCTTCTACTGCGTCAGCAACAGCGGTAATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGCCGACTTTAACCGC  >  NC_000913/2632712‑2633191

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