Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 508,616 A→C 5.7% intergenic (+57/+259) ybaQ → / ← copA putative DNA‑binding transcriptional regulator/copper transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913508,6160AC 5.7% 167.0 / 2.2 56intergenic (+57/+259)ybaQ/copAputative DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base C (2/2);  total (30/26)
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.74e-01

TGCATGTTCATCGTAATAGCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTAC  >  NC_000913/508318‑508915
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                           
tgcATGTTCATCGTAATAGCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCACTTCATCTCGTAGACCGGAGAAGGCGTTACCGCCcaatccggaatc                                                                                                                                                                                                                                                                                                           >  1:120926/1‑301 (MQ=255)
     tttCATCGTAATAGCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGGATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAga                                                                                                                                                                                                                                                                                                      <  1:20096/300‑2 (MQ=255)
           acccgacgctcttccgatcTGATTAATAACAATCGTAAACTCACTTCTGAGTTGGCATTTCGTCTGGCGAAAGTTTTTGATACCCCAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAgatgaagc                                                                                                                                                                                                                                                                                                <  1:432828/283‑8 (MQ=255)
            gTAATAGCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTg                                                                                                                                                                                                                                                                                               >  2:201296/1‑301 (MQ=255)
                  gCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCAACCGGCAACCGAGCACCATTGCCTGAt                                                                                                                                                                                                                                                                                         >  1:461858/1‑301 (MQ=255)
                           cACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGAAGCGACGCtt                                                                                                                                                                                                                                                                                >  2:183824/1‑301 (MQ=255)
                               gATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAgatgaagaagatcggaagagcacacgtc                                                                                                                                                                                                                                                                            >  2:433392/1‑274 (MQ=255)
                                aTCAATAACAATCGTAAACTCACTATTGGGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCg                                                                                                                                                                                                                                                                           <  1:73353/301‑1 (MQ=255)
                                   aaTAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCGTCATCTCGTAGGCCGGATAAGGCGTTCACACAGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCt                                                                                                                                                                                                                                                                        >  2:385249/1‑301 (MQ=255)
                                        cAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGCTCTTTGGGAAGTTGCAAACAACATGCGCACCCAGGAAGACTTGGGACGGCTTGCAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCACAAAAGGTCGCGTAAGGCCAACACATCCCTTCATCTCGTCAGCCGGACAAGGCTTTCACGCCCCATCCAGCCTCCGCCAACCATTGCCCGAtcagcaccctacccgccgcaacc                                                                                                                                                                                                                                                                   >  1:373367/1‑279 (MQ=255)
                                        cAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCa                                                                                                                                                                                                                                                                   >  2:442683/1‑301 (MQ=255)
                                                    cacttctttgaTGGCATTTCGTCTGGCGAAAGTTGTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACt                                                                                                                                                                                                                                                       <  1:201053/292‑1 (MQ=255)
                                                       tactGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                                                                                                                                                                                                                    >  2:639625/1‑301 (MQ=255)
                                                        tctGAGATTGCATTTCGTTTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgtg                                                                                                                                                                                                                                                   <  2:546262/300‑1 (MQ=255)
                                                          ggagaTGGCATTTCGTCTTGCGAAAGTTTTGTATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGcc                                                                                                                                                                                                                                                 <  1:488106/300‑1 (MQ=255)
                                                                  cATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCCGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCCCGCCGTGAAGAGCGTGCAACAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTCGGCCGGATAAGGCGTTCACGCCGCATCCGTCATCCGAGCCCCCTTGCCTGATGCGACGCTTACCCGTCTTATCCGGCCTACCAACTTATTCCTGAcccct                                                                                                                                                                                                                                         >  1:152706/1‑299 (MQ=255)
                                                                         tCTGTCGAAAGTTTTTGATACCACTGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCgg                                                                                                                                                                                                                                  <  2:82186/301‑1 (MQ=255)
                                                                           tGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGGTCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTATGCCTGAACCGTAGGCCAGaa                                                                                                                                                                                                                                >  2:166242/1‑300 (MQ=255)
                                                                             gCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa                                                                                                                                                                                                                              >  2:56411/1‑301 (MQ=255)
                                                                                     gaggtttgtccgttctgATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAAGATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCCGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTa                                                                                                                                                                                                                      <  1:497714/285‑1 (MQ=255)
                                                                                       ttGATACCTCAGTCGATTTTTGGCTAACCCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                                                                                                    <  2:462422/301‑1 (MQ=255)
                                                                                         gATACCACAGTCGATGTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGAGTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                                                                                  <  2:152857/301‑1 (MQ=255)
                                                                                             ccACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcaa                                                                                                                                                                                                              >  2:265907/1‑300 (MQ=255)
                                                                                              ccctgTCGATTTTTGGCTAAACCTCCGTGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcacc                                                                                                                                                                                                             <  1:183583/297‑3 (MQ=255)
                                                                                              aacaGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatc                                                                                                                                                                                                             <  1:56343/300‑1 (MQ=255)
                                                                                                   tCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccggca                                                                                                                                                                                                        >  2:290180/1‑301 (MQ=255)
                                                                                                     gATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCCGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTaagatcggaagagcaca                                                                                                                                                                                                      >  2:498361/1‑285 (MQ=255)
                                                                                                      aTTTTTGGGTAAACGTCCAGGGGGCGGTTGATCTTTGGGAAGTTGAAAACAACAGGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccggcatcc                                                                                                                                                                                                     <  1:289778/301‑1 (MQ=255)
                                                                                                               tAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACACACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCCTCCCAGcccaat                                                                                                                                                                                            >  1:447633/1‑296 (MQ=255)
                                                                                                                    cTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAAGCCGCATCCGGCATCCGAGAAACATTAccc                                                                                                                                                                                       >  2:331451/1‑300 (MQ=255)
                                                                                                                           gggcggTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGGGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGa                                                                                                                                                                                <  1:442119/300‑1 (MQ=255)
                                                                                                                                      tcttGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGt                                                                                                                                                                     <  2:77954/299‑1 (MQ=255)
                                                                                                                                         tGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCATCCGAGCCACATTGCCCGATACGACGCTTACGCGtaga                                                                                                                                                                  >  1:77864/1‑298 (MQ=255)
                                                                                                                                          gggAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTa                                                                                                                                                                 >  2:527268/1‑301 (MQ=255)
                                                                                                                                                 tGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCCGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAACAGGCGTTTACGCCGCATACGGCATCCGAGCACCAATGCCTGATCCGACGCTTAACCCTCTTAACAGccc                                                                                                                                                          >  1:49334/1‑301 (MQ=37)
                                                                                                                                                                 ctccAGGAATAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTg                                                                                                                                          <  2:121050/299‑1 (MQ=255)
                                                                                                                                                                     aGGAAGTATTGGGACGGTGTTAATCTGTGGCGGAATGTTTGGCACGCTGTGAAGATCGGGCAAAAAAGGTCGCGTAAGGCTAAAACATCGCTTCATCGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACGTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAAcc                                                                                                                                      <  1:166046/301‑1 (MQ=25)
                                                                                                                                                                             ttGGGACGTATTGAAACAGTGGCTGATTATTTGGCACGCCGTGAAGAGCGTGAAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCg                                                                                                                              <  2:49402/301‑1 (MQ=255)
                                                                                                                                                                                               gTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                            <  2:448197/301‑1 (MQ=255)
                                                                                                                                                                                                    tGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTATCCCGGATAAACCGACTACGCCCCATcc                                                                                                       >  1:110241/1‑301 (MQ=255)
                                                                                                                                                                                                         atTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCACCCCGCATCCGAGCACCATTGCCTGATGCGACGCTTCAGCCTCCTATCACGCCTACAAACCTGTGCCTGAACCCCCAGCAGGAAAAAGCGATAACGCAGCACCCGGCAg                                                                                                  >  1:353678/1‑301 (MQ=25)
                                                                                                                                                                                                            ggggaCGCCGTTAAGGGCGGGCAAAAAAGGTCGCGTAAGGCATAAAAATCGCTTCTGCTCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCATCCGAGCCCCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCa                                                                                               <  1:301499/297‑1 (MQ=25)
                                                                                                                                                                                                               cACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTTCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAACCCTGGGACTTAACCCTCAGCCGGaaaaagagatacacacgcataccggaaccacag                                                                                            >  1:412733/1‑272 (MQ=25)
                                                                                                                                                                                                                           agCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGGCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGATTACGCCGCATCCGGCaaccatccctcacaacacaa                                                                                >  1:146402/1‑282 (MQ=255)
                                                                                                                                                                                                                            gCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAGAAGGCGTTTACGCCGCATACGGCATCCGAGCACCATTGACTGATGCGACGCTTACGCGTCCTAACAGGCCCAACAACCTGTACCTGAAACGTAGACCGGATAAGGCGTTTACGCCCACACCGGCCGTAATGCGCCGATCCCacaa                                                                               >  1:353133/1‑298 (MQ=25)
                                                                                                                                                                                                                                   aaaacGGTTGCGTAAGGCCAATACGTCGAGTTATCTCGTAGGCCGGTGAAGGCGTTCACGCCACTGCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTGATCAGGCCTACAAACGTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCAc                                                                        <  1:384753/301‑1 (MQ=25)
                                                                                                                                                                                                                                        ggTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTaa                                                                   >  2:553779/1‑301 (MQ=255)
                                                                                                                                                                                                                                         gTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTaaa                                                                  >  2:477452/1‑301 (MQ=255)
                                                                                                                                                                                                                                            gcgTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCAACCCTAAaaca                                                               >  2:583843/1‑301 (MQ=255)
                                                                                                                                                                                                                                                     cAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCg                                                      <  2:146528/301‑1 (MQ=255)
                                                                                                                                                                                                                                                      aaatcgtcggttgatctgtgaGGCGGGATAAGGCGTGCAGTCCTCATCGGGCAGCCGGGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCAGCCGGCATCCGAGCACGATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGAGAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGc                                                     <  1:638748/281‑1 (MQ=37)
                                                                                                                                                                                                                                                           aTCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAatga                                                >  2:330891/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                    tctcGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCATCCGATCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTAt                                       <  2:353596/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGCCGCTTACGCGTCTTATCAGGCCTACCAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCCGCCATGCGTCGACGCCCAATGCCCCACCCTAAaccacccccacccccatcacgccccactgcc                                    >  1:243070/1‑274 (MQ=25)
                                                                                                                                                                                                                                                                          gggCCGGATAAGGCGTTCATGCCGCATCCGGCATCCGAGGACCTGTGCCTGATGCGACTCTTAGCCGTCTTATCAGGCCTACAAACTTGGGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTc                                 <  1:476851/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cGGATTAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCTTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGtt                             <  2:413296/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                cgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGCTAAGGCGTTTACGCCGCCTCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCCACTCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGACGCCCACCCGCAATCACATACTTATTCCACCAGTTTAAACCCCAACAACCgct           >  1:393405/1‑299 (MQ=25)
                                                                                                                                                                                                                                                                                                   atccggcatcggGGCACCATTGCCGGTGGCGACGTTTACCCTTCTGATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTgg        <  1:265639/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                         catccgTGCACCATTGCCTGATGCGACGTTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATGGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTAc  <  2:373863/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                           
TGCATGTTCATCGTAATAGCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTAC  >  NC_000913/508318‑508915

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