Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,946,588 G→C 10.0% P312A (CCG→GCG) ‡ yijP ← phosphoethanolamine transferase CptA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,946,5880GC10.0% 173.8 / 9.8 60P312A (CCG→GCG) ‡yijPphosphoethanolamine transferase CptA
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base C (4/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GTTGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTGGGCGTG  >  NC_002947/2946445‑2946736
                                                                                                                                               |                                                                                                                                                      
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)
                                                                                                                                               ggTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTGGGCgtg  <  2:4603/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
GTTGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGA‑GGAAACGGTCCGGGTTC‑TGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTGGGCGTG  >  NC_002947/2946445‑2946736

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