Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,739,682 T→C 18.6% E129G (GAG→GGG)  tyrA ← fused chorismate mutase T/prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,739,6820TC18.6% 135.8 / 5.3 65E129G (GAG→GGG) tyrAfused chorismate mutase T/prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (30/22);  new base C (4/9);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACA  >  NC_000913/2739447‑2739910
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                    
cAGGCTACCCGTGGCCGGCCCGAACATCGGGGGTAGCCCCAGCTCCGGGCCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCGTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCa                                                                                                                                                                                                                                   <  2:198657/239‑1 (MQ=255)
          cttgccGGACCCAACATCGGGTGTGGCCCCAGCCCCGGACCATCATGCCCCCCCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCAc                                                                                                                                                                                                                         <  1:485678/235‑1 (MQ=255)
                   ccGAACAGCGGGTGTAGCGCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAAc                                                                                                                                                                                                                 <  2:200004/238‑1 (MQ=255)
                       aCATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGc                                                                                                                                                                                                                           <  2:601600/224‑1 (MQ=255)
                       aCATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCgag                                                                                                                                                                                                            <  2:508696/239‑1 (MQ=255)
                       aCAGCGGGTGTAGCCCAAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCgag                                                                                                                                                                                                            <  2:442897/239‑1 (MQ=255)
                        aagcggtggtAGCCCCCGCCCCGGACACTCAGGCGCCCCCAGCATGGCCTGTAATGGCCCATTTTCCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCCCCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGGCAGGTTGCTCCAGAATCCGCACCTGATAACCCgaga                                                                                                                                                                                                           <  1:259015/231‑1 (MQ=255)
                                 tacccccACCCCCGGCCCCTCATGCCCCACCAGCTTGGCCTGTAATGGCCAATTTTTCCTTGATGCCAGATCGCCCAGAATACAATCTTTCGTGAAAGGCGG‑TAATTTGCCAATACCTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCa                                                                                                                                                                                                  <  1:88127/236‑1 (MQ=255)
                                 tAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATAACCATTCCGGCATCGGCACCCATAACAGCCCCCCGATCCCAGTCCTTTTTCTCCAAAAACCCCCCCTGATAAacccgctccttttaac                                                                                                                                                                                                  >  1:601600/1‑226 (MQ=255)
                                 tAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCAACCTAA‑CAATCACCATTCCGGCACCGGCAACAACATCCGCCGCTCGACCCCAGTCATTTTTCTCCCTAATCCCCACCTCTAAACCCCAGACGGGCacc                                                                                                                                                                                                   >  1:508696/1‑236 (MQ=255)
                                   gCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCCGTAACGTGGATTGGCCCACTAA‑CAATCACCATTCCGGCATCGGCACCAATATCAGCCGCTCGATCCCAGTCACGTTGCCCCAGAATCCGCCCCTGCTAAccccggatggccacca                                                                                                                                                                                                  >  2:541106/1‑224 (MQ=255)
                                        aGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCCCCAGAATCCGCACCTGATAACCCGAAAGGGGCAGCAtctttcc                                                                                                                                                                                           >  2:728925/1‑236 (MQ=255)
                                         gCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGTTCAGTAACGGGGATTGGCCAACTAA‑CAAACAACATTCCGGCATCGGCAAACAAATAAGCCCCCCGATCCCCACAACTTTGTTCCAAAACCCACACCCCACAACCCCGGAGGGTCAACATCTTCTCg                                                                                                                                                                                          >  1:804965/1‑239 (MQ=255)
                                             cggaccctcgaggccccccagattGGCCTTTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACa                                                                                                                                                                                      <  1:429948/216‑1 (MQ=255)
                                             cGGACCATCATGCGCCCCCAGCGTGGCCTGGAATGTCCCATTTTTCACTGATGCCAGAACGACCAGAATACAAACTTTGTGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACa                                                                                                                                                                                      <  1:108774/239‑1 (MQ=255)
                                               gACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTCA‑CAATCACCATTCCGGCATCGGCAACAAATTCAGCCGCTCGATCCCAGCCATGTTGCTCCCGAATCCGCACCTGGTAACCCGCGAGGGCCACCAcctttccggaaagg                                                                                                                                                                                    >  2:86471/1‑229 (MQ=255)
                                               gACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCACCTTCTCGAACAg                                                                                                                                                                                     >  2:150073/1‑238 (MQ=255)
                                                    tcatGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCCTCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCACCCCCTCATAACCCCAGAGGGTAAGCAACTTCTCGAACAAGCGTcc                                                                                                                                                                               >  1:721743/1‑239 (MQ=255)
                                                    tcatGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAAAATCAACCGCCCCATCCCAGCCCTGTTTCTCCACAAACCGCCCCCGATAACCCGAGAGGGTCCACATCTTCTCaaaaaaggacc                                                                                                                                                                                >  1:149401/1‑230 (MQ=255)
                                                         cgcCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTTGCACACTCA‑CAATCACCATTCCGGCATCGGCAACAATATTAGCCGCTCGATCCCAGTCATGTTGCCCCCGAATCCCCACCTGATAACCCGCGAGGGTCAGCATCTTCTCGACCAGGCGTCCCa                                                                                                                                                                             >  2:701346/1‑236 (MQ=255)
                                                            caccaGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGt                                                                                                                                                                                                       <  2:537111/207‑1 (MQ=255)
                                                                     gcccttattgtccccTTTTTTCCATGTTTCCAGTTCGACCAAAATAAACTCTTTCGTTAAAGGGGG‑TAATTTCCCAAAAACTTGCTCGTTACCTGGGTTTGGCCCCCTAA‑CAATCCCCATTCCGCCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATTTTGCTCCAAAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAATAGGCCTCCCATCTAaccgccacagc                                                                                                                                                              <  1:458443/227‑1 (MQ=255)
                                                                      gCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAAATTGCTCAGAAACGTGGATTGGCCCACTAA‑CAAACACAATTCCGGCATTGGCAACAATATAAGACGCCCGATCCCAGTCATGTTGTTCCCAAATCCGAACCTGGTAAACCGTGAGGGCAAGCATCTTCTACAACAGTCGGCCCCCCTGaacgccaccgcc                                                                                                                                                             >  2:66613/1‑239 (MQ=255)
                                                                            aagggCCCATTTTTCACGGAGGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCCCCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGAt                                                                                                                                                       <  2:477131/236‑1 (MQ=255)
                                                                            aaTGGCCCATTTTTCACTGATGCCAGATCGCCCAAAATAAAATCTTGCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGAAACGTGGATTGGCACACTAA‑CAATCCCCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGAt                                                                                                                                                       <  2:94817/239‑1 (MQ=255)
                                                                            aaTGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCCCCCGAATCCCCACCTGTTAACCCGACAAGGTCCGCAACTTTTCGAACAGGCGTCCCACCCTAcccccaccccccacagc                                                                                                                                                       >  1:734188/1‑236 (MQ=255)
                                                                                 cccATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTTCGGCATCGGCAAACATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGGCAACCACTTCTCGAAAAGGCGGCCCATCCGACCCCCCCCCCCGACGATaaac                                                                                                                                                   >  1:751510/1‑236 (MQ=255)
                                                                                 cccATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGACGTCCCACCTGACCGCCCCCGCCACCGAGAacc                                                                                                                                                   >  2:310678/1‑238 (MQ=255)
                                                                                        ttCACTGATCCCAGACCGACAAGAATCCAATCTTTCGGT‑AAGGCGGTTATTTTGCCAATAACTTGGTCAGTACCTGGGATTGGCACACTAA‑CAATCCCCTTCCCGGACTCGGCAAAAAGACCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACg                                                                                                                                           <  1:701346/239‑1 (MQ=255)
                                                                                              tgtGCCAAATCCACCAAAATAACATTTTTCGGTAAAGGGGG‑TAATTTGCAAAAAACTTGCTAAGTAACTTGGTTTGGCACCATAA‑CAATCCCCATCCCGGCACCGGCAACAATACCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGa                                                                                                                                     <  1:150073/237‑1 (MQ=255)
                                                                                              gctGCCAAATCGCCCAAAAAAAAATCTTTGGGAAAGGGCGT‑ATATTTCCAAATACTTGGCCCATAAACG‑GGGTTGGCCCCCCAACCAACCCCCTTCCCGGCATCGGAAAAAACATCAGCCCCTCGACCCCGGTCATGTTGCCCCAAAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTCCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGa                                                                                                                                     <  1:804204/237‑1 (MQ=255)
                                                                                                     gATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAAAAGGCGTCCCATCTGCCCGCCACCGCCGACGATAACACACGGACGCAGTGACGGACaa                                                                                                                              >  2:328114/1‑239 (MQ=255)
                                                                                                        ctcacaaaaTACATTCTTTGTGAAAAGGGCG‑TAATTTTCCAAATACGTGCCTAATAACGGGGTTTGGCACTCTAA‑CACACCCCATCCCGGCCTCGACAACAATCCCCGCCGGCCGATCCCATTCATGTTTCTCCAGAACCCGCCCCGAATAACCGGGGAGGGTCAGCATCTTCTCGAACAGGCTTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAgt                                                                                                                           <  1:681195/232‑1 (MQ=255)
                                                                                                        cGACCAGAATCCAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCCCACTAA‑CAATCCCCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAgt                                                                                                                           <  2:688886/239‑1 (MQ=255)
                                                                                                                         ttCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTTATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGGGACGGACAAAAGGTTTTAAATCCTTTTTc                                                                                                          >  2:718527/1‑239 (MQ=255)
                                                                                                                         ttCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGCGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGGGACGGACAAAGGGTTTTAAATCCtttt                                                                                                            >  2:547122/1‑236 (MQ=255)
                                                                                                                         ttCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGGTTTAAATCCTTTGTc                                                                                                          >  2:19056/1‑239 (MQ=255)
                                                                                                                         ttCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGAAGCAGTGCCGGACAAAGGGTTTTAAATCCTTTGTc                                                                                                          >  2:576237/1‑239 (MQ=255)
                                                                                                                         ttCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCCCCGACGATAACCACCGGACGCCGTGACGGCCCAAGGGTTTTAAATCCtttt                                                                                                            >  2:540385/1‑236 (MQ=255)
                                                                                                                         ttCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGGCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACCGTAACCACCGGACCCAGGGACGGCCAAAGGGTTTTTAATCCTTTTTc                                                                                                          >  2:395450/1‑239 (MQ=255)
                                                                                                                                  aggagttattttGCCAATACCTTGCCCAGTAACTTGGTTTGGCACACTAA‑CAATCCCCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTCCTCCAGAATCCGCCCCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCAc                                                                                                  <  2:819057/230‑1 (MQ=255)
                                                                                                                                         aaTTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATTCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGGCCCATCTGACCGCCCCCGCCTGCGAAAACCACCGGACGCAGGGCCGGACAAAGGGTTTTAAAACCTTTGTCGTTTTTCCTGGAGTa                                                                                           >  1:880572/1‑239 (MQ=255)
                                                                                                                                                             aGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGGGTTTTAAATCCTTTGTCGGTTTCACTGGAGTAAAATTCACCCATAAcccac                                                                        >  1:772335/1‑236 (MQ=255)
                                                                                                                                                               tAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGc                                                                     >  2:871097/1‑239 (MQ=255)
                                                                                                                                                                aaagTGGATTGGCCCCCAAA‑CAATCCCCATTCCGGCATCGGCAACAATACCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGAc                                                                       <  2:149401/233‑1 (MQ=255)
                                                                                                                                                                aaCGTGGATTGGCACACTAA‑CAATCACCATCCCGGAATCGGCAACAATATCAGCCGCTCGATCCCAGTCATTTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGCCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGAc                                                                       <  2:721743/236‑1 (MQ=255)
                                                                                                                                                                   gtgttttGGCCCCCAAA‑CACTCCCCATCCCGACATCGCCAACAATCACTGCCCCCCGACCCCGGTCATGTGGCCCAAAAACCCCCCCCTGATACCCCGAGAGGGTCAGCATCTACTCGAACAGGCGTCCCATCTGACCGCCTCCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCaaaa                                                                 <  1:196698/234‑1 (MQ=255)
                                                                                                                                                                    tGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCCGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTAACTGGAGGAAGATTCAAGAATAACCCGACGAAAAAc                                                                >  1:297212/1‑239 (MQ=255)
                                                                                                                                                                            cccacTCA‑CCATCCCCCTTCCGGCCTCGGCAAAAATTCCAGCCGCCCGCTCCCAGTAATGTTGCCCAAGAACCCGCACCTGATACCCCGAGAGGGTCAGCATCTTCCCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCaa                                                        <  1:395450/237‑1 (MQ=255)
                                                                                                                                                                             ccccTAA‑CAATCCCCTTTCCGCCTGCGGCAAAAATACCCGCCGCGCGATCCCAGTCATGTTGCCCCAGAATCCCCCCCTGATACCCCGAGAGGGTCATCATCTTCTCGAACAGGCGCCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAAt                                                       <  1:728925/236‑1 (MQ=255)
                                                                                                                                                                                cTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGCCGTTAACCACCGGACGCAGGGACGGACAAAGTGTTTTAAATCCCTTGTCGTTTTCACTGGAGTAAAATTCACGCACCACCCGCCGCAAAAAATCCCCaaaaag                                                    >  1:351270/1‑235 (MQ=255)
                                                                                                                                                                                   a‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACCAAGCGGTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCCCGCATCAACCCACCCAAAACAACCCCAAtaaaaca                                                 >  1:155886/1‑233 (MQ=255)
                                                                                                                                                                                      aaTCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGGGTTTTAAATCCTTTGTCGTTTTCCCTGGAGTAAGATTAACGCCTAACCCGACGCAAAACATCCTCAAACaacacgg                                               >  1:619202/1‑233 (MQ=255)
                                                                                                                                                                                          aCCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGACTCACCCCACGCAAAACATCCTCAATCAGACCTGGCgg                                           >  1:135259/1‑239 (MQ=255)
                                                                                                                                                                                            ctttCCGGCACGGGCAACAATATCAGCCGCTCGATCCCATTCAGGTTGCCCCAGAACCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTa                                         <  2:351270/237‑1 (MQ=255)
                                                                                                                                                                                            cATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGCCAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGCGTAAGATTCACGCATCACCCGACGCAAAAAAACCTCAATCAGaccggagggaa                                         >  1:337782/1‑229 (MQ=255)
                                                                                                                                                                                            cATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGCCGCCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGGGACGGACAAAGGGTTTTAAATCCTTTGTCGTTTTCACAGGAGTAAGATTCAAGACACACACCACCCAAAACAACCCCATCCAGATCGGGCggaa                                         >  1:717487/1‑237 (MQ=255)
                                                                                                                                                                                            cATTCCGGCAACGGCAACAATACCCGCCGCCCGATCCCAAGCATGTTGCCCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTa                                         <  2:804965/239‑1 (MQ=255)
                                                                                                                                                                                              ttCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGGGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGCCCCCAAACATCCCCAATCAAATCTGGCGGGac                                        >  1:772023/1‑238 (MQ=255)
                                                                                                                                                                                                  ggCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCCCCTTAGAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCa                                   <  2:832329/239‑1 (MQ=255)
                                                                                                                                                                                                    cctgggCAAAATTATCAGCCGCTCGATCCCAGTCATGTTTCCCCAAAATCCGCCCCTGATACCCCGAGGGGGTCCGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCaga                                 <  1:208359/235‑1 (MQ=255)
                                                                                                                                                                                                                          tCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGCCGGACAAAGGGTTTTAAATCCTTTGTCGTTTTCACTGGAGTACGATTCAGCCATCACCCGACGCAAACCATTCCCAATCAACTCTGGCGGTAGACACAAAGCTTCCCCCTCTGCacaaccag           >  1:385263/1‑236 (MQ=255)
                                                                                                                                                                                                                            gaacccAGGCATGGTGCCCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCCCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAg         <  2:535444/236‑1 (MQ=255)
                                                                                                                                                                                                                                cctatCATGTTGCCCCAGAACCGCCCCCTGATAACCCGGAGGGGTCAGCATCTTCTCGACCAGGCGTCCCATCGGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGcta     <  2:332655/235‑3 (MQ=255)
                                                                                                                                                                                                                                   gTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCCCCTCCCGACCCCAGGCCacaa  >  2:294585/1‑236 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                    
CAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGG‑TAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAA‑CAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACA  >  NC_000913/2739447‑2739910

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