Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,882,088 G→A 8.4% P59P (CCG→CCA)  dmlA → D‑malate oxidase, NAD‑dependent; putative tartrate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,882,0880GA 8.4% 136.9 / 2.4 49P59P (CCG→CCA) dmlAD‑malate oxidase, NAD‑dependent; putative tartrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (24/19);  new base A (2/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

GGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTACGGGGATTAATTCCTTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAACACCGAAGGCGAATATTCCTCGCTCGGCGGTAGAGTGAATGAAGGTACAGAGCATGAAGTCGTCATTCAGGAA  >  NC_000913/1881796‑1882388
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                            
ggTTAATTTTTCATGATTCGTTTTTTATGTTTCACCAGTTACGGGGATTAATTCCTTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTg                                                                                                                                                                                                                                                                                                      <  1:116765/301‑1 (MQ=255)
      ctattCATGATTCGTGTTTTATGTTTCACCAGTTACGGGGATTAATTCCTTTTCAGTCAATTATAGGATGGCGATGTTGTCAATTTTGATGGTCAGGAAGTGAGCACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTACTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGAGGATTCGCCTTGAGTTCTGAGCCACTGCAGTGGGCTCACCGCGACTATACCAGCCAACCAGGTAAAATGAACCAGCACCACCGGCATgc                                                                                                                                                                                                                                                                                                >  1:495694/2‑300 (MQ=255)
           cATGATTCGTGTTTTATGTTTCACCAGTTACGGGGATTAATTCCTTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGAAGACTGGCATGAGAAAc                                                                                                                                                                                                                                                                                           >  1:149999/1‑301 (MQ=255)
            ggtcgtgtgctcttccgatctaCCAGTTACGGGGATTAATTCCGTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGGGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCAGTCGTCCGGCAtcattttac                                                                                                                                                                                                                                                                                          <  2:493350/280‑9 (MQ=255)
             tGATTCGTTTTTTATGGTTCACCAGTTACGGGGATTAATTCCTTTTCAGTAAATTATAGGATGGTGATGTTGTCAATTTTGATGGGCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCAGTCGTCCGGCAtcattttacc                                                                                                                                                                                                                                                                                         <  2:516607/301‑10 (MQ=255)
                      ttttATGTTTCACCATTTACGGGGATTAATTCCTTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTc                                                                                                                                                                                                                                                                                <  2:392778/301‑1 (MQ=255)
                         tATGTTTCACCAGTTACGGGGATTAATTCCTTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGAc                                                                                                                                                                                                                                                                             >  2:26248/1‑301 (MQ=255)
                           tGTTTCACCAGTTACGGGGATTAATTCCTTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAACTGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGc                                                                                                                                                                                                                                                                           >  2:46479/1‑301 (MQ=255)
                                 aCCAGTTACGGGGATTAATTCCTTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGCAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGCGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCCTCACGGTAAAATGATGCCAGTCGTCCGGCAtcattttacaaatcgcacaaccgtcgcgta                                                                                                                                                                                                                                                                     >  1:492704/1‑272 (MQ=255)
                                            ggATTAATTCCTTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACCACTGGCATGAGCAACTTAACCGTTTCCACGCCAACCATTTTGGTAccc                                                                                                                                                                                                                                                          >  1:609786/1‑300 (MQ=255)
                                                    tCCTTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTgg                                                                                                                                                                                                                                                  <  2:150150/301‑1 (MQ=255)
                                                     ccTTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGc                                                                                                                                                                                                                                                 >  2:152183/1‑301 (MQ=255)
                                                                    ttaGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCgg                                                                                                                                                                                                                                  <  2:573702/300‑1 (MQ=255)
                                                                    ataGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCAGATGGCCGGATACCAGTCCgg                                                                                                                                                                                                                                  >  1:73366/1‑301 (MQ=255)
                                                                             tGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCAGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGcataattttaccgtgatggcagtaa                                                                                                                                                                                                                         >  2:251635/1‑277 (MQ=255)
                                                                                      cAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTggg                                                                                                                                                                                                                <  2:203686/301‑1 (MQ=255)
                                                                                      cAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTggg                                                                                                                                                                                                                >  2:120208/1‑301 (MQ=255)
                                                                                              cttcccctacacgacggtcttccgatctaaaaCGATGCGGATTGTGGCGATCCCGGGGGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCgctg                                                                                                                                                                                                        <  1:401753/273‑1 (MQ=255)
                                                                                                      ggaGTGGGAACCCAAGGGGGAAAGCGATGGGTGTTGCTGCGATCCGGGGGGACGGGATTGGAAAAGAGGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAAtt                                                                                                                                                                                                <  1:46426/299‑1 (MQ=255)
                                                                                                         ttGAGAACCCAATGATGAAAACGATGCGTATTGCGGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTtcg                                                                                                                                                                                             <  2:205347/300‑1 (MQ=255)
                                                                                                           gagaACCCCATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGGTACCGTTCCGGCACATATTTCGTTGTGGGTTTCTCTTCTTACATTtcgtc                                                                                                                                                                                           >  1:541849/1‑301 (MQ=255)
                                                                                                              ggttcagacgtgtgctcttccgattttGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTg                                                                                                                                                                                        <  2:196851/276‑1 (MQ=255)
                                                                                                                        tgaAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGAcaa                                                                                                                                                                              >  2:578173/1‑299 (MQ=255)
                                                                                                                          aaaaCGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGCCTGGCCTGAGCAACTTAGCCGTTTCGACGCCATCTCTTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCCTTATTCGTTTTTGGGTTCGCTGCTGAAATTTAGTCCTGAATTCGACCCGt                                                                                                                                                                            >  1:381291/1‑301 (MQ=255)
                                                                                                                          aaaaCGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCgctgagatcggaagagcacaagtctgaactcc                                                                                                                                                                            >  2:402267/1‑273 (MQ=255)
                                                                                                                           aaaCGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACTGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTa                                                                                                                                                                           <  2:89782/301‑1 (MQ=255)
                                                                                                                           aaaCGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCTGATACCGTTCCTGACCATCTTTCGTTGTTGGGTTCGCTGATGAAATTTCCTCGTTAATTCGAcacgtt                                                                                                                                                                           >  1:247111/1‑296 (MQ=255)
                                                                                                                                        tgctgcGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGagatcggaagagcgtcgtgtagggaa                                                                                                                                                              >  1:196608/1‑276 (MQ=255)
                                                                                                                                             ggATCCCGGGGGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGCTTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGt                                                                                                                                                         <  1:26197/300‑1 (MQ=255)
                                                                                                                                             cGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGt                                                                                                                                                         <  2:265198/301‑1 (MQ=255)
                                                                                                                                                   cGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATCATTCGTTGTGGGGTTCGCTGCTGAAATTTCGTAGTGACTTCGACCAGTACGTCAACCATCACCCTGTTCgctt                                                                                                                                                   >  1:594545/1‑298 (MQ=255)
                                                                                                                                                             gATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGc                                                                                                                                         <  2:381787/301‑1 (MQ=255)
                                                                                                                                                                          gTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTCTTTTGGTGCCGTCGGCTTGCCGGATACCGTTCCGGACCATATTTCGTTGCTGGGTTCGCTGCTGCAATTTCGTCCTGAATTCGACCAGTACGTCAACCTGACCCCGGTTCGTCACTTTCCTGGCGTTCCCTGCCCGc                                                                                                                            >  1:503255/1‑301 (MQ=255)
                                                                                                                                                                                tctGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCgg                                                                                                                      <  2:595306/300‑1 (MQ=255)
                                                                                                                                                                                  tGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGa                                                                                                                    <  2:610644/301‑1 (MQ=255)
                                                                                                                                                                                   gAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCCGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCCGGCGGGaa                                                                                                                   >  1:630177/1‑301 (MQ=255)
                                                                                                                                                                                             gcgTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTgg                                                                                                         >  2:258808/1‑301 (MQ=255)
                                                                                                                                                                                                                    cgcTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCCCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGACTGCCCACACCGATTTTTACGtaggca                                                                                  >  1:578428/1‑296 (MQ=255)
                                                                                                                                                                                                                                gtctttggTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACGGGCATGAGCAACTTAGCCGTTTCGACGCCATCTAGTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAACACCg                                                                      <  1:152032/294‑1 (MQ=255)
                                                                                                                                                                                                                                               aaaaTGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAACACCGAAGGCGAATATTCct                                                       <  1:120084/300‑1 (MQ=255)
                                                                                                                                                                                                                                                            cGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAACACCGAAGGCGAATATTCCACGATCGGCGGTag                                          >  2:558394/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                   ggAGTATTACAGCTATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAACACCGAAGGCGAATATTCCTCGCTCGGCGGTAGAGtgaat                                   <  2:391408/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAACACCGAAGGCGAATATTCCTCGCTCGGCGGTAGAGtgaaa                                   >  2:129691/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                 ttCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAACACCGAAGGCGAATATTCCTCGCTCGGCGGTAGAGTGAATGAAGGTACAGAGCa                     <  2:631054/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       gTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTACTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGACAACACCGAAGGCGAATATTACTCGCTCGGCGGTAGAGTGAATGAAGGTACAGCGCAAGAAgt               >  1:147143/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          tgctcatgCCAGTCGTCCGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTGCGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAAACCCGAAGGCGAATATTCCTCGACCGGCGGGCAAGTGAAGGAAGTTAAAAAGCATTAAatcat            >  1:31434/6‑299 (MQ=255)
                                                                                                                                                                                                                                                                                           aaTGATGCCTGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAACACCGAAGGCGAATATTCCTCGCTCGGCGGTAGAGTGAATGAAGGTACAGAGCATGAAgtcgtc           <  1:258540/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 gCCTGACGATTGGTATGATCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAACACCGAAGGCGAATATTCCTCGCTCGGCGGTAGAGTGAATGAAGGTACAGAGCATGAAGTCGTCATTCAg     <  2:73456/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   cGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAACACACCGAAGGCGAATCTTCCTCGCTCGGCGGTAGAGTGAATGAAGGTAAAGAGCATGAAGTAGTCATTCAGGa   >  1:577233/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                    ggACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAACACCGAAGGCGAATATTCCTCGCTCGGCGGTAGAGTGAATGAAGGTACAGAGCATGAAGTCGTCATTCAGGaa  >  2:80525/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                            
GGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTACGGGGATTAATTCCTTTTCAGTCAATTATAGGATGGTGATGTTGTCAATTTTGATGGTCAGGAAGTGAGAACCCAATGATGAAAACGATGCGTATTGCTGCGATCCCGGGAGACGGGATTGGCAAAGAAGTCCTTCCTGAAGGGATTCGCGTGTTACAGGCTGCCGCTGAGCGCTGGGGCTTCGCCTTGAGTTTTGAGCAAATGGAGTGGGCGAGCTGCGAGTATTACAGCCATCACGGTAAAATGATGCCGGACGACTGGCATGAGCAACTTAGCCGTTTCGACGCCATCTATTTTGGTGCCGTCGGCTGGCCGGATACCGTTCCGGACCATATTTCGTTGTGGGGTTCGCTGCTGAAATTTCGTCGTGAATTCGACCAGTACGTCAACCTGCGCCCGGTTCGTCTCTTTCCTGGCGTTCCCTGCCCGCTGGCGGGAAAACAGCCTGGCGACATCGATTTTTACGTGGTCAGGGAAAACACCGAAGGCGAATATTCCTCGCTCGGCGGTAGAGTGAATGAAGGTACAGAGCATGAAGTCGTCATTCAGGAA  >  NC_000913/1881796‑1882388

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