Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,439,309 G→A 7.0% P16L (CCT→CTT)  rluE ← 23S rRNA pseudouridine(2457) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,439,3090GA 7.0% 168.5 / 3.1 57P16L (CCT→CTT) rluE23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base G (22/31);  new base A (2/2);  total (24/33)
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.89e-01

GTATACGCCCGGAATGTCGATGAAGTCCTTGAGCGTCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAGCTCGGTGCGCAGCAGATCGAGGGCAGCAATGCCGATCTT  >  NC_002947/5439161‑5439457
                                                                                                                                                    |                                                                                                                                                    
gTATACGCCCGGAACGTCGATGAAGTCCTTGAGCGTCGGGCGCCCTGCCCCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTCAACAGGATCACGCGCGGCTCGGGGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAg                                                                                                                                                      <  6:33332/149‑1 (MQ=255)
       cccGGAATGTCGATGAAGTCCTTGAGCGTCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCgcgc                                                                                                                                               <  2:61904/149‑1 (MQ=255)
            aaTGTCGATGAAGTCCTTGAGCGTCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGATCGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTggg                                                                                                                                          <  2:248445/149‑1 (MQ=255)
             aTGTCGATGAAGTCGTTGCGCGTCGCGCGCCCTTCCCCGGCGCTCACCTGGGTCAGCACGTCTACCGGTTTGTTCAACAGTATCCGGCGCGGCGCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTgggg                                                                                                                                         <  4:119523/149‑1 (MQ=255)
              tGTCGATGAAGTCCTTGAGCGTCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTggggg                                                                                                                                        <  6:374572/149‑1 (MQ=255)
                  gacgAAGGCCTTGCGCGGCGCGCGCCCTTCACCGTCGCTGAACGGGGTCAGCACGTCGAACGGTTGGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGggcggc                                                                                                                                    <  6:314289/146‑1 (MQ=255)
                  gATGAAGTCCTTGAGCGTCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGggcggc                                                                                                                                    >  8:83128/1‑149 (MQ=255)
                  gATGAAGTCCTTGAGCGTCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGggcgg                                                                                                                                     <  6:191502/148‑1 (MQ=255)
                       aGTCCTTGAGCGTCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGgcggcgcggc                                                                                                                               >  8:336990/1‑149 (MQ=255)
                         tcgttGCGCGCCGCGCGCCCTTCCCCGTCGGTGACGTGGGTCCGCACGTCTACCGGTTTGTTGAACAGGATCTGGCGGGGCTCGGGGGGTGGGGCGTTGGCCGGGCGACGGGGGCCGGCGGCTGGCCGCGCCTGGGGGCGGCGCGGCtt                                                                                                                             <  4:171741/146‑1 (MQ=255)
                              gAGCGTCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGtt                                                                                                                          <  7:115849/147‑1 (MQ=255)
                              gAGCGTCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTAc                                                                                                                        >  8:203360/1‑149 (MQ=255)
                                   tCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTTAACAGGATCCGGCGCGGCACGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGt                                                                                                                   <  1:117214/149‑1 (MQ=255)
                                   tCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGt                                                                                                                   >  7:154982/1‑149 (MQ=255)
                                     gcgcgcCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTgg                                                                                                                 <  7:336990/149‑1 (MQ=255)
                                      cgcgcCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCtt                                                                                                                             <  3:298858/136‑1 (MQ=255)
                                      cgcgcCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCtt                                                                                                                             >  4:298858/1‑136 (MQ=255)
                                       gcgcCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTgg                                                                                                                 >  6:143211/1‑147 (MQ=255)
                                         gcCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGt                                                                                                             >  5:257112/1‑149 (MQ=255)
                                          ctcTTCCCGGTTGCTTAACTGGGTCAGAAGGTCTAAGGTTTTTTTCAACTGGATCAGGCGCGGCGCGGGGGGTGGGGGTTTGGCCTGGCGACGCGGGCGGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGt                                                                                                             <  2:323761/146‑1 (MQ=255)
                                          cccTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTg                                                                                                            <  4:289079/149‑1 (MQ=255)
                                          cccTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTg                                                                                                            <  4:283975/149‑1 (MQ=255)
                                             ttCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGag                                                                                                          >  1:323761/1‑147 (MQ=255)
                                             ctccccGTCGCTCACCTGGGGCCGCCCGTCCACCGGTTTGTTCACCCGGACCCGGCGCGGCGCGGCGGGGGGCGCCTTGGCCGGGCGCCGCGGCCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACa                                                                                                         <  8:154982/145‑1 (MQ=255)
                                                aCCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATgg                                                                                                      <  5:143211/149‑1 (MQ=255)
                                                        tGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGGGTGTTACGGGGTGGCAggg                                                                                                            >  7:148468/1‑133 (MQ=255)
                                                        tGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAg                                                                                              >  7:167838/1‑149 (MQ=255)
                                                            cTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCGGGGGGCGGCGCGGGATGTTACGGGGTGGCAGTGACATGGATCCGCAGCgg                                                                                           >  5:34552/1‑148 (MQ=255)
                                                             tGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGaa                                                                                         >  3:123195/1‑149 (MQ=255)
                                                                gTCAGCAGGTCTAACGGTTTGTTGAACAGGATCCGGCGCGGCTCGGCGGGGGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAAcgg                                                                                      <  3:151745/149‑1 (MQ=255)
                                                                 tCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACggg                                                                                                                     <  4:245629/117‑1 (MQ=255)
                                                                 tCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACggg                                                                                                                     >  3:245629/1‑117 (MQ=255)
                                                                        gTCGAAGGGTTTGTTGAACAGGAGCCGGCGGGGCGCGGCGGGGGGCGCCTTGGCGGGGCGACGCGGGCGGGGGGCTGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCAt                                                                              <  6:34552/149‑1 (MQ=255)
                                                                          cGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCGTGTTACGGGGTGGCAGTGACATGGCTCCTCAGCGGAACGGCGGCTCagc                                                                             >  4:67999/1‑146 (MQ=255)
                                                                           gAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTggg                                                                                                                                          <  3:174863/86‑1 (MQ=255)
                                                                           gAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTggg                                                                                                                                          >  4:174863/1‑86 (MQ=255)
                                                                               ggTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCt                                                                       <  3:67999/149‑1 (MQ=255)
                                                                                 tttGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTgc                                                                     <  4:123195/149‑1 (MQ=255)
                                                                                     tggacccggATCCGGCGCGGCGCGGCGGGTGGCGCCTGGGCCGGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGGGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCa                                                                  <  2:232144/141‑1 (MQ=255)
                                                                                         aCAGGATCCGGCGCGGCTCGGCGGGGGGCGCCTGGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTa                                                             <  4:132681/149‑1 (MQ=255)
                                                                                                  ggCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAAGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTg                                                     >  4:149469/1‑148 (MQ=255)
                                                                                                  ggCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAAGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTg                                                     >  4:155014/1‑148 (MQ=255)
                                                                                                          tCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGtt                                            <  8:167838/149‑1 (MQ=255)
                                                                                                             gcggGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGa                                                        <  7:110527/134‑1 (MQ=255)
                                                                                                             gcggGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGa                                                        >  8:110527/1‑134 (MQ=255)
                                                                                                             gcggGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAg                                         <  7:83128/149‑1 (MQ=255)
                                                                                                                    gcgcCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGt                                                      <  2:190614/129‑1 (MQ=255)
                                                                                                                    gcgcCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGt                                                      >  1:190614/1‑129 (MQ=255)
                                                                                                                      gcCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCGGGGGGCGGCGCGGCGTGTTACGGGGTGGGAGTGACATGGATCCTCAGCGGCACGGCGGCGCATAGAAGCTGCGCAGTTTTCGCGAGTGCAGCGAGTTGAGCTCGGTgcg                                >  4:251548/1‑149 (MQ=255)
                                                                                                                             ccTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATACTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAGCTCGGTGCgcag                             >  6:327266/1‑145 (MQ=255)
                                                                                                                                 gCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCa                                                                  <  8:1141/104‑1 (MQ=255)
                                                                                                                                  cGACGCGGGCCGGCGGCAAGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAGCTCGGTGCGCAGCAGATCGAg                    <  3:155014/149‑1 (MQ=255)
                                                                                                                                  cGACGCGGGCCGGCGGCAAGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAGCTCGGTGCGCAGCAGATCGAg                    <  3:149469/149‑1 (MQ=255)
                                                                                                                                      gcgGGCCGGCGGCAGGCCGCGCCTGGGGGCNGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAGCTCGGTGCGCAGCAGATCGAggg                  >  1:31914/1‑147 (MQ=255)
                                                                                                                                        gggCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAGCTCGGTGCGCAGCAGATCGAGGgcag               >  7:23787/1‑148 (MQ=255)
                                                                                                                                           ccGGCGGCCGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAGCTCGGTGCGCAGCAGATCGAGGGCAGCAAt           <  3:251548/149‑1 (MQ=255)
                                                                                                                                              gcggcAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAGCTCGGTGCGCAGCAGATCGAGGgca                <  5:366970/141‑1 (MQ=255)
                                                                                                                                              gcggcAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAGCTCGGTGCGCAGCAGATCGAGGgca                >  6:366970/1‑141 (MQ=255)
                                                                                                                                                    ggCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAGCTCGGTGCGCAGCAGATCGAGGGCAGCAATGCCGATCtt  >  4:270586/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTATACGCCCGGAATGTCGATGAAGTCCTTGAGCGTCGCGCGCCCTTCACCGTCGCTGAACTGGGTCAGCACGTCGAACGGTTTGTTGAACAGGATCAGGCGCGGCTCGGCGGGTGGCGCCTTGGCCTGGCGACGCGGGCCGGCGGCAGGCCGCGCCTGGGGGCGGCGCGGCTTGTTACGGGGTGGCAGTGACATGGATCCTCAGCGGAACGGCGGCTCATCGAAGCTGCGCAGTTTACGCGAGTGCAGCGAGTTGAGCTCGGTGCGCAGCAGATCGAGGGCAGCAATGCCGATCTT  >  NC_002947/5439161‑5439457

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: ‑