Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,141,385 T→A 8.4% L374* (TTG→TAG) ‡ 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,3850TA 8.4% 207.1 / 7.8 71L374* (TTG→TAG) ‡PP_1897DNA internalization‑related competence protein ComEC/Rec2
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base A (2/4);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

CTGGGTTCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTG  >  NC_002947/2141238‑2141525
                                                                                                                                                    |                                                                                                                                            
cTGGGTTCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTg                                                                                                                                             <  1:69804/149‑1 (MQ=255)
cTGGGTTCTGGCTGTCGTTAGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTg                                                                                                                                             >  3:151551/1‑149 (MQ=255)
     ttCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGt                                                                                                                                        <  8:15041/149‑1 (MQ=255)
          gCTGTCGTTCGGAAGCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCAAGCCGGTggcca                                                                                                                                   <  1:312164/149‑5 (MQ=255)
          gCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCTCACCGGCAAGCCGGTggc                                                                                                                                     <  8:148010/148‑3 (MQ=255)
          gCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGa                                                                                                                                      <  6:78911/146‑1 (MQ=255)
          gCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGa                                                                                                                                      <  2:362350/146‑1 (MQ=255)
          gCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGcc                                                                                                                                   >  2:181804/1‑149 (MQ=255)
          gCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGc                                                                                                                                    >  4:333278/1‑148 (MQ=255)
          gCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCAAGCCGGTggc                                                                                                                                     <  5:336777/147‑3 (MQ=255)
             gTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCt                                                                                                                                  >  3:139135/1‑147 (MQ=255)
             gTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCt                                                                                                                                  >  7:290727/1‑147 (MQ=255)
              tCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCa                                                                                                                                <  3:223879/148‑1 (MQ=255)
                    gccgccGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTgcgc                                                                                                                          <  1:59418/148‑1 (MQ=255)
                         cgTTGCCACGCTCATCTACTGCTTCAGTTCCCGGTTGGGTGGTTGGCTGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTgg                                                                                                                    <  6:333915/149‑1 (MQ=255)
                               cACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAAcc                                                                                                              >  5:250255/1‑149 (MQ=255)
                                      aTCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTtgtg                                                                                                         <  6:188352/147‑1 (MQ=255)
                                      aTCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTgg                                                                                                        <  2:119070/148‑1 (MQ=255)
                                       tCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCg                                                                                                      >  3:241346/1‑149 (MQ=255)
                                           cTGCTTCACTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGc                                                                                                  >  2:113538/1‑149 (MQ=255)
                                            tGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGcc                                                                                                 >  7:228263/1‑149 (MQ=255)
                                            tGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGATTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGc                                                                                                  >  3:90749/1‑148 (MQ=255)
                                              cttgattGCCCTGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGc                                                                                                  <  4:90749/140‑1 (MQ=255)
                                              cTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCTCACCGGCAAGCCGGTggc                                                                                                                                     <  6:44065/112‑3 (MQ=255)
                                              cTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCTCACCGGCAAGCCGGTggc                                                                                                                                     >  5:44065/1‑110 (MQ=255)
                                              cTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGt                                                                                               >  3:125954/1‑149 (MQ=255)
                                              cTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGt                                                                                               <  1:47239/149‑1 (MQ=255)
                                               ttCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTg                                                                                              >  7:232548/1‑149 (MQ=255)
                                               ttCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTg                                                                                              >  3:158225/1‑149 (MQ=255)
                                                 cAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTgg                                                                                             <  2:119451/148‑1 (MQ=255)
                                                  aGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCAGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGt                                                                                           >  2:61291/1‑149 (MQ=255)
                                                    tGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCa                                                                                         <  3:333278/149‑1 (MQ=255)
                                                      cccGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGc                                                                                       >  7:246733/1‑149 (MQ=255)
                                                           ttGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGc                                                                                  <  2:337196/149‑1 (MQ=255)
                                                               gtggtTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTc                                                                              >  1:249968/1‑149 (MQ=255)
                                                                  gtTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCtt                                                                            <  4:158225/148‑1 (MQ=255)
                                                                  gtTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCACTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCtt                                                                            <  4:139135/148‑1 (MQ=255)
                                                                          cccTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGc                                                                   >  3:285577/1‑149 (MQ=255)
                                                                           ccTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCg                                                                  >  3:114940/1‑149 (MQ=255)
                                                                              ggCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTg                                                               >  3:300715/1‑149 (MQ=255)
                                                                               gCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGc                                                              <  1:113538/149‑1 (MQ=255)
                                                                                acatccTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCt                                                             <  6:154214/145‑1 (MQ=255)
                                                                                  gtccTGGAGAGTAGAGCAATGGGTGATGGCGGTCGGTTTTCTTCTGGTGGTGCTGTCCACCGTCTTGCCGGTGAGCCTCAGTGCGCAGCTGGCCAAGCTTGTGGCGGTGGCGTTGGTAAGCCTGGCGGTCTTGCCGTTGGAGTTGCTgg                                                           <  4:300715/147‑1 (MQ=255)
                                                                                  ggCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTgg                                                           <  4:114940/149‑1 (MQ=255)
                                                                                      tGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCAc                                                       >  2:151856/1‑149 (MQ=255)
                                                                                          cacGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTg                                                   <  8:290727/149‑1 (MQ=255)
                                                                                             gAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATg                                                <  2:304548/149‑1 (MQ=255)
                                                                                                cgcAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGtt                                              <  4:151551/148‑1 (MQ=255)
                                                                                                       ggTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTgg                                      <  4:241346/149‑1 (MQ=255)
                                                                                                           aTCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCgg                                   <  2:241734/148‑1 (MQ=255)
                                                                                                            tCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGgt                                 <  1:151856/149‑1 (MQ=255)
                                                                                                            tCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGgt                                 >  6:94435/1‑149 (MQ=255)
                                                                                                             cGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGgtg                                >  5:98863/1‑149 (MQ=255)
                                                                                                             cGCCGTCGGCTTGCGACCGGTGTTGCTGGCCACCGGCTTGCCGGTGCGACTCAGTGCTCCGCTGGCCACCCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGAACGCTGATGTTGCCGCTGGCGGgtg                                >  1:206385/1‑149 (MQ=255)
                                                                                                               ccGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGtt                              <  8:246733/149‑1 (MQ=255)
                                                                                                                     gCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAg                        <  6:98863/149‑1 (MQ=255)
                                                                                                                      cTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAgg                       <  4:285577/149‑1 (MQ=255)
                                                                                                                       ttGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGc                      >  3:255317/1‑149 (MQ=255)
                                                                                                                         gCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTg                                                                           >  2:370612/1‑94 (MQ=255)
                                                                                                                         gCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTg                                                                           <  1:370612/94‑1 (MQ=255)
                                                                                                                          cTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGt                    <  1:181804/148‑1 (MQ=255)
                                                                                                                            aCCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTAGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGc                 >  1:303301/1‑149 (MQ=255)
                                                                                                                                gtgtTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTctg             <  1:61291/149‑1 (MQ=255)
                                                                                                                                      cTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAg                        >  3:316631/1‑132 (MQ=255)
                                                                                                                                      cTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAg                        <  4:316631/132‑1 (MQ=255)
                                                                                                                                          ccACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGt   >  4:320281/1‑149 (MQ=255)
                                                                                                                                           tcACCGGCAAGCCGGTGGGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGt   >  6:366318/2‑149 (MQ=255)
                                                                                                                                            cACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTg  <  6:250255/149‑1 (MQ=255)
                                                                                                                                            cACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTg  >  6:130599/1‑149 (MQ=255)
                                                                                                                                               cGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTg                                             <  7:191615/103‑1 (MQ=255)
                                                                                                                                               cGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTg                                             >  8:191615/1‑103 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CTGGGTTCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTG  >  NC_002947/2141238‑2141525

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: