Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 508,809 G→A 7.9% intergenic (+250/+66) 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,8090GA 7.9% 140.8 / 4.3 51intergenic (+250/+66)ybaQ/copAputative DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base G (20/27);  new base A (2/2);  total (22/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 = 2.63e-01

TGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCG  >  NC_000913/508522‑509071
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                      
tGGCACGCCGTAAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCCACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGCTGCGAAGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCCGAACAGGCGTTTAcgccg                                                                                                                                                                                                                                                                        >  2:230826/1‑288 (MQ=255)
 ggCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGCGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                                                                                                                                                      >  2:121397/1‑289 (MQ=255)
    aCGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCAGGATAAGGCGTTTACGCCGCATcc                                                                                                                                                                                                                                                                   >  1:222514/1‑289 (MQ=255)
      gCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                                                                                                                                                 <  2:130635/289‑1 (MQ=255)
       acgggAAGAGCACGGAAAAAAGGTCGCGTAGGGCCAAAACATCGCATCGACTCGTAGGCCGGATAAGGGGTTCACGCCGCCTCCGGTATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCGGAACCGTAGGCGGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                                                                                                                                                 <  2:117129/284‑1 (MQ=25)
                gCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                                                                                                                                                 >  1:34639/1‑279 (MQ=255)
                  gTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTc                                                                                                                                                                                                                                                     >  1:108276/1‑289 (MQ=255)
                    gCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCGACATACGTCGa                                                                                                                                                                                                                                                   >  2:298963/1‑289 (MQ=255)
                                gcgTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATgcgc                                                                                                                                                                                                                                       <  2:52288/289‑1 (MQ=255)
                                 cgTAAGGCCAAAACTGCGGTTGAGCGAGGATGCGGGATAAGGCGGTCGCGCCGCATTCGGCAGCCGAGCACCAGTGCCTGATGCGACGCTTACCCGTCGTATCAGGCGTACAAACTTGTGCCTGAGCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCATCATTGCCGGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCGGCATCCGGCAGTCATGCGTCGATGCCAAATgcgc                                                                                                                                                                                                                                       <  2:276613/288‑1 (MQ=21)
                                                 cGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCa                                                                                                                                                                                                                      >  1:75540/1‑289 (MQ=255)
                                                      cATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCACATCCGGCAGTCATGCGTCGATGCCAAATgcgc                                                                                                                                                                                                                                       >  2:306915/1‑267 (MQ=255)
                                                      cATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCACATCCGGCAGTCATGCGTCGATGCCAAATgcgc                                                                                                                                                                                                                                       <  1:306915/267‑1 (MQ=255)
                                                       aTCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCa                                                                                                                                                                                                                <  2:310418/289‑1 (MQ=255)
                                                                gCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGt                                                                                                                                                                                                         >  2:332353/1‑287 (MQ=255)
                                                                  cGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCACATCCGGCAGTCATGCg                                                                                                                                                                                                                                                       >  2:4025/1‑239 (MQ=255)
                                                                  cGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCACATCCGGCAATCATGCg                                                                                                                                                                                                                                                       <  1:4025/239‑1 (MQ=255)
                                                                               cACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGACGCCAAATGCGCCCCCCTAAAGCAGCACAACCGCAATGACGTCCTTATTCCTTCGGttt                                                                                                                                                                                        >  2:196149/1‑289 (MQ=255)
                                                                                                cggAGCCCGATGGCCTGATGCAAGGCTTACCCGTGTTACCAGGCCTACAAACTGGTGCCTGAACCGTAGGCCGGGTAAGGCGTATACGCCGCATCCGGCATCCGGGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGt                                                                                                                                                                       <  2:207206/287‑1 (MQ=25)
                                                                                                ccgGGCCCCATTGCCTGATGCGACGCTTACCCTTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGt                                                                                                                                                                       <  2:222514/289‑1 (MQ=255)
                                                                                                          ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGCCGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATTCGGCCGTCATGCGTCGATGCAAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGAAGTACTTATTCCTTCGGTTTAAACCGCACAACCCGGTTGGCGTTAc                                                                                                                                                              >  2:281376/1‑288 (MQ=255)
                                                                                                          ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGCACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATTCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAACAACCGGTTgggcgtac                                                                                                                                                              >  2:383627/1‑282 (MQ=255)
                                                                                                               tGATGCGCCGCTTACCCGTCTTATCAGGCCTACAACCTGGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAACCCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTa                                                                                                                                                        <  1:383627/289‑1 (MQ=255)
                                                                                                                                gtgttATCAGGCCTACAAACTTGTGCCAGAACCGTGGGCCGGATAAGGCGTCTACGCCGCACCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAACCCTGGGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAgc                                                                                                                                       <  1:196149/286‑1 (MQ=255)
                                                                                                                                                    ttGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAAACGTAGGCCGTATAAGGCGACTACGCCGCATCCGGCATTCACTCATCGACGCCAAATGCGCCACACTATAGCAGCGCAACCGACATGATGTACTTATTACTTCGGTTTAAACAGCAGCAACCGGTGTGCGCTACTCACTACGGTAATCGAAGAGCGCGCCATTGTCgacacggcac                                                                                                                    >  2:248339/1‑279 (MQ=21)
                                                                                                                                                         ccTGAACGCTAGTACGTATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACTCGTCTTATCAGGCCGACAAACCTGTGCCTGAACCGTATGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGCATGCGTCGATGCCAAAGGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTAcc                                                                                                               <  1:75670/288‑1 (MQ=37)
                                                                                                                                                           tGAACCGTAGACCGGATAAGGCGTTTACGCCGCCTCCGGCATCCGATCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTTAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCggg                                                                                                            <  1:322229/289‑1 (MQ=255)
                                                                                                                                                                  tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAAGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCa                                                                                                     <  1:129821/289‑1 (MQ=255)
                                                                                                                                                                   aGGCTGGATAAGGCGTTTAGGCCGCATCCGGCATCCGGGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCa                                                                                                     <  2:278427/288‑1 (MQ=255)
                                                                                                                                                                   aGGCAGGATAAGTGGTTACAGACGTATCCGTGATCCGAGAACCAGAGCTTAATGGGCGGCTTTCGCGTCATATCAGGCCGAAGACCCAGTGCCTGAACCTTAGGCCGGCTAAGTCATTCACGCCGCATACGGCAGTCATGCGTCGATGCCAAATGCGCCATCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCa                                                                                                     <  2:136721/288‑1 (MQ=255)
                                                                                                                                                                         gATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCCCTCCGGTAATCGACGAGAGGGCCATTTCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTcc                                                                                              >  2:375669/1‑289 (MQ=255)
                                                                                                                                                                                cttttACGCCGCATCCGGCATCCGAGCACCATTGCATGATGGGAGGCTTACGCGGTTTATCAGGCCTACACACCTGTGCCTGAACCGGAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGAGCCTCCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGAAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGAAACTACCGGGTTAAGCAGTGTTCCAGTGAAc                                                                                       <  1:395478/287‑1 (MQ=38)
                                                                                                                                                                                   ttACGCCGCATCCGGCATCCTAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGc                                                                                    <  1:251398/289‑1 (MQ=255)
                                                                                                                                                                                         cgcATCCGGCATCCGGGCACCATTGCCTGATGCGACGCTTACGCGTCTTGTCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACaa                                                                               <  2:191478/288‑1 (MQ=255)
                                                                                                                                                                                                        aGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTa                                                                                                                 <  1:273390/239‑1 (MQ=255)
                                                                                                                                                                                                        aGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTa                                                                                                                 >  2:273390/1‑239 (MQ=255)
                                                                                                                                                                                                          cACCATTGCCTGATGGCACCCTTACCCGTCTCAGCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGAAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCCTCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCgg                                                             <  1:411516/289‑1 (MQ=255)
                                                                                                                                                                                                            ccATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAACGCGTTTACGCCGCATCCGGCAGTCATGAGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCCATGATGTACTTATTCCGTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGTTTAAGCAGTGTTACAGTGAACGGCCACAAAATACCGCAGGCGACCGGaa                                                           >  1:283603/1‑289 (MQ=255)
                                                                                                                                                                                                                    tGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATa                                                   <  1:332353/289‑1 (MQ=255)
                                                                                                                                                                                                                    tGAGGCGAGGCTTAAGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATa                                                   <  2:283603/289‑1 (MQ=255)
                                                                                                                                                                                                                            cGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACTAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAACATACCGGCGGCGACCGGACTACCGATACTGTTGTa                                           >  1:281081/1‑289 (MQ=255)
                                                                                                                                                                                                                            cGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTAcc                                                                                                               >  2:29522/1‑221 (MQ=255)
                                                                                                                                                                                                                            cGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTAcc                                                                                                               <  1:29522/221‑1 (MQ=255)
                                                                                                                                                                                                                            cGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTa                                           >  1:243186/1‑289 (MQ=255)
                                                                                                                                                                                                                                       cTTATCAGGCCTACAAACCTAGGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCAc                                <  1:281376/289‑1 (MQ=255)
                                                                                                                                                                                                                                               gCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAgg                        >  2:388143/1‑289 (MQ=255)
                                                                                                                                                                                                                                               gCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAgg                        >  2:215290/1‑289 (MQ=255)
                                                                                                                                                                                                                                                   aCAAACCTGTGCCGGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCt                    <  2:75540/289‑1 (MQ=255)
                                                                                                                                                                                                                                                    cAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCAc                                <  2:346016/276‑1 (MQ=255)
                                                                                                                                                                                                                                                    cAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCAc                                >  1:346016/1‑276 (MQ=255)
                                                                                                                                                                                                                                                                tGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTtgtg       <  1:121397/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cGTAGGCCGGATAAGCCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTGGGCGTTACTCACTACGGTAATCGACGAGCGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCg  <  2:108276/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                      
TGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCG  >  NC_000913/508522‑509071

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