Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,618,312 T→C G28G (GGT→GGC deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,618,3120TC100.0% 258.9 / NA 79G28G (GGT→GGCdeoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/39);  total (40/39)

GGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCG  >  NC_000913/4618066‑4618558
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
gggccACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTc                                                                                                                                                                                                                                                    >  2:526620/5‑251 (MQ=255)
   ccACGGTGACGTTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTcgcg                                                                                                                                                                                                                                                 >  2:356997/2‑251 (MQ=255)
                     cgcCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAATGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAg                                                                                                                                                                                                                               <  1:54772/251‑1 (MQ=255)
                     cgcCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGa                                                                                                                                                                                                                                 >  1:540477/1‑249 (MQ=255)
                       ccAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAggg                                                                                                                                                                                                                             <  1:311341/251‑1 (MQ=255)
                          gcagcTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTGTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAg                                                                                                                                                                                                                          <  2:103663/251‑1 (MQ=255)
                           aagcTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGc                                                                                                                                                                                                                            >  2:340089/2‑248 (MQ=255)
                               tactaAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTgc                                                                                                                                                                                                                     >  1:7497/1‑251 (MQ=255)
                                actaAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTgcc                                                                                                                                                                                                                    <  1:169085/251‑1 (MQ=255)
                                   aagtaagATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAg                                                                                                                                                                                                                          <  1:456035/242‑1 (MQ=255)
                                   aagtaagATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAg                                                                                                                                                                                                                          >  2:456035/1‑242 (MQ=255)
                                      taagATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTAGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTc                                                                                                                                                                                                              >  2:128833/1‑251 (MQ=255)
                                       aagATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTc                                                                                                                                                                                                              <  2:480386/250‑1 (MQ=255)
                                           ggCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGGGTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGAt                                                                                                                                                                                                         <  2:89771/250‑1 (MQ=255)
                                            gCTTTACGCCTGATTCGCTGCGCTTATCAGTCCTACGATACGTATCGACCCTTAGGCCGGATAAGGCGTAGTCTCATCCTGCAAAAGCCGCCTCATTCTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAATCTTGATGGTCATGCACTGAGCGATTAAGAAATTCGTTTCTTTATCAACTGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCTCCCTcgctat                                                                                                                                                                                                         <  1:572150/250‑4 (MQ=255)
                                             cTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGTGTTGAAGAACTTCGTTTCTTTATCAACGGCATTCGCGACACCCCTATCTCCGAAGGGCAGATTGCCGCCATCGCGATgc                                                                                                                                                                                                       >  2:171309/1‑250 (MQ=255)
                                             cTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATg                                                                                                                                                                                                        >  1:347265/1‑250 (MQ=255)
                                               ttACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGAc                                                                                                                                                                                                      <  2:431959/250‑1 (MQ=255)
                                                 cctccTGATGCGCTGCGCTTATCAGTCCTACGAGACGTATCTACCCGTAGGCAGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCATGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATTCTCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCAt                                                                                                                                                                                                   <  1:7666/248‑1 (MQ=255)
                                                 aCGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCAt                                                                                                                                                                                                   <  2:581418/251‑1 (MQ=255)
                                                       gATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTttct                                                                                                                                                                                              <  1:327207/250‑1 (MQ=255)
                                                         tgcgctgcgcTTATCAGACCTACGGGGCGTAGCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATTCTCTTTGCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGAGGAAGAAATTCGTTTCGGTATCGACTGCATTCGCTACAACTCTAGCGCCGAAGGGCAGATTGCCGCCCTCGCGAGGACTATTttcttc                                                                                                                                                                                            <  2:563645/250‑1 (MQ=255)
                                                           cgctgcgctTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCAc                                                                                                                                                                                         <  1:204679/251‑1 (MQ=255)
                                                                    tATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGAc                                                                                                                                                                                 >  1:212237/1‑250 (MQ=255)
                                                                    gaTCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGATCATTTTCTTCCACGATATGa                                                                                                                                                                                  <  1:553922/248‑1 (MQ=255)
                                                                        aGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGc                                                                                                                                                                            >  1:71654/1‑251 (MQ=255)
                                                                        aGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCAGTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGc                                                                                                                                                                            <  2:439707/251‑1 (MQ=255)
                                                                           ccTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTg                                                                                                                                                                         >  1:556260/1‑251 (MQ=255)
                                                                              aCGAGACGTATCTACCCGTAGGCCGTATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATTCCTGAGc                                                                                                                                                                      <  1:442636/251‑1 (MQ=255)
                                                                              aCGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAg                                                                                                                                                                       <  2:81864/250‑1 (MQ=255)
                                                                                                     cGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATg                                                                                                                                               <  1:183946/251‑1 (MQ=255)
                                                                                                                    cGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAAcc                                                                                                                                >  1:584352/1‑251 (MQ=255)
                                                                                                                       aTCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGCCAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGtt                                                                                                                             >  1:247987/1‑251 (MQ=255)
                                                                                                                              aaaaGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATCTGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACt                                                                                                                      >  2:360611/1‑251 (MQ=255)
                                                                                                                               aaaGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACCAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGACATTCGTTTCTTTATCAACGGCATTAGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTg                                                                                                                     >  1:258595/1‑251 (MQ=255)
                                                                                                                               aaaGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACAATATCTCCGTAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGACCCGTTCTCGACTg                                                                                                                     >  1:410833/1‑251 (MQ=255)
                                                                                                                                       cTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGcc                                                                                                             >  1:451690/1‑251 (MQ=255)
                                                                                                                                            aCTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGTAAAAGCCTGCAt                                                                                                        <  2:359966/251‑1 (MQ=255)
                                                                                                                                                     cTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGc                                                                                               >  2:487111/1‑251 (MQ=255)
                                                                                                                                                       cGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGccc                                                                                             <  1:356997/251‑1 (MQ=255)
                                                                                                                                                        gggAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCg                                                                                            >  1:9598/1‑251 (MQ=255)
                                                                                                                                                         ggAGGTTACCTTGATTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTTGATAAGCCTTCATCTGAATGGCCCGa                                                                                           <  2:382885/251‑1 (MQ=255)
                                                                                                                                                          gAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCAACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCCGAATGGCCCGAt                                                                                          <  1:363135/251‑1 (MQ=255)
                                                                                                                                                            ggTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGAttg                                                                                        >  2:409533/1‑251 (MQ=255)
                                                                                                                                                             gTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGAttgt                                                                                       >  2:609919/1‑251 (MQ=255)
                                                                                                                                                                aCCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCAGCTGAATGGCCCGAttgttg                                                                                     <  1:485791/250‑1 (MQ=255)
                                                                                                                                                                  cTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGAt                                                                                   >  1:544011/1‑250 (MQ=255)
                                                                                                                                                                        tcgcGCACAAGAAATTATTCGTAAAAAACGTGATGGGCATGCACTGAGCGATGAAGAAATTCGTTTCTTTAGCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAacat                                                                             <  1:344456/247‑2 (MQ=255)
                                                                                                                                                                           cGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGGCTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATACACATTccc                                                                         >  1:385898/1‑250 (MQ=255)
                                                                                                                                                                             cacaAGAAATTATTCGTCAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGACAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCAct                                                                       >  2:33630/1‑250 (MQ=255)
                                                                                                                                                                                 aGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCGTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCgt                                                                                                                                                                    <  2:227157/154‑1 (MQ=255)
                                                                                                                                                                                     attattCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGt                                                               >  1:204571/1‑251 (MQ=255)
                                                                                                                                                                                           cGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAACTTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATAACAATGCCTg                                                                                                                                                                         >  1:624240/1‑139 (MQ=255)
                                                                                                                                                                                           cGTAAAAAACGTGAGGGTTATGCACGGAGCGATGAAGAAATTCGTTTCTTTATCATCGGCATTCGCGACAACACTATCGCCGAAGGGGAGATTGCCGTCCTCGCGATGACCATGTTCTCCCACTATATGACAATGCCTg                                                                                                                                                                         <  2:624240/139‑1 (MQ=255)
                                                                                                                                                                                            gTAAACAACGTGCTGGTCATGCACTGAGCGATTAAGAACTTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGACGTGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGCTATGACAATGCCTGAGCGTGTCTCGCTGTCCCTGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAACAGCCTGCATCTGACTGGCCCGATTTTTGATCACCATTCCACTGGTGGCGTCGGCGAt                                                        >  2:355238/1‑251 (MQ=255)
                                                                                                                                                                                            gTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGCCTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGAt                                                        >  2:250603/1‑251 (MQ=255)
                                                                                                                                                                                             tAAAAAACGTGATGGCCATGCACTGAGAGATGAAGAAATTCGTTTATTTAACAACGGCATTCGCGACCACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCATGTCTCGCTGCCCATCGCGATGCGAGCTTCAGGAACCGTTCTCGACTGGAAAAGCCTGCAACCGAATGTCACGATTGTTGATAAACATTcccctgggtgcctcggcgatg                                                       >  2:536594/1‑232 (MQ=255)
                                                                                                                                                                                              aaaaaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGAtgt                                                      <  1:303528/251‑1 (MQ=255)
                                                                                                                                                                                                aaaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACTCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGGTGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGtt                                                                                                                             <  2:403960/178‑1 (MQ=255)
                                                                                                                                                                                                aaaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGGTGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGtt                                                                                                                             >  1:403960/1‑178 (MQ=255)
                                                                                                                                                                                                  aaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACt                                                  >  1:163051/1‑251 (MQ=255)
                                                                                                                                                                                                   aCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACtt                                                 >  1:159040/1‑251 (MQ=255)
                                                                                                                                                                                                       gATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGACTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGACAAGCCTGCATCTGAATGGCCCGATTGTTGATCAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCt                                             >  1:421556/1‑251 (MQ=255)
                                                                                                                                                                                                           gTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATg                                         >  1:390690/1‑251 (MQ=255)
                                                                                                                                                                                                             cATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGtt                                       <  1:128833/251‑1 (MQ=255)
                                                                                                                                                                                                               tgcacTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTg                                                                                                                     >  2:163773/5‑171 (MQ=255)
                                                                                                                                                                                                               tgcacTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTg                                                                                                                     <  1:163773/167‑1 (MQ=255)
                                                                                                                                                                                                                gcacTGAGCGATGAAGAAAATCGTTTCTTTATCAACGGCATTCGCGACAACCCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGAGATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTggg                                    <  2:247987/248‑1 (MQ=255)
                                                                                                                                                                                                                     gAGCGATGAAGAAAGTCGTTTATTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTAATCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGcc                                 <  2:495803/249‑1 (MQ=255)
                                                                                                                                                                                                                       gCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATg                             <  1:349342/251‑1 (MQ=255)
                                                                                                                                                                                                                            gaagaaATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCg                         <  1:250603/250‑1 (MQ=255)
                                                                                                                                                                                                                                aaaTTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGAt                                                                                                                                                <  2:365098/127‑1 (MQ=255)
                                                                                                                                                                                                                                aaaTTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGAt                                                                                                                                                >  1:365098/1‑127 (MQ=255)
                                                                                                                                                                                                                                       tttCTTTATCAACGGCATTAGAGCCAACACTATCTACGAAGGGCAGATTGCCGCCCTCGCGCTGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTAGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCCGCCTgcggcg             >  2:223892/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cTTTATCCACGGCATTCGCGACAACACTATCTCCGTCGGGCAGATTGCCGTCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATGTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTAGCTGATGTTGGGGCCGCTGGTCGCAGCCTGCGGCGGCt          >  2:61376/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ttATCAACGGCATTCGCGACAACCCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCta         <  2:451690/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaCGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGAt                                                                                                                                                >  1:620306/1‑110 (MQ=255)
                                                                                                                                                                                                                                                 aaCGGCATTCGCGACAACACTATCGCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGAt                                                                                                                                                <  2:620306/110‑1 (MQ=255)
                                                                                                                                                                                                                                                  aCGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCg  <  1:33630/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
GGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCG  >  NC_000913/4618066‑4618558

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