Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,761,425:1 +T 5.3% coding (326/1623 nt) pvdP ← pyoverdine biosynthesis‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,761,4251.T 5.3% 905.6 / 40.4 208coding (326/1623 nt)pvdPpyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base . (97/100);  new base T (8/3);  total (105/103)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

GCCCCGCATTGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGCCAGATG  >  NC_002947/4761276‑4761574
                                                                                                                                                      |                                                                                                                                                     
gccCCGCATTGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTg                                                                                                                                                        <  2:3220808/150‑1 (MQ=255)
   ccGCATTGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCt                                                                                                                                                    >  2:4985067/1‑150 (MQ=255)
    cGCATTGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTc                                                                                                                                                   <  2:5282344/150‑1 (MQ=255)
        ttGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCa                                                                                                                                               <  1:3872653/150‑1 (MQ=255)
         tGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCAc                                                                                                                                              <  2:4218375/150‑1 (MQ=255)
                   cccccACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCaa                                                                                                                                    <  2:1697865/150‑1 (MQ=255)
                    ccccACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAg                                                                                                                                   >  2:2315359/1‑150 (MQ=255)
                    ccccACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAg                                                                                                                                   >  2:204174/1‑150 (MQ=255)
                    ccccACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAg                                                                                                                                   >  1:935040/1‑150 (MQ=255)
                    ccccACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAg                                                                                                                                   <  2:4578312/150‑1 (MQ=255)
                     cccACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAgc                                                                                                                                  <  2:4008616/150‑1 (MQ=255)
                     cccACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAgc                                                                                                                                  >  2:2171243/1‑150 (MQ=255)
                      ccACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAgcg                                                                                                                                 >  2:1535247/1‑150 (MQ=255)
                      ccACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAgcg                                                                                                                                 >  2:569322/1‑150 (MQ=255)
                          acTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTcg                                                                                                                             <  1:4195978/150‑1 (MQ=255)
                              cATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTcggcgg                                                                                                                         <  1:569322/150‑1 (MQ=255)
                                 ccAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTAt                                                                                                                      >  1:3135465/1‑150 (MQ=255)
                                  cAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATc                                                                                                                     <  1:4792645/150‑1 (MQ=255)
                                    gCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGa                                                                                                                   <  2:4161768/150‑1 (MQ=255)
                                     cACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCCATATCGTGGCTCGGCACCGCTGCGCAAGCGCACGGCGGTATCGAg                                                                                                                  >  1:2282759/1‑150 (MQ=255)
                                       cGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGt                                                                                                                >  2:2149820/1‑150 (MQ=255)
                                        gagaTGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTa                                                                                                               <  1:379751/150‑1 (MQ=255)
                                           aTGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAAcc                                                                                                            >  2:1046157/1‑150 (MQ=255)
                                             ggCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACcgc                                                                                                          <  2:4717596/150‑1 (MQ=255)
                                             ggCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACcgc                                                                                                          <  1:4906849/150‑1 (MQ=255)
                                             ggCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACcgc                                                                                                          >  1:2163435/1‑150 (MQ=255)
                                              gCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACcgcg                                                                                                         >  2:344147/1‑150 (MQ=255)
                                               cACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGa                                                                                                        <  2:2803644/150‑1 (MQ=255)
                                               cACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGa                                                                                                        >  1:2226518/1‑150 (MQ=255)
                                                aCTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGaa                                                                                                       >  2:1305674/1‑150 (MQ=255)
                                                 cTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAAt                                                                                                      >  2:4323632/1‑150 (MQ=255)
                                                    gTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAAtgct                                                                                                   >  1:1460745/1‑150 (MQ=255)
                                                           ggCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGccc                                                                                            >  2:2350642/1‑150 (MQ=255)
                                                               tGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCt                                                                                        <  2:2868497/150‑1 (MQ=255)
                                                               tGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCt                                                                                        <  1:4985067/150‑1 (MQ=255)
                                                                 cGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTgc                                                                                      <  1:300770/150‑1 (MQ=255)
                                                                           cGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGa                                                                            <  2:3529338/150‑1 (MQ=255)
                                                                                     gAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTc                                                                  >  1:2197226/1‑150 (MQ=255)
                                                                                     gAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTc                                                                  >  1:4506168/1‑150 (MQ=255)
                                                                                     gAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTc                                                                  >  1:2012448/1‑150 (MQ=255)
                                                                                     gAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGAGGAACAATTc                                                                  >  1:3401043/1‑150 (MQ=255)
                                                                                       aCAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCa                                                                <  2:4453475/150‑1 (MQ=255)
                                                                                       aCAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCa                                                                <  1:2076372/150‑1 (MQ=255)
                                                                                        cAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCac                                                               >  2:141401/1‑149 (MQ=255)
                                                                                         aaGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGc                                                              <  2:5100898/150‑1 (MQ=255)
                                                                                         aaGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGc                                                              <  1:1547593/150‑1 (MQ=255)
                                                                                           gcgcCAGCGTAACCGTCCCGGGTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCTAGCGCTCGGCGGTATCGAGGTAAGCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGt                                                            <  1:3300528/150‑1 (MQ=255)
                                                                                           gcgcCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGt                                                            >  2:3237919/1‑150 (MQ=255)
                                                                                           gcgcCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGt                                                            <  2:2258447/150‑1 (MQ=255)
                                                                                             gcCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGtct                                                          >  2:1018425/1‑150 (MQ=255)
                                                                                               cAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCCCCATCGATGAACAATTCCAGCGTCTCg                                                        >  1:5245384/1‑150 (MQ=255)
                                                                                               cAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCg                                                        >  2:4112835/1‑150 (MQ=255)
                                                                                                 gCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGcc                                                      <  2:2702463/150‑1 (MQ=255)
                                                                                                 gCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGAGGAACAATTCCAGCGTCTCGcc                                                      >  2:2996046/1‑150 (MQ=255)
                                                                                                 gCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGGGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGcc                                                      <  1:4151840/150‑1 (MQ=255)
                                                                                                  cGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGccc                                                     >  1:3526753/1‑150 (MQ=255)
                                                                                                   gTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGCTATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCt                                                    <  2:1958277/150‑1 (MQ=255)
                                                                                                   gTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCt                                                    <  1:375523/150‑1 (MQ=255)
                                                                                                     aaCCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGc                                                  <  1:2897525/150‑1 (MQ=255)
                                                                                                      aCCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATc                                                                              <  1:1386417/121‑1 (MQ=255)
                                                                                                       ccGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCaa                                                <  2:2608983/150‑1 (MQ=255)
                                                                                                          tCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCCAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCAGCAAGcc                                             >  1:1601968/1‑150 (MQ=255)
                                                                                                          tCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGcc                                             >  2:1996609/1‑150 (MQ=255)
                                                                                                           cccGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGccc                                            <  2:2485911/150‑1 (MQ=255)
                                                                                                            ccGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCg                                           >  2:2424815/1‑150 (MQ=255)
                                                                                                            ccGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCg                                           <  2:1934577/150‑1 (MQ=255)
                                                                                                              ggCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGc                                         <  2:1571777/150‑1 (MQ=255)
                                                                                                               gCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCa                                        <  1:1780236/150‑1 (MQ=255)
                                                                                                                cTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCaa                                       <  2:34055/150‑1 (MQ=255)
                                                                                                                 ttGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAg                                      >  1:770220/1‑150 (MQ=255)
                                                                                                                 ttGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAg                                      >  2:3935217/1‑150 (MQ=255)
                                                                                                                  tGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGc                                     <  2:1434117/150‑1 (MQ=255)
                                                                                                                    tCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCg                                   >  2:1665738/1‑150 (MQ=255)
                                                                                                                      ggcaggcaAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGc                                 >  1:2813216/1‑150 (MQ=255)
                                                                                                                       gcaggcaAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGcc                                <  1:1176516/150‑1 (MQ=255)
                                                                                                                          ggcaAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCTCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTgc                             <  1:1979307/150‑1 (MQ=255)
                                                                                                                              aGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGccc                         >  1:3980850/1‑150 (MQ=255)
                                                                                                                               ggtgCCGAGCCACGATATCGGGTGCTCGGCGCCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGcccc                        >  1:2393653/4‑150 (MQ=255)
                                                                                                                               ggtgCCGAGCCACGATATCGGGTGCTCGGCGCCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGcccc                        >  2:3648868/4‑150 (MQ=255)
                                                                                                                               gTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGcccc                        <  1:344147/150‑1 (MQ=255)
                                                                                                                               gTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGcccc                        <  2:5186116/150‑1 (MQ=255)
                                                                                                                                tCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGccccc                       <  1:2228744/150‑1 (MQ=255)
                                                                                                                                  gCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACCTCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCtt                     <  1:4112835/150‑1 (MQ=255)
                                                                                                                                      gagccacGATATCGGGTGCTCGGCGCCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGaa                  >  2:2172382/7‑150 (MQ=255)
                                                                                                                                      gagccacGATATCGGGTGCTCGGCGCCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGaa                  <  1:3648868/144‑1 (MQ=255)
                                                                                                                                        cACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGc                                 <  2:1249988/132‑1 (MQ=255)
                                                                                                                                        cACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGc                                 >  1:1249988/1‑132 (MQ=255)
                                                                                                                                           ccGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGccc                                            <  2:2928975/118‑1 (MQ=255)
                                                                                                                                           ccGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGccc                                            >  1:2928975/1‑118 (MQ=255)
                                                                                                                                           ccGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCa            >  2:1420392/1‑150 (MQ=255)
                                                                                                                                           ccGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCa            >  1:4705012/1‑150 (MQ=255)
                                                                                                                                            cGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAg           >  2:4682213/1‑150 (MQ=255)
                                                                                                                                            cGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAg           >  1:4547702/1‑150 (MQ=255)
                                                                                                                                            cGATATCGGGTGCTAGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCa            >  1:5331351/1‑150 (MQ=255)
                                                                                                                                                atCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGcc       <  1:639749/150‑1 (MQ=255)
                                                                                                                                                atCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGcc       <  1:1857895/150‑1 (MQ=255)
                                                                                                                                                 tCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCACCATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGCCa      >  1:2650773/1‑150 (MQ=255)
                                                                                                                                                   gTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGGAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGCCAGa    <  2:840423/150‑1 (MQ=255)
                                                                                                                                                     ggCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGCCAGATg  <  2:3135465/150‑1 (MQ=255)
                                                                                                                                                     ggCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGCCAGATg  <  1:1046157/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCCCCGCATTGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCGGCAGGCAAGTCGCCGAGCACCCGATATCGTGGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGCCAGATG  >  NC_002947/4761276‑4761574

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