Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,141,384 T→A 6.1% L374M (TTG→ATG)  PP_1897 → DNA internalization‑related competence protein ComEC/Rec2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,141,3840TA 6.1% 197.1 / 3.1 65L374M (TTG→ATG) PP_1897DNA internalization‑related competence protein ComEC/Rec2
Reads supporting (aligned to +/- strand):  ref base T (34/27);  new base A (2/2);  total (36/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 = 6.73e-01

GGGTTCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCT  >  NC_002947/2141240‑2141532
                                                                                                                                                 |                                                                                                                                                     
gggTTCTGGCTGTCGTTCGCCGCCGTGGCCACCCTCAACTACTGATTCAGTGCCCGGTTGGGTGGTTCGCGGCCCTGGCAGGCCTGTACACCAGCGCAATGGGTGATCACCTTCGGCTTGCTACCGGTGTTGCTGGCAAACGGCTTGcc                                                                                                                                                   >  8:250085/1‑149 (MQ=255)
  gTTCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCgg                                                                                                                                                 <  1:41122/149‑1 (MQ=255)
   ttCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGt                                                                                                                                                <  5:363279/149‑1 (MQ=255)
        gCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGa                                                                                                                                              <  7:185795/146‑1 (MQ=255)
        gCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGcc                                                                                                                                           >  8:314789/1‑149 (MQ=255)
        gCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGc                                                                                                                                            <  7:170962/148‑1 (MQ=255)
               ttCGCCGCCCTGGCCACGCTCAACTACAGCTTCAGTGCCCGGTTGGCTGGTTGGCGGCCCTGGCAGGACTGGACACGAGCGCAATGGGTGATCGAAGTCGCCTTGATACCCGAGATGATCGCCACCGGCTTGACGGTGAGCCTCAGTgc                                                                                                                                    >  4:48832/1‑149 (MQ=255)
                  gccgccGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGa                                                                                                                                              <  1:262734/136‑1 (MQ=255)
                  gccgccGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGa                                                                                                                                              >  2:262734/1‑136 (MQ=255)
                   ccgccgTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGcgccg                                                                                                                                >  7:123670/1‑149 (MQ=255)
                            ccACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGc                                                                                                                           >  1:191933/1‑145 (MQ=255)
                            ccACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGc                                                                                                                           <  2:191933/145‑1 (MQ=255)
                             cACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCAAGCCGGTggcca                                                                                                                                           <  7:308498/128‑5 (MQ=255)
                             cACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCAAGCCGGTggcca                                                                                                                                           >  8:308498/1‑124 (MQ=255)
                             cACGCTCATCTACTGCTTCAGTGCCCGGTGGGTTGGTTGGCGGCCCTGTCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTgcgc                                                                                                                                  <  6:169360/137‑1 (MQ=255)
                             cACGCTCATCTACTGCTTCAGTGACCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCTGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTgcgc                                                                                                                                  >  5:169360/1‑137 (MQ=255)
                                    aTCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTg                                                                                                                             >  3:70138/1‑135 (MQ=255)
                                    aTCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTg                                                                                                                             <  4:70138/135‑1 (MQ=255)
                                     tCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCg                                                                                                              >  3:1544/1‑149 (MQ=255)
                                        aCTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTg                                                                                                           >  1:351974/1‑149 (MQ=255)
                                         cTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGt                                                                                                            <  2:50989/147‑1 (MQ=255)
                                              tCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCTCACCGGCAAGCCGGTggcc                                                                                                                                            <  7:311118/111‑4 (MQ=255)
                                              tCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCTCACCGGCAAGCCGGTggcc                                                                                                                                            >  8:311118/1‑108 (MQ=255)
                                                  tGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGGCAACCTTGTGGGGGTGCCGTGGGt                                                                                                  >  8:183533/1‑147 (MQ=255)
                                                  tGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGt                                                                                                  >  3:194348/1‑147 (MQ=255)
                                                      cGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCt                                                                                            >  5:187222/1‑149 (MQ=255)
                                                           gggTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCg                                                                                        <  2:351974/148‑1 (MQ=255)
                                                                ggtGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCtt                                                                                   <  7:250085/146‑1 (MQ=255)
                                                                     cggcCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGtt                                                                             >  5:288592/1‑149 (MQ=255)
                                                                       gcCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTgg                                                                           >  7:319309/1‑149 (MQ=255)
                                                                             gCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGc                                                                     <  4:194348/149‑1 (MQ=255)
                                                                               aGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTg                                                                   <  7:314789/149‑1 (MQ=255)
                                                                                   cTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGc                                                                                         <  5:258515/123‑1 (MQ=255)
                                                                                   cTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGc                                                                                         >  6:258515/1‑123 (MQ=255)
                                                                                   cTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTgg                                                                           >  2:18442/1‑137 (MQ=255)
                                                                                   cTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTgg                                                                           <  1:18442/137‑1 (MQ=255)
                                                                                      gACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGc                                                            >  2:291212/1‑149 (MQ=255)
                                                                                          cGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGAt                                                        <  2:132278/149‑1 (MQ=255)
                                                                                           gAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTtgtg                                                                                                                 <  4:229301/92‑1 (MQ=255)
                                                                                           gAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTtgtg                                                                                                                 >  3:229301/1‑92 (MQ=255)
                                                                                           gAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATg                                                       >  1:86950/1‑149 (MQ=255)
                                                                                            aGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGt                                                      <  6:288592/149‑1 (MQ=255)
                                                                                              cgcAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTg                                                    >  8:331540/1‑149 (MQ=255)
                                                                                               gcAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGc                                                   <  7:331540/149‑1 (MQ=255)
                                                                                                   tGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCt                                               >  7:274210/1‑149 (MQ=255)
                                                                                                   tGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCt                                               >  4:80527/1‑149 (MQ=255)
                                                                                                           cGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGc                                                                     <  4:156703/119‑1 (MQ=255)
                                                                                                           cGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGc                                                                     >  3:156703/1‑119 (MQ=255)
                                                                                                           cGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGgt                                        <  8:274210/148‑1 (MQ=255)
                                                                                                                  ggCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGa                                <  8:123670/149‑1 (MQ=255)
                                                                                                                   gCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAg                               <  6:187222/149‑1 (MQ=255)
                                                                                                                    cTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAgg                              >  4:149987/1‑149 (MQ=255)
                                                                                                                     ttGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGc                             >  2:35999/1‑149 (MQ=255)
                                                                                                                         tACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTg                         >  3:202239/1‑149 (MQ=255)
                                                                                                                               tgtTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTctgg                   <  8:319309/149‑1 (MQ=255)
                                                                                                                                   gCTGGCCACCGACTTGCCGGTGAGCCTCAGTGCGCCGCTGGTCAACCTTGTGGCGGTGCCGTGGGTCAGCCTg                                                                                           <  6:220710/73‑1 (MQ=255)
                                                                                                                                   gCTGGCCACCGACTTGCCGGTGAGCCTCAGTGCGCCGCTGGTCAACCTTGTGGCGGTGCCGTGGGTCAGCCTg                                                                                           >  5:220710/1‑73 (MQ=255)
                                                                                                                                       gCCACCGGCTTGCCGGTGAGCCTCAGTGGGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAgg           >  5:316007/1‑149 (MQ=255)
                                                                                                                                        ccACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGt          <  4:202239/149‑1 (MQ=255)
                                                                                                                                          cACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTg         <  3:330838/149‑1 (MQ=255)
                                                                                                                                          cACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTg         >  6:40520/1‑149 (MQ=255)
                                                                                                                                           aCCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGGTCAGCCTGGCGGGATTGCCGTTGGCGTTGATGGGCCCGATGTTGTTGCCGATGGCGGGTGTTGGCGAGGCGTTGGTGTGGATGGCAGGTg         >  1:278002/1‑149 (MQ=255)
                                                                                                                                           aCCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTgg        >  7:172880/1‑149 (MQ=255)
                                                                                                                                             cGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGcc      >  3:97525/1‑149 (MQ=255)
                                                                                                                                             cGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGcc      >  6:9850/1‑149 (MQ=255)
                                                                                                                                                 ttGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCTTTGCTGGGCAGGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCctgct  <  3:48832/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GGGTTCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCT  >  NC_002947/2141240‑2141532

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: