Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,112,011 T→C 5.5% D134G (GAC→GGC)  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,0110TC 5.5% 220.7 / 2.1 73D134G (GAC→GGC) rluB23S rRNA pseudouridylate synthase B
Reads supporting (aligned to +/- strand):  ref base T (35/34);  new base C (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GTGAACTTGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTCGCCGTCGGGTTTGTT  >  NC_002947/5111867‑5112159
                                                                                                                                                |                                                                                                                                                    
gTGAACTTGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGTATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCggt                                                                                                                                                  >  8:219187/1‑149 (MQ=255)
 tGAACTTGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCggtg                                                                                                                                                 >  5:141075/1‑149 (MQ=255)
   aaCTTGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCggtggt                                                                                                                                               <  7:82500/149‑1 (MQ=255)
      ttGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCgg                                                                                                                                                   <  5:163209/142‑1 (MQ=255)
      ttGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCgg                                                                                                                                                   >  6:163209/1‑142 (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                                                                            <  1:48454/141‑1 (MQ=255)
                                                                              acacGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGc                                                                            >  2:48454/1‑141 (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                                                                      <  2:217877/109‑1 (MQ=255)
                                                                                                                    gTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGc                                                                      >  1:217877/1‑109 (MQ=255)
                                                                                                                     tGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACg                             >  1:152282/1‑149 (MQ=255)
                                                                                                                       cATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACggg                           <  3:59135/149‑1 (MQ=255)
                                                                                                                        aTCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGt                          >  7:183820/1‑149 (MQ=255)
                                                                                                                           aGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCgg                                      <  2:84771/134‑1 (MQ=255)
                                                                                                                           aGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCgg                                      >  1:84771/1‑134 (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                            >  4:169768/1‑135 (MQ=255)
                                                                                                                                    gCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACgg                            <  3:169768/135‑1 (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)
                                                                                                                                                |                                                                                                                                                    
GTGAACTTGGCCGGGCCGTCCTCCAGCATTACGCCGGCTTTCAGGCGCTCGATCATCTCGTCGTCGACTTCACCACGCACACGTACCGCGTACTCACGGTCCATCTCGTAGGATGGGTGCATCAGGCGGTTGGCCAGTTCACCGTCGGTGGTGAACAGCAGCAAGCCAGTGGTGTTGATGTCGAGGCGGCCGATGTTGATCCAGCGGCCTTCTTTTGGCCGTGGCAAGCGGTCGAACACGGTCGGGCGCCCTTCCGGGTCGTCACGGGTGCAGATTTCGCCGTCGGGTTTGTT  >  NC_002947/5111867‑5112159

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: