Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,907,573 T→G 69.1% 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,5730TG69.1% 49.7 / 30.1 39K7Q (AAA→CAA) yobFDUF2527 family heat‑induced protein
Reads supporting (aligned to +/- strand):  ref base T (5/7);  new base G (12/15);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

CAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAG  >  NC_000913/1907285‑1907855
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                          
cAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGcc                                                                                                                                                                                                                                                                                <  2:461757/301‑1 (MQ=255)
  aCGTTCTTACCTTCGGCCAGAGTTTGTAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCAGCAGCCTGAGTCATGAAGCCAAAACCTTTAGACTCGTTGCACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATCACAGAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCcac                                                                                                                                                                                                                                                                              <  1:197579/301‑1 (MQ=255)
               cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGaaa                                                                                                                                                                                                                                                                 >  2:224417/1‑300 (MQ=255)
                 gCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTg                                                                                                                                                                                                                                                               >  2:546528/1‑301 (MQ=255)
                          tttcaGCCATTTCCCTTGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGtt                                                                                                                                                                                                                                                      <  2:249399/297‑1 (MQ=255)
                                 cATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAAccca                                                                                                                                                                                                                                               >  2:193418/1‑300 (MQ=255)
                                   ttACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCAGAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTAAAAATGCAACACATAAACAGAACTGTACCTCGTATAACCCGaa                                                                                                                                                                                                                                             >  1:356404/1‑301 (MQ=255)
                                      ggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTAGTTTAACCCGAAATc                                                                                                                                                                                                                                          <  2:467554‑M1/176‑1 (MQ=255)
                                                   agaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATAc                                                                                                                                                                                                                             <  2:108087/301‑1 (MQ=255)
                                                    ttggagttcagacgtgtgctcttccgatctaGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACa                                                                                                                                                                                                                            <  2:443880/271‑1 (MQ=255)
                                                            cGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCCTTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTGGACTGAAAATGCCAAACATAAACAGAACTGAACCTCATTTAACACAAAATATGTTATCACAAACAACGTTaaa                                                                                                                                                                                                                    >  1:555016/1‑300 (MQ=255)
                                                             gAACTCATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTATACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAAtt                                                                                                                                                                                                                   <  1:561561/301‑1 (MQ=255)
                                                                              cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTa                                                                                                                                                                                                  >  2:210008/1‑301 (MQ=255)
                                                                                  aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACaaaatcggaagaacgacgtgtagggaaaaag                                                                                                                                                                                              >  1:443316/1‑272 (MQ=255)
                                                                                    ccGGAGTAATGAAGCCAAAACC‑TTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGt                                                                                                                                                                                           >  1:49625/1‑301 (MQ=255)
                                                                                             tGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTTACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAATCAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCa                                                                                                                                                                                   <  2:555730/301‑1 (MQ=255)
                                                                                                aGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATTTCTCGATCAAAc                                                                                                                                                                                >  1:77405/1‑301 (MQ=255)
                                                                                                 gCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACg                                                                                                                                                                               <  1:545855/301‑1 (MQ=255)
                                                                                                              aaCTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGAt                                                                                                                                                                  <  1:193175/300‑1 (MQ=255)
                                                                                                                    ttGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCt                                                                                                                                                            <  2:49693/301‑1 (MQ=255)
                                                                                                                            caaacatgatctcatatcagggacttgttcgcaccttccTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCAGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCAATTTTTAAACATGTCTCGATCAAACGCAAAAATCTCCATACGTCTTCCACttt                                                                                                                                                    <  1:258322‑M1/262‑1 (MQ=255)
                                                                                                                                  ttGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGc                                                                                                                                              <  2:310108/300‑1 (MQ=255)
                                                                                                                                                              aTTCCTTTGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGat                                                                                                                  <  1:209765/301‑1 (MQ=255)
                                                                                                                                                                                                                cgacgctcttccgatctcACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGt                                                                <  1:190934/284‑1 (MQ=255)
                                                                                                                                                                                                                  caTTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGTAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTAt                                                              >  1:75245/1‑301 (MQ=255)
                                                                                                                                                                                                                   aTTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTAGCGTAtt                                                             >  1:629939/1‑301 (MQ=255)
                                                                                                                                                                                                                                 gatgctcttccgatctggCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTc                                               <  1:399129/285‑1 (MQ=255)
                                                                                                                                                                                                                                 cACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTagatcggaagagcacaa                                               >  2:191177/1‑285 (MQ=255)
                                                                                                                                                                                                                                 cACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTc                                               <  2:75335/301‑1 (MQ=255)
                                                                                                                                                                                                                                   cccacacgaccctcttccgatctgtgggaTTAAACGTCAACGTGTTTACTTAGTACTTCTTGACGGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATAACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCaa                                             <  1:336605/273‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggCAGAGAAGCGGTATTCAGCGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTACAgttgtt                               >  2:143346/1‑301 (MQ=255)
                                                                                                                                                                                                                                                 ggCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCagatcgaaagagcaca                               >  2:399643/1‑286 (MQ=255)
                                                                                                                                                                                                                                                    tcttccgatctATTCAACGTCAACGTCTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCa                            <  1:144029/291‑1 (MQ=255)
                                                                                                                                                                                                                                                          gCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGAGTATTTTCTAACATAGTTCaaagatcggaagagcacacgtcaga                      >  2:337010/1‑278 (MQ=255)
                                                                                                                                                                                                                                                             gTATTCAACGTCAACGTGTTTACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTgg                   <  2:630816/301‑1 (MQ=255)
                                                                                                                                                                                                                                                               aTTCAACGTCAACGTCTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCaagatcggaaga                 >  2:144155/1‑290 (MQ=255)
                                                                                                                                                                                                                                                                  caacgtcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCAAGTGgaagaa              >  2:192116/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                         aacgtGTTTACTCATGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGttatta       <  1:191873/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ttttACTCAGGACTTCTTGACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAg  <  2:67306/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                          
CAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAGATAAAACTGACACATTACTGCATGAGGCACCAATATAAGGCTCGGCAGAGAAGCGGTATTCAACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTTCAGTTGTTCCAGGCCAGTGGGAGAAGTTATTACATAG  >  NC_000913/1907285‑1907855

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