Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,298,604 G→A 5.2% D506D (GAC→GAT sbcC ← exonuclease subunit SbcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,298,6040GA 5.2% 300.2 / 3.2 96D506D (GAC→GATsbcCexonuclease subunit SbcC
Reads supporting (aligned to +/- strand):  ref base G (46/45);  new base A (3/2);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

GCGTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGTCAGCC  >  NC_002947/2298463‑2298749
                                                                                                                                             |                                                                                                                                                  
gCGTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGAt                                                                                                                                               <  3:247621/147‑1 (MQ=255)
gCGTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCt                                                                                                                                             >  2:251215/1‑149 (MQ=255)
gCGTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCt                                                                                                                                             <  6:211797/149‑1 (MQ=255)
  gTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGc                                                                                                                                           >  5:298745/1‑149 (MQ=255)
      cagcaTGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGATCGACCTGCTCgg                                                                                                                                       <  3:2390/149‑1 (MQ=255)
       agcaTGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGc                                                                                                                                      <  4:36020/149‑1 (MQ=255)
          aTGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAg                                                                                                                                   >  4:169377/1‑149 (MQ=255)
             gCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGc                                                                                                                                >  5:327067/1‑149 (MQ=255)
                  ctgctgcCGTTCTCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGGTGTCCTTCATCATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTg                                                                                                                             <  6:262710/147‑1 (MQ=255)
                  ctgctgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTggg                                                                                                                           <  6:219462/149‑1 (MQ=255)
                   tgctgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTc                                                                                                                                         <  4:212652/134‑1 (MQ=255)
                   tgctgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTc                                                                                                                                         >  3:212652/1‑134 (MQ=255)
                     ctgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGtt                                                                                                                        <  3:169377/149‑1 (MQ=255)
                      tgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTc                                                                                                                       >  3:196547/1‑149 (MQ=255)
                      tgcCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTACAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTc                                                                                                                       >  6:64383/1‑149 (MQ=255)
                        ccGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCgg                                                                                                                     >  4:300882/1‑149 (MQ=255)
                        ccGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCg                                                                                                                      <  6:307274/148‑1 (MQ=255)
                          gTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGcc                                                                                                                   >  3:77400/1‑149 (MQ=255)
                          gTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTAACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATTCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGcc                                                                                                                   >  3:283661/1‑149 (MQ=255)
                              gcgcAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTgg                                                                                                               <  6:52334/149‑1 (MQ=255)
                               cgcAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTACAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGc                                                                                                              <  5:64383/149‑1 (MQ=255)
                                 cAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCtt                                                                                                            <  6:298745/149‑1 (MQ=255)
                                  aGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCtt                                                                                                            >  3:281964/1‑148 (MQ=255)
                                  aGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCtt                                                                                                            >  3:286595/1‑148 (MQ=255)
                                        cccAGCTGCCGGCGCAGCGCCGGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGCTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGaa                                                                                                     <  6:95821/148‑1 (MQ=255)
                                              tgctgGTGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAg                                                                                               >  8:107516/1‑149 (MQ=255)
                                                ctgGCGCAGGTCCTGTTCACGGCCTTGCAGGCGTGACATTTCCTCTACGGCGCGATGCTGTTTGCGGTTGTCATGCAACATGCCGCCGTGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAgg                                                                                              <  4:283661/148‑1 (MQ=255)
                                                ctgGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGt                                                                                            >  4:14871/1‑149 (MQ=255)
                                                    cgcAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTa                                                                                        <  4:71093/149‑1 (MQ=255)
                                                       aGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGc                                                                                     >  5:104226/1‑149 (MQ=255)
                                                          tCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGaacaac                                                                                                 <  7:26180/135‑1 (MQ=255)
                                                          tCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGaacaac                                                                                                 >  8:26180/1‑135 (MQ=255)
                                                          tCCTGTACACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGt                                                                                         <  3:95460/142‑1 (MQ=255)
                                                          tCCTGTACACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGt                                                                                         >  4:95460/1‑142 (MQ=255)
                                                            cTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTcgc                                                                                >  4:48248/1‑149 (MQ=255)
                                                             tGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGaa                                                                                                     <  6:211550/128‑1 (MQ=255)
                                                             tGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGaa                                                                                                     >  5:211550/1‑128 (MQ=255)
                                                                  aCGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTg                                                                                                                             <  7:26424/99‑1 (MQ=255)
                                                                  aCGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTg                                                                                                                             >  8:26424/1‑99 (MQ=255)
                                                                   cGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGc                                                                         <  6:104226/149‑1 (MQ=255)
                                                                      ccGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAgg                                                                                                                                  <  1:323683/90‑1 (MQ=255)
                                                                      ccGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAgg                                                                                                                                  >  2:323683/1‑90 (MQ=255)
                                                                       cGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGc                                                                     <  4:77400/149‑1 (MQ=255)
                                                                       cGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGCCGTCGCGCTCGGCCTg                                                                      >  2:232795/1‑148 (MQ=255)
                                                                                gTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAgc                                                                 <  7:91686/144‑1 (MQ=255)
                                                                                gTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAgc                                                                 >  8:91686/1‑144 (MQ=255)
                                                                                   aCATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTg                                                         >  5:219437/1‑149 (MQ=255)
                                                                                   aCATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTg                                                         >  3:179021/1‑149 (MQ=255)
                                                                                    cATTTCCTCTACGGCGCGCTGCTGTTTGCCGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGTCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTg                                                         <  4:196547/148‑1 (MQ=255)
                                                                                         ccTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCctg                                                    <  3:14871/148‑1 (MQ=255)
                                                                                           tctACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCctggct                                                 >  7:42144/1‑149 (MQ=255)
                                                                                              aCGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCtggctggc                                              >  7:200796/1‑149 (MQ=255)
                                                                                               cGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTggcc                                             >  6:108957/1‑149 (MQ=255)
                                                                                                  cgcgcTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTggcctgg                                          <  3:48248/149‑1 (MQ=255)
                                                                                                     gctgctGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGCTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTggcctgggc                                        <  1:251215/149‑1 (MQ=255)
                                                                                                     gctgctGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTggcctgggcc                                       >  7:120909/1‑149 (MQ=255)
                                                                                                      ctgctgTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTggcctgggcct                                      >  6:85078/1‑149 (MQ=255)
                                                                                                       tgctgTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTggcctgggcctg                                     >  4:116987/1‑149 (MQ=255)
                                                                                                        gctgTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTggcctgggcctg                                     >  6:287819/1‑148 (MQ=255)
                                                                                                                cGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGTTCTGCGCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCggg                            <  3:301530/149‑1 (MQ=255)
                                                                                                                cGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCgg                             <  7:159550/148‑1 (MQ=255)
                                                                                                                 ggTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCTTGGGCCTGCAGGCCggg                            >  6:267447/1‑148 (MQ=255)
                                                                                                                 ggTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGt                           <  6:219437/149‑1 (MQ=255)
                                                                                                                 ggTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGt                           <  4:320265/149‑1 (MQ=255)
                                                                                                                  gTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTa                          <  5:74493/149‑1 (MQ=255)
                                                                                                                   ttGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAg                         <  4:286595/149‑1 (MQ=255)
                                                                                                                   ttGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAg                         <  4:281964/149‑1 (MQ=255)
                                                                                                                   ttGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAg                         <  4:151138/149‑1 (MQ=255)
                                                                                                                     gTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGct                       >  3:52659/1‑149 (MQ=255)
                                                                                                                      tCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGctc                      >  5:16006/1‑149 (MQ=255)
                                                                                                                        cTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGctct                     <  8:42144/148‑1 (MQ=255)
                                                                                                                        cTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGctct                     >  4:69184/1‑148 (MQ=255)
                                                                                                                        cTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCtctt                    >  6:184726/1‑149 (MQ=255)
                                                                                                                         tGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGccaatgc                          <  3:111098/142‑6 (MQ=255)
                                                                                                                         tGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGccaatgc                          >  4:111098/1‑137 (MQ=255)
                                                                                                                         tGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCtcttc                   >  1:47373/1‑149 (MQ=255)
                                                                                                                         tGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCtcttc                   >  2:278764/1‑149 (MQ=255)
                                                                                                                         tGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCtcttc                   >  8:286727/1‑149 (MQ=255)
                                                                                                                           cAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCtcttct                  <  5:85078/148‑1 (MQ=255)
                                                                                                                            aaCATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCg                >  5:247845/1‑149 (MQ=255)
                                                                                                                             aCATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCgg               <  5:108957/149‑1 (MQ=255)
                                                                                                                              cATGCCGCCGCGCAGGTCGATCTGCTCGGCCAGGGCCGGGGGTTCGGCCTGGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCCTTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTGCTCGGc              >  7:122165/1‑149 (MQ=255)
                                                                                                                               aTGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAg                         >  2:47413/1‑137 (MQ=255)
                                                                                                                               aTGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAg                         <  1:47413/137‑1 (MQ=255)
                                                                                                                                 ctggccgAGCAGATCGACCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCt           >  1:120138/4‑149 (MQ=255)
                                                                                                                                  ccgccgAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCtt          <  3:300882/149‑1 (MQ=255)
                                                                                                                                  ccgccgAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCtt          >  1:305547/1‑149 (MQ=255)
                                                                                                                                   cgccgAGCAGGTCGATCTGCTCGGCGAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGGCGAGCTCGTACTTGAATCGCGCATTGGCCTGTCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCACTTg         <  5:228254/149‑1 (MQ=255)
                                                                                                                                   cgccgAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTg         <  2:22306/149‑1 (MQ=255)
                                                                                                                                     ccgAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGt       <  4:133870/149‑1 (MQ=255)
                                                                                                                                     ccgAGCAGATCGACCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGt       >  8:30475/1‑149 (MQ=255)
                                                                                                                                      cgAGCAGATCGACCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTtctc                 >  2:217345/1‑138 (MQ=255)
                                                                                                                                      cgAGCAGATCGACCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTtctc                 <  1:217345/138‑1 (MQ=255)
                                                                                                                                       gAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGTCa     <  3:69184/149‑1 (MQ=255)
                                                                                                                                        agcagcTCGATCTGCTCGGCCAGGGCCAGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGTCAg    >  5:66661/1‑149 (MQ=255)
                                                                                                                                         gCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGTCAGc   >  1:84801/1‑149 (MQ=255)
                                                                                                                                          cAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGTCAGcc  <  2:47373/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCGTTGCAGCATGGCCTGCTGCTGCCGTTCGCGCAGGGCATCCAGCTGCTGGCGCAGCTCCTGTTCACGGCCGTGCAGGCGTGACATTTCCTCTACGGCGCGCTGCTGTTTGCGGTTGTCCTGCAACATGCCGCCGAGCAGGTCGATCTGCTCGGCCAGGGCCTGGGGTTCGGCCTTGGCTTCACGGAACAACAGGTCGTAGGCGTCGCGCTCGGCCTGCAAGCGCGCATTGGCCTGGCTGGCCTGGGCCTGCAGGCCGGGTAGCTCTTCTCGGCCCTTGGTCAGCC  >  NC_002947/2298463‑2298749

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: