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.

CGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGC  >  NC_000913/1553780‑1554235
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
cgcAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCCGAAAt                                                                                                                                                                                                                            >  1:1356779/1‑239 (MQ=255)
cgcAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGGTTGCCCGaaaa                                                                                                                                                                                                                            >  2:889983/1‑238 (MQ=255)
cgcAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTCGGGTTGCCAGaaaa                                                                                                                                                                                                                            >  1:2000714/1‑238 (MQ=255)
               gTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATgg                                                                                                                                                                                                                >  1:1000433/1‑236 (MQ=255)
                 aaCGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGTGCGGGGATCGGGTTGCCCGGAAATATCGtaa                                                                                                                                                                                                                    >  2:1616876/1‑228 (MQ=255)
                 aaCGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGAAGTCGCGGGATTAGGCCTGGCAAAAATTATCGTTAATgagccgt                                                                                                                                                                                                           >  1:2337900/1‑233 (MQ=255)
                              tgatCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCAGAAAATATCCTCAAAGGCCTCTTGCAGGGCCTcc                                                                                                                                                                                               >  1:1709858/1‑237 (MQ=255)
                                   cTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGCTTGTCCGAAATTATCTTCAAAGGGCTGTTTCCGGGCGTCGGCagag                                                                                                                                                                                         >  2:59911/1‑239 (MQ=255)
                                          acaacTGCAGCTTTCATGTTTCCTCTTTTGGGGATATGTTTCTGCATAGAGGGGGGGGAAAAGGGAATTTATTCGCCACTACCCAGTATGGATGAGATCTAAAAAGGGGAGAGGAAAATACCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGGGGTCCGGCGGTAGTCGCGGTATCCGGCTTGCCAAAAATTATCGTCAATGCCCTGTTGCAGGGTTTCGGCAGAGTTTTTCa                                                                                                                                                                                  <  2:1511054/239‑1 (MQ=255)
                                             acTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATTGTCCATTGCCCTTTGCCGGGGTTCGGCAGGGGgttttaac                                                                                                                                                                                >  1:1511054/1‑234 (MQ=255)
                                                 cAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCCTGCCCGAAAATATCTGCAAAGGCCTATTGCCGGGGTGCGGCGGCGGATTTTAGCGCCAc                                                                                                                                                                           >  2:1713398/1‑239 (MQ=255)
                                                 cAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTCGGCCTGCCCGAAATTATCGGCAATGGCCTGGTGCAGGGGTTTGGCAGAGGGTTTTACCGCCAc                                                                                                                                                                           >  1:1204504/1‑239 (MQ=255)
                                                   gCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCTGCTTGCCAGAAAATATCGTCAAAGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACccccccc                                                                                                                                                                          >  1:1862242/1‑235 (MQ=255)
                                                      ttCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAAt                                                                                                                                                                                                                            >  1:2299312/1‑185 (MQ=255)
                                                                ctcctTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAATATCGTCAATGGCCTGTTTCAGGGCTTCGGCAGAGGGTTTTACCGCCCCGCCTTGCTGCTgcg                                                                                                                                                             >  1:1441072/1‑238 (MQ=255)
                                                                     tttCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCCGAAATTATCGTCAAAGGCCTGTTggtgg                                                                                                                                                                                                     >  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)
                                                                                                                                                           gggAAATAGCCCGTTAGCCTTCACTACCGGGGGCAGGTTTAGAGGGAGTAGCGGCGGTAGTCGGGGTATCCGGCTGGCCAGAAATTATCGTCAAGGGCCTGTTGCAGGGTTTCGGCAGAGGTTTTCCCCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCa                                                                  <  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)
                                                                                                                                                                   gCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGTATTCGGCGTGGCAGAAATTATCGGCAATGGGCTGTTGCCGGGGCTCGGCAGAGGGTTTTACCCCCCCGCCCTGCTTCTGCGGCACTTTTCCCACCGCCAAAGCCATTGCCCGGGGGGCCTTTTGGATTTTTTTTCGTTCCGGGCGGAACCTACCTTCGCCGGTGCGCCCCACTGGTGGAt                                                         >  2:1323098/1‑239 (MQ=255)
                                                                                                                                                                         tAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTCGCTTTCCCGAAAATATCGTCAATGGCCCGTTGCAGGGCTTCGGCCGAGGTTTTTCCCCCCCCGCCCTTCTGCCGCGCCCTTTTTCCAACCCCCAACGCAATTTCGCGGGGGACCTTTTGTAAATCTTTTAATTCCGGCAGTACCCAACCCTTCCCCTTCAGACCCACATGTGAATACTgc                                                    >  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)
                                                                                                                                                                                        ggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCCGAAAATATTGTCAAAGGCCTGTTGCAGGGCCTTGGCCGAGGGTTTTACCGCCCCGCCTTTCTGCTGCGCCCTTTTTCCCACCCCCAAAGCCATTTCGCGGGGGGCTTTTTGAATATTTTTCAATTTCCGCCGGACCCAACCTTCCCCGTTTCGGCCCCAAGGTGGATACTTCCCCACCGGTTc                                         >  2:738544/1‑234 (MQ=255)
                                                                                                                                                                                        ggCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGCTTGCCCGAAAATATATGCAATGGGCTGTTGGAGGGCTTTGGCAGAGGGTTTTACCGCCACGCCCTGCTGCTGCGCCCCTTTTCCAACCCCAAACACAAATTGGCGGGGGACTTTTTGAATATCTTTCAGTTTCGGCGGTAACAAACCCTTCCCGTGTCGCACCACACGTGGATAATGGGCCCGCGTGTCACtt                                     >  2:1221716/1‑238 (MQ=255)
                                                                                                                                                                                          cgcAGGCTTAGATGGAGGTACGGCGGTGGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTTTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTgcg                                                                                                                                                             <  1:1849535/117‑1 (MQ=255)
                                                                                                                                                                                          cgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTTTCGTCAATGGCCTGTTGCAGGGCCTCGGCAGAGGGTTTCACCGCCACGCCTtgctg                                                                                                                                                                  >  1:432727/1‑112 (MQ=255)
                                                                                                                                                                                           gaaGGCTTAGATGGGGGACCGGGGGTAGGCGGGGTTTTCGGTTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGTTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGcc                                  <  1:530852/236‑1 (MQ=255)
                                                                                                                                                                                               gCTTAGATGGAGGTGCGGGGGTAGTCGGGGTTTTCGGTTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCa                                                                  <  1:917746/202‑1 (MQ=255)
                                                                                                                                                                                                    gATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTTCCCGAAATTATCGTCAATGGCCTGTTGGAGGGCTTTGGCCGGGGGTTTTACCGCCACCCCTTGCTGCTGCGCCCTTTTTcc                                                                                                                                                   >  2:1849535/1‑117 (MQ=255)
                                                                                                                                                                                                         aGGTACGGCGGTAGTCGCGGGAATCGGGCTGCCAGGAATTATCGTCAATGGGCTGTTTCAGGGCTTCTGCAGAGGGTTTCACCCCCACCCCTTTCTGCCGCGCCACTTTTCCCACCCCCAACGCCATTGCGCGGGGGAATTTCTGGATATTTTTCAGTTTCGGCCGGACCAAACCTTTCCCCTTTAGTACCAAAAGTGAATa                                                        >  2:917746/1‑202 (MQ=255)
                                                                                                                                                                                                           gtccgggggTAGTCGCGGTTTTCGGCTTGCCAAAAATTATCGTCAAGGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCAt                  <  1:1616876/235‑1 (MQ=255)
                                                                                                                                                                                                               ggcggTAGTCGCGGGATTTGGCCTGCCCGAAATTATCGTCCATGGCCTGGTGCAGGGGTTCGGCAGAGGGTTTTACCGCCACGCCCTGCTGCTGCGCCCTTTTGCCAACCGCAAACGCAATTTCGCCGGGGGCCTTCTGAATATCTTTCAGGTCCCGCAGTATCATAACTTCGCCCGTTAGCAACAATTGGTAATACTGCGCCCACGTTTTCCTTGCCGACATCAGCAGCCCACCgggg             >  2:508381/1‑237 (MQ=255)
                                                                                                                                                                                                                           ggTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGAcgc  >  1:1214981/1‑238 (MQ=255)
                                                                                                                                                                                                                           ggTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTACCGACATCAGCATCTCATCGGTGATACGTGAcgc  >  1:927556/1‑238 (MQ=255)
                                                                                                                                                                                                                             tATTCGGTTTGCCAAAAATTATCGTAGATGGCCTTTTGCGGGCCTGCGGCAGAGGTTTTCCCCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGa                                <  1:2024283/206‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
CGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTACCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCATCTCATCGGTGATACGTGACGC  >  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: