Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,322,365 T→A 7.0% E296V (GAG→GTG)  mrcB ← transglycosylase/transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,322,3650TA 7.0% 165.6 / 3.3 57E296V (GAG→GTG) mrcBtransglycosylase/transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (26/27);  new base A (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CTGCCTGGCTGGCCACGCCCTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTCGCGTTTGTCGTAGTGCAGTTCGAGCAGGACCGCCATCA  >  NC_002947/5322219‑5322512
                                                                                                                                                  |                                                                                                                                                   
cTGCCTGGCTGGCCACGCCCTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCt                                                                                                                                                     <  7:270656/147‑1 (MQ=255)
 tGCCTGGCTGGCCACGCCCTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGCCGCCGCGGGGTGTAGTAGGACGGGCCTATGCCCATGCCTACCAGCAGGGCGAGAAGCTGCAGTTTCAGcgcg                                                                                                                                                   >  3:215910/1‑148 (MQ=255)
   ccTGGCTGGCCACGCCCTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGa                                                                                                                                                 >  1:57994/1‑148 (MQ=255)
            ccACGCCCTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCggctgg                                                                                                                                       >  8:337804/1‑149 (MQ=255)
            ccACGCCCTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCggctgg                                                                                                                                       <  2:158045/149‑1 (MQ=255)
             cACGCCCTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCggctggc                                                                                                                                      <  6:89160/149‑1 (MQ=255)
             cACGCCCTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCggctggc                                                                                                                                      >  6:61808/1‑149 (MQ=255)
                gCCCTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCggctggct                                                                                                                                     >  8:260232/1‑147 (MQ=255)
                gCCCTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCggctggct                                                                                                                                     <  7:260232/147‑1 (MQ=255)
                 cccTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGgctggctg                                                                                                                                    >  1:134021/1‑147 (MQ=255)
                   cTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTgaag                                                                                                                                 >  3:118673/1‑148 (MQ=255)
                     gCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTgaag                                                                                                                                 >  5:360463/1‑146 (MQ=255)
                      cTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTgaagaa                                                                                                                               <  7:337656/147‑1 (MQ=255)
                         ggCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACt                                                                                                                          <  6:233964/149‑1 (MQ=255)
                            caccaGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGc                                                                                                                       >  8:182050/1‑149 (MQ=255)
                            caccaGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGc                                                                                                                       >  3:299519/1‑149 (MQ=255)
                              ccaGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTa                                                                                                                     <  4:299519/149‑1 (MQ=255)
                                aGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGc                                                                                                                   <  5:61808/149‑1 (MQ=255)
                                  gTCGAGCACCAGGTTGCGGCGAGCCAGGGCCATCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGAGGATGGCTGAAGAAGAACTGGCTAGc                                                                                                                   >  5:329738/1‑147 (MQ=255)
                                     gAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGAc                                                                                                              <  7:337804/149‑1 (MQ=255)
                                          ccAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAg                                                                                                         >  6:38739/1‑149 (MQ=255)
                                           cAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGc                                                                                                        >  6:222042/1‑149 (MQ=255)
                                              gTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGt                                                                                                     <  5:222042/149‑1 (MQ=255)
                                               ttGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTg                                                                                                    >  1:114222/1‑149 (MQ=255)
                                               ttGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTg                                                                                                    <  5:38739/149‑1 (MQ=255)
                                                            gggCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGa                                                                                       <  5:54027/149‑1 (MQ=255)
                                                            gggCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGa                                                                                       <  4:118673/149‑1 (MQ=255)
                                                                ccGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGTTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGt                                                                                   >  3:52697/1‑149 (MQ=255)
                                                                  gCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTcc                                                                                 >  1:213579/1‑149 (MQ=255)
                                                                      ggggTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTg                                                                               <  7:336132/147‑1 (MQ=255)
                                                                            aCGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCACGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCcgacg                                                                        >  6:170040/1‑148 (MQ=255)
                                                                              ggcgagggTTGTAGTAGGACGGGCCGTTGACCATGCCTACCAGCAGGGCTATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGTCTAGCCAGACCAAATCCTTTCACCGCGCGCTGACCGTCCTGGCCTACGaag                                                                     <  6:329738/144‑2 (MQ=255)
                                                                              gccgcGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCACGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAac                                                                     <  5:170040/149‑1 (MQ=255)
                                                                                      ttGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGtt                                                             <  2:114222/149‑1 (MQ=255)
                                                                                      ttGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGtt                                                             <  4:301667/149‑1 (MQ=255)
                                                                                           gtagGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTg                                                                               <  3:335317/126‑1 (MQ=255)
                                                                                           gtagGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTg                                                                               >  4:335317/1‑126 (MQ=255)
                                                                                           gtagGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGt                                                        <  2:57994/149‑1 (MQ=255)
                                                                                            tagGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTa                                                       <  7:182050/149‑1 (MQ=255)
                                                                                               gACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGc                                                    >  3:43978/1‑149 (MQ=255)
                                                                                                    gCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTc                                                 >  6:156931/1‑147 (MQ=255)
                                                                                                    gCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTCGCCGACGAACACCTCGTTGAGGTAGGCCTc                                                 >  3:275930/1‑147 (MQ=255)
                                                                                                          gACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGATGAACACCTCGTTGAGGTAGGCCTCAAGGAtt                                          <  2:134021/148‑1 (MQ=255)
                                                                                                           aCCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTc                                        >  8:288080/1‑149 (MQ=255)
                                                                                                           aCCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTc                                        >  2:280640/1‑149 (MQ=255)
                                                                                                           aCCATGCCAACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTc                                        <  7:288080/149‑1 (MQ=255)
                                                                                                             cATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTcgc                                      >  3:74588/1‑149 (MQ=255)
                                                                                                                gCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGt                                                                                                     >  2:213010/1‑83 (MQ=255)
                                                                                                                gCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTAGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGt                                                                                                     <  1:213010/83‑1 (MQ=255)
                                                                                                                 ccTACCAGCCTGCCTATGTCATGCTGTTTGAGCTCGGACAGCGGCTGGCGGAAGACGACCTGGCTAGCCAGACCAAAGGCGTGCCCCGCGCGCTGACGTTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTCGCGttt                                  <  4:215910/149‑1 (MQ=255)
                                                                                                                     ccAGGAGTGCGATATGATGCAGTTTCTGCTCGGAGAGCTGCTGGCTGAATAAGAACTGGCTAGCCAGACCTAACGCTTGCACCGCGCGCTGACCGTCATGGCCGACGAAGACTTCGTTTAGGTAGGCCTCAAGGATTTCGCGTTTGTCg                              <  6:360463/149‑1 (MQ=255)
                                                                                                                     ccAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTCGCGTTTGTc                               <  4:43978/148‑1 (MQ=255)
                                                                                                                      cagcagGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTCGCGTTTGTCg                              >  3:284686/1‑148 (MQ=255)
                                                                                                                      cagcagGGCGATCTGATGCAGTTTCAGCACGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCCCCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTCGCGTTTGTCg                              >  7:188956/1‑148 (MQ=255)
                                                                                                                        gcagGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTCGCGTTTGTCGTAg                           >  1:210653/1‑149 (MQ=255)
                                                                                                                          agGGCGATCTGATGCAGTTTCAGCTCTGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTCGCGTTTGTCGTAGt                          <  4:284686/148‑1 (MQ=255)
                                                                                                                                     aTGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTCGCGTTTGTCGTAGTGCAGTTCGAGCa              <  2:213579/149‑1 (MQ=255)
                                                                                                                                             tCAGCACGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTCGCGTTTGTCGTAGTGCAGTTCGAGCAGGACCGc       <  8:188956/148‑1 (MQ=255)
                                                                                                                                                 cTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTCGCGTTTGTCGTCGTGCAGTGAGAGCAGGACGGCcatca  >  8:165832/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGCCTGGCTGGCCACGCCCTGCTCGGCCACCAGGTCGAGCACCAGGTTGCGGCGAGCCAGGGCCCGCTCGGGGTGACGCCGCGGGTTGTAGTAGGACGGGCCTTTGACCATGCCTACCAGCAGGGCGATCTGATGCAGTTTCAGCTCGGACAGCGGCTGGCTGAAGAAGAACTGGCTAGCCAGACCAAAGCCGTGCACCGCGCGCTGACCGTCCTGGCCGACGAACACCTCGTTGAGGTAGGCCTCAAGGATTTCGCGTTTGTCGTAGTGCAGTTCGAGCAGGACCGCCATCA  >  NC_002947/5322219‑5322512

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