Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,554,0000TG39.7% 104.2 / 19.5 68Y557S (TAC→TCC) maeAmalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme
Reads supporting (aligned to +/- strand):  ref base T (9/32);  new base G (27/0);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.86e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCG‑TTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGC  >  NC_000913/1553780‑1554235
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
cgcAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCCGAAAt                                                                                                                                                                                                                            >  1:1356779/1‑239 (MQ=255)
cgcAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGGTTGCCCGaaaa                                                                                                                                                                                                                            >  2:889983/1‑238 (MQ=255)
cgcAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTCGGGTTGCCAGaaaa                                                                                                                                                                                                                            >  1:2000714/1‑238 (MQ=255)
               gTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATgg                                                                                                                                                                                                                >  1:1000433/1‑236 (MQ=255)
                 aaCGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGTGCGGGGATCGGGTTGCCCGGAAATATCGtaa                                                                                                                                                                                                                    >  2:1616876/1‑228 (MQ=255)
                 aaCGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGAAGTCGCGGGATTAGGCCTGGCAAAAATTATCGTTAATgagccgt                                                                                                                                                                                                           >  1:2337900/1‑233 (MQ=255)
                              tgatCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCAGAAAATATCCTCAAAGGCCTCTTGCAGGGCCTcc                                                                                                                                                                                               >  1:1709858/1‑237 (MQ=255)
                                   cTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGCTTGTCCGAAATTATCTTCAAAGGGCTGTTTCCGGGCGTCGGCagag                                                                                                                                                                                         >  2:59911/1‑239 (MQ=255)
                                          acaacTGCAGCTTTCATGTTTCCTCTTTTGGGGAT‑ATGTTTCTGCATAGAGGGGGGGGAAAAGGGAATTTATTCGCCACTACCCAGTATGGATGAGATCTAAAAAGGGGAGAGGAAAATACCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGGGGTCCGGCGGTAGTCGCGGTATCCGGCTTGCCAAAAATTATCGTCAATGCCCTGTTGCAGGGTTTCGGCAGAGTTTTTCa                                                                                                                                                                                  <  2:1511054/239‑1 (MQ=255)
                                             acTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATTGTCCATTGCCCTTTGCCGGGGTTCGGCAGGGGgttttaac                                                                                                                                                                                >  1:1511054/1‑234 (MQ=255)
                                                 cAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCCTGCCCGAAAATATCTGCAAAGGCCTATTGCCGGGGTGCGGCGGCGGATTTTAGCGCCAc                                                                                                                                                                           >  2:1713398/1‑239 (MQ=255)
                                                 cAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTCGGCCTGCCCGAAATTATCGGCAATGGCCTGGTGCAGGGGTTTGGCAGAGGGTTTTACCGCCAc                                                                                                                                                                           >  1:1204504/1‑239 (MQ=255)
                                                   gCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCTGCTTGCCAGAAAATATCGTCAAAGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACccccccc                                                                                                                                                                          >  1:1862242/1‑235 (MQ=255)
                                                      ttCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAAt                                                                                                                                                                                                                            >  1:2299312/1‑185 (MQ=255)
                                                                ctcctTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAATATCGTCAATGGCCTGTTTCAGGGCTTCGGCAGAGGGTTTTACCGCCCCGCCTTGCTGCTgcg                                                                                                                                                             >  1:1441072/1‑238 (MQ=255)
                                                                     tttCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCCGAAATTATCGTCAAAGGCCTGTTggtgg                                                                                                                                                                                                     >  2:653927/1‑189 (MQ=255)
                                                                                          aGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCCGGGCTTCGGCAGAGGATTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAAtt                                                                                                                                   >  1:278594/1‑239 (MQ=255)
                                                                                          aGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGGTTTTGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCgg                                                                                                                                                                                              >  2:797745/1‑180 (MQ=255)
                                                                                                                            cccAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCTGCTTTCCAGAAATTATCGTCAATGGGCTGTTTCAgg                                                                                                                                                                                                     >  2:512559/1‑139 (MQ=255)
                                                                                                                                      aTGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGAtt                                                                                                                                                                                                                                           >  1:1413766/1‑91 (MQ=255)
                                                                                                                                           aTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGCACGGCGGCGGTCGCGGTATTCGGCGTGCCAGAAATAATTGTCAACGGGCTGtt                                                                                                                                                                                                          >  2:155892/1‑119 (MQ=255)
                                                                                                                                             cTGAAAAAGGGAGAGGGAAATAGCCCGATAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGGGGGATTTgg                                                                                                                                                                                                                                        >  1:333081/1‑87 (MQ=255)
                                                                                                                                                aaaaaGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCCGAAAATATCGTCCAAGGCCCTTTGCAGGGCTTCGGCAGAGGGTTTTACCGCCACGCCCTGCTGCTGCGCCCTTTTTCCAACCGCACACGCCATTGCGCGGGGGACTTTTTGAATATCTTTTTGTTCCGGCCGTAACATACCtct                                                                                 >  2:530852/1‑233 (MQ=255)
                                                                                                                                                        agagGGAAATAGCCCGGTAGCCTTCCCTCCCGGGCGCAGGCTTAGATGGAGGTGCGGCGGTAGTCGCGGTTTTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGTTTCGGCAGAGGTTTTCACCGcc                                                                                                                                                                             <  2:1062441/135‑1 (MQ=255)
                                                                                                                                                        agagGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGcc                                                                                                                                                                                                              <  1:1300812/102‑1 (MQ=255)
                                                                                                                                                          agGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGc                                                                                                                                                                                                               <  2:1809855/99‑1 (MQ=255)
                                                                                                                                                          agGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGc                                                                                                                                                                                                               >  1:1809855/1‑99 (MQ=255)
                                                                                                                                                           gggAAATAGCCCGTTAGCCTTCACTACCGGGGGCAGGTTTAGAGGGAGTAGCGGCGGTAGTCGGGGTATCCGGCTGGCCAGAAATTATCGTCAAGGGCCTGTTGCAGGGTTTCGGCAGAGGTTTTCCCCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCG‑TTCAGCACCa                                                                  <  2:1441072/238‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTCCCGGGCGCAGGCTTAGAGGGGGGTACGGGGGTAGGCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGTTTCGGCAGAGGTTTTCACCGCCa                                                                                                                                                                            <  2:1000433/131‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCt                                                                                                                                                                                                                                      <  2:2299312/73‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCt                                                                                                                                                                                                                                      <  1:1756494/73‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGcc                                                                                                                                                                                                                                  <  1:1837419/77‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAAtt                                                                                                                                                                                                                           <  1:2328362/84‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTAt                                                                                                                                                                                                                         <  1:2073772/86‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTAt                                                                                                                                                                                                                         <  2:1356779/86‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCa                                                                                                                                                                                                                    <  1:155892/91‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATg                                                                                                                                                                                                                 <  2:958286/94‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGc                                                                                                                                                                                                               <  1:653927/96‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGc                                                                                                                                                                                                               <  1:512559/96‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGc                                                                                                                                                                                                               <  1:119063/96‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGc                                                                                                                                                                                                        <  1:797745/103‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGc                                                                                                                                                                                                        <  1:889983/103‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGTTTCGGCAGAGGtttt                                                                                                                                                                                    <  1:1095421/123‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCtt                                                                                                                                                                                                 <  1:306754/110‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCa                                                                                                                                                                                            <  1:59911/115‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCgg                                                                                                                             <  1:1262330/178‑1 (MQ=255)
                                                                                                                                                             gAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTAAAGGGCCTGTTGCAGGGTTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGgag                                                                                                                          <  1:239593/181‑1 (MQ=255)
                                                                                                                                                             aaaatgccccGGTAGCCTTCCATCCGGGGCGCAGGTTAAGAGGGGGGTACGGGGTTGGTCGGGGTATTCGGCTTGCAAAAAATTATCTACAATGCCCTGTTGCGGGGCTTCGGCAGGGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCgg                                                                                                <  2:1204504/200‑1 (MQ=255)
                                                                                                                                                             aaaaaAGCCCGGTGGCCTTCCCTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTTTTGCAGGGCTTCGGCAGAGGTTTTCCCCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCt                                                                                                                   <  1:977426/187‑1 (MQ=255)
                                                                                                                                                                  aGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTCGGGTTGCCAGAAATTATTGTCAATGGCCTGTTGCAGGGCTTCGGGAGAGGGTTTTACCGCCACGCCTtgct                                                                                                                                                                   >  1:1062441/1‑135 (MQ=255)
                                                                                                                                                                  aGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTCGGCTTGCCAGAAAATATCGTCAATGGGCTGTTGCAggg                                                                                                                                                                                                    >  2:1300812/1‑102 (MQ=255)
                                                                                                                                                                   gCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGTATTCGGCGTGGCAGAAATTATCGGCAATGGGCTGTTGCCGGGGCTCGGCAGAGGGTTTTACCCCCCCGCCCTGCTTCTGCGGCACTTTTCCCACCGCCAAAGCCATTGCCCGGGGGGCCTTTTGGATTTTTTTTCGTTCCGGGCGGAACCTACCTTCGCCG‑GTGCGCCCCACTGGTGGAt                                                         >  2:1323098/1‑239 (MQ=255)
                                                                                                                                                                         tAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTCGCTTTCCCGAAAATATCGTCAATGGCCCGTTGCAGGGCTTCGGCCGAGGTTTTTCCCCCCCCGCCCTTCTGCCGCGCCCTTTTTCCAACCCCCAACGCAATTTCGCGGGGGACCTTTTGTAAATCTTTTAATTCCGGCAGTACCCAACCCTTCCCC‑TTCAGACCCACATGTGAATACTgc                                                    >  2:2227622/1‑238 (MQ=255)
                                                                                                                                                                                cACTCCCGGGCGCAGACTTAGATGGAGGTGGGGCGGTGGTCGCGTTATTCGGCTTGCCAAAATTTATCGCCAAGGGCCTGTGGCGGGGTTTCGGCAGAGGTTTTCACCGCCa                                                                                                                                                                            <  2:432727/112‑1 (MQ=255)
                                                                                                                                                                                        ggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGc                                                                                                                                                                                                               <  2:2237641/69‑1 (MQ=255)
                                                                                                                                                                                        ggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGc                                                                                                                                                                                                               >  1:2237641/1‑69 (MQ=255)
                                                                                                                                                                                        ggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCCGAAAATATTGTCAAAGGCCTGTTGCAGGGCCTTGGCCGAGGGTTTTACCGCCCCGCCTTTCTGCTGCGCCCTTTTTCCCACCCCCAAAGCCATTTCGCGGGGGGCTTTTTGAATATTTTTCAATTTCCGCCGGACCCAACCTTCCCCGTTTCGGCCCCAA‑GGTGGATACTTCCCCACCGGTTc                                         >  2:738544/1‑234 (MQ=255)
                                                                                                                                                                                        ggCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGCTTGCCCGAAAATATATGCAATGGGCTGTTGGAGGGCTTTGGCAGAGGGTTTTACCGCCACGCCCTGCTGCTGCGCCCCTTTTCCAACCCCAAACACAAATTGGCGGGGGACTTTTTGAATATCTTTCAGTTTCGGCGGTAACAAACCCTTCCCG‑TGTCGCACCACACGTGGATAATGGGCCCGCGTGTCACtt                                     >  2:1221716/1‑238 (MQ=255)
                                                                                                                                                                                          cgcAGGCTTAGATGGAGGTACGGCGGTGGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTTTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTgcg                                                                                                                                                             <  1:1849535/117‑1 (MQ=255)
                                                                                                                                                                                          cgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTTTCGTCAATGGCCTGTTGCAGGGCCTCGGCAGAGGGTTTCACCGCCACGCCTtgctg                                                                                                                                                                  >  1:432727/1‑112 (MQ=255)
                                                                                                                                                                                           gaaGGCTTAGATGGGGGACCGGGGGTAGGCGGGGTTTTCGGTTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGTTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCG‑TTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGcc                                  <  1:530852/236‑1 (MQ=255)
                                                                                                                                                                                               gCTTAGATGGAGGTGCGGGGGTAGTCGGGGTTTTCGGTTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCG‑TTCAGCACCa                                                                  <  1:917746/202‑1 (MQ=255)
                                                                                                                                                                                                    gATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTTCCCGAAATTATCGTCAATGGCCTGTTGGAGGGCTTTGGCCGGGGGTTTTACCGCCACCCCTTGCTGCTGCGCCCTTTTTcc                                                                                                                                                   >  2:1849535/1‑117 (MQ=255)
                                                                                                                                                                                                         aGGTACGGCGGTAGTCGCGGGAATCGGGCTGCCAGGAATTATCGTCAATGGGCTGTTTCAGGGCTTCTGCAGAGGGTTTCACCCCCACCCCTTTCTGCCGCGCCACTTTTCCCACCCCCAACGCCATTGCGCGGGGGAATTTCTGGATATTTTTCAGTTTCGGCCGGACCAAACCTTTCCCC‑TTTAGTACCAAAAGTGAATa                                                        >  2:917746/1‑202 (MQ=255)
                                                                                                                                                                                                           gtccgggggTAGTCGCGGTTTTCGGCTTGCCAAAAATTATCGTCAAGGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCG‑TTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCAt                  <  1:1616876/235‑1 (MQ=255)
                                                                                                                                                                                                               ggcggTAGTCGCGGGATTTGGCCTGCCCGAAATTATCGTCCATGGCCTGGTGCAGGGGTTCGGCAGAGGGTTTTACCGCCACGCCCTGCTGCTGCGCCCTTTTGCCAACCGCAAACGCAATTTCGCCGGGGGCCTTCTGAATATCTTTCAGGTCCCGCAGTATCATAACTTCGCCC‑GTTAGCAACAATTGGTAATACTGCGCCCACGTTTTCCTTGCCGACATCAGCAGCCCACCgggg             >  2:508381/1‑237 (MQ=255)
                                                                                                                                                                                                                           ggTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCG‑TTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGAcgc  >  1:1214981/1‑238 (MQ=255)
                                                                                                                                                                                                                           ggTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCG‑TTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTACCGACATCAGCATCTCATCGGTGATACGTGAcgc  >  1:927556/1‑238 (MQ=255)
                                                                                                                                                                                                                             tATTCGGTTTGCCAAAAATTATCGTAGATGGCCTTTTGCGGGCCTGCGGCAGAGGTTTTCCCCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCG‑TTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGa                                <  1:2024283/206‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
CGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATG‑TTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCG‑TTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGC  >  NC_000913/1553780‑1554235

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