Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,763,418 C→G 5.6% V454L (GTG→CTG)  pvdP ← pyoverdine biosynthesis‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,763,4180CG 5.6% 1587.9 / 57.4 466V454L (GTG→CTG) pvdPpyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base C (205/235);  new base G (15/11);  total (220/246)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 470 total aligned reads displayed.

CATGGTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGCGCC  >  NC_002947/4763269‑4763566
                                                                                                                                                     |                                                                                                                                                    
  tGGTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGa                                                                                                                                                    >  1:6405638/1‑150 (MQ=255)
   ggTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGAc                                                                                                                                                   <  1:183580/150‑1 (MQ=255)
    gTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACt                                                                                                                                                  >  2:4974096/1‑150 (MQ=255)
          ggATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCa                                                                                                                                            >  2:5198524/1‑150 (MQ=255)
          ggATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCa                                                                                                                                            >  2:6271143/1‑150 (MQ=255)
              ccATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGt                                                                                                                                        <  2:6110428/150‑1 (MQ=255)
               cATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTg                                                                                                                                       <  2:5831196/150‑1 (MQ=255)
               cATCTCCAAGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTg                                                                                                                                       <  1:1693976/150‑1 (MQ=255)
                aTCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTgg                                                                                                                                      >  1:356402/1‑150 (MQ=255)
                 tctcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGa                                                                                                                                     <  1:3858314/150‑1 (MQ=255)
                  ctcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGaa                                                                                                                                    <  2:5746472/150‑1 (MQ=255)
                   tcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAAt                                                                                                                                   >  2:5908279/1‑150 (MQ=255)
                   tcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAAt                                                                                                                                   >  2:3109419/1‑150 (MQ=255)
                   tcCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGGGGAAt                                                                                                                                   >  1:5552360/1‑150 (MQ=255)
                    ccACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATc                                                                                                                                  >  2:3368035/1‑150 (MQ=255)
                    ccACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATc                                                                                                                                  <  1:890737/150‑1 (MQ=255)
                     cACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCg                                                                                                                                 >  1:5162558/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)
                         cTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCa                                                                                                                             <  1:6522145/150‑1 (MQ=255)
                          tACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCAt                                                                                                                            >  1:3874545/1‑150 (MQ=255)
                          tACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCAt                                                                                                                            >  2:5981503/1‑150 (MQ=255)
                          tACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCAt                                                                                                                            >  1:5994583/1‑150 (MQ=255)
                             gCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATggg                                                                                                                         >  2:6141682/1‑150 (MQ=255)
                               ccTGGTCGCAGCCCCGCAACGTAACTGCAGCCGAGCGTGTCAGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCAGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGgcg                                                                                                                       >  2:2682981/1‑150 (MQ=255)
                                cTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGgcgc                                                                                                                      <  1:4318996/150‑1 (MQ=255)
                                   ggCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGa                                                                                                                   <  2:3519207/150‑1 (MQ=255)
                                    gcgcAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATgg                                                                                                                          >  2:430423/1‑142 (MQ=255)
                                     cgcAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAAc                                                                                                                 <  2:3813974/150‑1 (MQ=255)
                                       cAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCa                                                                                                               <  1:292018/150‑1 (MQ=255)
                                        agtcccGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCAc                                                                                                                                                      >  1:6190392/4‑110 (MQ=255)
                                        aGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAg                                                                                                              >  2:1123068/1‑150 (MQ=255)
                                          ccccGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTc                                                                                                            >  1:6435445/1‑150 (MQ=255)
                                             cGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTc                                                                                                         <  1:2603108/150‑1 (MQ=255)
                                             cGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTc                                                                                                         <  1:4461304/150‑1 (MQ=255)
                                             cGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTc                                                                                                         <  2:6560243/150‑1 (MQ=255)
                                              gCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCg                                                                                                        >  2:100931/1‑150 (MQ=255)
                                                     tACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgt                                                                                                 <  1:753366/150‑1 (MQ=255)
                                                        tGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGaa                                                                                                                                    <  2:1201935/112‑1 (MQ=255)
                                                        tGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgt                                                                                              >  1:6042781/1‑150 (MQ=255)
                                                        tGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgt                                                                                              >  1:5292860/1‑150 (MQ=255)
                                                        tGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgt                                                                                              >  2:3331657/1‑150 (MQ=255)
                                                        tGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgt                                                                                              <  1:3641230/150‑1 (MQ=255)
                                                         gCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgtc                                                                                             >  1:4401511/1‑150 (MQ=255)
                                                         gCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgtc                                                                                             <  2:2413528/150‑1 (MQ=255)
                                                          cAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGtcgtcg                                                                                            >  2:105311/1‑150 (MQ=255)
                                                           aGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGa                                                                                           >  2:3959906/1‑150 (MQ=255)
                                                           aGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGa                                                                                           <  2:6425248/150‑1 (MQ=255)
                                                           aGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGa                                                                                           <  2:819917/150‑1 (MQ=255)
                                                             ccGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATc                                                                                         <  2:4185664/150‑1 (MQ=255)
                                                             ccGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATc                                                                                         >  2:5114201/1‑150 (MQ=255)
                                                               gtgcgtgTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCa                                                                                       <  2:1605963/150‑1 (MQ=255)
                                                                   gtgtCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCa                                                                                   >  1:4190179/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)
                                                                            cccAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAgcgc                                                                          <  1:2776247/150‑1 (MQ=255)
                                                                                 ccAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGaa                                                                     <  1:6331638/150‑1 (MQ=255)
                                                                                  cAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAac                                                                    >  1:3507340/1‑150 (MQ=255)
                                                                                   aaGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAaca                                                                   >  1:265331/1‑150 (MQ=255)
                                                                                      gATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGACGCGCCAGAACACAg                                                                >  2:2303212/1‑150 (MQ=255)
                                                                                      gATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAg                                                                >  2:1948761/1‑150 (MQ=255)
                                                                                        tCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAggg                                                              <  2:1813705/150‑1 (MQ=255)
                                                                                                cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTg                                                      >  1:3141944/1‑150 (MQ=255)
                                                                                                cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTg                                                      >  1:3960766/1‑150 (MQ=255)
                                                                                                cTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTg                                                      >  1:4665390/1‑150 (MQ=255)
                                                                                                 tCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTgg                                                     >  2:1351824/1‑150 (MQ=255)
                                                                                                  ccACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTggg                                                    >  1:4507557/1‑150 (MQ=255)
                                                                                                   cACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAAcaca                                                                 >  2:674783/1‑136 (MQ=255)
                                                                                                   cACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGACGATCCAGCCATGGAAGCGCCAGAACACAGGGGTGACGTGGGa                                                   >  1:2198846/1‑150 (MQ=255)
                                                                                                     gttcgcgcCGGGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAag                                                 >  1:4826494/5‑150 (MQ=255)
                                                                                                         gcgcCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGaa                                             <  2:2064342/150‑1 (MQ=255)
                                                                                                         gcgcCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGaa                                             >  2:1829776/1‑150 (MQ=255)
                                                                                                           gcCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACg                                           <  2:2635711/150‑1 (MQ=255)
                                                                                                              ccGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGt                                        >  1:3711920/1‑150 (MQ=255)
                                                                                                               cGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGACGAGAACGGGTc                                       >  2:481930/1‑150 (MQ=255)
                                                                                                                     gAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCgag                                 >  1:3988837/1‑150 (MQ=255)
                                                                                                                      aaCCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCgaga                                >  2:2134819/1‑150 (MQ=255)
                                                                                                                        ccctGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGaaa                              <  1:6279272/147‑1 (MQ=255)
                                                                                                                        ccATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGaaa                              <  1:5213792/150‑1 (MQ=255)
                                                                                                                        ccATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGaaa                              <  1:2022020/150‑1 (MQ=255)
                                                                                                                        ccATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGaaa                              <  1:1264649/150‑1 (MQ=255)
                                                                                                                            ggCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCg                          >  2:2838360/1‑150 (MQ=255)
                                                                                                                               aCGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATAGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTc                       >  1:2306824/1‑150 (MQ=255)
                                                                                                                                cGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCt                      <  1:3362280/150‑1 (MQ=255)
                                                                                                                                  ccATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgc                    <  2:4890735/150‑1 (MQ=255)
                                                                                                                                    aTTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcgc                  >  1:239542/1‑150 (MQ=255)
                                                                                                                                      tGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTgcgcgc                <  2:3631346/150‑1 (MQ=255)
                                                                                                                                        aCCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCa              <  1:5670038/150‑1 (MQ=255)
                                                                                                                                          ctcctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATa            >  1:6748769/1‑150 (MQ=255)
                                                                                                                                          ctcctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATa            >  1:6585518/1‑150 (MQ=255)
                                                                                                                                          ctcctGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATa            >  1:326625/1‑150 (MQ=255)
                                                                                                                                          accctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATa            >  2:165461/3‑150 (MQ=255)
                                                                                                                                            cctGGCGCAGTCGTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGgtcgtc                                                                                             <  1:1197648/67‑1 (MQ=37)
                                                                                                                                             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)
                                                                                                                                                     |                                                                                                                                                    
CATGGTTTCCGGATCCATCTCCACGCTACGCCCTGGGCGCAGCCCCGCCACGTTACTGCAGCCGTGCGTGTCCGGCCCCAGCCAAGGATCGGCGACCTCCACCCAGCGCCCCGGCGCGAACCATGGCACGCCATTGACCTCCTGGCGCACGACTGCGCCAGGGTGGAATCGTTCATGGGCGCGGAACCAGTCCTCGACACGGTCGTCGATCCAGCCATGGAAGCGCCAGAACACAGGGTTGACGTGGGAAGAGAACGGGTCCCCGAGAAAGTCGTTCTGCGCGCCATACCAGCGCGCC  >  NC_002947/4763269‑4763566

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: