Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,763,419 G→T 5.6% V453V (GTC→GTA) ‡ pvdP ← pyoverdine biosynthesis‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,763,4190GT 5.6% 1586.1 / 53.8 466V453V (GTC→GTA) ‡pvdPpyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base G (204/235);  new base T (15/11);  total (220/246)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 470 total aligned reads displayed.

TGGTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGCGCCGC  >  NC_002947/4763271‑4763568
                                                                                                                                                    |                                                                                                                                                     
 ggTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGAc                                                                                                                                                     <  1:183580/150‑1 (MQ=255)
  gTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACt                                                                                                                                                    >  2:4974096/1‑150 (MQ=255)
        ggATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCa                                                                                                                                              >  2:5198524/1‑150 (MQ=255)
        ggATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCa                                                                                                                                              >  1:3816809/1‑150 (MQ=255)
            ccATCTCCCCGCTAAGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCTAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGt                                                                                                                                          <  2:2090841/150‑1 (MQ=255)
            ccATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGt                                                                                                                                          <  2:849191/150‑1 (MQ=255)
            ccATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGt                                                                                                                                          <  2:6110428/150‑1 (MQ=255)
            ccATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCCCCTGGCGCAGTCGTGCGCCAggag                                                                                                                                          <  1:894569/150‑3 (MQ=255)
             cATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTg                                                                                                                                         <  2:5831196/150‑1 (MQ=255)
             cATCTCCAAGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTg                                                                                                                                         <  1:1693976/150‑1 (MQ=255)
              aTCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCCAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTgg                                                                                                                                        >  1:4209754/1‑150 (MQ=255)
              aTCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTgg                                                                                                                                        >  2:6282431/1‑150 (MQ=255)
                ctcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGaa                                                                                                                                      <  2:5746472/150‑1 (MQ=255)
                 tcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAAt                                                                                                                                     >  2:6841900/1‑150 (MQ=255)
                 tcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAAt                                                                                                                                     >  2:3109419/1‑150 (MQ=255)
                 tcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAAt                                                                                                                                     <  2:1955295/150‑1 (MQ=255)
                 tcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAAt                                                                                                                                     >  2:5908279/1‑150 (MQ=255)
                 tcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGGGGAAt                                                                                                                                     >  1:5552360/1‑150 (MQ=255)
                  ccACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATc                                                                                                                                    >  2:3368035/1‑150 (MQ=255)
                  ccACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATc                                                                                                                                    <  1:6222048/150‑1 (MQ=255)
                   cACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCg                                                                                                                                   >  1:5162558/1‑150 (MQ=255)
                   cACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCg                                                                                                                                   >  1:3875571/1‑150 (MQ=255)
                    aCGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGACTCGt                                                                                                                                  >  1:5419183/1‑150 (MQ=255)
                    aCGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGt                                                                                                                                  >  2:289998/1‑150 (MQ=255)
                      gCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTg                                                                                                                                         >  1:609647/1‑141 (MQ=255)
                      gCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTc                                                                                                                                <  2:5620660/150‑1 (MQ=255)
                        tACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCAt                                                                                                                              >  1:5994583/1‑150 (MQ=255)
                        tACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCAt                                                                                                                              >  2:5981503/1‑150 (MQ=255)
                           gCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATggg                                                                                                                           >  2:6141682/1‑150 (MQ=255)
                           gCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATggg                                                                                                                           >  1:6510894/1‑150 (MQ=255)
                           gCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATggg                                                                                                                           <  1:4341244/150‑1 (MQ=255)
                             ccTGGTCGCAGCCCCGCAACGTAACTGCAGCCGAGCGTGTCAGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCAGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGgcg                                                                                                                         >  2:2682981/1‑150 (MQ=255)
                              cTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGgcgc                                                                                                                        <  1:6343917/150‑1 (MQ=255)
                              cTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGGTCATGGgcgc                                                                                                                        >  2:1174412/1‑150 (MQ=255)
                                 ggCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGa                                                                                                                     <  2:3519207/150‑1 (MQ=255)
                                  gcgcAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATgg                                                                                                                            >  2:430423/1‑142 (MQ=255)
                                   cgcAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAAc                                                                                                                   <  2:3813974/150‑1 (MQ=255)
                                   cgcAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAAc                                                                                                                   <  1:39711/150‑1 (MQ=255)
                                     cAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCa                                                                                                                 <  1:292018/150‑1 (MQ=255)
                                     cAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCa                                                                                                                 <  1:4415609/150‑1 (MQ=255)
                                           cGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCCGGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTc                                                                                                           <  1:6281579/150‑1 (MQ=255)
                                               aCGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGaca                                                                                                       <  2:6577072/150‑1 (MQ=255)
                                                   tACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgt                                                                                                   <  1:753441/150‑1 (MQ=255)
                                                    aCTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtc                                                                                                  <  2:185531/150‑1 (MQ=255)
                                                      tGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgt                                                                                                >  2:3331657/1‑150 (MQ=255)
                                                      tGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgt                                                                                                >  1:5292860/1‑150 (MQ=255)
                                                      tGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgt                                                                                                >  1:6042781/1‑150 (MQ=255)
                                                        cAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGtcgtcg                                                                                              >  1:26019/1‑150 (MQ=255)
                                                         aGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGa                                                                                             >  2:3959906/1‑150 (MQ=255)
                                                           ccGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATc                                                                                           <  2:4185664/150‑1 (MQ=255)
                                                           ccGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATc                                                                                           <  2:1816756/150‑1 (MQ=255)
                                                           ccGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATc                                                                                           >  2:5114201/1‑150 (MQ=255)
                                                               gcgtgTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGc                                                                                       >  2:242017/1‑150 (MQ=255)
                                                               gcgtgTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCACGGAACCAGTCCTCGACACGGTCGTCGATCCAGc                                                                                       >  1:6528382/1‑150 (MQ=255)
                                                                  tgtCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCAt                                                                                    <  1:1415482/150‑1 (MQ=255)
                                                                     ccGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATc                                                                                                                                    >  1:5811661/1‑99 (MQ=255)
                                                                         ccccAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAgcg                                                                             <  1:4974096/150‑1 (MQ=255)
                                                                             aGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGACCCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAg                                                                         >  2:1610152/1‑150 (MQ=255)
                                                                             aGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAg                                                                         >  2:1371578/1‑150 (MQ=255)
                                                                             aGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAg                                                                         >  1:3645374/1‑150 (MQ=255)
                                                                             aGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAg                                                                         >  1:6611464/1‑150 (MQ=255)
                                                                               ccAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGaa                                                                       >  1:4518617/1‑150 (MQ=255)
                                                                               acgtggATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGaa                                                                       <  1:2750307/146‑1 (MQ=255)
                                                                                cAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAac                                                                      >  1:3507340/1‑150 (MQ=255)
                                                                                 aaGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGACGCGCCCGaccc                                                                     >  2:1214859/1‑147 (MQ=255)
                                                                                 aaGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAaca                                                                     <  2:2305213/150‑1 (MQ=255)
                                                                                 aaGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAaca                                                                     >  1:6413487/1‑150 (MQ=255)
                                                                                      tCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAggg                                                                <  1:2393583/150‑1 (MQ=255)
                                                                                        ggCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGtt                                                              <  1:401054/150‑1 (MQ=255)
                                                                                        ggCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGtt                                                              <  1:1956543/150‑1 (MQ=255)
                                                                                              cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTg                                                        >  1:4665390/1‑150 (MQ=255)
                                                                                              cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTg                                                        >  1:3141944/1‑150 (MQ=255)
                                                                                                ccACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACTGGGTTGACGTggg                                                      <  2:6440659/150‑1 (MQ=255)
                                                                                                ccACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTggg                                                      <  1:2645241/150‑1 (MQ=255)
                                                                                                   gttcgcgcCGGGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAag                                                   >  1:4826494/5‑150 (MQ=255)
                                                                                                         gcCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACg                                             <  2:2635711/150‑1 (MQ=255)
                                                                                                          ccccGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTAGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACgg                                            <  2:2189436/150‑1 (MQ=255)
                                                                                                             cGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGACGAGAACGGGTc                                         >  2:481930/1‑150 (MQ=255)
                                                                                                              ggCGCGAACCAGGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTcc                                        <  2:1023905/150‑1 (MQ=255)
                                                                                                                      ccATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGaaa                                <  1:3610111/150‑1 (MQ=255)
                                                                                                                      ccATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGaaa                                <  1:5213792/150‑1 (MQ=255)
                                                                                                                      ccATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGaaa                                <  1:2167676/150‑1 (MQ=255)
                                                                                                                      ccATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGaaa                                <  2:1075650/150‑1 (MQ=255)
                                                                                                                      ccATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGaaa                                <  1:1264649/150‑1 (MQ=255)
                                                                                                                      ccATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGaaa                                <  1:4491760/150‑1 (MQ=255)
                                                                                                                          ggCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCg                            >  2:2838360/1‑150 (MQ=255)
                                                                                                                            cACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGtt                          <  1:3865406/150‑1 (MQ=255)
                                                                                                                              cGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCt                        <  1:3362280/150‑1 (MQ=255)
                                                                                                                                ccATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgc                      <  2:4890735/150‑1 (MQ=255)
                                                                                                                                    tGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcgcgc                  <  2:3631346/150‑1 (MQ=255)
                                                                                                                                      aCCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCa                <  1:5670038/150‑1 (MQ=255)
                                                                                                                                        ctcctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATa              >  1:326625/1‑150 (MQ=255)
                                                                                                                                        accctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATa              >  2:165597/3‑150 (MQ=255)
                                                                                                                                          cctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgtc                                                                                               <  1:1197648/67‑1 (MQ=37)
                                                                                                                                           ctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCa           >  1:2273781/1‑150 (MQ=255)
                                                                                                                                           ctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCa           <  1:104516/150‑1 (MQ=255)
                                                                                                                                              gcgcACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAgcg        >  2:408487/1‑150 (MQ=255)
                                                                                                                                               cgcACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAgcgc       <  2:2455936/150‑1 (MQ=255)
                                                                                                                                                  aCGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGcgcc    >  2:4634680/1‑150 (MQ=255)
                                                                                                                                                  aCGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGcgcc    <  1:6841900/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGCGCCGC  >  NC_002947/4763271‑4763568

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