Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,507,249 C→T D30N (GAC→AAC)  yhfT ← uncharacterized protein YhfT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,507,2490CT83.0% 114.8 / 17.6 53D30N (GAC→AAC) yhfTuncharacterized protein YhfT
Reads supporting (aligned to +/- strand):  ref base C (3/6);  new base T (23/21);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

AAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTCATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATGCCTGGTTTGGCAATGGTGCAACTTTT  >  NC_000913/3507006‑3507491
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
aaGTAAGGATCAACACGCCCAAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGCTACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATTCCGTTATGacaa                                                                                                                                                                                                                                             >  2:426271/1‑247 (MQ=255)
    aaGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCAAAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACa                                                                                                                                                                                                                                          <  1:342485/250‑1 (MQ=255)
    aaGGATCAACACGCCCCAGATAGCGACAAGACCAAACGCCATCCGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCCGGCGGGTGCCGTTATGAAAAACa                                                                                                                                                                                                                                          >  2:71954/1‑250 (MQ=255)
     aGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGATCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGCCAGGAACGGGCGGATGCCGTTCTGAAAAAAAgc                                                                                                                                                                                                                                        >  1:477627/1‑251 (MQ=255)
           aaCACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAgccgcct                                                                                                                                                                                                                                   >  2:283834/1‑249 (MQ=255)
            acacGCCCCAGATAGCGCCAAAACAAAAATCCATCAGGCTGTTTATCGCCGGGACGCCGCGGATACCGGTAGTAAGAAAGAGTAACCAGTCGTGGAGCAGCCCGGTTTTAAGGGTAAAAGACATCCCCACCGGGACCACAAAACCAATGCTCAGACCAAAAGCGTTACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAg                                                                                                                                                                                                                                 <  2:304720/251‑1 (MQ=255)
                gCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGAt                                                                                                                                                                                                                             >  2:137335/1‑251 (MQ=255)
                  cccAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTCATGAAAAACAGCCGCCGAGCGATgc                                                                                                                                                                                                                           <  2:243303/251‑1 (MQ=255)
                                                  cTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGGAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGAACGGGCCGATGCCGGTATGAAAAACAGCCGGCGAGACATACGCCAGAAGCGATGTCATACCCGTCaagaacgc                                                                                                                                                                                           >  2:91440/1‑244 (MQ=255)
                                                      ttATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCCTACCCGTCAGGCACGccacc                                                                                                                                                                                       >  2:153977/1‑251 (MQ=255)
                                                             cAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTCATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATaa                                                                                                                                                                                <  1:413957/251‑1 (MQ=255)
                                                              aggaggCCGAGGATATCGGTAGGAAGAAAAAGTACCCATGCGTTGAGCCGCCCGGTTTTCAGGGTAAAAGAGATCCCCCCCGAGGCCACAAAACCAATGCTCAGACAAAAAGCGATCCTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTCATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAAt                                                                                                                                                                               <  2:330667/251‑1 (MQ=255)
                                                                gACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCt                                                                                                                                                                             >  1:23684/1‑251 (MQ=255)
                                                                         gATATCGGTAGGAAGAAAGAGTAACCATGCGTTAAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCTATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACa                                                                                                                                                                    <  2:161417/251‑1 (MQ=255)
                                                                            ataggTAGGAAGAAAGAGTAACCCTGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCCATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCCGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGACGCCGAGCGATGCGCCAGAAGCGATGTACTACCCGTCAGGCACGCCACCCCGATAATCTGAAtataa                                                                                                                                                                     >  2:301846/4‑246 (MQ=255)
                                                                              cGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATcc                                                                                                                                                               >  1:265436/1‑251 (MQ=255)
                                                                                  aGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTCATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCAt                                                                                                                                                             >  2:264401/1‑249 (MQ=255)
                                                                                    gaagaaAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGTTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTGCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATTGCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTc                                                                                                                                                         <  2:483428/251‑1 (MQ=255)
                                                                                        aaagaGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTTTGCAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATCACCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCNTGc                                                                                                                                                     >  2:438934/1‑251 (MQ=255)
                                                                                                 ccATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTAttt                                                                                                                                            >  1:491552/1‑251 (MQ=255)
                                                                                                       gTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCCTCGTCGTGCCTTTTATTTCAGATg                                                                                                                                      >  2:101791/1‑251 (MQ=255)
                                                                                                          gAGCAGCCCGTTTTTCAGGTAAAAAGAGTTCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGtt                                                                                                                                   <  2:13058/251‑1 (MQ=255)
                                                                                                               gcCCGGTTTTCAGGGTAAAAGAGATCCCCCCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAg                                                                                                                             <  1:426271/251‑1 (MQ=255)
                                                                                                               ccccGGTTTGAAGGGAAAAAGAGATCCCCCCCGAGGCCACAAACCCAATGCTCAGACCAAAAGCGATCCCCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCa                                                                                                                              <  2:248935/250‑1 (MQ=255)
                                                                                                                cccGGTTTTCAGGTGAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAg                                                                                                                             <  1:377726/251‑1 (MQ=255)
                                                                                                                  cGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAg                                                                                                                                                                   >  1:315964/1‑211 (MQ=255)
                                                                                                                   ggTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGCCCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCAt                                                                                                                          <  1:440825/251‑1 (MQ=255)
                                                                                                                      tttCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCg                                                                                                                        <  1:301846/250‑1 (MQ=255)
                                                                                                                                                 ggCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGCTCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCCCGATAATCTGAATATACAGATCCATAGCCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCCACGTGCTCAACgg                                                                                             >  2:409648/1‑250 (MQ=255)
                                                                                                                                                           ccAATGCTCAGACCAAAAGCGATACCCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAAc                                                                                  <  1:368361/251‑1 (MQ=255)
                                                                                                                                                             aaTGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGc                                                                                >  1:200106/1‑251 (MQ=255)
                                                                                                                                                                gCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTCATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCAc                                                                             <  1:450650/251‑1 (MQ=255)
                                                                                                                                                                     aacaaaaaGCGTTCCTCCCCGCCTCGCGACGGTTCACATAGCCGTCAATCAGTTGGGGCAGGATCGGGTGGAGTCCGTTATGAAAAACAGCCGCCGAGAGATGCGCCAGAAGTGCTGTCATACCCTGCAGGCACGCCACCACGATAATCTGAATATACAGATCCATATGCGTGCCTGTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACTTGCGCAACGGAAAGGCCAAACGCCACtcggc                                                                        <  2:95820/247‑5 (MQ=255)
                                                                                                                                                                      aCCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTcc                                                                       <  1:283834/251‑1 (MQ=255)
                                                                                                                                                                        cAAAAGAGATACTCCCCGCCTCGCGACGGGTAATATAGCCGTAAATCAGTTGCGGCAGGAGCGGGCGGATGCCGTTATGAAAAACAGCCGCCTAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCAGAGTCGTGCCTTTTATTTCAGAGGGTTAATCAGCATCGGAATCGCGTGTTCGACGGGCTCAACGGAAAGGCCAAACGACACTTTTCCtt                                                                     <  2:401211/251‑1 (MQ=255)
                                                                                                                                                                               gATACTCCCCGCCTCGCGACGGTTCCTATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGCAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTATGATGGTTAATCAGCATCGGAATTGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGAt                                                              >  2:194414/1‑251 (MQ=255)
                                                                                                                                                                               gATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGAt                                                              <  1:101791/251‑1 (MQ=255)
                                                                                                                                                                               gATACTCCCCGCCGCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGag                                                              <  2:310224/251‑2 (MQ=255)
                                                                                                                                                                                aTACTCCCCGCCTCGCTACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGTGCGGATGCCGTCATGAAAAACAGCCGCCGAGCGATGCGCCATAAGCGATGTCATACCCGTCAGGCACTCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTGCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACTGGCTCAACGGAAAGTCCATTCGCCGCTTTTCCTTCTGCGATc                                                             <  2:312830/251‑1 (MQ=255)
                                                                                                                                                                                aTACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATc                                                             >  1:107759/1‑251 (MQ=255)
                                                                                                                                                                                     ccccGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTTTGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCAtttt                                                        <  1:438934/251‑1 (MQ=255)
                                                                                                                                                                                     ccccGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTTTTCGACGTGCTCAACGGAAAGGCCAAACGCCTCTTTTCCTTCAGCGATCAtt                                                          <  1:301659/249‑1 (MQ=255)
                                                                                                                                                                                                      tCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATTAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCtt                                       >  2:344975/1‑251 (MQ=255)
                                                                                                                                                                                                                 tCAATCAGTTGCGGCAGGATCTGGCGGATGCCGTCATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAAt                             <  1:79579/250‑1 (MQ=255)
                                                                                                                                                                                                                 tCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAAt                             <  1:71954/250‑1 (MQ=255)
                                                                                                                                                                                                                 tCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGACAAAAGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATg                            >  1:295799/1‑251 (MQ=255)
                                                                                                                                                                                                                   aaTCAGTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGTTCGTCTTTCGCTTTAATGcc                          >  1:374810/1‑251 (MQ=255)
                                                                                                                                                                                                                        gTTGCGGCAGGATCGGGCGGATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATAACCGTCAGGCACGCCACCACGATAATCTGCATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATAGGAATCGGGTGTTCGACGTGCTCCACGGAAAAGCCAACAGCCACTTTTCCTTAAGCGATCATTTTGGCGGTATGCTCGTCTTTCGCTTTAATGCCTGGtt                     >  2:378401/1‑251 (MQ=255)
                                                                                                                                                                                                                               cAGGATCGGGAGGATGCCGTCATGAAAAACAGCCGCCCAGCGGTGCGCCAGAAGCGATGTCATACCCGTAAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAGAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATGCCTGGTTTGCCAAt              >  2:62713/1‑251 (MQ=255)
                                                                                                                                                                                                                                  taTCGGGCGGATGCCGTCATGAAAAACAGCCGCCGCGCGATGCGCCAGAAGCCATGTCATACCCGTCATGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGCCGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATGCCTGGTTTGGCAATGGt           >  2:301068/2‑251 (MQ=255)
                                                                                                                                                                                                                                          ggATGCCGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATcg                                                                                                                  >  2:305934/1‑140 (MQ=255)
                                                                                                                                                                                                                                          ggATGCCGGTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATcg                                                                                                                  <  1:305934/140‑1 (MQ=255)
                                                                                                                                                                                                                                           gATGCAGTTATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGAGCCATAGTCGTGCCTGTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGAGATGCTCGTCTTTCGCTTTAATGCCTGGTTTGGCAATGGTGCAACtttt  <  2:52094/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
AAGTAAGGATCAACACGCCCCAGATAGCGCCAAGACCAAACGCCATCAGGCTGTTTATCGCCAGGACGCCGAGGATATCGGTAGGAAGAAAGAGTAACCATGCGTTGAGCAGCCCGGTTTTCAGGGTAAAAGAGATCCCCACCGAGGCCACAAAACCAATGCTCAGACCAAAAGCGATACTCCCCGCCTCGCGACGGTTCATATAGCCTTCAATCAGTTGCGGCAGGATCGGGCGGATGCCGTCATGAAAAACAGCCGCCGAGCGATGCGCCAGAAGCGATGTCATACCCGTCAGGCACGCCACCACGATAATCTGAATATACAGATCCATAGTCGTGCCTTTTATTTCAGATGGTTAATCAGCATCGGAATCGCGTGTTCGACGTGCTCAACGGAAAGGCCAAACGCCACTTTTCCTTCAGCGATCATTTTGGCGATATGCTCGTCTTTCGCTTTAATGCCTGGTTTGGCAATGGTGCAACTTTT  >  NC_000913/3507006‑3507491

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: