Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,424,646 C→A 8.0% noncoding (2138/2904 nt) rrlD ← 23S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,424,6460CA 8.0% 358.3 / 7.1 146noncoding (2138/2904 nt)rrlD23S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base C (65/67);  new base A (2/11);  total (67/79)
Fisher's exact test for biased strand distribution p-value = 2.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

CTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGAC  >  NC_000913/3424505‑3424792
                                                                                                                                             |                                                                                                                                                  
ctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTa                                                                                                                                             <  3:136639/149‑1 (MQ=31)
ctccGTTACTCTTTAGGAGGAGACCGCACCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGTGTCAACGTTAGAACATAAAACAATAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCt                                                                                                                                              >  4:219869/1‑148 (MQ=2)
 tccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCt                                                                                                                                              >  7:323906/1‑147 (MQ=11)
   cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTAGCCAACAGACACTGTCCGCTACCCGGATTTCGGGTCAACGTTAGAAAATCAAACATTTAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTaa                                                                                                                                            <  2:372686/147‑1 (MQ=11)
   cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACGAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGc                                                                                                                                          <  8:275002/149‑1 (MQ=17)
   cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTaa                                                                                                                                            <  1:334151/147‑1 (MQ=32)
   cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTaa                                                                                                                                            <  7:255327/147‑1 (MQ=32)
        ctctTTCGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTcc                                                                                                                                      <  4:240350/148‑1 (MQ=17)
        ctctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTccc                                                                                                                                     <  4:312196/149‑1 (MQ=32)
         tctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGGGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTcccc                                                                                                                                    >  3:211902/1‑148 (MQ=17)
           tttAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAATGGTGGTATTTCAAGGTAGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAcc                                                                                                                                  >  8:136694/1‑149 (MQ=17)
           tttAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAcc                                                                                                                                  >  6:177909/1‑149 (MQ=32)
             tAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAc                                                                                                                                   <  2:278153/146‑1 (MQ=32)
              aggaggAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTcc                                                                                                                                      >  7:359668/1‑142 (MQ=32)
              aggaggAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTcc                                                                                                                                      <  8:359668/142‑1 (MQ=32)
                gaggagACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATcc                                                                                                                             >  3:350787/1‑149 (MQ=32)
                   gagaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTcc                                                                                                                                      >  3:89425/1‑137 (MQ=32)
                   gagaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTcc                                                                                                                                      <  4:89425/137‑1 (MQ=32)
                      aCCGCCCCAGTCAAACTACCCCCCAGACACTGTCCGCAACCCGGATTACGGGTCAACGGTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTAcaca                                                                                                                       >  5:367348/1‑149 (MQ=17)
                        cGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGTTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATc                                                                                                                     <  1:162161/149‑1 (MQ=17)
                          ccccAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACCCATc                                                                                                                     >  1:293713/1‑147 (MQ=17)
                           cccAGTCAAACTACCCATCAGACACTGTCAGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCAACACATCAAg                                                                                                                  >  5:185163/1‑149 (MQ=12)
                           cccAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAg                                                                                                                  >  6:48546/1‑149 (MQ=32)
                            ccAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAgg                                                                                                                 <  5:10559/149‑1 (MQ=32)
                                  aaaCTACCCTCCAGACAATGTCAGCAACGCTGATTACGGGTCAACGTTAGAACATCAAACATTAAAGTGTGGTATTTCAAGGTCGTCTCCATGCATACTGGAGTCCACAATTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCaa                                                                                                            <  3:309565/148‑1 (MQ=11)
                                  aaaCTACCCACCAGACACTGTCCGCAACGCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTTGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAAt                                                                                                           <  7:350235/149‑1 (MQ=17)
                                  aaaCTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCATATTTTTTTTAACAATTTTATGCGTCCAACCTTATCTACACATCAAGGCTCAAt                                                                                                           <  2:259369/149‑1 (MQ=11)
                                  aaaCTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAg                                                                                                                                           <  4:393199/117‑1 (MQ=32)
                                  aaaCTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAg                                                                                                                                           >  3:393199/1‑117 (MQ=32)
                                  aaaCTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAAt                                                                                                           >  5:304541/1‑149 (MQ=32)
                                     cTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGtt                                                                                                        <  7:187320/149‑1 (MQ=32)
                                       aCCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCGCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCAACCTATCCTACACATCAAGGCTCAATGTTCa                                                                                                      <  8:413630/149‑1 (MQ=17)
                                        cccACCAGACACTGTCCGCAACGCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTTGTATTTCAAGGGAGGCTCCATGCAGACTGGCGTCCAAAATTCTAAGCCTCCCACCTATCCTACACATGAAGGCTCAATGTTCa                                                                                                      <  5:292929/148‑1 (MQ=11)
                                         ccaccaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCa                                                                                                                                    <  6:22857/117‑1 (MQ=32)
                                         ccaccaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCa                                                                                                                                    >  5:22857/1‑117 (MQ=32)
                                         ccaccaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCACACCTATCCTACACATCAAGGCTCAATGTTCAgt                                                                                                    <  3:169836/149‑1 (MQ=17)
                                          caccaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAgtg                                                                                                   >  2:161956/1‑149 (MQ=32)
                                                acacTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGc                                                                                             >  4:51205/1‑149 (MQ=32)
                                                      tCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTa                                                                                                                                             <  5:314469/95‑1 (MQ=32)
                                                      tCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTa                                                                                                                                             >  6:314469/1‑95 (MQ=32)
                                                         gCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTACGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTa                                                                                      >  5:91631/1‑147 (MQ=17)
                                                         gCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTa                                                                                      >  5:7106/1‑147 (MQ=32)
                                                               cGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGt                                                                                                         >  2:314291/1‑122 (MQ=32)
                                                                ggATTACGGGTCAACGTTAGAACATCAAACATTACAGTGTGGTATTTCAAGGTCGGCGCCAGGCAGACTGGCGGCCACACTTCTAAGCCTCCCACCTATCCTACACTTCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCAc                                                                             <  3:4984/149‑1 (MQ=11)
                                                                ggATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCAc                                                                             <  3:233659/149‑1 (MQ=32)
                                                                   ttACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGGCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACggg                                                                          <  5:390345/149‑1 (MQ=17)
                                                                     ccGGGTCAACGTTAGACCATCAAACATTAAAGGGTGGTATTGCAAGGTCGGCTCCATGCAGAATGTCGTCAACACTTCTAAGCCTCCAACCTATCCTAAACATCAATGCTCAATGTTAAGTGTCAAGCTATAGTAAAGGTTCACggg                                                                          <  6:91631/146‑1 (MQ=11)
                                                                     aCGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATTCAGACTGGAGTCCACACTTCTAAGCCTACCACCTATCCTACACATCAAGGCTCAATTTTCAGTGTCACGCTATAGTAAAGGTTCACgggg                                                                         <  3:144285/148‑1 (MQ=11)
                                                                      cGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTc                                                                       >  1:114629/1‑149 (MQ=32)
                                                                         gTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCttt                                                                    >  2:275110/1‑149 (MQ=32)
                                                                          tCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTc                                                                   >  2:419061/1‑149 (MQ=32)
                                                                            aaCGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCg                                                                 >  1:101996/1‑149 (MQ=32)
                                                                               gTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCa                                                                                                      <  5:276609/109‑1 (MQ=32)
                                                                               gTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCt                                                              >  1:91750/1‑149 (MQ=32)
                                                                               gTTAGAACATCAAACATTAAAGGCTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTAACATACACATCAAGGCGAAATTTTCAGTGACAAGCTATAGTAAAGGTTCACGGGGGCTTTCCGGCt                                                              >  1:300778/1‑149 (MQ=11)
                                                                                   gAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGcc                                                          >  7:121769/1‑149 (MQ=32)
                                                                                   gAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGcc                                                          >  5:233783/1‑149 (MQ=32)
                                                                                   gAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCAACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGcc                                                          >  1:319157/1‑149 (MQ=17)
                                                                                    aaCATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCcg                                                         >  4:338234/1‑149 (MQ=32)
                                                                                       aTCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGtt                                                                                                        >  5:80554/1‑99 (MQ=33)
                                                                                       aTCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCAATCATACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCgg                                                      >  4:135396/1‑149 (MQ=17)
                                                                                        tCAAACATTAAAGGGTGGTATTTCAAGTTCGTCTCCATGCAGATTTTTTTTTTAATTTTTATGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCt                                                              <  6:415012/140‑1 (MQ=2)
                                                                                            aCATTAAAGGGTGGTATTTCAAGGTCGGCTACATGCAGTATTTATTTTAAATTTTTATGCGTCAAACGTATTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTaca                                                 <  8:214545/149‑1 (MQ=11)
                                                                                             cATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTGCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTacac                                                <  5:48546/149‑1 (MQ=17)
                                                                                              aTTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTc                                                               >  1:409212/1‑133 (MQ=32)
                                                                                              aTTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACt                                               <  5:224201/149‑1 (MQ=32)
                                                                                              aTTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCAAAAATTCTAAGCCTCCCACCTATACTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTc                                                               <  2:409212/133‑1 (MQ=255)
                                                                                                 aaaGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGc                                             <  3:51205/148‑1 (MQ=32)
                                                                                                  aaGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCa                                            <  2:91750/148‑1 (MQ=32)
                                                                                                  aaGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCAt                                           >  4:234464/1‑149 (MQ=32)
                                                                                                     ggtggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCta                                                                                           >  6:275153/1‑98 (MQ=33)
                                                                                                     ggtggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCtt                                        <  1:419061/149‑1 (MQ=32)
                                                                                                       tggtATTTCAAGGTCGGCTCCATGCAGACTGGCTTCCAAACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAATCTATAGTAAAGGTTAACGGGGTATTTCAGACTTTCCTCGTGTACACTGCATAttaa                                      >  4:215822/1‑147 (MQ=11)
                                                                                                       tggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTca                                      <  1:275110/149‑1 (MQ=32)
                                                                                                       tggtATTTCAAGGTCGGCTCCATGCAGAATGGCCTCCACACTTCTAAGCCTCCAACCTATCCTAAACATCAAGGCTAAATGAACATTGTCCAGCTATATTAAATGTAAAAGGGGTCTATCCGTCTTGCCGAGTGTACACTGCATCTGca                                      >  4:210466/1‑149 (MQ=11)
                                                                                                        ggtATTTCAAGGTCGGCTCCATGCAGATTTTTTTTTATAATTTTATGCGTCCAACGTTTTCTACACATCAAGGATCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTca                                      <  6:185163/148‑1 (MQ=11)
                                                                                                         gtATTTCAAGGTCGGCTCCAAGCAGACTGGCGTCCACACTTCTAAGCCACCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGAGCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTcaca                                    <  1:133357/149‑1 (MQ=255)
                                                                                                           aTTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                                  >  7:12438/1‑149 (MQ=32)
                                                                                                            tttCAAGGTCGGCTCCATGTAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                                  >  7:9027/1‑148 (MQ=17)
                                                                                                            tttCAAGGGCGGCTCCATGCAGTCTGGGGTCGACAGTTCTAAGCCTCCCACGTATCCTACACATCAATGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                                  <  4:200372/148‑1 (MQ=11)
                                                                                                             ttCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCATTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGGCTGGCCGCGCGTACACTGCATCTTCACAGc                                  >  5:56422/1‑147 (MQ=255)
                                                                                                             ttCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                                  >  2:184560/1‑147 (MQ=32)
                                                                                                              tCAAGGTCGGCTCCATGCAGACTGGAGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGa                                <  6:304541/148‑1 (MQ=17)
                                                                                                                aaGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGTCGCGGGTACACTGCATCTTCACAGCGAGtt                             <  2:101996/149‑1 (MQ=17)
                                                                                                                aaGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCAAGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGt                              >  6:353586/1‑148 (MQ=17)
                                                                                                                    tCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCTAGTTCAAt                         >  2:6398/1‑149 (MQ=17)
                                                                                                                    tCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAAt                         >  3:369869/1‑149 (MQ=32)
                                                                                                                      ggCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAAttt                       >  7:166638/1‑149 (MQ=32)
                                                                                                                       gCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCCCGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAAttt                       >  7:385241/1‑148 (MQ=17)
                                                                                                                        cTCCATGCAGACTGGCGGCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCa                     <  1:161956/149‑1 (MQ=17)
                                                                                                                           gaTGGAGAGTGGCGTGCTCACTTCTAAGACTGCCTCCTATCCTACAGATCAAGGATCAATGATCAGTGTCAAGCTATAGTAAAGGTTCACTGGGTCTTTCCGTCTTGCCTCGGGTACACTGCATCTTCTCAGCGAGTTCAATTTCACt                   <  1:254934/147‑1 (MQ=11)
                                                                                                                            aTGCAGACTGTCGTCCAAACTTATAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGa                 <  4:47823/149‑1 (MQ=12)
                                                                                                                            aTGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACt                   >  6:67366/1‑147 (MQ=32)
                                                                                                                              gCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGt               >  4:366294/1‑149 (MQ=32)
                                                                                                                               cAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCAAGTTCAATTTCACTTAGt               <  7:136694/148‑1 (MQ=17)
                                                                                                                               cAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGATATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGt               <  3:226329/148‑1 (MQ=17)
                                                                                                                                   cTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCgg          <  6:7106/149‑1 (MQ=32)
                                                                                                                                        gTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGGCTCGGGTGga     <  6:294862/149‑1 (MQ=17)
                                                                                                                                         tCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGgag    <  5:13785/149‑1 (MQ=32)
                                                                                                                                           cacaCTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGAc  <  3:338234/149‑1 (MQ=32)
                                                                                                                                             |                                                                                                                                                  
CTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGAC  >  NC_000913/3424505‑3424792

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