Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,020,784 T→A 6.7% intergenic (‑50/+184) PP_5679 ← / ← PP_4421 hypothetical protein/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,020,7840TA 6.7% 730.5 / 29.6 225intergenic (‑50/+184)PP_5679/PP_4421hypothetical protein/aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (98/112);  new base A (7/8);  total (105/120)
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.35e-01

Only 100 of 226 total aligned reads displayed.

GTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGG  >  NC_002947/5020638‑5020932
                                                                                                                                                  |                                                                                                                                                     
   gCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCgg                                                                                                                                                 >  1:4960899/1‑150 (MQ=255)
      gCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGgcgc                                                                                                                                              >  2:1653718/1‑150 (MQ=255)
       cTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGcgcg                                                                                                                                             <  1:5048626/150‑1 (MQ=255)
        ttGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCgg                                                                                                                                            <  2:2144996/150‑1 (MQ=255)
         tGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCGCCGGAGACCGGTGCggg                                                                                                                                           <  2:4895824/150‑1 (MQ=255)
            agagTGCCTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCGCCGGAGACCGGTGCGGGCtt                                                                                                                                        <  1:3890686/150‑1 (MQ=255)
              agTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCtttt                                                                                                                                      <  1:3309961/150‑1 (MQ=255)
               gTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCACACCGGTCTCCGGCGCGGGCTTTTa                                                                                                                                     <  1:2253442/150‑1 (MQ=255)
                    tgctgGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTc                                                                                                                                <  1:5364834/150‑1 (MQ=255)
                     gctgGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTcc                                                                                                                               <  2:4781965/150‑1 (MQ=255)
                       tgGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCgg                                                                                                                             >  2:3673498/1‑150 (MQ=255)
                         ggCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGc                                                                                                                           <  1:1206042/150‑1 (MQ=255)
                         ggCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGc                                                                                                                           >  2:4543864/1‑150 (MQ=255)
                           cACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCATTTACTTTCCGGGCCt                                                                                                                         <  2:2057279/150‑1 (MQ=255)
                            aCGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCtt                                                                                                                        <  1:4256651/150‑1 (MQ=255)
                              gAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTgg                                                                                                                      >  2:2892254/1‑150 (MQ=255)
                                 gCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACa                                                                                                                   <  2:3795694/150‑1 (MQ=255)
                                 gCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACa                                                                                                                   <  2:4495704/150‑1 (MQ=255)
                                     ctcGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCGTCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTAt                                                                                                               <  2:1419728/150‑1 (MQ=255)
                                     ctcGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTAt                                                                                                               >  2:398436/1‑150 (MQ=255)
                                      tcGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATc                                                                                                              <  2:1393659/150‑1 (MQ=255)
                                      tcGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATc                                                                                                              <  2:289836/150‑1 (MQ=255)
                                           aGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGAcc                                                                                                         <  1:1972352/150‑1 (MQ=255)
                                              ccTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGa                                                                                                      <  1:2146253/150‑1 (MQ=255)
                                              ccTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCACGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGa                                                                                                      >  2:2913568/1‑150 (MQ=255)
                                                tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGaac                                                                                                    >  1:1260267/1‑149 (MQ=255)
                                                tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGaaa                                                                                                    <  1:5466400/150‑1 (MQ=255)
                                                tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGaaa                                                                                                    >  2:3829771/1‑150 (MQ=255)
                                                tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGaaa                                                                                                    <  2:2254933/150‑1 (MQ=255)
                                                tGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGaaa                                                                                                    <  2:4256518/150‑1 (MQ=255)
                                                    cGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTAc                                                                                                <  2:575381/150‑1 (MQ=255)
                                                     gaaaaGGTAGCCACTGCTTATTTGATCGGAGCCAGCTCCGAAATACTAGCGCTCTCGTATGAGTCCACGATGAGATGAATAAGCCCGCAACTGTCTCCGGCGCGGGCTTTAACTTTCCGCGCCTTGGACAGTATCTGACCCGAAAGTACg                                                                                               >  1:3218482/5‑150 (MQ=255)
                                                          ggTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCtt                                                                                          <  1:1788342/150‑1 (MQ=255)
                                                           gTAGCCACTGCTGCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTCTCTGACCCGACAGTACGTCCTtg                                                                                         >  2:1023442/1‑150 (MQ=255)
                                                            tAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTtgt                                                                                        <  1:3689025/150‑1 (MQ=255)
                                                            tAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTtgt                                                                                        <  2:1319214/150‑1 (MQ=255)
                                                            tAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTtgt                                                                                        <  2:3492177/150‑1 (MQ=255)
                                                            tAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGACCTtgt                                                                                        >  1:2839746/1‑150 (MQ=255)
                                                               ccACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGaa                                                                                     <  2:1839240/150‑1 (MQ=255)
                                                               ccACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGaa                                                                                     <  1:4627511/150‑1 (MQ=255)
                                                                cACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTGTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAAc                                                                                    >  1:4904491/1‑150 (MQ=255)
                                                                cACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAAc                                                                                    <  1:398348/150‑1 (MQ=255)
                                                                cACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAAc                                                                                    >  1:3187007/1‑150 (MQ=255)
                                                                cACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAAc                                                                                    >  1:1736586/1‑150 (MQ=255)
                                                                  cTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtg                                                                                  >  1:3330125/1‑150 (MQ=255)
                                                                    gCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACtgat                                                                                <  1:1415320/150‑1 (MQ=255)
                                                                         tttGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAg                                                                           >  1:551660/1‑150 (MQ=255)
                                                                          ttGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGc                                                                          >  1:1034699/1‑150 (MQ=255)
                                                                            gATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCttt                                                                                                                                       >  1:266944/1‑87 (MQ=255)
                                                                            gATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGt                                                                        <  1:1005223/150‑1 (MQ=255)
                                                                             aTCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTc                                                                       >  2:3621394/1‑150 (MQ=255)
                                                                             aTCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTc                                                                       <  2:4411976/150‑1 (MQ=255)
                                                                             aTCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTc                                                                       <  2:2615672/150‑1 (MQ=255)
                                                                             aTCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTc                                                                       <  1:3646518/150‑1 (MQ=255)
                                                                              tCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCt                                                                      >  2:2614035/1‑150 (MQ=255)
                                                                              tCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCt                                                                      >  1:2565744/1‑150 (MQ=255)
                                                                               cGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCtt                                                                     <  2:881114/150‑1 (MQ=255)
                                                                               cGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCtt                                                                     <  1:856165/150‑1 (MQ=255)
                                                                               cGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCtt                                                                     <  1:3555582/150‑1 (MQ=255)
                                                                                 gAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTgg                                                                   >  2:3901866/1‑150 (MQ=255)
                                                                                  aGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGc                                                                  <  1:4810304/150‑1 (MQ=255)
                                                                                     cAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCttt                                                               <  2:4321298/150‑1 (MQ=255)
                                                                                             tCATCCTAGCGCTCTCGTTTGATTACACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGTGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGTCTTTGATTGTGa                                                       <  2:3218482/150‑1 (MQ=255)
                                                                                             tCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGa                                                       >  2:3375329/1‑150 (MQ=255)
                                                                                               aTCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATa                                                     <  2:4680016/150‑1 (MQ=255)
                                                                                                tCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAg                                                    >  1:3426278/1‑150 (MQ=255)
                                                                                                  cTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCa                                                  >  2:1553289/1‑150 (MQ=255)
                                                                                                   tAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCAc                                                 >  1:1566870/1‑150 (MQ=255)
                                                                                                    aGCGCTCTCGATAGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCAATTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACt                                                <  2:5132263/150‑1 (MQ=255)
                                                                                                     gcgcTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTg                                                                    >  2:2806055/1‑129 (MQ=255)
                                                                                                     gcgcTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTg                                                                    <  1:2806055/129‑1 (MQ=255)
                                                                                                         tctcGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaa                                          >  1:5194232/1‑150 (MQ=255)
                                                                                                         tctcGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaa                                          >  2:5399448/1‑150 (MQ=255)
                                                                                                           tcGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaa                                        >  2:660512/1‑150 (MQ=255)
                                                                                                           tcGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaa                                        <  1:4115351/150‑1 (MQ=255)
                                                                                                             gTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaaa                                      <  2:1037895/150‑1 (MQ=255)
                                                                                                               ttGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaaaaa                                    >  1:2628181/1‑150 (MQ=255)
                                                                                                                  aGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAg                                 >  1:3507077/1‑150 (MQ=255)
                                                                                                                     ccACGATGATATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGccc                              <  1:3452106/150‑1 (MQ=255)
                                                                                                                     ccACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACATTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGccc                              <  2:4960899/150‑1 (MQ=255)
                                                                                                                     ccACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaaaaaaaagcc                              >  1:1451726/1‑146 (MQ=255)
                                                                                                                      cACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAATAAAATAGCCCg                             >  2:937559/1‑150 (MQ=255)
                                                                                                                      cACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCg                             >  1:4522077/1‑150 (MQ=255)
                                                                                                                               atgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGc                    >  1:140739/1‑150 (MQ=255)
                                                                                                                                 gtaTAAGCCCGCACCGGTCTCGGGCGCGGGCTTTTACTTTCCGGGCCTTGGACTGTATCTGACCCGAAATTACGTCCTTGTGAAGTGTTGACAGCGTCTTGGCTTTGATTGTGATTGCAGTACAAAAATAAAATGCCCGCCTAAAGGCgg                  <  1:696008/148‑1 (MQ=255)
                                                                                                                                 gaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATCGCACTACAAAAAAAAAAAAGCCCGCCTAAAGGCg                   <  1:439228/150‑1 (MQ=255)
                                                                                                                                 gaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCgg                  >  1:5304849/1‑150 (MQ=255)
                                                                                                                                 gaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCggg                 >  1:4577982/1‑150 (MQ=255)
                                                                                                                                  aaTAAGCCCGCACCGGTCTCCGGCGCGGTCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCTAAAGTACGTCCTTGTGAACTGATGACAGATTATTGTCTTTTATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCggg                 <  1:1104696/150‑1 (MQ=255)
                                                                                                                                  aaTAAGCCCGCACCGGTATCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCggg                 <  2:412116/150‑1 (MQ=255)
                                                                                                                                     aaGCCCGCGCCGGAGACCGGTGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACa                                            <  1:4651817/120‑1 (MQ=255)
                                                                                                                                     aaGCCCGCGCCGGAGACCGGTGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACa                                            >  2:4651814/1‑120 (MQ=255)
                                                                                                                                       gCCCGCGCCGGAGACCGGTGCGGGCTTTTACTTTCCGGGCCTTGGACCGTAGCTGACCCGAAAGGACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCGAACGGCGGGCGTa             >  2:3560139/1‑150 (MQ=255)
                                                                                                                                       gCCCGCGCCGGAGACCGGTGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGa                                                       >  2:1628986/1‑108 (MQ=39)
                                                                                                                                         ccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTaaaa          <  1:3829771/150‑1 (MQ=255)
                                                                                                                                             aCCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGaa      <  2:5541362/150‑1 (MQ=255)
                                                                                                                                              ccGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGaa      <  1:3673498/150‑1 (MQ=255)
                                                                                                                                               cGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAaga    <  2:3347589/150‑1 (MQ=255)
                                                                                                                                               cGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACCGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAaga    <  1:2896178/150‑1 (MQ=255)
                                                                                                                                               cGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGAACGTCCTGGTGAACTGGTGACAGCGTCGTGGCTTTGATTGTGAGAGCGCTACCACAAAAAATAGCCCGCCGACAGGCGGGCATAACACGGagg    >  2:5418221/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTAGCCGCTTGCAGAGTGACTGCTGGGCACGAAGCTTCTCGGCAGCCCTGGACGTATAGGTAGCCACTGCTTCTTTGATCGGAGCCAGCTCGGTCATCCTAGCGCTCTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGG  >  NC_002947/5020638‑5020932

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