Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 327,608:1 +G 8.3% coding (1126/1473 nt) betB ← betaine aldehyde dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913327,6081.G 8.3% 185.9 / 12.2 48coding (1126/1473 nt)betBbetaine aldehyde dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base . (23/21);  new base G (2/2);  total (25/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 = 9.66e-01

CGTCATCACGCCGTTCTCGCGACCAATGCCGGAGTGTTTGTAGCCGCCAACGGGCATCTCTGCCGGGGATTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTG  >  NC_000913/327327‑327894
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                              
cGTCATCACGCCGTTCTCGCGACCAATGCCGGAGTGTTTGTAGCCGCCAACGGGCATCTCTGCCGGGGACTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGaa                                                                                                                                                                                                                                                                                          <  2:279493/288‑1 (MQ=255)
             ttCTCGCGACCAATGCCGGATTGTTTGTATCCGCCAACGGGCATCTCTGCCGGGGATTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGCTCCGCTGTCACGCTGCCCGCCGCCAGGCCGTAGACGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGAAATC‑GGTGAACACTGTCGGTGCa                                                                                                                                                                                                                                                                             >  1:183464/1‑288 (MQ=255)
                 cgggACCAAGGCCGGGGTGGTTGTAGCGGCCAACGGGCATCTCTGCCGGGGATTCGCCCCAGGTGTTGTTCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCa                                                                                                                                                                                                                                                                        <  2:240306/286‑1 (MQ=255)
                    taCCCATGCCGGGGTGTTTGTAGCCGCCAACGGGCATCTCTGCCGGGGATTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATg                                                                                                                                                                                                                                                                      <  2:295665/287‑1 (MQ=255)
                        aaTGCCGGAGTGTTTGTAGCCGCCAACGGGCATCTCTGCCGGGGATTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCg                                                                                                                                                                                                                                                                 <  1:61864/289‑1 (MQ=255)
                                      tGTAGCCGCCAACGGGCATCTCTGCCGGGGATTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCAt                                                                                                                                                                                                                                                    >  1:305483/1‑288 (MQ=255)
                                          gccgccAACGGGCATCTCTGCCGGGGATTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCAt                                                                                                                                                                                                                                                    <  2:345878/284‑1 (MQ=255)
                                          gccgccAACGGGCATCTCTGCCGGGGATTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCAt                                                                                                                                                                                                                                                    >  1:345878/1‑284 (MQ=255)
                                                          cctGCCGGGGATTCAGCCCAGGTGTTGTTCCAGCAAACACCCGCTTCCAGCTGATGAAGGACGCACTGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCGGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATcgccgc                                                                                                                                                                                                                               <  1:252948/288‑1 (MQ=255)
                                                            tGCCGGGGATTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATcgccgccg                                                                                                                                                                                                                             <  2:288158/289‑1 (MQ=255)
                                                              ccGGGGATTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGTAGCTCTCTGCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCAAAGCCATCGCCTTTCAGTACATcgccaccgc                                                                                                                                                                                                                            >  2:82131/1‑288 (MQ=255)
                                                              ccGGGGATTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATcgccgccgc                                                                                                                                                                                                                            >  2:367655/1‑288 (MQ=255)
                                                                                tgtTGATCCAGTAAATACCCTCTGCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGcc                                                                                                                                                                                                            <  1:303330/286‑1 (MQ=255)
                                                                                    gATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTC‑CGCACGATGGTCATATCTTCGCTGCAATC‑GGTGAACACTGTCGGGGCAACCCATGCGCCGTTATCGAACCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCCCGCGCCTTCCTCtt                                                                                                                                                                                                    >  2:288377/1‑289 (MQ=255)
                                                                                         aGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTg                                                                                                                                                                                                  >  1:53253/1‑287 (MQ=255)
                                                                                            aaaTACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCt                                                                                                                                                                                               <  1:367655/287‑1 (MQ=255)
                                                                                                aCCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCg                                                                                                                                                                                         >  1:318022/1‑289 (MQ=255)
                                                                                                             gagagaatgaCGCGGTGCGCGCTGTGCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGCCTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCTCTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATCGGGTGAACACTGTCCGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCgcag                                                                                                                                                                              <  1:315440/284‑1 (MQ=255)
                                                                                                                tgaatgaCGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTCTCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCCCTTTGCCTTTGGCGACATAGCGCAGCACGt                                                                                                                                                                         >  2:99126/1‑289 (MQ=255)
                                                                                                                             gcgcgcGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATCGGGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCgcag                                                                                                                                                                              >  2:378516/1‑272 (MQ=255)
                                                                                                                             gcgcgcGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATCGGGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCgcag                                                                                                                                                                              <  1:378516/272‑1 (MQ=255)
                                                                                                                             gcgcgcGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATCGGGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCCCACAGTACGCCCGCGCCTTCCTCTTTGCCTTTGGCGATCTCGCCCCGCACGGTATcccctgtc                                                                                                                                                              >  2:315440/1‑281 (MQ=255)
                                                                                                                                               cTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGcc                                                                                                                                          >  1:327472/1‑289 (MQ=255)
                                                                                                                                                           cggccgccAGGCGGTGTTCGGTATGGTTAGCGCGGCGAATGGCTGCGTCGTCGGACTCGTAGGTCAAAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCTTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTGCCTCTTTGCCTTTGGCGATATAGTGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTg                                                                                                                              <  2:307251/287‑1 (MQ=255)
                                                                                                                                                              ccgccAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCACTC‑GGTGAACACTGTCGGTGCACCCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACCACGCCGCCGCCCAGTACGCCCGCGCCTTCCCCTTTTCCTTTCGCGCTCTAGCCCCGCCCGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGCaggtcgcttcggg                                                                                                                           >  2:2352/1‑278 (MQ=255)
                                                                                                                                                                         tAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGACGTTAGTTTGCGGATCGaaaagg                                                                                                                 >  1:142630/1‑286 (MQ=255)
                                                                                                                                                                             cGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGc                                                                                                             <  2:183464/288‑1 (MQ=255)
                                                                                                                                                                                tATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCCCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTTCCTTTGGCGATATAGCGCAGCACGTTATCGCGCTGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGGAAACGGCGCCcg                                                                                                          >  2:63852/1‑288 (MQ=255)
                                                                                                                                                                                 aTCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCcgc                                                                                                         >  1:256185/1‑288 (MQ=255)
                                                                                                                                                                                                                gACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGACGTTAGTTTGCGGATCGAAAACGTCTCCCGCCCGAATGCGCCCAACGCGCCCCAGAATTTTCt                                                                         >  2:133692/1‑289 (MQ=255)
                                                                                                                                                                                                                  gtcagaGGTCAGAATGGGAATCACGGGCCCGAAGATCTCTTCAGGCAGCAGGGTCATAGCGTCGCTGCACTC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCTCGCGCGTTCCTCTTGGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTg                                                                        <  2:6409/283‑1 (MQ=255)
                                                                                                                                                                                                                          tCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATg                                                                <  2:327472/288‑1 (MQ=255)
                                                                                                                                                                                                                           cAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGGTCGAAAACGTCGCCCGCGCGAATGCCCCCAACCCGCGCCAGAATTTTCTGCTCAAATCCg                                                              >  2:390705/1‑289 (MQ=255)
                                                                                                                                                                                                                                        aCTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCACTC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATTTAGCGCAGCACGATATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGCTCGAAAACGTCGCCCGCGCGAATGCGCCCAACCCGCGCCAGAATTTTCTGCTCAACTGCTGCTTTGCCTTTc                                                  >  1:77161/1‑288 (MQ=255)
                                                                                                                                                                                                                                          tGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCCCAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGACGCTAGGTTGCGGATCGAAAACGGCGCCCGCGCGGATGCGCTCAACCGCCGCCAGAAATTTCCGCTCACACGCGGCGTTGGATTTCGc                                                >  2:152845/1‑288 (MQ=255)
                                                                                                                                                                                                                                                 aaGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGa                                        >  2:87143/1‑289 (MQ=255)
                                                                                                                                                                                                                                                           cacgcacgATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGt                              >  2:343520/1‑289 (MQ=255)
                                                                                                                                                                                                                                                             cccacgATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGc                            <  2:53253/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                acgGTGGTCATATCGTCGCTGCAATC‑TGTGAACACTGTCGGTTCAACCCAGGCGCCGTTATCGAAGCCTTCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGAGTGCCAt                         <  1:2352/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                acgATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCAt                         <  2:318022/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                acgATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCAt                         <  1:343520/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ggTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTa                    <  1:87143/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                cgCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCCCAACGCGCGCCAGCATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGAc          >  2:392039/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                cgCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGATATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGAc          <  1:287465/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                cgCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGATATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCCCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGGCCCCTTGGTACACACCTGAc          >  2:287465/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                 gCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTa                    >  2:335542/1‑277 (MQ=255)
                                                                                                                                                                                                                                                                                 gCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTa                    <  1:335542/277‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       aTC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCGGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTg  <  1:373286/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                              
CGTCATCACGCCGTTCTCGCGACCAATGCCGGAGTGTTTGTAGCCGCCAACGGGCATCTCTGCCGGGGATTCGCCCCAGGTGTTGATCCAGCAAATACCCGCTTCCAGCTGATGAATGACGCGATGCGCGCGGTTCAGGTCCGCTGTCACGATGCCCGCCGCCAGGCCGTAGTCGGTATCGTTAGCGCGGCGAATGACTTCGTCTTCCGACTCGTAGGTCAGAATGGACATCACTGGCCCGAAGATCTCTTCACGCACGATGGTCATATCGTCGCTGCAATC‑GGTGAACACTGTCGGTGCAACCCATGCGCCGTTATCGAAGCCATCGCCTTTCAGTACATCGCCGCCGCACAGTACGCGCGCGCCTTCCTCTTTGCCTTTGGCGATATAGCGCAGCACGTTATCGCGATGCGGGAAGCTGACCAGCGGGCCGAAGTTAGTTTGCGGATCGAAAACGTCGCCCGCGCGAATGCGCTCAACGCGCGCCAGAATTTTCTGCTCAAATGCGGCTTTGCATTTCGCCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTG  >  NC_000913/327327‑327894

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