Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,206,172 T→C 5.8% T29A (ACG→GCG)  ttdR ← transcriptional activator of ttdABT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,1720TC 5.8% 263.4 / 3.7 81T29A (ACG→GCG) ttdRtranscriptional activator of ttdABT
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base C (2/3);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

GGTAATAGCTGGCGCAATATGGCTGCGCCCAAAACCAAAACTGCAACCAATACGAATCATCCCCTCCGGGCGCGTTTTGATTTGCGTGACGTCATCGACCAGTCGCTGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGG  >  NC_000913/3205909‑3206412
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
ggTAATAGCTGGCGAAATATGGCTGCGCCCAAAACCAAAACTGCAACCAATACGAATCATCCCCTCCGGGCGCGTTTTGATTTGCGTGACGTCATCGACCAGTCGCTGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTgcgg                                                                                                                                                                                                                                       <  2:484367/275‑1 (MQ=255)
           gcgcAATATGGCTGCGCCCAAAACCAAAACTGCAACCAATACGAATCATCCCCTCCGGGCGCGTTTTGATTTGCGTGACGTCATCGACCAGTCGCTGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCACCGTTGTTGCCAGCGTATTTTCACGGATTTGGATCCGCTTGGTGACAATTGCAGGCGTTTNNNCACGCGCTCCCGCCGCGGCGCCaaacaca                                                                                                                                                                                                                            >  3:149151/1‑271 (MQ=255)
            cgcAATATGGCTGCGCCCAAAACCAAAACTGCAACCAATACGAATCATCCCCTCCGGGCGCGTTTTGATTTGCGTGACGTCATCGACCAGTCGCTGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTcc                                                                                                                                                                                                                           >  2:272468/1‑275 (MQ=255)
            cgcAATATGGCTGCGCCCAAAACCAAAACTGCAACCAATACGAATCATCCCCTCCGGGCGCGTTTTGATTTGCGTGACGTCATCGACCAGTCGCTGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTAAGCGCCACGCCGCGAGCGGAGCGGTTCACCAGCTTTGTTGCCAGGGGCTTTTCAAGGATTTGGACCCGCTGGGGGACAAATCAAGGGGTTGTCCCAAAGGTACCCGCTGCGGCGCTaaaatccc                                                                                                                                                                                                                           >  4:389439/1‑271 (MQ=255)
               aaTAGGGCTGCGCCCAAAACCAAAACTGCAACCAATACGAATCATCCCCTCCGGGCGCGTTTTGATTTGCGTGACGTCATCGACCAGTCGCTGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTgc                                                                                                                                                                                                                                         <  1:149151/258‑1 (MQ=255)
                           cccAAAACCAAAACTGCAACCAATACGAATCATCCCCTCCGGGCGCGTTTTGATTTGCGTGACGTCATCGACCAGTCGCTGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAttt                                                                                                                                                                                                            <  2:426686/275‑1 (MQ=255)
                                     aaacTGAAACCAAGACGACGCACCCCGNNGGGGTGGCTTTTGAGTGCCGGGAAGTCATGNACCAATCGTAGAAAGTGGGTTAGGATTTCCAACGCATGTTCGTAACAGCGTTGGCCAATTTCAGTCAGCGCCACNCCTCGAGCGGAGCGATTCAGAAGCGTTGTTGACAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTGCTACCAGt                                                                                                                                                                                                    <  3:70878/273‑1 (MQ=255)
                                           cAACCAATACGAATCATCCCCTCCGGGCGCGTTTTGATTTGCGTGACGTCATCGACCAGTCGCTGATATTGGGTGAGGATTTCAAACGCATGTTCGTAACAACCTTCGCCACTTTCAGTCAGCGCCACGGCGCGAGCGTAGCGGTTCAGCCGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGACCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCCTCGCCCCTGCGCCCCTAAAACTCCCCCTGTGAACGATTTCTCCCACCACCTgaca                                                                                                                                                                                            >  4:19268/1‑272 (MQ=255)
                                                          atcCCCTCCGGGCGCGTTTTGATTTGCGTGACATCATCGACCAGTCGCTGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCACGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGCACGATTTCTACCAGTACCTGCCAGCCTTTGGTTAAAgg                                                                                                                                                                             >  1:2337/1‑275 (MQ=255)
                                                            ccccTCCGGGCGCGTTTTGATTTGCGTGACGTCATCGACCAGTCGCTGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGcc                                                                                                                                                                           >  2:276542/1‑275 (MQ=255)
                                                                       gcgTTTTGATTTGCGTGACGTCATCGACCAGTCGCTGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGc                                                                                                                                                                <  2:389444/275‑1 (MQ=255)
                                                                                    cGTGACGTCATTGACCAGTCGCCGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGACCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGTAAGTCTTTGGCTAAAGGCCAGCTATTCAGcctgtctcttata                                                                                                                                                    >  2:191220/1‑262 (MQ=255)
                                                                                                                    tGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTc                                                                                                                                                                                                                                                 >  5:281800/1‑149 (MQ=255)
                                                                                                                     gAGGATTTCCAACGCATGGTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGc                                                                                                                     <  2:19268/272‑1 (MQ=255)
                                                                                                                            tCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCtg                                                                                                                                                                                                                      >  2:503867/1‑168 (MQ=255)
                                                                                                                               aaCGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTgcggc                                                                                                                                                                                                                                      <  7:264674/149‑1 (MQ=255)
                                                                                                                                aCGCATGTTCGTAACNACGTNGGCCACTTTCAGTCAGCGCCACGNCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGaa                                                                                                                                                                                                                  <  4:503862/168‑1 (MQ=255)
                                                                                                                                 cGCATGTTCGTAACAACGTTTGCCACTTTCAGTCAGAGCCACGCCGCGAGCGGAGCGGTTCAGCAGAGTTGTTGCCAGCGTCTTTTCAAGGATTTGGATCCGCTTGGTGACATATGCAGTCGTCTGCCCAAGCGTCGCCGCTACGGCGCTAAAACTCACGCAGTGAACGTTTTCTGCACGTACCTGCAGCTCTTTAGCTCAAGGCCAGCTAATATGCATCGTTAATTTCCCCCCTTTGTCATAAGCGGCGTTACATTGGCAAAACATTTCCCttg                                                                                                      >  4:286452/1‑274 (MQ=255)
                                                                                                                                    aTGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttg                                                                                                                                                                                                                                 <  5:370270/149‑3 (MQ=255)
                                                                                                                                     tGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttgg                                                                                                                                                                                                                                <  7:200017/149‑4 (MQ=255)
                                                                                                                                           tAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTccc                                                                                                                                                                                                                          >  5:390656/1‑149 (MQ=255)
                                                                                                                                           tAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTccc                                                                                                                                                                                                                          >  8:439750/1‑149 (MQ=255)
                                                                                                                                             acaacGTTGGCCACTTTCAGTCAGGGCCACGCCGCGGGCGGAGCGGTTAAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGgcg                                                                                          <  2:9024/275‑1 (MQ=255)
                                                                                                                                             acaacGTTGGCAACTTTGAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGgcg                                                                                          <  2:43037/275‑1 (MQ=255)
                                                                                                                                                 cGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCtgtg                                                                                                                                                                                                                    >  6:248682/1‑149 (MQ=255)
                                                                                                                                                   ttGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGaa                                                                                                                                                                                                                  <  7:439750/149‑1 (MQ=255)
                                                                                                                                                    tGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAAc                                                                                                                                                                                                                 >  8:31156/1‑149 (MQ=255)
                                                                                                                                                    tGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAAc                                                                                                                                                                                                                 >  8:270882/1‑149 (MQ=255)
                                                                                                                                                        cACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAtt                                                                                                                                                                                                             >  5:231369/1‑149 (MQ=255)
                                                                                                                                                        cACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAtt                                                                                                                                                                                                             <  8:458571/149‑1 (MQ=255)
                                                                                                                                                          cTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTc                                                                                                                                                                                                           <  8:280158/149‑1 (MQ=255)
                                                                                                                                                                     cgcCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTc                                                                      >  2:284326/1‑271 (MQ=255)
                                                                                                                                                                      gcCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTg                                                                                                                                                                                               >  5:350782/1‑149 (MQ=255)
                                                                                                                                                                          cgccgcGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttg                                                                                                                                                                                                                                 >  7:118148/1‑109 (MQ=255)
                                                                                                                                                                          cgccgcGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttg                                                                                                                                                                                                                                 <  8:118148/111‑3 (MQ=255)
                                                                                                                                                                           gccgcGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGt                                                                                                                                                                                          <  6:390656/149‑1 (MQ=255)
                                                                                                                                                                               cgAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                                                                                                                      <  7:31156/149‑1 (MQ=255)
                                                                                                                                                                               cgAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                                                                                                                      <  7:270882/149‑1 (MQ=255)
                                                                                                                                                                                gagcggagcgGTTCAGCAGCGTTGGTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAGATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGaa                                                       <  1:8648/275‑1 (MQ=255)
                                                                                                                                                                                 agcggagcggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTgg                                                                                                                                                                                    <  5:248682/149‑1 (MQ=255)
                                                                                                                                                                                  gcggagcggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTgg                                                                                                                                                                                    >  5:171198/1‑148 (MQ=255)
                                                                                                                                                                                  gcggagcggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTgg                                                                                                                                                                                    <  6:171198/148‑1 (MQ=255)
                                                                                                                                                                                                     ttgttgCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTaaaa                                                                                                                                                                                                                               >  2:525542/1‑86 (MQ=255)
                                                                                                                                                                                                     ttgttgCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAgg                                                                                                                                                                             >  1:6591/1‑136 (MQ=255)
                                                                                                                                                                                                     ttgttgCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTAAGCATCGTTAATTATCCGCAGTTGTGATAAGCCCAGTGTATTTCGCAAAACATTCCCATTCCCTGCTCTGGCGAAATATTCATTCTTATTTTCGCGTTGATTAACGAACTTGTAATTACATATTCACgg                                  >  2:148712/1‑275 (MQ=255)
                                                                                                                                                                                                     ttgtggCAAGTGGATTTTCAAGGATTTGGATCCGCTTGGGGCAAAATGCAGCTGTCTGCCCAAGGGTAGCCGCGTCGGCGCTAAAACTCCCGCTGGGAACGATTGCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACg                                   <  4:272464/274‑1 (MQ=255)
                                                                                                                                                                                                         tgCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTcc                                                                                                                                                                                                                           <  4:525537/86‑1 (MQ=255)
                                                                                                                                                                                                         tgCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAg                                                                                                                                                                         <  3:6591/136‑1 (MQ=255)
                                                                                                                                                                                                             aGCGTGTTTTTAAGGATGTGGATCCGCTTGTGGACAAATGCAGGCGTCTTCCCAGGCGTCGCCGCTGCTGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCaaa                           <  4:90971/274‑1 (MQ=255)
                                                                                                                                                                                                               cGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCTCAAGCGTCGCAGCTGCGGCGCTAAAACTACCGCTGTGAACGATTTCTACCAGTAACTGCAAGTCTTTGGCTAAAGGCCAGATATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCATGCTCTGGCGCAATATTAAGTCTTATTTTCCCGTTGATTCACTCAAATTTAATTCACTATTCATGGTTGGaaaaaa                        >  4:5079/1‑275 (MQ=255)
                                                                                                                                                                                                               cGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCAATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCTTTGATTCACGAACTGTTAATTTCATATTCAAGGTTGCtaaaaa                        >  4:476163/1‑275 (MQ=255)
                                                                                                                                                                                                                gTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTa                                                                                                                                                     >  6:429395/1‑149 (MQ=255)
                                                                                                                                                                                                                 tATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATCTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTAt                                                                                                                                                    <  6:388212/149‑1 (MQ=255)
                                                                                                                                                                                                                        aaGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGt                                                                                                                                             >  6:296768/1‑149 (MQ=255)
                                                                                                                                                                                                                         aGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGtt                                                                                                                                            >  8:329704/1‑149 (MQ=255)
                                                                                                                                                                                                                         aGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGAt                                                              >  2:239260/1‑227 (MQ=255)
                                                                                                                                                                                                                           gATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTaa                                                                                                                                                        >  6:146069/1‑135 (MQ=255)
                                                                                                                                                                                                                           gATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTaa                                                                                                                                                        <  5:146069/135‑1 (MQ=255)
                                                                                                                                                                                                                            aTTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTtgtg                                                                                                                                         >  7:187973/1‑149 (MQ=255)
                                                                                                                                                                                                                             tttGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGa                                                                                                                                        >  7:188980/1‑149 (MQ=255)
                                                                                                                                                                                                                             tttGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCAc                                                          <  4:239256/227‑1 (MQ=255)
                                                                                                                                                                                                                                 gATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCGGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTa      <  2:300182/275‑1 (MQ=255)
                                                                                                                                                                                                                                   tCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGtt                                                                                                                                                          >  7:179256/1‑125 (MQ=255)
                                                                                                                                                                                                                                   tCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGtt                                                                                                                                                          <  8:179256/125‑1 (MQ=255)
                                                                                                                                                                                                                                     cGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATTTTCACGGTGGCAAAAAATTTAAAACCACATTTTCAGTgc  >  4:376267/1‑274 (MQ=255)
                                                                                                                                                                                                                                     cGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGATAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTgg  <  2:451636/275‑1 (MQ=255)
                                                                                                                                                                                                                                          ggTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTa                                                                                                                           >  7:311745/1‑149 (MQ=255)
                                                                                                                                                                                                                                            tGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTAtt                                                                                                                         >  5:356966/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aCAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTc                                                                                                                       <  7:329704/149‑1 (MQ=255)
                                                                                                                                                                                                                                               cAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACGCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGtt                                                                 <  1:7890/202‑1 (MQ=255)
                                                                                                                                                                                                                                                      aGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                                                                                                                      >  8:464018/1‑78 (MQ=255)
                                                                                                                                                                                                                                                      aGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCttt                                                                                                                                                                                      <  7:464018/78‑1 (MQ=255)
                                                                                                                                                                                                                                                      aGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACa                                                                                                               >  7:450534/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        gCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACAtt                                                                                                             <  6:231369/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          gTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGc                                                                                                           <  8:450534/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           tCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGcc                                                                                                          <  8:187973/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                ccAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTAtt                                                                                                                                                                    >  6:326268/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                ccAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTAtt                                                                                                                                                                    <  5:326268/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAATGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTccc                                                                                                    <  8:311745/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGc                                                                                                           >  8:243367/1‑137 (MQ=255)
                                                                                                                                                                                                                                                                      gTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGc                                                                                                           <  7:243367/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
GGTAATAGCTGGCGCAATATGGCTGCGCCCAAAACCAAAACTGCAACCAATACGAATCATCCCCTCCGGGCGCGTTTTGATTTGCGTGACGTCATCGACCAGTCGCTGATATTGGGTGAGGATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGG  >  NC_000913/3205909‑3206412

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: