Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,463,587 G→C 8.6% L272V (CTG→GTG)  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,5870GC 8.6% 280.1 / 14.5 93L272V (CTG→GTG) PP_3078ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (44/41);  new base C (4/4);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

CGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATT  >  NC_002947/3463440‑3463735
                                                                                                                                                    |                                                                                                                                                    
cGTTGTAGCTCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGAGGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGt                                                                                                                                                    >  2:1739286/1‑150 (MQ=255)
 gTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTc                                                                                                                                                   >  1:880613/1‑150 (MQ=255)
 gTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGA‑TGGCGAGCAGGTc                                                                                                                                                   >  2:784075/1‑150 (MQ=255)
         ccAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAg                                                                                                                                                      <  2:1680604/139‑1 (MQ=255)
         ccAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAg                                                                                                                                                      >  1:1680604/1‑139 (MQ=255)
          cAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGTGGGCACGTGACCTGcccg                                                                                                                                          >  2:1669514/1‑149 (MQ=255)
          cAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCcca                                                                                                                                          >  2:1485267/1‑150 (MQ=255)
          cAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCATAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCcca                                                                                                                                          >  1:362322/1‑150 (MQ=255)
           aaGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAgg                                                                                                                                                     >  2:1915420/1‑138 (MQ=255)
           aaGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAgg                                                                                                                                                     <  1:1915420/138‑1 (MQ=255)
           aaGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCccac                                                                                                                                         >  2:972639/1‑150 (MQ=255)
            aGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCccacc                                                                                                                                        <  1:906226/150‑1 (MQ=255)
             ggTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGTGGGCACGTGACCTGCccgcca                                                                                                                                       <  1:1669514/150‑1 (MQ=255)
              gtgtAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCATGGTGGGCACGTGACCTGCccgcca                                                                                                                                       <  1:548289/150‑1 (MQ=255)
                   ggCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCgg                                                                                                                                 >  2:918363/1‑150 (MQ=255)
                    gCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTCCCTGGCGCGCCACTTCGGTGTCGCCGGAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCggg                                                                                                                                <  2:1895847/150‑1 (MQ=255)
                           cGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGtt                                                                                                                         <  2:880613/150‑1 (MQ=255)
                            ggAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGGCGCGTTc                                                                                                                        >  2:2056019/1‑150 (MQ=255)
                             gagaACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCg                                                                                                                       >  2:423759/1‑150 (MQ=255)
                               gaACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGaa                                                                                                                     <  1:288036/150‑1 (MQ=255)
                                aaCTCGCGGTTGTAGTCGTCGCCACCGCATTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAg                                                                                                                    <  1:125883/150‑1 (MQ=255)
                                 aCTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAgg                                                                                                                   <  2:152443/150‑1 (MQ=255)
                                   tCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGTGGGCACGTGACCTGCccgcca                                                                                                                                       <  2:411899/128‑1 (MQ=255)
                                   tCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGTGGGCACGTGACCTGCccgcca                                                                                                                                       >  1:411899/1‑128 (MQ=255)
                                    cgcgGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGa                                                                                                                <  1:1485267/150‑1 (MQ=255)
                                        gTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGAtgct                                                                                                            >  2:2289624/1‑150 (MQ=255)
                                        gTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGATATGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGAtgct                                                                                                            >  1:633433/1‑150 (MQ=255)
                                               tCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGt                                                                                                                  <  2:1435923/137‑1 (MQ=255)
                                               tCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGt                                                                                                                  >  1:1435923/1‑137 (MQ=255)
                                               tCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGAt                                                                                                               <  1:336870/140‑1 (MQ=255)
                                               tCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGAt                                                                                                               >  2:336870/1‑140 (MQ=255)
                                                 gTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCTACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGtt                                                                                                   <  1:120295/150‑1 (MQ=255)
                                                 gTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGtt                                                                                                   <  2:1144680/150‑1 (MQ=255)
                                                  tAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTg                                                                                                  >  1:366610/1‑150 (MQ=255)
                                                    gCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTc                                                                                                >  1:395755/1‑150 (MQ=255)
                                                      cACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGt                                                                                                                                                    <  2:75296/96‑1 (MQ=255)
                                                      cACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGt                                                                                                                                                    >  1:75296/1‑96 (MQ=255)
                                                       aCCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGAt                                                                                             >  2:2117688/1‑150 (MQ=255)
                                                       aCCGCCTTTGAGTAACTGGCGCGCCACTTACTTGTCGTCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGAt                                                                                             <  2:633433/150‑1 (MQ=255)
                                                        ccGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATg                                                                                            >  2:417838/1‑150 (MQ=255)
                                                            cTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTg                                                                                                                 <  1:1363104/125‑1 (MQ=255)
                                                            cTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTg                                                                                                                 >  2:1363104/1‑125 (MQ=255)
                                                            cTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGt                                                                                                 <  2:2000826/141‑1 (MQ=255)
                                                            cTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGt                                                                                                 >  1:2000826/1‑141 (MQ=255)
                                                             tttGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGCAGGTGATGCTGCGGCCGTTGTCGATgcgcac                                                                                       >  1:1699685/1‑148 (MQ=255)
                                                                  gTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGc                                                                                  >  2:2190391/1‑150 (MQ=255)
                                                                   tACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCCGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCg                                                                                 >  2:754243/1‑150 (MQ=255)
                                                                   tACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCCCGGTCGAAGGTGATGCTGCGGCCGTTGTCCATGCGCACCATGCg                                                                                 >  2:541056/1‑150 (MQ=255)
                                                                      cTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAGGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTa                                                                              >  2:544205/1‑150 (MQ=255)
                                                                      cTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGAt                                                                                             <  1:505505/135‑1 (MQ=255)
                                                                      cTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGAt                                                                                             >  2:505505/1‑135 (MQ=255)
                                                                         atcgagtaCTTCGGCGACACCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTCGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATgg                                                                           >  2:1751028/8‑150 (MQ=255)
                                                                          cgcgcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGc                                                                          <  1:423759/150‑1 (MQ=255)
                                                                             gcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCg                                                                       >  2:1113118/1‑150 (MQ=255)
                                                                              ccACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGAC‑TGGCGAGCAGGTCACGTGCCCACCATGCCGAGTCTCTTTCTAAGTGGATGCTGCGGCCGTTGTCGATGCGCCCAATACGGTATGGCCCGc                                                                      >  2:274579/1‑150 (MQ=255)
                                                                                 cTTCGGCGACACCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGGGGTATGGCCCGCTGc                                                                   >  2:1680535/1‑150 (MQ=255)
                                                                                     ggTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGAc                                                               <  1:972639/150‑1 (MQ=255)
                                                                                       tgtCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTg                                                             <  2:362322/150‑1 (MQ=255)
                                                                                        gtCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTgg                                                            <  2:1254020/150‑1 (MQ=255)
                                                                                         tCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGc                                                           >  1:2082262/1‑150 (MQ=255)
                                                                                           gCCGAAGTACTCGATGCGCAGTTTGTCGAAATTATACCCACCACGGG‑TAGCAGGCAGGTCACGTGCCCACCATGCGGGGTCGCGTTCGCAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCtt                                                         >  1:2224234/1‑150 (MQ=255)
                                                                                           gCCGAAGTACTCGATGCACAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCtt                                                         >  2:886353/1‑150 (MQ=255)
                                                                                            ccGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTg                                                        <  1:1242094/150‑1 (MQ=255)
                                                                                            ccGAAGTACTCGATGCGCAGTTTGTCGAAACTGTCGAAGGCCTGGA‑AGAGCGCCTGGCCGTGTGCCCACCAGTCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTg                                                        >  1:1252953/1‑150 (MQ=255)
                                                                                                 gTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGaa                                                   <  1:2117688/150‑1 (MQ=255)
                                                                                                      cGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCgg                                              <  2:1049198/150‑1 (MQ=255)
                                                                                                           cgcAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCg                                         <  2:2082262/150‑1 (MQ=255)
                                                                                                              aGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAGGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTg                                      <  1:544205/150‑1 (MQ=255)
                                                                                                               gTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTgg                                     >  1:485015/1‑150 (MQ=255)
                                                                                                               gTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTgg                                     <  1:918363/150‑1 (MQ=255)
                                                                                                                   gTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGtt                                                                                                                         <  2:2027356/62‑1 (MQ=255)
                                                                                                                   gTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGtt                                                                                                                         >  1:2027356/1‑62 (MQ=255)
                                                                                                                    tCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGCACACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAg                                <  2:123497/150‑1 (MQ=255)
                                                                                                                      gAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAAt                                                                                    <  1:1963474/96‑1 (MQ=255)
                                                                                                                      gAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAAt                                                                                    >  2:1963474/1‑96 (MQ=255)
                                                                                                                      gAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCg                                               <  2:728876/133‑1 (MQ=255)
                                                                                                                      gAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCg                                               >  1:728876/1‑133 (MQ=255)
                                                                                                                          ttGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCgcgc                          <  2:1252953/150‑1 (MQ=255)
                                                                                                                              aGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGGTTGTTGAAGCCGGCAGAGTTGGCAAAGTGGCgcggct                      >  2:1035619/1‑150 (MQ=255)
                                                                                                                               gCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGctgctg                     >  2:892523/1‑150 (MQ=255)
                                                                                                                                  gCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAg                  <  2:366610/150‑1 (MQ=255)
                                                                                                                                  gCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATCCGGTATGGCACGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAg                  <  1:784075/150‑1 (MQ=255)
                                                                                                                                   ccACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTg                                      >  2:1260618/1‑129 (MQ=255)
                                                                                                                                   ccACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCAGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTg                                      <  1:1260618/129‑1 (MQ=255)
                                                                                                                                    cACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGt                >  1:836618/1‑150 (MQ=255)
                                                                                                                                     aCGGC‑TGGCCGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTAGGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTc               <  1:716247/150‑1 (MQ=255)
                                                                                                                                       ggC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAg                                                  <  2:1813891/113‑1 (MQ=255)
                                                                                                                                       ggC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAg                                                  >  1:1813891/1‑113 (MQ=255)
                                                                                                                                        gC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTg            <  2:2224234/150‑1 (MQ=255)
                                                                                                                                         c‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGATGCAGGTCGTGc           <  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)
                                                                                                                                                    |                                                                                                                                                    
CGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATT  >  NC_002947/3463440‑3463735

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: ‑