Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,946,587 G→C 9.8% P312R (CCG→CGG) ‡ yijP ← phosphoethanolamine transferase CptA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,946,5870GC 9.8% 172.6 / 8.5 60P312R (CCG→CGG) ‡yijPphosphoethanolamine transferase CptA
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base C (4/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GTTGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTG  >  NC_002947/2946445‑2946730
                                                                                                                                              |                                                                                                                                                 
gTTGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCg                                                                                                                                              <  5:187405/148‑1 (MQ=255)
 ttGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCt                                                                                                                                            <  3:71550/149‑1 (MQ=255)
     gcTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑g                                                                                                                                        >  4:278281/1‑148 (MQ=255)
       tGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGGGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGa                                                                                                                                      >  6:96966/1‑148 (MQ=255)
         gttgttCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAAc                                                                                                                                   >  5:164421/1‑149 (MQ=255)
                  gTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGt                                                                                                                          >  5:317056/1‑149 (MQ=255)
                  gTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGt                                                                                                                          >  3:38864/1‑149 (MQ=255)
                                   tttGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAAc                                                                                                        >  4:325738/1‑149 (MQ=255)
                                    ttGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGAGCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAAc                                                                                                        <  2:6781/148‑1 (MQ=255)
                                       ctgcGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGa                                                                                                    <  6:317056/149‑1 (MQ=255)
                                         gcGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑t                                                                                                                      <  7:5225/129‑1 (MQ=255)
                                         gcGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑t                                                                                                                      >  8:5225/1‑129 (MQ=255)
                                         gcGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAgg                                                                                                  <  6:92711/149‑1 (MQ=255)
                                          cGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAgg                                                                                                  <  3:278281/148‑1 (MQ=255)
                                                ggtgggGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATctgct                                                                                           <  2:340569/149‑1 (MQ=255)
                                                 gtggtgAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATctgctg                                                                                          >  5:247306/1‑149 (MQ=255)
                                                  tggtgAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTcc                                                                                                                              >  7:221639/1‑113 (MQ=255)
                                                  tggtgAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTcc                                                                                                                              <  8:221639/113‑1 (MQ=255)
                                                         cATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACtt                                                                                  >  3:362098/1‑149 (MQ=255)
                                                          aTGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATctgct                                                                                           <  1:125499/139‑1 (MQ=255)
                                                          aTGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATctgct                                                                                           >  2:125499/1‑139 (MQ=255)
                                                             gtgtTGCGCTTGGTCTAGGTCTGCTGGTTGGTGAGCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGAt                                                                              <  5:240565/149‑1 (MQ=255)
                                                             gtgtTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGAt                                                                              <  8:115655/149‑1 (MQ=255)
                                                                ttGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGgt                                                                           <  4:362098/149‑1 (MQ=255)
                                                                ttGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGgt                                                                           <  4:38864/149‑1 (MQ=255)
                                                                         gTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGgcg                                                                  >  3:112202/1‑149 (MQ=255)
                                                                            aaGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGATC‑TGCTCGTCGCCGAACGTGAGGATCTGCTCCATCACTTCGATGGTGTACGGGCgcgg                                                               >  5:203161/1‑149 (MQ=255)
                                                                                tCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGa                                                            >  6:300648/1‑148 (MQ=255)
                                                                                    ctgGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACTGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATa                                                          <  5:300648/146‑1 (MQ=255)
                                                                                    ctgGTTGGTGATCCAGAAGCTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATa                                                          <  2:135301/146‑1 (MQ=255)
                                                                                            tgttCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGTTCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCgg                                               <  5:329713/146‑1 (MQ=255)
                                                                                            tGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCgg                                               <  6:247306/149‑1 (MQ=255)
                                                                                              aTCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGaa                                             >  6:144670/1‑149 (MQ=255)
                                                                                                 cAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAacac                                          >  6:214732/1‑149 (MQ=255)
                                                                                                       gTCTTCTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGACGCTGA‑GGAAACGATCCGACTTC‑TGCTCATCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACCGGCGCGGCGATACCACGTCGCGGAACACGCTCAg                                    >  5:267686/1‑149 (MQ=255)
                                                                                                          ttGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATc                                                                                               >  5:150333/1‑87 (MQ=255)
                                                                                                          ttGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATc                                                                                               <  6:150333/87‑1 (MQ=255)
                                                                                                                    cctgTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTg                       >  3:73373/1‑149 (MQ=255)
                                                                                                                    cctaTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTg                       >  8:98860/1‑149 (MQ=255)
                                                                                                                     ctgTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTgg                      <  4:73373/149‑1 (MQ=255)
                                                                                                                       gTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTCGTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTggc                     >  1:291507/1‑149 (MQ=255)
                                                                                                                       gTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTggcg                    >  7:157978/1‑149 (MQ=255)
                                                                                                                        tttCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTggcgg                   >  1:33580/1‑149 (MQ=255)
                                                                                                                           cagcagGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCa                >  8:357445/1‑149 (MQ=255)
                                                                                                                           cagcagGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATAACACGTCGCGGAACACGGTCAGGCTCTTGTCGCTggcggc                  >  3:75080/1‑147 (MQ=255)
                                                                                                                            agcagGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACg                                         >  8:365831/1‑123 (MQ=255)
                                                                                                                            agcagGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACg                                         <  7:365831/123‑1 (MQ=255)
                                                                                                                             gcagGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAgc              <  7:357445/149‑1 (MQ=255)
                                                                                                                              cagGTTGATCAGCGACGGGTCGCTGAGGGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAgc              >  1:4603/1‑149 (MQ=255)
                                                                                                                               cgtttcctCAGCGACCCGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAgcgc            >  3:355624/8‑149 (MQ=255)
                                                                                                                                  ttGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTc         <  2:6326/149‑1 (MQ=255)
                                                                                                                                  ttGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTc         <  8:157978/149‑1 (MQ=255)
                                                                                                                                    gATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCg        <  4:75080/148‑1 (MQ=255)
                                                                                                                                      tCAGCGACCCGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTggc                     >  3:313877/1‑133 (MQ=255)
                                                                                                                                      tCAGCGACCCGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTggc                     <  4:313877/133‑1 (MQ=255)
                                                                                                                                      tCAGCGACCCGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTggc                     <  4:352209/133‑1 (MQ=255)
                                                                                                                                      tCAGCGACCCGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGAACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTggc                     >  3:352209/1‑133 (MQ=255)
                                                                                                                                        aGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGtt   <  5:96966/149‑1 (MQ=255)
                                                                                                                                        aGCGACCCGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTTCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGtt   >  8:59332/1‑149 (MQ=255)
                                                                                                                                         gCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTg  >  6:313479/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GTTGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTG  >  NC_002947/2946445‑2946730

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