Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,424,661 C→A 9.2% noncoding (2123/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,6610CA 9.2% 307.7 / 9.8 134noncoding (2123/2904 nt)rrlD23S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base C (65/55);  new base A (2/12);  total (67/67)
Fisher's exact test for biased strand distribution p-value = 8.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

CTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGACAGCCTGGC  >  NC_000913/3424513‑3424800
                                                                                                                                                    |                                                                                                                                           
 tctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGGGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTcccc                                                                                                                                            >  3:211902/1‑148 (MQ=17)
   tttAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAATGGTGGTATTTCAAGGTAGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAcc                                                                                                                                          >  8:136694/1‑149 (MQ=17)
   tttAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAcc                                                                                                                                          >  2:38001/1‑149 (MQ=32)
   tttAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAcc                                                                                                                                          >  6:177909/1‑149 (MQ=32)
       ggaggaGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATc                                                                                                                                      <  2:271844/149‑1 (MQ=32)
        gaggagACCGCCCCAGTCATACTTCCCACAAGTCACTGTCCGCAACCCGGATTACGTGTCAACGGTAGAACATCAATCATTAAAGGGTGGTATTTCAAGGGCGGCTCCATGCAGACTGGCGGCCACACTTCTAAGCCTCCCACCTATcc                                                                                                                                     <  2:155089/149‑1 (MQ=11)
           gagaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACCTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTac                                                                                                                                  <  3:111373/149‑1 (MQ=17)
              aCCGCCCCAGTCAAACTACCCCCCAGACACTGTCCGCAACCCGGATTACGGGTCAACGGTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTAcaca                                                                                                                               >  5:367348/1‑149 (MQ=17)
                cGCCCCAGTCAACCTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCCCCTATCCAACACAt                                                                                                                              >  4:144285/1‑148 (MQ=12)
                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)
                         cAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCTACCACACTACTAACCCTCCCACCTATCCTACACATCAAGGCTCa                                                                                                                     >  4:226329/1‑148 (MQ=255)
                          aaaCTACCCTCCAGACAATGTCAGCAACGCTGATTACGGGTCAACGTTAGAACATCAAACATTAAAGTGTGGTATTTCAAGGTCGTCTCCATGCATACTGGAGTCCACAATTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCaa                                                                                                                    <  3:309565/148‑1 (MQ=11)
                          aaaCTACCCACCAGACACTGTCCGCAACGCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTTGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAAt                                                                                                                   <  7:350235/149‑1 (MQ=17)
                          aaaCTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTTGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCaa                                                                                                                    <  7:182526/148‑1 (MQ=17)
                          aaaCTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCATATTTTTTTTAACAATTTTATGCGTCCAACCTTATCTACACATCAAGGCTCAAt                                                                                                                   <  2:259369/149‑1 (MQ=11)
                          aaaCTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAAt                                                                                                                   <  1:235514/149‑1 (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)
                                  caccaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCCCCTATCCTACACATCAAGGCTCAATGTTCAgtg                                                                                                           >  5:294862/1‑149 (MQ=17)
                                  caccaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAgtg                                                                                                           >  2:161956/1‑149 (MQ=32)
                                        acacTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGc                                                                                                     >  4:51205/1‑149 (MQ=32)
                                         cacTGTCAGAAACCCTGATTACGGGTCAACGTTAGAACATAAAACATTAAAGGGTTGTATTTCAAGGTCTGCTCCATTCAGATTTTTTTTTTAATTTTTATGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCt                                                                                                    <  6:367348/149‑1 (MQ=2)
                                            tGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGAGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCta                                                                                                   >  5:155484/1‑147 (MQ=17)
                                              tCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCtata                                                                                                 <  7:357173/147‑1 (MQ=32)
                                                 gCAACCCGGATTACGGGTCAACGTTAGAACATCACACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTCGCGTCCACACTTCTAAGCCTCCCACCAATCCTACACATCAAGGCTCAATGTACAGTGTCAAGCAATAGTaa                                                                                             >  4:4984/1‑148 (MQ=11)
                                                 gCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTACGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTa                                                                                              >  5:91631/1‑147 (MQ=17)
                                                       cGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGt                                                                                                                 <  1:314291/122‑1 (MQ=32)
                                                       cGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGt                                                                                                                 >  2:314291/1‑122 (MQ=32)
                                                        ggATTACGGGTCAACGTTAGAACATCAAACATTACAGTGTGGTATTTCAAGGTCGGCGCCAGGCAGACTGGCGGCCACACTTCTAAGCCTCCCACCTATCCTACACTTCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCAc                                                                                     <  3:4984/149‑1 (MQ=11)
                                                        ggATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCAc                                                                                     <  3:233659/149‑1 (MQ=32)
                                                           ttACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGATTTTTTTTTAAATTTTTATGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACggg                                                                                  <  6:310873/149‑1 (MQ=11)
                                                             ccGGGTCAACGTTAGACCATCAAACATTAAAGGGTGGTATTGCAAGGTCGGCTCCATGCAGAATGTCGTCAACACTTCTAAGCCTCCAACCTATCCTAAACATCAATGCTCAATGTTAAGTGTCAAGCTATAGTAAAGGTTCACggg                                                                                  <  6:91631/146‑1 (MQ=11)
                                                              cGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTc                                                                               >  1:114629/1‑149 (MQ=32)
                                                               gggTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCt                                                                              >  1:46198/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                                                                                                              >  6:276609/1‑109 (MQ=32)
                                                                           gAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGcc                                                                  >  5:233783/1‑149 (MQ=32)
                                                                           gAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCAACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGcc                                                                  >  1:319157/1‑149 (MQ=17)
                                                                            aaCATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCcg                                                                 >  4:338234/1‑149 (MQ=32)
                                                                              cATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTAAACATCAAGGCTCAATGTGCAGTGGCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCcgcg                                                               >  5:80278/1‑149 (MQ=12)
                                                                               aTCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGtt                                                                                                                <  6:80554/99‑1 (MQ=33)
                                                                               aTCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGtt                                                                                                                >  5:80554/1‑99 (MQ=33)
                                                                                tCAAACATTAAAGGGTGGTATTTCAAGTTCGTCTCCATGCAGATTTTTTTTTTAATTTTTATGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCt                                                                      <  6:415012/140‑1 (MQ=2)
                                                                                tCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCATCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCggg                                                             <  1:285332/149‑1 (MQ=17)
                                                                                tCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTAAAGCTATAGTAAAGGTTCACGGGGTCTTTCAGTCt                                                                      >  5:415012/1‑140 (MQ=17)
                                                                                    aCATTAAAGGGTGGTATTTCAAGGTCGGCTACATGCAGTATTTATTTTAAATTTTTATGCGTCAAACGTATTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTaca                                                         <  8:214545/149‑1 (MQ=11)
                                                                                     cATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTGCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTacac                                                        <  5:48546/149‑1 (MQ=17)
                                                                                      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)
                                                                                             ggtggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCta                                                                                                   <  5:275153/98‑1 (MQ=33)
                                                                                             ggtggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCtt                                                <  1:419061/149‑1 (MQ=32)
                                                                                              gtggtATTTCAAGGTCGGCTCCATGCAGATTTTTTTTTTAATTTTTATGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTc                                               <  4:350787/149‑1 (MQ=11)
                                                                                               tggtATTTCAAGGTCGGCTCCATGCAGACTGGCTTCCAAACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAATCTATAGTAAAGGTTAACGGGGTATTTCAGACTTTCCTCGTGTACACTGCATAttaa                                              >  4:215822/1‑147 (MQ=11)
                                                                                               tggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTca                                              <  1:275110/149‑1 (MQ=32)
                                                                                                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)
                                                                                                    tttCAAGGTCGGCTCCATGCAGATTTTTTTTTTAATTTTTATGCGTCCAACGTTTTCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTc                                               <  4:285388/143‑1 (MQ=11)
                                                                                                    tttCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATACTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGAGGGTACACTGCATCTTc                                               >  3:285388/1‑143 (MQ=17)
                                                                                                    tttCAAGGGCGGCTCCATGCAGTCTGGGGTCGACAGTTCTAAGCCTCCCACGTATCCTACACATCAATGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                                          <  4:200372/148‑1 (MQ=11)
                                                                                                     ttCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCATTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGGCTGGCCGCGCGTACACTGCATCTTCACAGc                                          >  5:56422/1‑147 (MQ=255)
                                                                                                      tCAAGGTCGGCTCCATGCAGACTGGAGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGa                                        <  6:304541/148‑1 (MQ=17)
                                                                                                        aaGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGTCGCGGGTACACTGCATCTTCACAGCGAGtt                                     <  2:101996/149‑1 (MQ=17)
                                                                                                        aaGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCAAGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGt                                      >  6:353586/1‑148 (MQ=17)
                                                                                                         aGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTc                                    <  1:38001/149‑1 (MQ=32)
                                                                                                           gTCGGCTCCATGCAGACTGGCGTCCACAATTCTTAGCCTCCCACGTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCaa                                  <  4:211902/149‑1 (MQ=12)
                                                                                                            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)
                                                                                                                   cATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACt                           <  8:166638/148‑1 (MQ=32)
                                                                                                                      gCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGt                       >  4:366294/1‑149 (MQ=32)
                                                                                                                       cAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGATATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGt                       <  3:226329/148‑1 (MQ=17)
                                                                                                                          aCTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCg                   >  7:124476/1‑149 (MQ=32)
                                                                                                                           cTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCgg                  <  6:7106/149‑1 (MQ=32)
                                                                                                                              gcttCCACAATTTTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGGGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGTGTTCAATTTCACTGAGTCTCGGGTg               <  2:300778/146‑1 (MQ=11)
                                                                                                                                gTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGga             <  5:353586/149‑1 (MQ=32)
                                                                                                                                gTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGga             <  3:354631/149‑1 (MQ=32)
                                                                                                                                 tCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTg               >  3:63449/1‑146 (MQ=32)
                                                                                                                                   cacaCTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGAc          <  3:338234/149‑1 (MQ=32)
                                                                                                                                      acTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                                          <  3:228014/114‑1 (MQ=32)
                                                                                                                                        ttCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGa                                        <  4:200233/114‑1 (MQ=32)
                                                                                                                                        ttCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGa                                        >  3:200233/1‑114 (MQ=32)
                                                                                                                                           tAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGACAGCCTGGc  >  5:336636/1‑149 (MQ=11)
                                                                                                                                                    |                                                                                                                                           
CTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGAGACAGCCTGGC  >  NC_000913/3424513‑3424800

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: