Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,463,593 G→C 8.5% R270G (CGT→GGT)  PP_3078 ← ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,463,5930GC 8.5% 282.5 / 9.9 93R270G (CGT→GGT) PP_3078ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (43/42);  new base C (4/4);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCGGCA  >  NC_002947/3463450‑3463742
                                                                                                                                                |                                                                                                                                                     
cAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGTGGGCACGTGACCTGcccg                                                                                                                                                 >  2:1669514/1‑149 (MQ=255)
cAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCcca                                                                                                                                                 >  2:1485267/1‑150 (MQ=255)
cAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCATAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCcca                                                                                                                                                 >  1:362322/1‑150 (MQ=255)
 aaGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCccac                                                                                                                                                >  2:972639/1‑150 (MQ=255)
  aGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCccacc                                                                                                                                               <  1:906226/150‑1 (MQ=255)
   ggTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGTGGGCACGTGACCTGCccgcca                                                                                                                                              <  1:1669514/150‑1 (MQ=255)
    gtgtAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCATGGTGGGCACGTGACCTGCccgcca                                                                                                                                              <  1:548289/150‑1 (MQ=255)
         ggCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCgg                                                                                                                                        >  2:918363/1‑150 (MQ=255)
          gCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTCCCTGGCGCGCCACTTCGGTGTCGCCGGAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCggg                                                                                                                                       <  2:1895847/150‑1 (MQ=255)
                 cGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGtt                                                                                                                                <  2:880613/150‑1 (MQ=255)
                  ggAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGGCGCGTTc                                                                                                                               >  2:2056019/1‑150 (MQ=255)
                   gagaACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCg                                                                                                                              >  2:423759/1‑150 (MQ=255)
                     gaACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGaa                                                                                                                            <  1:288036/150‑1 (MQ=255)
                      aaCTCGCGGTTGTAGTCGTCGCCACCGCATTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAg                                                                                                                           <  1:125883/150‑1 (MQ=255)
                       aCTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAgg                                                                                                                          <  2:152443/150‑1 (MQ=255)
                         tCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGTGGGCACGTGACCTGCccgcca                                                                                                                                              >  1:411899/1‑128 (MQ=255)
                         tCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGTGGGCACGTGACCTGCccgcca                                                                                                                                              <  2:411899/128‑1 (MQ=255)
                          cgcgGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGa                                                                                                                       <  1:1485267/150‑1 (MQ=255)
                              gTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGAtgct                                                                                                                   >  2:2289624/1‑150 (MQ=255)
                              gTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGATATGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGAtgct                                                                                                                   >  1:633433/1‑150 (MQ=255)
                                     tCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGt                                                                                                                         <  2:1435923/137‑1 (MQ=255)
                                     tCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGt                                                                                                                         >  1:1435923/1‑137 (MQ=255)
                                     tCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGAt                                                                                                                      >  2:336870/1‑140 (MQ=255)
                                     tCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGAt                                                                                                                      <  1:336870/140‑1 (MQ=255)
                                       gTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCTACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGtt                                                                                                          <  1:120295/150‑1 (MQ=255)
                                       gTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGtt                                                                                                          <  2:1144680/150‑1 (MQ=255)
                                        tAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTg                                                                                                         >  1:366610/1‑150 (MQ=255)
                                          gCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTc                                                                                                       >  1:395755/1‑150 (MQ=255)
                                             aCCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGAt                                                                                                    >  2:2117688/1‑150 (MQ=255)
                                             aCCGCCTTTGAGTAACTGGCGCGCCACTTACTTGTCGTCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGAt                                                                                                    <  2:633433/150‑1 (MQ=255)
                                              ccGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATg                                                                                                   >  2:417838/1‑150 (MQ=255)
                                                  cTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTg                                                                                                                        <  1:1363104/125‑1 (MQ=255)
                                                  cTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTg                                                                                                                        >  2:1363104/1‑125 (MQ=255)
                                                  cTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGt                                                                                                        <  2:2000826/141‑1 (MQ=255)
                                                  cTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGt                                                                                                        >  1:2000826/1‑141 (MQ=255)
                                                   tttGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGCAGGTGATGCTGCGGCCGTTGTCGATgcgcac                                                                                              >  1:1699685/1‑148 (MQ=255)
                                                        gTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGc                                                                                         >  2:2190391/1‑150 (MQ=255)
                                                         tACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCCGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCg                                                                                        >  2:754243/1‑150 (MQ=255)
                                                         tACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCCCGGTCGAAGGTGATGCTGCGGCCGTTGTCCATGCGCACCATGCg                                                                                        >  2:541056/1‑150 (MQ=255)
                                                            cTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAGGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTa                                                                                     >  2:544205/1‑150 (MQ=255)
                                                            cTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGAt                                                                                                    <  1:505505/135‑1 (MQ=255)
                                                            cTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGAt                                                                                                    >  2:505505/1‑135 (MQ=255)
                                                               atcgagtaCTTCGGCGACACCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTCGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATgg                                                                                  >  2:1751028/8‑150 (MQ=255)
                                                                cgcgcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGc                                                                                 <  1:423759/150‑1 (MQ=255)
                                                                   gcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCg                                                                              >  2:1113118/1‑150 (MQ=255)
                                                                    ccACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGACTGGCGAGCAGGTCACGTGCCCACCATGCCGAGTCTCTTTCTAAGTGGATGCTGCGGCCGTTGTCGATGCGCCCAATACGGTATGGCCCGc                                                                             >  2:274579/1‑150 (MQ=255)
                                                                       cTTCGGCGACACCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGGGGTATGGCCCGCTGc                                                                          >  2:1680535/1‑150 (MQ=255)
                                                                           ggTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGAc                                                                      <  1:972639/150‑1 (MQ=255)
                                                                             tgtCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTg                                                                    <  2:362322/150‑1 (MQ=255)
                                                                              gtCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTgg                                                                   <  2:1254020/150‑1 (MQ=255)
                                                                               tCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGc                                                                  >  1:2082262/1‑150 (MQ=255)
                                                                                 gCCGAAGTACTCGATGCGCAGTTTGTCGAAATTATACCCACCACGGGTAGCAGGCAGGTCACGTGCCCACCATGCGGGGTCGCGTTCGCAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCtt                                                                >  1:2224234/1‑150 (MQ=255)
                                                                                 gCCGAAGTACTCGATGCACAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCtt                                                                >  2:886353/1‑150 (MQ=255)
                                                                                  ccGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTg                                                               <  1:1242094/150‑1 (MQ=255)
                                                                                  ccGAAGTACTCGATGCGCAGTTTGTCGAAACTGTCGAAGGCCTGGAAGAGCGCCTGGCCGTGTGCCCACCAGTCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTg                                                               >  1:1252953/1‑150 (MQ=255)
                                                                                       gTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGaa                                                          <  1:2117688/150‑1 (MQ=255)
                                                                                            cGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCgg                                                     <  2:1049198/150‑1 (MQ=255)
                                                                                                 cgcAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCg                                                <  2:2082262/150‑1 (MQ=255)
                                                                                                    aGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAGGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTg                                             <  1:544205/150‑1 (MQ=255)
                                                                                                     gTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTgg                                            >  1:485015/1‑150 (MQ=255)
                                                                                                     gTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTgg                                            <  1:918363/150‑1 (MQ=255)
                                                                                                         gTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGtt                                                                                                                                >  1:2027356/1‑62 (MQ=255)
                                                                                                         gTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGtt                                                                                                                                <  2:2027356/62‑1 (MQ=255)
                                                                                                          tCGAAATTGTAGCGGCCACGGCTGGCGGGCAGCACACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAg                                       <  2:123497/150‑1 (MQ=255)
                                                                                                            gAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAAt                                                                                           <  1:1963474/96‑1 (MQ=255)
                                                                                                            gAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAAt                                                                                           >  2:1963474/1‑96 (MQ=255)
                                                                                                            gAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCg                                                      >  1:728876/1‑133 (MQ=255)
                                                                                                            gAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCg                                                      <  2:728876/133‑1 (MQ=255)
                                                                                                                ttGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCgcgc                                 <  2:1252953/150‑1 (MQ=255)
                                                                                                                    aGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGGTTGTTGAAGCCGGCAGAGTTGGCAAAGTGGCgcggct                             >  2:1035619/1‑150 (MQ=255)
                                                                                                                     gCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGctgctg                            >  2:892523/1‑150 (MQ=255)
                                                                                                                        gCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAg                         <  2:366610/150‑1 (MQ=255)
                                                                                                                        gCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATCCGGTATGGCACGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAg                         <  1:784075/150‑1 (MQ=255)
                                                                                                                         ccACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTg                                             >  2:1260618/1‑129 (MQ=255)
                                                                                                                         ccACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCAGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTg                                             <  1:1260618/129‑1 (MQ=255)
                                                                                                                          cACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGt                       >  1:836618/1‑150 (MQ=255)
                                                                                                                           aCGGCTGGCCGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTAGGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTc                      <  1:716247/150‑1 (MQ=255)
                                                                                                                             ggCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAg                                                         >  1:1813891/1‑113 (MQ=255)
                                                                                                                             ggCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAg                                                         <  2:1813891/113‑1 (MQ=255)
                                                                                                                              gCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTg                   <  2:2224234/150‑1 (MQ=255)
                                                                                                                               cTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGATGCAGGTCGTGc                  <  1:2289624/150‑1 (MQ=255)
                                                                                                                                 tggtggGCACGTGACCTGCCCGCCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCt                 >  1:1825199/1‑150 (MQ=255)
                                                                                                                                 tGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCt                 >  2:1649600/1‑150 (MQ=255)
                                                                                                                                   gCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCg               <  1:2190391/150‑1 (MQ=255)
                                                                                                                                    tgggCACGTGACCTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGtt                                              >  2:741340/2‑118 (MQ=255)
                                                                                                                                    tgggCACGTGACCTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGtt                                              <  1:741340/117‑1 (MQ=255)
                                                                                                                                    cGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCgg              <  1:1035619/150‑1 (MQ=255)
                                                                                                                                         aGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTAGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCAtt         <  1:1751028/150‑1 (MQ=255)
                                                                                                                                            tCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGTACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTAcc      >  1:1813894/1‑150 (MQ=255)
                                                                                                                                             cACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCg     <  1:1113118/150‑1 (MQ=255)
                                                                                                                                              aCGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCg                                                                                        <  1:2235524/66‑1 (MQ=255)
                                                                                                                                              aCGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCg                                                                                        >  2:2235524/1‑66 (MQ=255)
                                                                                                                                              aCGTGCCCACCATGCCGGGTCGCGCTCGACGGGGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTGGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCgg    >  1:573840/1‑150 (MQ=255)
                                                                                                                                                gTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCGGCa  >  2:2265504/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCTGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCGGCA  >  NC_002947/3463450‑3463742

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