Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,946,587 G→C 8.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 8.8% 187.6 / 5.7 64P312R (CCG→CGG) ‡yijPphosphoethanolamine transferase CptA
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base C (4/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

TTGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTG  >  NC_002947/2946446‑2946730
                                                                                                                                             |                                                                                                                                               
ttGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCt                                                                                                                                          >  1:57404/1‑149 (MQ=255)
   cgcTGGGTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACCCGTCGCTGat                                                                                                                                       <  4:106209/149‑2 (MQ=255)
        gttgttCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGaaa                                                                                                                                   <  1:423153/148‑1 (MQ=255)
        gttgttCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGaaa                                                                                                                                   <  1:429841/148‑1 (MQ=255)
        gttgttCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGaaa                                                                                                                                   <  2:316791/148‑1 (MQ=255)
            ttCAGGTACACCGGTGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGt                                                                                                                               <  1:161799/148‑1 (MQ=255)
            ttCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGt                                                                                                                               <  4:203466/148‑1 (MQ=255)
             tCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTcc                                                                                                                             >  8:335167/1‑149 (MQ=255)
               aGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCTTGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCgg                                                                                                                           <  4:324227/149‑1 (MQ=255)
                     acCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTg                                                                                                                     <  1:318609/149‑1 (MQ=255)
                           gcgTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCtcgtc                                                                                                               <  6:163875/149‑1 (MQ=255)
                             ggatGTTTGCTGCGAGACGGTGGTGAGCAGTGTGTTGCGCTTGGGCCAGGGCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGc                                                                                                             <  4:247000/146‑1 (MQ=255)
                                  tttGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAAc                                                                                                        >  8:314591/1‑149 (MQ=255)
                                         cGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGa                                                                                                 >  2:213218/1‑149 (MQ=255)
                                           agaAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATc                                                                                               >  4:205244/1‑149 (MQ=255)
                                           agaAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATc                                                                                               <  7:314591/149‑1 (MQ=255)
                                           agaAGGTGGTGAGCATGGTGCTGCGCTTGGTCAAGGTCGGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATc                                                                                               <  4:291479/149‑1 (MQ=255)
                                             aaGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATctg                                                                                             >  3:380877/1‑149 (MQ=255)
                                                   gtgAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCAt                                                                                       >  3:269564/1‑149 (MQ=255)
                                                      aGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCAc                                                                                    <  3:205244/149‑1 (MQ=255)
                                                        cATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACtt                                                                                  <  7:335167/149‑1 (MQ=255)
                                                              gtTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATgg                                                                            >  8:22165/1‑149 (MQ=255)
                                                               ttacgcTTGATCAAGCTCTGCTCGTTGGAGATCCACAAGGTCTTGTAGCCTGCCTGTTTCATAAGGTTGATCAGCGACCCGTCGCTGAGCAAACGGTCCGGGTTCTCCTCGTCGCCGACCCTCAGGAACTGCTGCATCACTACGAtctt                                                                           >  8:95417/4‑146 (MQ=255)
                                                                       ggTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCCCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGgc                                                                   >  8:218569/1‑149 (MQ=255)
                                                                              gTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGa                                                            >  3:374665/1‑149 (MQ=255)
                                                                              gTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCCGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGa                                                            >  5:163986/1‑149 (MQ=255)
                                                                                    tggttggtGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCAc                                                      >  1:341178/1‑149 (MQ=255)
                                                                                      gttggtGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGt                                                    >  3:211991/1‑149 (MQ=255)
                                                                                       ttggtGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAgg                                                                                                  <  3:402874/102‑1 (MQ=255)
                                                                                       ttggtGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAgg                                                                                                  >  4:402874/1‑102 (MQ=255)
                                                                                         ggtGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTcg                                                  <  2:57404/148‑1 (MQ=255)
                                                                                         ggtGATCCAGAAGGTCTTGTAGCCAGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTAGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTcgc                                                 >  2:181533/1‑149 (MQ=255)
                                                                                          gtGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTcgcg                                                <  7:22165/149‑1 (MQ=255)
                                                                                           tGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTcgcg                                                >  4:309294/1‑148 (MQ=255)
                                                                                                  gAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAAc                                                                                                        >  6:363768/1‑85 (MQ=255)
                                                                                                  gAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAAc                                                                                                        <  5:363768/85‑1 (MQ=255)
                                                                                                  gAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACgg                                        >  7:211919/1‑149 (MQ=255)
                                                                                                  gAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACgg                                        <  1:67281/149‑1 (MQ=255)
                                                                                                    aGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTc                                      <  8:279900/149‑1 (MQ=255)
                                                                                                       tCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAgg                                   <  4:269564/149‑1 (MQ=255)
                                                                                                         ttGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGc                                  >  2:177436/1‑148 (MQ=255)
                                                                                                           gTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCTTCGCCGAACGTGAGGATCTGCTTCATCACTTCGATGGTGTATGGTCTCGGCGATACCACGTCGCGGAACAGTGTCAGtctc                                <  4:380877/148‑1 (MQ=255)
                                                                                                           gTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGctct                               >  7:411603/1‑149 (MQ=255)
                                                                                                           gTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGctc                                <  4:374665/148‑1 (MQ=255)
                                                                                                           gTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGc                                  <  1:213218/146‑1 (MQ=255)
                                                                                                            tAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGctct                               <  2:341178/148‑1 (MQ=255)
                                                                                                                ctgcctgTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCg                          >  8:369604/1‑149 (MQ=255)
                                                                                                                   cctgTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCt                        <  8:411603/148‑1 (MQ=255)
                                                                                                                   cctgTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGACGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTg                       >  5:220880/1‑149 (MQ=255)
                                                                                                                    ctgTTTCAGCAGGTTGATCAGCGACGGTTCGCTGAGGAAACGGTCCGGGTGCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTAGGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTgg                      <  7:95417/149‑1 (MQ=255)
                                                                                                                       tttCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTggcgg                   >  8:391757/1‑149 (MQ=255)
                                                                                                                           agcagGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGcc                 <  6:351578/147‑1 (MQ=255)
                                                                                                                           agcagGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGcc                 >  5:351578/1‑147 (MQ=255)
                                                                                                                            gcagGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAgc              <  8:211919/149‑1 (MQ=255)
                                                                                                                            gcagGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATAACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGAGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAgc              >  5:72391/1‑149 (MQ=255)
                                                                                                                               ggTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGACCACGGTCAGGCTCTTGTCGCTGGCGGCCAGcgcg           >  4:298680/1‑149 (MQ=255)
                                                                                                                                gTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGCACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGt          >  2:117206/1‑149 (MQ=255)
                                                                                                                                    ctCAGCGACCCGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAg      >  8:167735/2‑149 (MQ=255)
                                                                                                                                    ctCAGCGACCCGTCGCTGAGGAAACGGTCCGGGGTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGt                                       <  4:60122/115‑1 (MQ=255)
                                                                                                                                    ctCAGCGACCCGTCGCTGAGGAAACGGACCGGGTTCTGCTCGTCGCCGACCGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACgct                                       >  3:60122/2‑114 (MQ=255)
                                                                                                                                    aTCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTAGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGTCGCGGCGATACCACGTCCCGGATCACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGGCGac      >  8:390366/1‑148 (MQ=255)
                                                                                                                                     tCAGCGACCCGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAgg     >  7:61148/1‑149 (MQ=255)
                                                                                                                                        gCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTg  >  8:54124/1‑149 (MQ=255)
                                                                                                                                        gCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCCCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTg  <  7:218569/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
TTGCGCTGGTTGTTCAGGTACACCGGCGCGTCTGTTTGCTGCGAGAAGGTGGTGAGCATGGTGTTGCGCTTGGTCAAGGTCTGCTGGTTGGTGATCCAGAAGGTCTTGTAGCCTGCCTGTTTCAGCAGGTTGATCAGCGACGGGTCGCTGAGGAAACGGTCCGGGTTCTGCTCGTCGCCGAACGTGAGGATCTGCTGCATCACTTCGATGGTGTACGGGCGCGGCGATACCACGTCGCGGAACACGGTCAGGCTCTTGTCGCTGGCGGCCAGCGCGTCGAGGTTG  >  NC_002947/2946446‑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: ‑