Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,474,747 C→T 6.7% intergenic (+30/‑115) bdcR → / → tabA transcriptional repressor for divergent bdcA/biofilm modulator regulated by toxins; DUF386 family protein, cupin superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,7470CT 6.7% 181.2 / 2.4 60intergenic (+30/‑115)bdcR/tabAtranscriptional repressor for divergent bdcA/biofilm modulator regulated by toxins; DUF386 family protein, cupin superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base T (2/2);  total (28/32)
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.25e-01

TGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTTCAATATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGT  >  NC_000913/4474464‑4475027
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                        
tGGTACTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTTCAATATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGtt                                                                                                                                                                                                                                                                           <  1:578348/299‑1 (MQ=255)
       tGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTTCAATATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATaa                                                                                                                                                                                                                                                                    <  1:440527/299‑1 (MQ=255)
        gAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTTCAATATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAg                                                                                                                                                                                                                                                                   >  1:527618/1‑299 (MQ=255)
          agttATGCATAGTCATGATCCACAAGCGCGTGATATTGCCGTTCAATATTAGCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGAt                                                                                                                                                                                                                                                                                <  2:308020/281‑1 (MQ=255)
                 cATAGTCATGATCCACAAGCGCGTTATATTGCCGTTCAATAGTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGGGACTGATTGTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTtgtg                                                                                                                                                                                                                                                                        <  2:83001/285‑1 (MQ=255)
                   tAGTCATGATCCACAAGCGCGTGATATTGCCGTTCAATATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTc                                                                                                                                                                                                                                                         >  1:386586/1‑298 (MQ=255)
                      tCATGATCCACAAGCGCGTGATATTGCCGTTCAATATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATCGAACAACACTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTAGTCTTGAGGAGTGATGCTGGCCTTGACCGGAAACGCGGGAAAGCGCCTGCCGAGAAGGTGTGCTAAg                                                                                                                                                                                                                                                                   >  2:40216/1‑285 (MQ=255)
                          gATCCACAAGCGCGTGATATTGCCGTTCAATATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGATa                                                                                                                                                                                                                                                 >  1:587720/1‑299 (MQ=255)
                            tCCACAAGCGCGTGATATTGCCGTTCAATATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAg                                                                                                                                                                                                                                               <  1:12074/299‑1 (MQ=255)
                                      cgTGATATTGCCGTTCAATATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACCATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCAGATACGGTGTGATAAGACAATACTTCAcgcga                                                                                                                                                                                                                                                   >  2:8410/1‑284 (MQ=255)
                                          ttatTGCCGTTCAATATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGGGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAg                                                                                                                                                                                                                                               <  2:3527/284‑1 (MQ=255)
                                                   ctCAATATTATCACGCCGCAGAAACGACCATTTATTACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCtttt                                                                                                                                                                                                                        <  1:15474/298‑1 (MQ=255)
                                                      aaTATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGAGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCAGAAAGGCGGGAAAGCGCCTGCCGAACAGTTGTGATAAGCCAATAATTCAAGAGTTAAggatcgaggattg                                                                                                                                                                                                                                   >  2:12753/1‑274 (MQ=255)
                                                       atatTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTg                                                                                                                                                                                                                                  >  2:373660/1‑285 (MQ=255)
                                                          ttatCACGCCGAAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTTATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATaa                                                                                                                                                                                                                 <  1:540900/299‑1 (MQ=255)
                                                                        aaaCGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTAAAGCGTTAAGGCTAGCGGAAATGATTATCTTTTCGTATAACGATAGAAATGaaa                                                                                                                                                                                                   >  1:180614/1‑299 (MQ=255)
                                                                                aTTTATGATTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGttt                                                                                                                                                                                           <  1:341274/299‑1 (MQ=255)
                                                                                 tttATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGtttt                                                                                                                                                                                          <  1:46755/299‑1 (MQ=255)
                                                                                  ttATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTAGGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAAt                                                                                                                                                                                                       <  2:387316/285‑1 (MQ=255)
                                                                                  gtATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTGCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTtta                                                                                                                                                                                         <  1:54878/298‑1 (MQ=255)
                                                                                            aTCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGt                                                                                                                                                                               >  1:570559/1‑299 (MQ=255)
                                                                                              cGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATTAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGttt                                                                                                                                                                                           <  2:113550/285‑1 (MQ=255)
                                                                                                         aTCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGGGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAgatgat                                                                                                                                                                   <  1:446647/298‑1 (MQ=255)
                                                                                                             acaACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCg                                                                                                                                                              <  1:476566/299‑1 (MQ=255)
                                                                                                             acaACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCg                                                                                                                                                              >  1:207607/1‑299 (MQ=255)
                                                                                                                     cacaATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAgatg                                                                                                                                                                     <  2:271292/284‑1 (MQ=255)
                                                                                                                                    aTTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTCTAACGATAGCAATGCAACGGTGTTTTAATTAAGGCGTGGAAAAGATGATCATCGGAAatatac                                                                                                                                                     >  2:9472/1‑283 (MQ=255)
                                                                                                                                    aTTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTACTTAAGGAGTGGAAAAGATGATAATCGGAAATATTCATAATCTTCAGCCg                                                                                                                                       >  1:68160/1‑299 (MQ=255)
                                                                                                                                       ttATGAGGACCGGGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAAGGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTgg                                                                                                                                    <  1:227385/299‑1 (MQ=255)
                                                                                                                                           taGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAc                                                                                                                                <  1:466500/298‑1 (MQ=255)
                                                                                                                                                    gATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGt                                                                                                                                      <  2:182702/284‑1 (MQ=255)
                                                                                                                                                     aTGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTgg                                                                                                                                    <  2:158984/285‑1 (MQ=255)
                                                                                                                                                        tCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCt                                                                                                                                  >  2:506738/1‑284 (MQ=255)
                                                                                                                                                          tGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAc                                                                                                                                >  2:431517/1‑284 (MQ=255)
                                                                                                                                                           gggCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCa                                                                                                                >  1:167369/1‑299 (MQ=255)
                                                                                                                                                             gCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTTATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAg                                                                                                               <  1:556795/298‑1 (MQ=255)
                                                                                                                                                                  cTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTAAGCCAGGCGAtt                                                                                                         >  1:244291/1‑299 (MQ=255)
                                                                                                                                                                    gCGAAGGCACGAGAGGGGCACTCTATAGAACAACGCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTa                                                                                                                     <  2:460206/285‑1 (MQ=255)
                                                                                                                                                                                gagGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAg                                                                                            <  1:9881/298‑1 (MQ=255)
                                                                                                                                                                                   gatctctctATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGcaca                                                                                        <  1:188195/294‑1 (MQ=255)
                                                                                                                                                                                      cACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGTGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATATTCAGCCGGGGCTAACCCAGGAGTTACGCAAGGCGATTGAGCat                                                                                                   >  2:228773/1‑285 (MQ=255)
                                                                                                                                                                                         tctATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATc                                                                                                >  2:165686/1‑285 (MQ=255)
                                                                                                                                                                                           tataGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGTGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGcagatcggaaga                                                                                >  1:119668/1‑289 (MQ=255)
                                                                                                                                                                                           tataGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACgg                                                                                <  1:494049/299‑1 (MQ=255)
                                                                                                                                                                                            ataGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGc                                                                               >  1:523244/1‑299 (MQ=255)
                                                                                                                                                                                             taGAACAACTCTGTGCAACAGCGGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGTGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAg                                                                                            <  2:104823/285‑1 (MQ=255)
                                                                                                                                                                                              aGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAg                                                                                            <  2:215076/284‑1 (MQ=255)
                                                                                                                                                                                                    acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATACAAACTAATCTTGAGGAGTGATGCGGGACTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACg                                                                                      <  2:445512/284‑1 (MQ=255)
                                                                                                                                                                                                          tgCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGc                                                                               >  2:63458/1‑285 (MQ=255)
                                                                                                                                                                                                            cAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCCAAGGGc                                                               >  1:116536/1‑299 (MQ=255)
                                                                                                                                                                                                                  cTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTAGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACg                                                                        >  2:185924/1‑284 (MQ=255)
                                                                                                                                                                                                                      aaTGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTCTAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAACGCACACGTTACGGCAGAAACGCCAAAGGGCAAGAACGata                                                     >  1:233270/1‑299 (MQ=255)
                                                                                                                                                                                                                                             cTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGTGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTttat                              <  1:259549/299‑1 (MQ=255)
                                                                                                                                                                                                                                             cTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTGCGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTtta                               >  1:596447/1‑298 (MQ=255)
                                                                                                                                                                                                                                                ttCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATc                                         <  2:456836/285‑1 (MQ=255)
                                                                                                                                                                                                                                                ttCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTttatctt                           >  1:500381/1‑299 (MQ=255)
                                                                                                                                                                                                                                                             aTGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAACGGGCAAGCACGATATCCAAGGCCATCGACTGTTttatct                            >  2:32239/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                        gACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTGCGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGaa                  <  2:519257/284‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ccGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATCTCGAAGGCAATCGACTGTTTTATCTTCTCTCGGAAGata              >  2:519874/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                                        aaGCGCCTTCCGATAATTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGAATATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTTGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGt  <  2:338887/284‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                        
TGGTTCTTGAAGGTATTCATAGTCATGATCCACAAGCGCGTGATATTGCCGTTCAATATTATCACGCCGCAGAAACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGT  >  NC_000913/4474464‑4475027

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