Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,996 T→G D70A (GAT→GCT)  thiE ← thiamine phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,9960TG81.7% 163.4 / 29.2 72D70A (GAT→GCT) thiEthiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (3/10);  new base G (29/29);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

CCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCC  >  NC_000913/4193752‑4194236
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
ccctGCCAGCTGTTCCAGCCCCGGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACCCATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGa                                                                                                                                                                                                                                            <  2:402350/248‑1 (MQ=255)
 cGTGCCAGCTGTCCCCGCCCCTGCGGTGCAAAAGGCTTCTTTTTGGTTTGCTGCGGGAACACCTGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGGCATGTGTCGAAACGCCCAGCCGCAGTCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAGCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTTATCCCCAGCCGCCAGTAATCGTTGAt                                                                                                                                                                                                                                           <  2:39113/251‑1 (MQ=255)
  gTGCCAGATGTTCAAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGAGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGCTGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATAGCAAGCCGCCCGTAAGCGttngt                                                                                                                                                                                                                                           >  2:442816/1‑247 (MQ=255)
     ccAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATctgtgttaaac                                                                                                                                                                                                                                       >  2:214472/1‑241 (MQ=255)
          tGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATcg                                                                                                                                                                                                                                  >  1:223473/1‑251 (MQ=255)
             tCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATcgcgc                                                                                                                                                                                                                               >  1:9638/1‑251 (MQ=255)
              ccAGCCGCTGCGGTGCAGAAGGCATCGGTTTGGTTTGCGTCGGGAACACAGGTCCCAGCGCGATATAAGAGGGGCGTGCCGCCAGCGCGACGTCGATTTCCATATCGTCATGGGGCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAAGGGACGCCATACGCCTGATGCTTGTTCGCCAGCCGCCAGTAAGCGTTGATAAACAATcgcgc                                                                                                                                                                                                                               <  2:317730/250‑1 (MQ=255)
                aGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGtt                                                                                                                                                                                                                            >  1:52994/1‑251 (MQ=255)
                aGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTAATCTCCAGACGCCAATAAGCGTTGTTAAACAATCGCGCGtt                                                                                                                                                                                                                            >  2:292566/1‑251 (MQ=255)
                      tGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcg                                                                                                                                                                                                                      >  1:124006/1‑251 (MQ=255)
                      tGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcg                                                                                                                                                                                                                      >  2:380991/1‑251 (MQ=255)
                        cGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcggc                                                                                                                                                                                                                    <  1:18754/251‑1 (MQ=255)
                           ggCAGAGGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATTTAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCggcggc                                                                                                                                                                                                                 <  2:420249/250‑1 (MQ=255)
                             cAGAAGGCATCTGGTTGGTTGGCGTCGGGACCAAATGTCCGAGGGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTAATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGGTCGGTGGCTTGCAAATCTTCCTGCCCCAAAGGGACGCCATACGCCTGATTCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGccc                                                                                                                                                                                                               <  2:213137/251‑1 (MQ=255)
                                      attttttggttgtGCGTCGAGAACAAATGTCACAGCGACAGATAAGAGGGGCGTCTTGCACGCGCGACGTCGTTTTCCATATCGTCATGTGTCGAATCGCCCATCCTCGGGCCGCCCGCGCGGATGTCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGTCGCCTTACGCCTGATGCTTGCTCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAAt                                                                                                                                                                                                       <  2:464497/247‑1 (MQ=255)
                                         gtttgttttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAAt                                                                                                                                                                                                       <  1:64566/247‑1 (MQ=255)
                                              gtttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCAc                                                                                                                                                                                              >  1:365879/1‑251 (MQ=255)
                                                ttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACg                                                                                                                                                                                             >  1:372162/1‑250 (MQ=255)
                                                            acacATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTcc                                                                                                                                                                                 <  1:311591/251‑1 (MQ=255)
                                                                    cccAGCGCGATATAAGAGGGGCGTGCTGCCAGCGGGACGTCGATTTCCATATCTTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGGGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCTCCTCTTc                                                                                                                                                                         <  1:352054/251‑1 (MQ=255)
                                                                         cgggATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAAGGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATcgc                                                                                                                                                                    <  1:363420/248‑1 (MQ=255)
                                                                            gATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGcc                                                                                                                                                                 >  1:376353/1‑251 (MQ=255)
                                                                                taAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATc                                                                                                                                                             >  1:5858/1‑251 (MQ=255)
                                                                                 aaGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCt                                                                                                                                                            >  1:391216/1‑251 (MQ=255)
                                                                                         ggTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGGGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATGGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCg                                                                                                                                                    <  1:131987/250‑1 (MQ=255)
                                                                                          gTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGt                                                                                                                                                   >  1:379139/1‑251 (MQ=255)
                                                                                                 gaGCGGGACGTCGATTTCCAAATCGTCATGTTGGGAAACGCCCAGCCGCGGGCCGGCCGGGCGTATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCGTGATGCTTGATCGCCAGCCGCCAGTAAGCGGTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGGTCTTTGATGCGTAGCTGga                                                                                                                                            <  2:65238/250‑1 (MQ=255)
                                                                                                   gcgcGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCAAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGgaga                                                                                                                                          <  2:372162/251‑1 (MQ=255)
                                                                                                     gcgACGTCGATTTCCATATCGTCATGTGTAGAAACGCCCCGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCAAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAgt                                                                                                                                        <  1:318265/251‑1 (MQ=255)
                                                                                                      cgACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAgt                                                                                                                                        >  2:11369/1‑250 (MQ=255)
                                                                                                            cGATTTCCATTTCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAgtacgtac                                                                                                                                  <  1:117298/250‑1 (MQ=255)
                                                                                                                    atatCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGAGGGCATTGAGATCGGTGGCTTGCAAAACTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATc                                                                                                                         <  2:434162/251‑1 (MQ=255)
                                                                                                                     tatCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCCGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGGACGTACGCCTGCATcc                                                                                                                        >  1:122085/1‑251 (MQ=255)
                                                                                                                            aTGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGAc                                                                                                                 <  1:327908/251‑1 (MQ=255)
                                                                                                                            aTGTGTCGAAACGCCAGGCCGAAGGCCGGCCGCGCGTATGGCATTGAGATCGTTGGCTTGAAAATCTTCCTGCCCAAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGAc                                                                                                                 <  2:419905/251‑1 (MQ=255)
                                                                                                                                 gcGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATGGAGATCGGTGGCTTGCAAATCTTCCTCCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGAGCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCg                                                                                                            <  1:292956/250‑1 (MQ=255)
                                                                                                                                     aaCGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTTCCGCCACGACATCGGCTTCCACCTCTTCCTCGCGCCGCTCTTTGATGCGTAGCTGGAGAGTACGCACGCCTGCATCCAACAGACGTTCGATcc                                                                                                        >  2:194809/1‑251 (MQ=255)
                                                                                                                                                gCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTg                                                                                             >  2:189106/1‑251 (MQ=255)
                                                                                                                                                         ccGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCaccacc                                                                                    <  1:329585/251‑1 (MQ=255)
                                                                                                                                                         ccGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGCACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCaccacc                                                                                    >  2:456109/1‑251 (MQ=255)
                                                                                                                                                         ccGCGCGAATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCaccacc                                                                                    <  2:365879/251‑1 (MQ=255)
                                                                                                                                                               ggATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTAc                                                                              >  1:204753/1‑251 (MQ=255)
                                                                                                                                                                   ggCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGt                                                                           <  1:11369/250‑1 (MQ=255)
                                                                                                                                                                      aTTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTg                                                                       >  1:448520/1‑251 (MQ=255)
                                                                                                                                                                         gaaaTCGGGGGCTTGCAAATCTCCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCCCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAAc                                                                    <  2:185814/248‑1 (MQ=255)
                                                                                                                                                                               ggTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGCTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGCTGCGTAGCTGGAGAGTACGTACGCCCGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACACTCCTGAACGaaa                                                                >  2:280305/1‑249 (MQ=255)
                                                                                                                                                                                gTGGCTTGCAAATCTTCCTGCCCCAAATGGTCGCCATATGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCCGCGCCATTGCCGCCACGACATCGGATTCCACCTCTTCATCGCGCCGATCTTTTATGCGTAGCCGGAGAGGACGTAAGCCTGCATCCAAATGAACTTCGGTCCACTGGACGCTGTTCACCCACGGGGAAAGTCCTGAACGAAAAgg                                                             >  2:248402/1‑251 (MQ=255)
                                                                                                                                                                                   gCTTGCAATTCTTCCGGCCCCAAATGGCCGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCATTTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTAc                                                          <  1:442816/251‑1 (MQ=255)
                                                                                                                                                                                         aaaTCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACaggagg                                                    <  1:36907/251‑1 (MQ=255)
                                                                                                                                                                                                ccTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGCTGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACCAAAAGGCACAGGAGGAAAATc                                              >  2:45083/1‑250 (MQ=255)
                                                                                                                                                                                                 cTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATc                                              >  1:244346/1‑249 (MQ=255)
                                                                                                                                                                                                   cccccAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGc                                          <  2:204753/250‑1 (MQ=255)
                                                                                                                                                                                                    ccccAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCt                                         >  1:420325/1‑251 (MQ=255)
                                                                                                                                                                                                    ccccAAATGGACGCCATACGCCTGATGCTTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCt                                         <  2:379139/251‑1 (MQ=255)
                                                                                                                                                                                                       cAAATGGACGCCATACGCCGGATGCTTGATCGCCAGCCGCCAGTAACCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGAt                                      <  1:93381/251‑1 (MQ=255)
                                                                                                                                                                                                        aaaTGGACGCCATCCGCCGGAGGTTTGAGCCCCACCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATa                                     <  2:52994/251‑1 (MQ=255)
                                                                                                                                                                                                        aaaTGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGAt                                      >  1:229379/1‑250 (MQ=255)
                                                                                                                                                                                                         aaTGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATAc                                    <  1:380991/251‑1 (MQ=255)
                                                                                                                                                                                                           gggACGCCATACGCCTGATGCTGGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACAt                                  <  2:122085/250‑1 (MQ=255)
                                                                                                                                                                                                             gACGCCATCCGCCTGATGCTTTATCCCAAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGTCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCa                                <  2:55776/251‑1 (MQ=255)
                                                                                                                                                                                                              aCGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCAc                               <  2:106922/251‑1 (MQ=255)
                                                                                                                                                                                                              aCGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCAc                               <  1:214472/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCt                            >  2:12061/1‑251 (MQ=255)
                                                                                                                                                                                                                      cGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTtcctc                       >  1:90918/1‑251 (MQ=255)
                                                                                                                                                                                                                       gCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTtcctcc                      >  1:106859/1‑251 (MQ=255)
                                                                                                                                                                                                                          ggATGCTGGATCGCCAGCAGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCtt                    <  2:233943/249‑1 (MQ=255)
                                                                                                                                                                                                                           gATGCTTGATCCCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTa                   <  2:361624/250‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCCATTGCCGCCACGACATCGGCTTCCACCTCTTAATCGCGCCGATCTTTGATGCGTAGCTGGAGAGCACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGCAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGCAGAtt       >  2:224095/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTtc     <  1:194809/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ccgccATTAAGCGTCGATAAACAATCGCGCGTTATGGCGGCGGCCCAGCCCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCGTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTcc  <  2:223473/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 gTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGCTCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTa           >  2:441691/1‑236 (MQ=255)
                                                                                                                                                                                                                                                 gTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTa           <  1:441691/236‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
CCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCC  >  NC_000913/4193752‑4194236

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: