Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,322,834 A→C 5.4% S140A (TCC→GCC)  mrcB ← transglycosylase/transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,322,8340AC 5.4% 301.7 / 7.0 93S140A (TCC→GCC) mrcBtransglycosylase/transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (44/44);  new base C (3/2);  total (47/46)
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.25e-01

CGTTGCCACCAAGGTTTCGAGCAGATATTTGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGTTGGCCGCGCTTTCG  >  NC_002947/5322691‑5322979
                                                                                                                                               |                                                                                                                                                 
cGTTGCCACCAAGGTTTCGAGCAGATATTTGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAg                                                                                                                                                <  3:152040/147‑1 (MQ=255)
cGTTGCCACCAAGGTTTCGAGCAGATATTTGGGTACCTGATCGATCTTGATCAGGATGCGGCCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAgg                                                                                                                                               <  3:326401/148‑1 (MQ=255)
       accaAGGTTTCGAGCAGATATTTGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAAc                                                                                                                                       >  3:15936/1‑149 (MQ=255)
        ccaAGGTTTCGAGCAGATATTTGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACg                                                                                                                                      >  6:305767/1‑149 (MQ=255)
         caAGGTTTCGAGCAGATATTTGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACg                                                                                                                                      >  7:146884/1‑148 (MQ=255)
              tttCGAGCAGATATTTGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCg                                                                                                                                >  6:20117/1‑149 (MQ=255)
                    gCAGATATTTGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGc                                                                                                                               >  5:326016/1‑144 (MQ=255)
                    gCAGATATTTGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGc                                                                                                                               <  6:326016/144‑1 (MQ=255)
                         tatTTGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAgcg                                                                                                                     <  5:20117/149‑1 (MQ=255)
                          atTTGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGccc                                                                                                                                            <  8:253992/125‑1 (MQ=255)
                           tttGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCa                                                                                                                   >  4:196285/1‑149 (MQ=255)
                           tttGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCAGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAgcgc                                                                                                                    <  6:280679/148‑1 (MQ=255)
                               ggTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCag                                                                                                                            <  5:45183/136‑1 (MQ=255)
                               ggTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCAcgc                                                                                                                <  6:307050/148‑1 (MQ=255)
                               ggTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACCTAAACGCCac                                                                                                                            >  6:45183/1‑135 (MQ=255)
                                gTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCAcg                                                                                                                 >  8:201979/1‑146 (MQ=255)
                                gTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCAcg                                                                                                                 <  7:201979/146‑1 (MQ=255)
                                 tACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGa                                                                                                             >  1:193788/1‑149 (MQ=255)
                                  aCCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGAc                                                                                                            >  1:147631/1‑149 (MQ=255)
                                  aCCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGAc                                                                                                            >  5:47236/1‑149 (MQ=255)
                                    cTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACg                                                                                                           >  1:292047/1‑148 (MQ=255)
                                        tcgatcTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTTCCGGTGCCAACGGCAAg                                                                                                                                                >  8:183368/1‑107 (MQ=255)
                                        tcgatcTTGATCAGGATGCGCTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTTCCGGTGCCAACGGCAAg                                                                                                                                                <  7:183368/107‑1 (MQ=255)
                                            tcTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGcc                                                                                                  >  3:344256/1‑149 (MQ=255)
                                              ttGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCgg                                                                                                >  3:298053/1‑149 (MQ=255)
                                                gATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGtt                                                                                              >  3:21471/1‑149 (MQ=255)
                                                 aTCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGtt                                                                                              <  4:298053/148‑1 (MQ=255)
                                                  tCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTTCCGGTGCCAACGGCAAGGcgc                                                                                                                                            >  4:353000/1‑99 (MQ=38)
                                                   cAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCa                                                                                           >  2:91877/1‑149 (MQ=255)
                                                   cAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCa                                                                                           >  7:306667/1‑149 (MQ=255)
                                                    aGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCa                                                                                           <  1:293677/148‑1 (MQ=255)
                                                    aGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCAt                                                                                          <  4:15936/149‑1 (MQ=255)
                                                           ggTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCg                                                                                   >  3:79608/1‑149 (MQ=255)
                                                               ttCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAg                                                                                >  8:272449/1‑148 (MQ=255)
                                                               ttCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAg                                                                                >  7:275227/1‑148 (MQ=255)
                                                                 cGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGaa                                                                              <  3:8385/148‑1 (MQ=255)
                                                                   agaTTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTg                                                                           >  6:100102/1‑149 (MQ=255)
                                                                   agaTTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTg                                                                           <  2:292047/149‑1 (MQ=255)
                                                                    gaTTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTgg                                                                          <  7:272449/149‑1 (MQ=255)
                                                                      tttttCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGaa                                                                        <  3:68030/149‑1 (MQ=255)
                                                                            gggTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGATCCATGCCCTCGTAGAACTGGAAGCCCCg                                                                  >  4:296330/1‑149 (MQ=255)
                                                                             ggTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACgg                                                                 <  6:47236/149‑1 (MQ=255)
                                                                             ggTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACgg                                                                 <  5:100102/149‑1 (MQ=255)
                                                                             ggTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACgg                                                                 <  8:275227/149‑1 (MQ=255)
                                                                             ggTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACgg                                                                 <  2:147631/149‑1 (MQ=255)
                                                                              gTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACggg                                                                >  8:329285/1‑149 (MQ=255)
                                                                                aGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGATAGGCCCGCAACGTAATCGCCAGAGAAGCGCGCGCGGCCGTACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTa                                                              >  5:394436/1‑149 (MQ=255)
                                                                                   tGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCATCGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTAtt                                                            <  8:146884/148‑1 (MQ=255)
                                                                                    gccgccGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTc                                                                                    <  5:78389/123‑1 (MQ=255)
                                                                                    gccgccGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTc                                                                                    >  6:78389/1‑123 (MQ=255)
                                                                                        ccgATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGt                                                       >  4:225015/1‑148 (MQ=255)
                                                                                        ccgATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGt                                                       <  3:225015/148‑1 (MQ=255)
                                                                                            tcatcaGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACg                                                  >  4:345139/1‑149 (MQ=255)
                                                                                              atcaGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAg                                                         <  2:385358/140‑1 (MQ=255)
                                                                                              atcaGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAg                                                         >  1:385358/1‑140 (MQ=255)
                                                                                              atcaGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACgg                                                 >  8:78626/1‑148 (MQ=255)
                                                                                               tcaGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTa                                               >  5:121255/1‑149 (MQ=255)
                                                                                                caGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTg                                                                                                      <  2:32590/93‑1 (MQ=255)
                                                                                                caGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTg                                                                                                      >  1:32590/1‑93 (MQ=255)
                                                                                                 aGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTAtt                                             <  7:329285/149‑1 (MQ=255)
                                                                                                  gCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTg                                            <  5:305767/149‑1 (MQ=255)
                                                                                                  gCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTg                                            >  1:305486/1‑149 (MQ=255)
                                                                                                  gCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTg                                            <  8:306667/149‑1 (MQ=255)
                                                                                                  gCGGATCCAGCCGCACCACGTCAAGGGCCTTGCCGTTTGCACCGGAAATGCCCGCAACTTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGTCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTg                                            <  5:395351/149‑1 (MQ=255)
                                                                                                    ggCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGcc                                          <  3:196285/149‑1 (MQ=255)
                                                                                                       tCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGtt                                       <  5:138902/149‑1 (MQ=255)
                                                                                                        ccAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGACAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGtt                                       >  6:17944/1‑148 (MQ=255)
                                                                                                         cAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGtt                                       <  3:345139/147‑1 (MQ=255)
                                                                                                          aGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTg                                                                           <  3:49694/110‑1 (MQ=255)
                                                                                                          aGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTg                                                                           >  4:49694/1‑110 (MQ=255)
                                                                                                          aGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTgac                                    <  2:305486/149‑1 (MQ=255)
                                                                                                           gCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGa                                                          >  1:184235/1‑126 (MQ=255)
                                                                                                           gCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGa                                                          <  2:184235/126‑1 (MQ=255)
                                                                                                           gCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTgacg                                   <  5:17944/149‑1 (MQ=255)
                                                                                                                accacGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACg                                                                  >  4:28694/1‑113 (MQ=255)
                                                                                                                accacGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACg                                                                  <  3:28694/113‑1 (MQ=255)
                                                                                                                  cacGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAg                            >  2:287303/1‑149 (MQ=255)
                                                                                                                   acGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAgc                           >  8:46226/1‑149 (MQ=255)
                                                                                                                     gTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAgcgc                         >  8:25170/1‑149 (MQ=255)
                                                                                                                         aGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCgg                     <  4:79608/149‑1 (MQ=255)
                                                                                                                          gcgcCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCggg                    >  7:104719/1‑149 (MQ=255)
                                                                                                                             ccTTTCCGGTGCCAACGGCAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCg                 >  5:241989/1‑149 (MQ=255)
                                                                                                                              cTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGt                <  1:91877/149‑1 (MQ=255)
                                                                                                                                    gTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGTTGGCcg          <  8:104719/149‑1 (MQ=255)
                                                                                                                                     ttGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAg                            <  6:66020/130‑1 (MQ=255)
                                                                                                                                     ttGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAg                            >  5:66020/1‑130 (MQ=255)
                                                                                                                                         cACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGTTGGCCGCGCtt     <  2:193788/149‑1 (MQ=255)
                                                                                                                                          aCCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGtt               <  4:102658/138‑1 (MQ=255)
                                                                                                                                          aCCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGtt               >  3:102658/1‑138 (MQ=255)
                                                                                                                                          aCCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGTTGGCCGCGCttt    >  2:373942/1‑149 (MQ=255)
                                                                                                                                           ccGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGTTGGCCGCGCTTTc   <  7:46226/149‑1 (MQ=255)
                                                                                                                                            cGGCAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGTTGGCCGCGCTTTcg  <  6:241989/149‑1 (MQ=255)
                                                                                                                                            cGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGTTGGCCGCGCTTTcg  >  4:105661/1‑149 (MQ=255)
                                                                                                                                               aaaGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGtt               <  3:351354/133‑1 (MQ=255)
                                                                                                                                               aaaGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGtt               >  4:351354/1‑133 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CGTTGCCACCAAGGTTTCGAGCAGATATTTGGGTACCTGATCGATCTTGATCAGGATGCGGTCTTCGAGATTTTTCGGGTAGATGCCGCCGATCATCAGCGGCTCCAGCCGCACCACGTCAAGCGCCTTGCCGTTGGCACCGGAAAGGCCCGCAACGTAATCGCCAGAGAAGCGCACGCGGACGAACTGCGCCGGTTCCATGCCCTCGTAGAACTGGAAGCCACGGGTATTGAGGTCGACGGTATTGCCGTTGACGGACGCAGCGCCCGGGCCGTTGGCCGCGCTTTCG  >  NC_002947/5322691‑5322979

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: