Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,112,015 T→C 5.8% T133A (ACC→GCC)  rluB ← 23S rRNA pseudouridylate synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,112,0150TC 5.8% 206.0 / 2.3 69T133A (ACC→GCC) rluB23S rRNA pseudouridylate synthase B
Reads supporting (aligned to +/- strand):  ref base T (32/33);  new base C (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GTGAACTTGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTCGCCGTCGGGTTTGTTGTA  >  NC_002947/5111867‑5112162
                                                                                                                                                    |                                                                                                                                                   
gTGAACTTGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGTATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCggt                                                                                                                                                     >  8:219187/1‑149 (MQ=255)
 tGAACTTGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCggtg                                                                                                                                                    >  5:141075/1‑149 (MQ=255)
   aaCTTGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCggtggt                                                                                                                                                  <  7:82500/149‑1 (MQ=255)
          ccgggccgTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAAcagc                                                                                                                                           <  1:234477/149‑1 (MQ=255)
          ccgggccgTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCACCGACGGTGAActgg                                                                                                                                           <  8:222310/149‑4 (MQ=255)
             ggccgTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAAcagcag                                                                                                                                         >  6:117899/1‑148 (MQ=255)
               ccgTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCaa                                                                                                                                      <  4:28792/149‑1 (MQ=255)
               ccatcctccAGCATTACGCCGGCTTTCAGGCGCGCGATCCTCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGCTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCaa                                                                                                                                      <  8:295853/146‑1 (MQ=255)
                cgTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCACCGACGGTGAAc                                                                                                                                              <  3:162348/140‑1 (MQ=255)
                cgTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCACCGACGGTGAAc                                                                                                                                              >  4:162348/1‑140 (MQ=255)
                 gTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGc                                                                                                                                    <  4:240938/149‑1 (MQ=255)
                     tccAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAgt                                                                                                                                >  7:118216/1‑149 (MQ=255)
                     tccAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAgt                                                                                                                                >  2:281209/1‑149 (MQ=255)
                       cAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAgtg                                                                                                                               <  1:281209/148‑1 (MQ=255)
                       cAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAgtg                                                                                                                               <  5:15397/148‑1 (MQ=255)
                                 ccGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACAGGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTTAACAGCAGCAAGCCAGTGGTGTTGATGTc                                                                                                                    <  6:287650/149‑1 (MQ=255)
                                 ccGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTc                                                                                                                    >  4:208614/1‑149 (MQ=255)
                                  cGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAAca                                                                                                                                             <  4:198279/123‑1 (MQ=255)
                                          aGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGcc                                                                                                           >  4:112775/1‑149 (MQ=255)
                                               cTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCTGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGt                                                                                                      <  1:212229/149‑1 (MQ=255)
                                               cTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGGTGATGTCGAGGCGGCCGATGt                                                                                                      <  8:118216/149‑1 (MQ=255)
                                                 cGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTg                                                                                                    <  3:208614/149‑1 (MQ=255)
                                                  gATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAgtggt                                                                                                                             >  7:86287/1‑123 (MQ=255)
                                                  gATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAgtggt                                                                                                                             <  8:86287/123‑1 (MQ=255)
                                                   atcatcTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGAt                                                                                                  >  5:13539/1‑149 (MQ=255)
                                                      atcTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGGTc                                                                                                 >  3:73214/1‑147 (MQ=255)
                                                           gtcgtcgACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAgtg                                                                                                                               >  3:95908/1‑112 (MQ=255)
                                                           gtcgtcgACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAgtg                                                                                                                               <  4:95908/112‑1 (MQ=255)
                                                               tcgACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCcttc                                                                                      <  7:219187/149‑1 (MQ=255)
                                                                   cTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCtttt                                                                                  <  5:117899/149‑1 (MQ=255)
                                                                    ttCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTg                                                                                 >  5:19058/1‑149 (MQ=255)
                                                                       accacGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCg                                                                                            <  2:293051/135‑1 (MQ=255)
                                                                       accacGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCg                                                                                            >  1:293051/1‑135 (MQ=255)
                                                                       accacGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGcc                                                                              >  6:258119/1‑149 (MQ=255)
                                                                         cacGCACACGTACCGCGTTCTCACGGTCAATATCGTAGCATGGGTGCATAAGTCGGTTGGCCAGTTAACCGTCGGTGGTGAAAAGCAGCAAGCCAGTGGTGTTGATGTCTAGGCGGCCGATGTTGATACAGCCGCCTTCATTTGGcc                                                                              >  8:211446/1‑147 (MQ=255)
                                                                         cacGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGt                                                                            <  6:19058/149‑1 (MQ=255)
                                                                            gCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGc                                                                         <  6:13539/149‑1 (MQ=255)
                                                                              acacGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGc                                                                               >  2:48454/1‑141 (MQ=255)
                                                                              acacGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGc                                                                               <  1:48454/141‑1 (MQ=255)
                                                                                  gTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCgg                                                                   >  8:3876/1‑149 (MQ=255)
                                                                                   tACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGt                                                                  <  3:112775/149‑1 (MQ=255)
                                                                                   tACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGt                                                                  >  1:10813/1‑149 (MQ=255)
                                                                                               aCGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTgg                                                                          <  5:219107/129‑1 (MQ=255)
                                                                                               aCGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTgg                                                                          >  6:219107/1‑129 (MQ=255)
                                                                                               aCGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCg                                                      >  4:59135/1‑149 (MQ=255)
                                                                                                      aTCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGcc                                                >  3:196046/1‑148 (MQ=255)
                                                                                                      aTCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCAAGCGGCCTGCTTTTGGCCGTGGCAAGCGGTCGCACACGGTCGGGCGccc                                               >  4:180698/1‑149 (MQ=255)
                                                                                                       tctcGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCt                                              >  1:129275/1‑149 (MQ=255)
                                                                                                           gTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTTGTGAACAGCAGCAAGCCAGTGGCGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTc                                            <  3:306883/147‑1 (MQ=255)
                                                                                                           gTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTcc                                           <  2:10813/148‑1 (MQ=255)
                                                                                                                    gTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGc                                                                         >  1:217877/1‑109 (MQ=255)
                                                                                                                    gTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGc                                                                         <  2:217877/109‑1 (MQ=255)
                                                                                                                     tGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACg                                >  1:152282/1‑149 (MQ=255)
                                                                                                                       cATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACggg                              <  3:59135/149‑1 (MQ=255)
                                                                                                                        aTCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGt                             >  7:183820/1‑149 (MQ=255)
                                                                                                                           aGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCgg                                         >  1:84771/1‑134 (MQ=255)
                                                                                                                           aGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCgg                                         <  2:84771/134‑1 (MQ=255)
                                                                                                                           aGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTg                            >  2:184358/1‑147 (MQ=255)
                                                                                                                                 ttGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGGTCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTc                    <  4:73214/149‑1 (MQ=255)
                                                                                                                                    tgctgTTCACCACCGACGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTCGcc                 >  2:42054/5‑149 (MQ=255)
                                                                                                                                    gCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTc                                                       >  7:85158/1‑111 (MQ=255)
                                                                                                                                    gCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTc                                                       <  8:85158/111‑1 (MQ=255)
                                                                                                                                    gCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACgg                               <  3:169768/135‑1 (MQ=255)
                                                                                                                                    gCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACgg                               >  4:169768/1‑135 (MQ=255)
                                                                                                                                     ccAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCgg                                         >  2:262792/1‑124 (MQ=255)
                                                                                                                                     ccAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCgg                                         <  1:262792/124‑1 (MQ=255)
                                                                                                                                       aGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATAAAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTCGCCGTc              >  8:32829/1‑149 (MQ=255)
                                                                                                                                            aCCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTCGCCGTCGGGtt         >  7:156440/1‑149 (MQ=255)
                                                                                                                                            aCCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTCGCCGTCGGGtt         >  3:13217/1‑149 (MQ=255)
                                                                                                                                             ccGTCGGTGGTGAACAGCATCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTCGCCGTCGGGttt        <  4:196046/149‑1 (MQ=255)
                                                                                                                                              cGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTCCGCCGTCGGGTttg       <  6:141075/149‑1 (MQ=255)
                                                                                                                                                tCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTCGCCGTCGGGTttgtt     <  5:258119/149‑1 (MQ=255)
                                                                                                                                                   gtggtgAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTCGCCGTCGGGTTTGTTGTa  >  8:96942/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGAACTTGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTCGCCGTCGGGTTTGTTGTA  >  NC_002947/5111867‑5112162

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