Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,553,762 C→T R180C (CGC→TGC)  hemF → coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,553,7620CT98.2% 202.8 / NA 57R180C (CGC→TGC) hemFcoproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/28);  minor base G (0/1);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGAT  >  NC_000913/2553470‑2554055
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                     
cGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTAt                                                                                                                                                                                                                                                                                                <  1:440633/300‑1 (MQ=255)
cGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTAt                                                                                                                                                                                                                                                                                                <  2:393827/300‑1 (MQ=255)
                       gCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGAcc                                                                                                                                                                                                                                                                         >  1:363518/1‑300 (MQ=255)
                       gCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCACGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGAcc                                                                                                                                                                                                                                                                         >  2:756192/1‑300 (MQ=255)
                          tttCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGCATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGa                                                                                                                                                                                                                                                                      >  2:19903/1‑300 (MQ=255)
                             cACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAaca                                                                                                                                                                                                                                                                   >  1:54134/1‑300 (MQ=255)
                             cACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTNGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAaca                                                                                                                                                                                                                                                                   >  1:564719/1‑300 (MQ=255)
                              aCTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCTAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAacac                                                                                                                                                                                                                                                                  <  1:410125/300‑1 (MQ=255)
                               cTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCAATTCTATGGTTTTGCAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCACGTTACAAAAAGTGGTGCGACGAATCCTTCTAACTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGATGTTTTTTGATGACCTGAacaca                                                                                                                                                                                                                                                                 >  2:389020/1‑299 (MQ=255)
                                 ggTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGcc                                                                                                                                                                                                                                                               >  1:796394/1‑300 (MQ=255)
                                  gTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGTGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCa                                                                                                                                                                                                                                                              >  2:435597/1‑300 (MQ=255)
                                        cATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTNTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAANACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTc                                                                                                                                                                                                                                                        >  2:638554/1‑300 (MQ=255)
                                                  aCCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGtt                                                                                                                                                                                                                                              >  2:673541/1‑300 (MQ=255)
                                                          gTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCttt                                                                                                                                                                                                                                      >  2:53221/1‑300 (MQ=255)
                                                             cccACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATg                                                                                                                                                                                                                                   >  2:472382/1‑300 (MQ=255)
                                                                    tccACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATTGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACTGCTGTTTTGCCTTTATGCAGGCgg                                                                                                                                                                                                                            <  2:73553/299‑1 (MQ=255)
                                                                               gtgCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGc                                                                                                                                                                                                                 <  1:498004/300‑1 (MQ=255)
                                                                                    gttgtttATTGCCGAAAAACCGTGTGCCGATCCCGTCTGTTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTCGGCGAAGACGTTTATCCCCTTTACAAAAAGTGGGGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCGTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTacac                                                                                                                                                                                                            <  2:77951/296‑1 (MQ=255)
                                                                                    gTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTTGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTacac                                                                                                                                                                                                            <  1:143334/300‑1 (MQ=255)
                                                                                    gTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTacac                                                                                                                                                                                                            <  1:20922/300‑1 (MQ=255)
                                                                                    attgtttgttGCCGAAAAACCGGGTGCCGATCCCGTCTGTGGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGGGCGGTATTGGCGGGCTGTTCTTTGATGATCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTacac                                                                                                                                                                                                            <  2:487234/292‑1 (MQ=255)
                                                                                                             gCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCNTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAg                                                                                                                                                                                   >  2:599962/1‑300 (MQ=255)
                                                                                                                                    cGGTGGTTTCGACGTAACCCCATTCTATAGTTTTGAAGAAGATGCTATTCACTGGCATCACACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACgg                                                                                                                                                            <  2:384962/300‑1 (MQ=255)
                                                                                                                                    cGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACgg                                                                                                                                                            <  2:213315/300‑1 (MQ=255)
                                                                                                                                    cGGTGGCTTCGACTTAACCCCATTCTATGGTTTGGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACgg                                                                                                                                                            <  2:145748/300‑1 (MQ=255)
                                                                                                                                      gTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGgcg                                                                                                                                                          <  1:780507/300‑1 (MQ=255)
                                                                                                                                               aCTTAACCCCATTCTATGGTTTTGACGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCCACGACCAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGCCAGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTACTTACCAATTGTCGAGCGAAGGAAAGCGATGGCCTAAGgcgagcgcgagc                                                                                                                                                 >  2:3304/1‑300 (MQ=255)
                                                                                                                                               aCTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTNTACCTCAAACATCGCAACGAACAGTGNGGTATTGGCGGGCTGTTCTTTGATGACCTGAACAGGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATGTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGgcgagcgcgagc                                                                                                                                                 >  2:527294/1‑300 (MQ=255)
                                                                                                                                                    aCCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAAt                                                                                                                                            >  1:164839/1‑300 (MQ=255)
                                                                                                                                                                  ttttGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTTTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATAtcg                                                                                                                              <  1:698949/300‑1 (MQ=255)
                                                                                                                                                                             aTGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGtt                                                                                                                   >  2:126138/1‑300 (MQ=255)
                                                                                                                                                                                 tATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGt                                                                                                               <  1:14680/300‑1 (MQ=255)
                                                                                                                                                                                   ttCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCg                                                                                                             <  1:12791/300‑1 (MQ=255)
                                                                                                                                                                                     cACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGGCGGAAAGCGATGGACTACGGCGAGCGCGAGCGCAATTTCCAGTCATATCGTCGCGGTCGTTCTGTCGAg                                                                                                           >  2:237775/1‑300 (MQ=255)
                                                                                                                                                                                          gCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCaa                                                                                                      >  1:106309/1‑300 (MQ=255)
                                                                                                                                                                                                  ccGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAAtctggtct                                                                                              >  2:497097/1‑300 (MQ=255)
                                                                                                                                                                                                    gCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAAtctggtctgg                                                                                            >  1:7635/1‑300 (MQ=255)
                                                                                                                                                                                                     cccGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTggg                                                                                           >  1:494172/1‑300 (MQ=255)
                                                                                                                                                                                                      ccGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGa                                                                                          <  1:267760/300‑1 (MQ=255)
                                                                                                                                                                                                       cGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGAt                                                                                         <  1:716660/300‑1 (MQ=255)
                                                                                                                                                                                                           aCCTGTGCCTGCCATTTGGCGGAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGTGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATcgcg                                                                                     <  2:730339/300‑1 (MQ=255)
                                                                                                                                                                                                           aCCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATcgcg                                                                                     <  2:730359/300‑1 (MQ=255)
                                                                                                                                                                                                                          ttGGCGAATACGTTTATCCCCGTTACAAAAAGTGGTTCTACGAAGACTTCTATCTCAAACCTCGCAACGAACACTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAATGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAATCGATGGCCTACGGCGAGCACGAGCGCAATTTCCAGTTACACCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGTTGCACGCTGTTTGGcc                                                                      <  2:8149/300‑1 (MQ=255)
                                                                                                                                                                                                                          ttGGCGAATACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGTGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCNGTAGGCAAAGGCTACACCGACGCTTANTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGcc                                                                      <  2:523806/300‑1 (MQ=255)
                                                                                                                                                                                                                                                        aGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATcc                                        <  1:373821/300‑1 (MQ=255)
                                                                                                                                                                                                                                                        aGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGCTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGNAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATNGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATcc                                        <  2:638331/300‑1 (MQ=255)
                                                                                                                                                                                                                                                         gtggtgCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCt                                       >  2:286443/1‑300 (MQ=255)
                                                                                                                                                                                                                                                             tgCGACGAATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATg                                   <  1:69383/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cgaATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAAt                              <  1:335087/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gaATACTTCTACCTCAAACATCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATg                             >  1:31679/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                               cTCAAACATCACACCGAACACTGCGGTATGGGCGTGCTGTTCGTTGATGACCGGAACACGCCAGATTGCGACCGCTGTTTGGCCTTGATGCAGGCGGTAGTCAAAGGCTACACCGACGCTTATTGACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCCGTTACATCGTCGCGGTCTTTATGTCGAGTTTAATCTGGTCTGGGGTCGCGGCACGCTGTTTGGCCTGTAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTggga                 <  2:298613/300‑3 (MQ=255)
                                                                                                                                                                                                                                                                                tCAAACATCGCAACGAACAGTGCGGTATTTGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTANACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTAc                <  1:535715/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       tCGCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTNGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGa         >  1:593074/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                         gcaccGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCACTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTAACAATTGTCGAGCGACGGAAAACGATGGCCTACGGCGAGCGCGAGCGCACTTTCCATTTATATCGTCGCGGTCGCTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGACGCACAGAGTTTCTCCTAATGTCAATACCGCCAATGATATGCTGGGAat       >  2:27747/5‑300 (MQ=255)
                                                                                                                                                                                                                                                                                         gCAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAat       >  1:423083/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                          cAACGAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAata      >  1:769986/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                              gAACAGTGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGAt  <  1:511589/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                     
CGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTAGGCAAAGGCTACACCGACGCTTATTTACCAATTGTCGAGCGACGGAAAGCGATGGCCTACGGCGAGCGCGAGCGCAATTTCCAGTTATATCGTCGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCTGTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGAT  >  NC_000913/2553470‑2554055

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