Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,695,023 G→A 6.4% R330Q (CGG→CAG)  wspC → biofilm formation methyltransferase WspC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,695,0230GA 6.4% 188.7 / 2.1 62R330Q (CGG→CAG) wspCbiofilm formation methyltransferase WspC
Reads supporting (aligned to +/- strand):  ref base G (25/33);  new base A (2/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATC  >  NC_002947/1694882‑1695169
                                                                                                                                             |                                                                                                                                                  
cgacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGcc                                                                                                                                             <  6:35106/149‑1 (MQ=255)
       gCCCGCACCGGTTGATCCGCAGCCCCGTGTGCCGTTGCCTCGTGCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGGTGCTGGCCGGCTTTGCCCGGCTGGCCAATGcc                                                                                                                                       <  5:14288/148‑1 (MQ=255)
        cccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCgg                                                                                                                                     >  2:113618/1‑149 (MQ=255)
        cccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCgg                                                                                                                                     >  2:48058/1‑149 (MQ=255)
        cccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                                      <  5:215800/148‑1 (MQ=255)
         ccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGt                                                                                                                                    <  4:57696/149‑1 (MQ=255)
             aCCGGTTTATCCGCAGCCCAGTGTGCCGTTNCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGc                                                                                                                                  >  1:25846/1‑147 (MQ=255)
             aCCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCa                                                                                                                                <  2:300325/149‑1 (MQ=255)
                  tttATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAg                                                                                                                           <  1:48058/149‑1 (MQ=255)
                      tCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAgg                                                                                                                       >  6:171588/1‑149 (MQ=255)
                        cGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGcc                                                                                                                     >  2:173530/1‑149 (MQ=255)
                                  tgtgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAAt                                                                                                           <  2:145972/149‑1 (MQ=255)
                                     gCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGc                                                                                                                                  <  1:247242/123‑1 (MQ=255)
                                     gCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGc                                                                                                                                  >  2:247242/1‑123 (MQ=255)
                                             cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGTGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                <  5:171588/149‑1 (MQ=255)
                                             cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAAc                                                                                                     >  5:12007/1‑144 (MQ=255)
                                             cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAAc                                                                                                     <  6:12007/144‑1 (MQ=255)
                                              tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCa                                                                                                                                <  4:216173/116‑1 (MQ=255)
                                              tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCa                                                                                                                                >  3:216173/1‑116 (MQ=255)
                                              tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                               >  6:181463/1‑149 (MQ=255)
                                                gTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGGCCAGCCGGGCAATGCCggc                                                                                                                                    <  2:57555/110‑2 (MQ=39)
                                                gTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGGCCAGCCGGGCAATGCCggc                                                                                                                                    >  1:57555/1‑109 (MQ=39)
                                                gTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                              <  1:173530/148‑1 (MQ=255)
                                                 ttCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                              >  7:208949/1‑147 (MQ=255)
                                                 ttCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGa                                                                                            <  8:233958/149‑1 (MQ=255)
                                                    gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                          <  5:195768/148‑1 (MQ=255)
                                                    gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                          <  3:296357/148‑1 (MQ=255)
                                                    gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                          >  6:256693/1‑148 (MQ=255)
                                                     cggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCa                                                                                        <  5:256693/149‑1 (MQ=255)
                                                            tgctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                          >  2:157933/1‑140 (MQ=255)
                                                            tgctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                          <  1:157933/140‑1 (MQ=255)
                                                            tgctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACccg                                                                                 >  2:309476/1‑149 (MQ=255)
                                                            tgctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACccg                                                                                 <  1:113618/149‑1 (MQ=255)
                                                              ctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACccgcc                                                                               >  5:314770/1‑149 (MQ=255)
                                                                 ccGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCAATGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTc                                                                            >  7:289968/1‑149 (MQ=255)
                                                                         ggccaggccaGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAgg                                                                    >  3:86892/1‑149 (MQ=255)
                                                                          gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                                                   <  6:314770/149‑1 (MQ=255)
                                                                           ccaggccagCGCGAGCCCGCTAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCAATGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                                                  <  8:289968/149‑1 (MQ=255)
                                                                            caggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAAt                                                                                                                                          >  6:81208/1‑76 (MQ=255)
                                                                            caggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAAt                                                                                                                                          <  5:81208/76‑1 (MQ=255)
                                                                                 cagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtactac                                                            <  8:208949/149‑1 (MQ=255)
                                                                                        aGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctgg                                                     <  1:215258/149‑1 (MQ=255)
                                                                                        aGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctgg                                                     <  4:86892/149‑1 (MQ=255)
                                                                                          cccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGt                                                   >  2:206813/1‑149 (MQ=255)
                                                                                            cgcgAGCACAGCCACGAGGGGGCCAGCGTGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGttt                                                 <  5:181871/149‑1 (MQ=255)
                                                                                              cgAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGc                                               <  4:244807/149‑1 (MQ=255)
                                                                                                  cacaGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                           <  5:181463/149‑1 (MQ=255)
                                                                                                        cACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAc                                     <  2:25846/149‑1 (MQ=255)
                                                                                                         aCGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAcc                                    >  4:231405/1‑149 (MQ=255)
                                                                                                             gggggCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAg                                >  1:36333/1‑149 (MQ=255)
                                                                                                                         ggtGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcagc                    <  8:28640/147‑1 (MQ=255)
                                                                                                                            gCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAc                                     <  8:120455/129‑1 (MQ=255)
                                                                                                                             cTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgc                >  1:197875/1‑149 (MQ=255)
                                                                                                                              tGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcg               >  1:50658/1‑149 (MQ=255)
                                                                                                                               ggCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc              <  2:197875/149‑1 (MQ=255)
                                                                                                                                gCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCt             <  2:36333/149‑1 (MQ=255)
                                                                                                                                 ccGGCATTGGCCAGCCGGGCAATGCCG‑‑GCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctgg                                                     >  2:136964/1‑106 (MQ=37)
                                                                                                                                 ccGGCATTGGCCAGCCGGGCAATGCCG‑‑GCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctgg                                                     <  1:136964/106‑1 (MQ=37)
                                                                                                                                  cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGa           <  3:231405/149‑1 (MQ=255)
                                                                                                                                    gCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGc         >  2:305508/1‑149 (MQ=255)
                                                                                                                                         gCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTa    >  1:301669/1‑149 (MQ=255)
                                                                                                                                           ccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATc  >  3:222294/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATC  >  NC_002947/1694882‑1695169

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