Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,907,573 T→G K7Q (AAA→CAA)  yobF ← DUF2527 family heat‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,907,5730TG96.1% 175.3 / ‑1.5 51K7Q (AAA→CAA) yobFDUF2527 family heat‑induced protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (29/20);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTG  >  NC_000913/1907275‑1907861
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                
aTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGAc                                                                                                                                                                                                                                                                                                <  4:410332/301‑1 (MQ=255)
aTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGAc                                                                                                                                                                                                                                                                                                >  4:2918/1‑301 (MQ=255)
             cGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCcac                                                                                                                                                                                                                                                                                    >  4:192373/1‑300 (MQ=255)
               ttCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCcacac                                                                                                                                                                                                                                                                                  >  3:355028/1‑300 (MQ=255)
                                 gTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTc                                                                                                                                                                                                                                                               >  4:94999/1‑301 (MQ=255)
                                                                            caTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACTTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGc                                                                                                                                                                                                                    >  3:403699/1‑301 (MQ=255)
                                                                              tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCaa                                                                                                                                                                                                                  >  2:55572/1‑301 (MQ=255)
                                                                                tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGc                                                                                                                                                                                                                <  2:112166/301‑1 (MQ=255)
                                                                                                ggAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGCCAAGCCATTTTTAAACATGtt                                                                                                                                                                                                >  3:18757/1‑300 (MQ=255)
                                                                                                   gTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAACCATTtctct                                                                                                                                                                                             >  3:278853/1‑300 (MQ=255)
                                                                                                     aaTGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGAt                                                                                                                                                                                           <  1:55640/301‑1 (MQ=255)
                                                                                                            ccAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGc                                                                                                                                                                                    <  4:349558/301‑1 (MQ=255)
                                                                                                                      tAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCg                                                                                                                                                                          >  3:415572/1‑301 (MQ=255)
                                                                                                                          cTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGAtcc                                                                                                                                                                      >  4:373363/1‑299 (MQ=255)
                                                                                                                                      cTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttt                                                                                                                                                          >  3:399583/1‑301 (MQ=255)
                                                                                                                                             tGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACATCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGca                                                                                                                                                   >  2:135454/1‑301 (MQ=255)
                                                                                                                                                        cTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCa                                                                                                                                         >  3:132454/1‑300 (MQ=255)
                                                                                                                                                           tGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCa                                                                                                                                                              >  4:140345/1‑276 (MQ=255)
                                                                                                                                                           tGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCa                                                                                                                                                              <  3:140344/276‑1 (MQ=255)
                                                                                                                                                             ccATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGa                                                                                                                                    >  4:220780/1‑300 (MQ=255)
                                                                                                                                                                        aTTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGat                                                                                                                        <  4:18757/301‑1 (MQ=255)
                                                                                                                                                                            cTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACCCTTATGCAACGGAATACGCGCCGatatct                                                                                                                    >  3:344619/1‑299 (MQ=255)
                                                                                                                                                                                  ttGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGAt                                                                                                                                                                        <  3:250753/243‑1 (MQ=255)
                                                                                                                                                                                  ttGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGAt                                                                                                                                                                        >  4:250757/1‑243 (MQ=255)
                                                                                                                                                                                           tcGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAAcc                                                                                                      <  4:132455/300‑1 (MQ=255)
                                                                                                                                                                                                        tAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACtttttt                                                                                                                                                       <  4:68530/238‑1 (MQ=255)
                                                                                                                                                                                                        tAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACtttttt                                                                                                                                                       >  3:68530/1‑238 (MQ=255)
                                                                                                                                                                                                          aCATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGTAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACATTGCCAGTTCAGGc                                                                                        >  3:277931/1‑299 (MQ=255)
                                                                                                                                                                                                               gATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTg                                                                                                                                                     <  3:221264/233‑1 (MQ=255)
                                                                                                                                                                                                               gATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTg                                                                                                                                                     >  4:221267/1‑233 (MQ=255)
                                                                                                                                                                                                                  aaaaCTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACAt                                                                                                                                                                                                                                     >  3:198785/1‑150 (MQ=255)
                                                                                                                                                                                                                  aaaaCTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACAt                                                                                                                                                                                                                                     <  4:198788/150‑1 (MQ=255)
                                                                                                                                                                                                                  aaaaCTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTaa                                                                              <  4:355035/301‑1 (MQ=255)
                                                                                                                                                                                                                       tGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCCGATCCTTAAACTgg                                                                         >  3:327239/1‑301 (MQ=255)
                                                                                                                                                                                                                         acacATTACTGCATGAGGAACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCg                                                                       <  3:2918/301‑1 (MQ=255)
                                                                                                                                                                                                                                    cATGAGGAACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAAc                                                            <  4:403709/301‑1 (MQ=255)
                                                                                                                                                                                                                                                     aGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTAGTTTCTAACATAGTTCAATAAGTTc                                            >  4:136205/1‑300 (MQ=255)
                                                                                                                                                                                                                                                      ggCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGa                                                                                                                         >  4:192741/1‑222 (MQ=255)
                                                                                                                                                                                                                                                      ggCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGa                                                                                                                         <  3:192739/222‑1 (MQ=255)
                                                                                                                                                                                                                                                          cGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTAATTTCTAACATAGTTCAATAAGTTCAgttgt                                      >  3:416642/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                agaAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAgg                                <  4:344625/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGcc                              >  3:310033/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                     cGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAg                            <  4:310038/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                             aacgtcaacgtGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAg                   <  4:399593/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   aacgtGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttatta             >  2:103650/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                   aacgtGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttatta             <  3:373356/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      gtgtTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACAt          <  4:415582/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          ttACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAgtt      >  3:458520/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                           tACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAgtgc     >  4:277883/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                            aCTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCCCGTTTAACCCGTAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGTACAAATCCCGATACGTCTTCCCCTTATTTTGCACACTTATGCAACGGAATACGCGCCGATCTATCATTGTGCTTAACTTTGCCTGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCACTAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAgtgcg    >  4:277534/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               cAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAgtgcgtg  >  2:33701/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                
ATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAGTGCGTG  >  NC_000913/1907275‑1907861

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