Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,418,845 A→C 5.1% intergenic (+234/‑197) yhdV → / → yhdW putative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,418,8450AC 5.1% 273.8 / 5.5 79intergenic (+234/‑197)yhdV/yhdWputative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (43/32);  new base C (2/2);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

TGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGA  >  NC_000913/3418580‑3419119
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                  
tGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGgcg                                                                                                                                                                                                                                                                           <  2:211735/275‑1 (MQ=255)
                    ctATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCa                                                                                                                                                                                                                                                       >  2:64294/1‑275 (MQ=255)
                              aaCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGgcg                                                                                                                                                                                                                                             <  2:347169/275‑1 (MQ=255)
                               aCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGgcgc                                                                                                                                                                                                                                            <  2:437240/275‑1 (MQ=255)
                                                                                    tAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTgcg                                                                                                                                                                                       >  2:53268/1‑275 (MQ=255)
                                                                                            aaTATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCgg                                                                                                                                                                                                                                                 >  3:12114/1‑209 (MQ=255)
                                                                                            aaTATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTttatta                                                                                                                                                                                                                                 >  4:361348/1‑225 (MQ=255)
                                                                                                ttGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGgcg                                                                                                                                                                                                                                             <  1:12114/209‑1 (MQ=255)
                                                                                                ttGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTc                                                                                                                                                                                                                             <  2:361349/225‑1 (MQ=255)
                                                                                                 tGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTTCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTTCGCAAAAGCTGGCac                                                                                                                                                                          >  2:425303/1‑275 (MQ=255)
                                                                                                 tGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCac                                                                                                                                                                          >  2:257588/1‑275 (MQ=255)
                                                                                                     atatTCTGCTGGGTAATTCCCTGCTGATTGAGCATGGTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACTCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCTTTGCGGGGCGCTGTTTTTATTATTTCCTGTCTacac                                                                                                                                                                                                                   <  1:116416/230‑1 (MQ=255)
                                                                                                          ctgctgGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACtttta                                                                                                                                                                  >  2:196908/1‑273 (MQ=255)
                                                                                                           tgctgGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCTCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACtttttc                                                                                                                                                                 >  2:229351/1‑272 (MQ=255)
                                                                                                                     ttCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTc                                                                                                                                                                                                                             >  2:232806/1‑204 (MQ=255)
                                                                                                                     ttCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATAATTATTCTGCTTAGTAAAACACGCNTAAATANTTTGCGCCCAAAtg                                                                                                                                                                                                                                                                                >  3:27081/1‑153 (MQ=255)
                                                                                                                         ctgctgATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGgc                                                                                                                                                                                                                                                                            <  1:27081/153‑1 (MQ=255)
                                                                                                                         ctgctgATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGt                                                                                                                                                                                                                         <  4:232805/204‑1 (MQ=255)
                                                                                                                               aTTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCCTGTTTCATTTTCGCACCATTGTGGGGCGCTGTTTTTATTCTTTCCTTCCTACACTCAGCTTTTTGCGACGCGAAATCGTCCGAAACAGCTGTCACTCCTTTTGCTTATAAGAGTATGGCCacgga                                                                                                                                            >  1:54886/1‑275 (MQ=255)
                                                                                                                               aTTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATANTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTTCTGTNTACACTCAGCTTATTCCGACGCGAAATAGTGCGCAAAATCTGGTCNCNCTTTTTCTCATAATAGGATGTCcacaaa                                                                                                                                            >  3:51551/1‑273 (MQ=255)
                                                                                                                                tttAGAATGTTAAATCCCGTTATTTCGCTGCTTTCGAAATAACGAATCAATTTAATGTTCACACTTCCCACCAGAAATATAATCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTGGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCacaga                                                                                                                                            <  4:64293/274‑1 (MQ=255)
                                                                                                                                         ttACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATCTACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAAACAGGTAAAAg                                                                                                                                  >  2:115379/1‑275 (MQ=255)
                                                                                                                                               cgcgTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTCCGCAAAAGCTGGCACTCCTTTTGCTTATAAGAGGATGGCCCCAGACAGCTCAAAGACGttt                                                                                                                            >  4:262203/1‑275 (MQ=255)
                                                                                                                                                 cgTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAag                                                                                                                                                          >  2:22651/1‑243 (MQ=255)
                                                                                                                                                   ttATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCAtttt                                                                                                                                                                                                                                                              >  4:431296/1‑141 (MQ=255)
                                                                                                                                                     aTTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGa                                                                                                                                                      <  4:22651/243‑1 (MQ=255)
                                                                                                                                                       ttCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCa                                                                                                                                                                                                                                                          <  2:431297/141‑1 (MQ=255)
                                                                                                                                                             gctTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTc                                                                                                                                                                                                                                                                   >  6:281287/1‑126 (MQ=255)
                                                                                                                                                             gctTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTc                                                                                                                                                                                                                                                                   <  5:281287/126‑1 (MQ=255)
                                                                                                                                                             gctTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGgcg                                                                                                                                                                                                                                             >  8:238208/1‑148 (MQ=255)
                                                                                                                                                                            aaTCAATTTAATTTTCCCCCTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGc                                                                                               <  4:115378/274‑1 (MQ=255)
                                                                                                                                                                              tCAATTTAATTTTCACACTTGCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAAtgtg                                                                                                                                                                                                                                                                              <  8:821848/98‑1 (MQ=255)
                                                                                                                                                                              tCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAAtgtg                                                                                                                                                                                                                                                                              >  7:821848/1‑98 (MQ=255)
                                                                                                                                                                                 atttaatttTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGt                                                                                                                                                                                                                         <  7:238208/148‑1 (MQ=255)
                                                                                                                                                                                        tttCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAAtgt                                                                                                                                                                                                                                                                               <  6:396193/87‑1 (MQ=255)
                                                                                                                                                                                        tttCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAAtgt                                                                                                                                                                                                                                                                               >  5:396193/1‑87 (MQ=255)
                                                                                                                                                                                         ttCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGgc                                                                                                                                                                                                                                              >  5:56454/1‑119 (MQ=255)
                                                                                                                                                                                         ttCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGgc                                                                                                                                                                                                                                              <  6:56454/119‑1 (MQ=255)
                                                                                                                                                                                                ttCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCACTGTTTTTATTATTTNCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGNATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTTCTATAACANANATTTTCCCCACACAGGAAAAATTaaaaa                                                                           >  3:50422/1‑271 (MQ=255)
                                                                                                                                                                                                 tCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGaaa                                                                          <  4:257587/275‑1 (MQ=255)
                                                                                                                                                                                                  cccACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTgg                                                                                                                                                                                                                                                                             >  7:981265/1‑79 (MQ=255)
                                                                                                                                                                                                  cccACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTgg                                                                                                                                                                                                                                                                             <  8:981265/79‑1 (MQ=255)
                                                                                                                                                                                                     accaGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTtta                                                                                                                                                                                                                                    <  8:1115779/117‑1 (MQ=255)
                                                                                                                                                                                                     accaGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTtta                                                                                                                                                                                                                                    >  7:1115779/1‑117 (MQ=255)
                                                                                                                                                                                                     accaGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAga                                                                      <  1:51551/275‑1 (MQ=255)
                                                                                                                                                                                                       ccgAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGttt                                                                                                                                                                                                                                                                    <  8:191683/81‑1 (MQ=255)
                                                                                                                                                                                                       caGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGttt                                                                                                                                                                                                                                                                    >  7:191683/1‑83 (MQ=255)
                                                                                                                                                                                                        aGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGtt                                                                                                                                                                                                                                        >  5:261486/1‑110 (MQ=255)
                                                                                                                                                                                                        aGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGtt                                                                                                                                                                                                                                        <  6:261486/110‑1 (MQ=255)
                                                                                                                                                                                                          aaaTATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAAtg                                                                                                                                                                                                                                                                                >  7:609465/1‑68 (MQ=255)
                                                                                                                                                                                                          aaaTATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAAtg                                                                                                                                                                                                                                                                                <  8:609465/68‑1 (MQ=255)
                                                                                                                                                                                                            atatACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCAtttt                                                                                                                                                                                                                                                              >  6:857732/1‑84 (MQ=255)
                                                                                                                                                                                                            atatACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCAtttt                                                                                                                                                                                                                                                              <  5:857732/84‑1 (MQ=255)
                                                                                                                                                                                                                   cACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGc                                                                                                                                                                                                                                                           >  7:531709/1‑80 (MQ=255)
                                                                                                                                                                                                                   cACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGc                                                                                                                                                                                                                                                           <  8:531709/80‑1 (MQ=255)
                                                                                                                                                                                                                    aCGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTTCGGGGCGCTGTTTTTATTATTTCCTGTCTACACCCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGTACTACTTTTGCTTATAAGAGTATGGCCACAGACGGGTAAAAGACGTTCCCCAAACGTCCTTTAACGATAATTTTCGACTCACAGGNNNCATACTGAAAAAGATGAAGATAGTcaca                                                        >  3:148565/1‑274 (MQ=255)
                                                                                                                                                                                                                     cGCTAACAATCATTATTCTGCTTAGTAGAACACGGGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGttttt                                                                                                                                                                                                                                     <  8:988744/100‑1 (MQ=255)
                                                                                                                                                                                                                     cGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGttttt                                                                                                                                                                                                                                     >  7:988744/1‑100 (MQ=255)
                                                                                                                                                                                                                          aCAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGttttt                                                                                                                                                                                                                                     <  5:678397/95‑1 (MQ=255)
                                                                                                                                                                                                                          aCAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGttttt                                                                                                                                                                                                                                     >  6:678397/1‑95 (MQ=255)
                                                                                                                                                                                                                                attattCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCCAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCTTCACAGGAAGCATTACTAAAAAGATGATGTTAGCCACATTGGCTcccgcc                                            >  1:35686/1‑274 (MQ=255)
                                                                                                                                                                                                                                attattCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGTATGGCCACAGACAGGGAAAAGACGTTTCCCCAACGTCCTATACCGATAATTTTCGCCTCACAGGAAGCATCATGAAAAAGATGATTATAGCCACACTGTCTgcctccg                                           >  1:115516/1‑274 (MQ=255)
                                                                                                                                                                                                                                      cTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGNCACTACTTTTGCTTATAAGAGGATGGCCACGGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCCNACAGGAAGCATTATGAAAAAAAATANAATAGCCACACTGGCTTCCGCCCGCGtgc                                     >  3:97308/1‑275 (MQ=255)
                                                                                                                                                                                                                                         cTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCact                                                                                                                                                                         >  4:189884/1‑140 (MQ=255)
                                                                                                                                                                                                                                         cTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGtgctg                                   >  2:271168/1‑274 (MQ=255)
                                                                                                                                                                                                                                         cTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGAANCACTTTTGCTTATCAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATGTTAGNCTCAACGGACGGACTATGAAAAAGCAGATANNAGCCAAACATGATGCCGACaacaggcggc                                  >  3:58365/1‑266 (MQ=255)
                                                                                                                                                                                                                                         cTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTACGAACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAAGTAAAAGACGTTTCCCAAACGTCCTATAAGCATAATTTTTGCCTAACAGGAAACATTACGAACAAGATGATGATAGCAACACTGGCTGCCGCCAGCGtgcagg                                  >  4:295596/1‑272 (MQ=255)
                                                                                                                                                                                                                                            aGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCAAAGGAAGCATTATGAAAAAGATGATGAAAACAACACTGGCTGCCGCCAGCGTGCTGCTTg                               >  4:63088/1‑275 (MQ=255)
                                                                                                                                                                                                                                             gTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACtt                                                                                                                                                                     <  2:189885/140‑1 (MQ=255)
                                                                                                                                                                                                                                             gTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCACTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGc                              <  1:50422/275‑1 (MQ=255)
                                                                                                                                                                                                                                                   cacGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACCACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGGGCCGCTTGCGGCTGc                        >  1:93705/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       cgTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGttt                                                                                                                                                                                                                                       <  7:648781/64‑1 (MQ=255)
                                                                                                                                                                                                                                                       cgTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGttt                                                                                                                                                                                                                                       >  8:648781/1‑64 (MQ=255)
                                                                                                                                                                                                                                                           attagtttGCGCTCAAATGTGGCGGATGTTTCATTTTCGCACCATTGCGGGTCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAgg                <  1:90682/270‑1 (MQ=255)
                                                                                                                                                                                                                                                                 atGCGCCACATTTTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCGAAAGCt                                                                                                                                                                               >  7:786892/2‑110 (MQ=255)
                                                                                                                                                                                                                                                                 atGCGCCACATTTTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCGAAAGCt                                                                                                                                                                               <  8:786892/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tGCGCCACATTTTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCtt                                                                                                                                                                                                            >  7:163828/1‑80 (MQ=255)
                                                                                                                                                                                                                                                                  tGCGCCACATTTTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCtt                                                                                                                                                                                                            <  8:163828/80‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCcctggc         >  4:324293/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                         aaaTGTGGCGCATGTTTTATTTTCGCACCATTGCGGGGCGCTGTTTTTATTAGTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGa  <  1:2822/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                  
TGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGTTATTTCGCTGCTTTCGAAATAACGCATCAATTTAATTTTCACACTTCCCACCAGAAATATACTCACGCTAACAATCATTATTCTGCTTAGTAAAACACGCGTAAATATTTTGCGCCAAAATGTGGCGCATGTTTCATTTTCGCACCATTGCGGGGCGCTGTTTTTATTATTTCCTGTCTACACTCAGCTTATTGCGACGCGAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGA  >  NC_000913/3418580‑3419119

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