Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,482,292 A→G 6.3% Q73R (CAG→CGG)  dusB → tRNA dihydrouridine synthase DusB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,482,2920AG 6.3% 182.6 / 2.5 63Q73R (CAG→CGG) dusBtRNA dihydrouridine synthase DusB
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base G (2/2);  total (31/32)
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.48e-01

GGACCAGCCTTTCCGTACACTTTGCCAGCGCCTGGGCGCCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAGGCTCTGCTTTATTGAGAGATG  >  NC_002947/5482146‑5482435
                                                                                                                                                  |                                                                                                                                                 
ggACCAGCCTTTCCGTACACTTTGCCAGCGCCTGGGCGCCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCGGa                                                                                                                                                 >  8:105579/1‑149 (MQ=255)
         tttCCGTACACTTTGCCAGCGCCTGGGCGCCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCcggcg                                                                                                                                        >  3:79244/1‑149 (MQ=255)
              gTACACTTTGCCAGCGCCTGGGCGCCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGAt                                                                                                                                   <  7:280965/149‑1 (MQ=255)
                 cacTTTGCCAGCGCCTGGGCGCCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATgc                                                                                                                                 <  7:243499/148‑1 (MQ=255)
                    tttGCCAGCGCCTGGGCGCCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCcgg                                                                                                                                          >  1:143491/1‑136 (MQ=255)
                    tttGCCAGCGCCTGGGCGCCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATgcg                                                                                                                                >  2:97992/1‑146 (MQ=255)
                          aGCGCCTGGGCGCCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGatgatg                                                                                                                       <  4:28602/149‑1 (MQ=255)
                           gcgcCTGGGCGCCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATgg                                                                                                                      <  7:277524/149‑1 (MQ=255)
                             gcCTGGGCGCCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGCGCCCCGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGGTGATGGCt                                                                                                                    >  8:322333/1‑149 (MQ=255)
                                     gcCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGcggcgg                                                                                                             >  7:107730/1‑148 (MQ=255)
                                        ggTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCaaa                                                                                                        <  8:96753/149‑1 (MQ=255)
                                         gTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAg                                                                                                       <  4:79244/149‑1 (MQ=255)
                                              gtggtgTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCTTCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCaa                                                                                                  >  3:46418/1‑149 (MQ=255)
                                               tggtgTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCGGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAAt                                                                                                 <  7:105579/149‑1 (MQ=255)
                                                 gtgtCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGt                                                                                              >  4:264021/1‑149 (MQ=255)
                                                  tgtCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTc                                                                                             >  5:334101/1‑149 (MQ=255)
                                                   gtCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCg                                                                                            <  3:264021/149‑1 (MQ=255)
                                                         aaTTGTTCGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGACGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCgg                                                                                      <  4:229281/149‑1 (MQ=255)
                                                             gTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGt                                                                                              >  4:291677/1‑137 (MQ=255)
                                                             gTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGt                                                                                              <  3:291677/137‑1 (MQ=255)
                                                                    gcGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGat                                                                            >  6:30882/1‑148 (MQ=255)
                                                                        cATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCACAAGGCCAATGTCGAGGCGGGTGCCCAGatcatc                                                                        <  1:48491/149‑1 (MQ=255)
                                                                         aTGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGa                                                                      >  1:20091/1‑149 (MQ=255)
                                                                         aTGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGa                                                                      >  1:25741/1‑149 (MQ=255)
                                                                             gCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATc                                                                  <  6:12547/149‑1 (MQ=255)
                                                                                tgtgGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGatcatc                                                                        <  6:35378/140‑1 (MQ=255)
                                                                                tgtgGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGatcatc                                                                        >  5:35378/1‑140 (MQ=255)
                                                                                  tgGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAAc                                                               <  3:323920/147‑1 (MQ=255)
                                                                                    gAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATgg                                                           >  1:172936/1‑149 (MQ=255)
                                                                                     aaCAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATgcg                                                                                                                                >  7:170569/1‑81 (MQ=255)
                                                                                     aaCAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATgcg                                                                                                                                <  8:170569/81‑1 (MQ=255)
                                                                                          ccGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTc                                                     <  8:107730/149‑1 (MQ=255)
                                                                                            gCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCg                                                   >  6:268352/1‑149 (MQ=255)
                                                                                                  cGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCaaa                                              >  1:299327/1‑148 (MQ=255)
                                                                                                     gCCTGTGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCaaaaa                                            >  3:131453/1‑147 (MQ=255)
                                                                                                          tgccgcATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTg                                     <  4:131453/148‑1 (MQ=255)
                                                                                                          cgccgcATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCt                                      <  5:268352/148‑1 (MQ=255)
                                                                                                          cgccgcATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTg                                     <  2:20091/149‑1 (MQ=255)
                                                                                                          cgccgcATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTg                                     <  2:25741/149‑1 (MQ=255)
                                                                                                           gccgcATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGc                                    <  5:30882/149‑1 (MQ=255)
                                                                                                            acgcATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGc                                    >  7:232133/2‑148 (MQ=255)
                                                                                                                aTCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATggg                                                          <  5:329395/122‑1 (MQ=255)
                                                                                                                aTCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATggg                                                          >  6:329395/1‑122 (MQ=255)
                                                                                                                  ccATGAAGGTGATCCCGAGCCACGCTCGGTGCGGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATc                                                                  <  6:191067/112‑1 (MQ=255)
                                                                                                                  ccATGAAGGTGATCCCGAGCCACGCTCGGTGCGGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATc                                                                  >  5:191067/1‑112 (MQ=255)
                                                                                                                       aaGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCa                                                                 >  7:9513/1‑108 (MQ=255)
                                                                                                                       aaGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCa                                                                 <  8:9513/108‑1 (MQ=255)
                                                                                                                        aGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAg                       <  2:299327/149‑1 (MQ=255)
                                                                                                                          gTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCg                                                                                       <  2:159532/83‑1 (MQ=255)
                                                                                                                          gTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCg                                                                                       >  1:159532/1‑83 (MQ=255)
                                                                                                                             aTCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAGGctct                  >  2:332799/1‑149 (MQ=255)
                                                                                                                               cccGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAGGCTCTGc                >  6:203430/1‑149 (MQ=255)
                                                                                                                                ccGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAGGCTCTGc                <  2:172936/148‑1 (MQ=255)
                                                                                                                                ccGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAGGCTCTGCt               >  3:300360/1‑149 (MQ=255)
                                                                                                                                   aGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCaaaaaaa                                          <  4:3352/119‑1 (MQ=255)
                                                                                                                                   aGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCaaaaaaa                                          >  3:3352/1‑119 (MQ=255)
                                                                                                                                     ccACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAGGCTCTGCTTTAt           >  1:107088/1‑149 (MQ=255)
                                                                                                                                       aCGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAGGCTCTGCTTTAttag        >  7:302724/1‑147 (MQ=255)
                                                                                                                                       aCGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAGGCTCTGCTTTAt           <  8:302724/146‑1 (MQ=255)
                                                                                                                                             ggTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAGGCTCTGCTTTATTGAGAGAt   <  4:300360/148‑1 (MQ=255)
                                                                                                                                             ggTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAGGCTCTGCTTTATTGAGAGATg  <  7:322333/149‑1 (MQ=255)
                                                                                                                                               tGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGatcatc                                                                        >  4:194366/1‑77 (MQ=255)
                                                                                                                                               tGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGatcatc                                                                        <  3:194366/77‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGACCAGCCTTTCCGTACACTTTGCCAGCGCCTGGGCGCCGGTATGGTGGTGTCGGAAATGGTTAGCAGCGACATGAGCCTGTGGAACAGCCGCAAGTCGAGCCTGCGCCGCATCCATGAAGGTGATCCCGAGCCACGCTCGGTGCAGATCGCCGGCGGCGATGCGCAGATGATGGCTGCGGCGGCAAAGGCCAATGTCGAGGCGGGTGCCCAGATCATCGATATCAACATGGGCTGTCCGGCAAAAAAAGTCTGCAACAAAGCCGCAGGCTCTGCTTTATTGAGAGATG  >  NC_002947/5482146‑5482435

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: