Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,497,475 C→A 9.1% E164* (GAG→TAG)  fryA ← putative PTS enzyme, Hpr component/enzyme I component/enzyme IIA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,497,4750CA 9.1% 940.0 / 44.3 273E164* (GAG→TAG) fryAputative PTS enzyme, Hpr component/enzyme I component/enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base C (134/114);  new base A (11/14);  total (145/128)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

GATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCC  >  CP009273/2497326‑2497623
                                                                                                                                                     |                                                                                                                                                    
gatgatCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTc                                                                                                                                                      >  2:3850526/1‑150 (MQ=255)
   gatCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGc                                                                                                                                                   >  1:3606932/1‑150 (MQ=255)
     tCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCaa                                                                                                                                                 <  1:1108640/150‑1 (MQ=255)
       gccgccCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGc                                                                                                                                               <  2:891034/150‑1 (MQ=255)
          gccCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGtt                                                                                                                                            <  1:3490151/150‑1 (MQ=255)
          gccCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGtt                                                                                                                                            <  1:4157281/150‑1 (MQ=255)
          gccCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGtt                                                                                                                                            <  1:909842/150‑1 (MQ=255)
               cAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAg                                                                                                                                       >  1:2767888/1‑150 (MQ=255)
               cAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAg                                                                                                                                       <  2:308714/150‑1 (MQ=255)
                 gCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGc                                                                                                                                     >  1:748108/1‑150 (MQ=255)
                   cncTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCtg                                                                                                                                   >  1:626575/1‑150 (MQ=255)
                   cccTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCtg                                                                                                                                   >  1:4201589/1‑150 (MQ=255)
                    ccTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCtgt                                                                                                                                  <  1:455982/150‑1 (MQ=255)
                    ccTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCtgt                                                                                                                                  >  1:2277469/1‑150 (MQ=255)
                         tGCTGTTCTGTCATCAGGCGACGGAGATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCACTTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTcc                                                                                                                             >  2:4990186/1‑150 (MQ=255)
                          gCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCa                                                                                                                            <  1:1412763/150‑1 (MQ=255)
                              ttCTGTCATCAGGCGACGGATATAGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCg                                                                                                                        >  1:1591479/1‑150 (MQ=255)
                               tCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCgg                                                                                                                       <  1:4263439/150‑1 (MQ=255)
                               tCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCgg                                                                                                                       <  1:1685758/150‑1 (MQ=255)
                                 tGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGt                                                                                                                     >  2:4360219/1‑150 (MQ=255)
                                  gTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTg                                                                                                                    >  2:3584957/1‑150 (MQ=255)
                                   tcatcaGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTgg                                                                                                                   >  2:3981554/1‑150 (MQ=255)
                                       caGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATc                                                                                                               >  2:1052433/1‑150 (MQ=255)
                                       caGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATc                                                                                                               >  1:184218/1‑150 (MQ=255)
                                           cGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTgc                                                                                                           <  1:782133/150‑1 (MQ=255)
                                              cGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCa                                                                                                        <  2:3141070/150‑1 (MQ=255)
                                              cGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCa                                                                                                        >  1:1044951/1‑150 (MQ=255)
                                               ggATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCAc                                                                                                       >  1:3125467/1‑150 (MQ=255)
                                                gATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACt                                                                                                      <  2:4659679/150‑1 (MQ=255)
                                                gATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACt                                                                                                      >  2:3235321/1‑150 (MQ=255)
                                                  tatTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTgg                                                                                                    <  2:263996/150‑1 (MQ=255)
                                                   atggCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCCTGCGCACTGGc                                                                                                   <  2:4391998/147‑1 (MQ=255)
                                                     tGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCgg                                                                                                 >  1:3672866/1‑150 (MQ=255)
                                                     tGCCTGCAAATTCATCATCCTGAATCAGCGAAAAATGTGCGCTGAGGATAGTTTTGCTTTCGCCGTCAAGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGTCTGTGCTCCAGACGGGTGGAATATTGCTCACTCGCgg                                                                                                 >  1:2523126/1‑150 (MQ=255)
                                                        cTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGa                                                                                              <  2:1566265/150‑1 (MQ=255)
                                                        cTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGa                                                                                              >  1:2242826/1‑150 (MQ=255)
                                                            aaaTTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGc                                                                                          >  1:4316108/1‑150 (MQ=255)
                                                              aTTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCg                                                                                       <  1:2212544/150‑1 (MQ=255)
                                                               ttCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTAGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCg                                                                                       <  1:3885267/150‑1 (MQ=255)
                                                                tcatcaTCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGa                                                                                      <  1:3561054/150‑1 (MQ=255)
                                                                     atcCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACt                                                                                 >  1:1299636/1‑150 (MQ=255)
                                                                     atcCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACt                                                                                 >  1:4838218/1‑150 (MQ=255)
                                                                        cTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTc                                                                              >  1:2650787/1‑150 (MQ=255)
                                                                         tGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCg                                                                             >  2:3503357/1‑150 (MQ=255)
                                                                         tGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGGTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCg                                                                             <  2:4913195/150‑1 (MQ=255)
                                                                          gAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGa                                                                            >  1:2693583/1‑150 (MQ=255)
                                                                           aaTCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAg                                                                           >  2:1844907/1‑150 (MQ=255)
                                                                             tCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGc                                                                         >  1:5134821/1‑150 (MQ=255)
                                                                              cAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCt                                                                        <  2:146444/150‑1 (MQ=255)
                                                                                 cGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTc                                                                     >  2:3429387/1‑150 (MQ=255)
                                                                                  gACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCg                                                                    >  2:4544809/1‑150 (MQ=255)
                                                                                   aCAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGc                                                                   >  2:75952/1‑150 (MQ=255)
                                                                                    cAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCt                                                                  <  2:4233213/150‑1 (MQ=255)
                                                                                     aaaTGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGTGTGGAATCTTGCGCACTGGCGGGGATTGCCCTATAACTGTCGAGGCTGTCGCtt                                                                 >  2:434076/1‑150 (MQ=255)
                                                                                     aaaTGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCtt                                                                 <  1:2838114/150‑1 (MQ=255)
                                                                                      aaTGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCttt                                                                >  1:5139494/1‑150 (MQ=255)
                                                                                      aaTGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCttt                                                                >  1:1246030/1‑150 (MQ=255)
                                                                                       aTGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAGTCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTg                                                               <  2:1966462/150‑1 (MQ=255)
                                                                                       aTGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTg                                                               >  2:2801933/1‑150 (MQ=255)
                                                                                       aTGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTg                                                               >  2:2248950/1‑150 (MQ=255)
                                                                                       aTGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACAGTCGAGGCTGTCGCTTTg                                                               >  1:2444775/1‑150 (MQ=255)
                                                                                         gggCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTAGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTa                                                             <  2:205840/150‑1 (MQ=255)
                                                                                           gcgcTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAAc                                                           >  2:934641/1‑150 (MQ=255)
                                                                                            cgcTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACa                                                          >  2:705189/1‑150 (MQ=255)
                                                                                             gcTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTAGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAg                                                         >  1:3329389/1‑150 (MQ=255)
                                                                                                  ggATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCa                                                    <  1:2964176/150‑1 (MQ=255)
                                                                                                  ggATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCa                                                    <  2:4904863/150‑1 (MQ=255)
                                                                                                     tAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCg                                                 <  1:832649/150‑1 (MQ=255)
                                                                                                      aGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGt                                                >  2:1695932/1‑150 (MQ=255)
                                                                                                       gTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTa                                               <  1:1203854/150‑1 (MQ=255)
                                                                                                          ttGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGGTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTAccc                                            <  2:2691598/150‑1 (MQ=255)
                                                                                                            gCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCac                                          <  2:160751/150‑1 (MQ=255)
                                                                                                              tttCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgc                                        >  2:4053574/1‑150 (MQ=255)
                                                                                                               ttCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgcc                                       <  2:4301922/150‑1 (MQ=255)
                                                                                                                 cGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgccga                                     >  1:1022568/1‑150 (MQ=255)
                                                                                                                  gCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGTGTACCCacgccgac                                    >  1:4350858/1‑150 (MQ=255)
                                                                                                                  gCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgccgac                                    <  1:3808298/150‑1 (MQ=255)
                                                                                                                    cGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgccgacgc                                  >  1:2208158/1‑150 (MQ=255)
                                                                                                                     gTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgccgacgcc                                 <  2:3615593/150‑1 (MQ=255)
                                                                                                                     gTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgccgacgcc                                 >  1:4998423/1‑150 (MQ=255)
                                                                                                                      tcacgctcacgcAATTGCTGGTTCAGTTGCTAGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgccgacgccg                                >  1:4640519/1‑150 (MQ=255)
                                                                                                                           ctcacgcAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGc                           >  2:1296349/1‑150 (MQ=255)
                                                                                                                            tcacgcAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGcc                          >  1:1783767/1‑150 (MQ=255)
                                                                                                                            tcacgcAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGcc                          <  2:327244/150‑1 (MQ=255)
                                                                                                                            tcacgcAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGcc                          >  2:4086553/1‑150 (MQ=255)
                                                                                                                            tcacgcAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGcc                          >  1:1784134/1‑150 (MQ=255)
                                                                                                                            tcacgcAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGcc                          <  1:2088438/150‑1 (MQ=255)
                                                                                                                            tcacgcAATTGCTGGTTCAGTTGCTAGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGcc                          >  1:3205596/1‑150 (MQ=255)
                                                                                                                             cacgcAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCa                         >  1:323196/1‑150 (MQ=255)
                                                                                                                             caagcAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCa                         >  1:215646/4‑150 (MQ=255)
                                                                                                                                gcAATTGCTGGGTCAGTTGCTCGGCTAGCGTTGCCCGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATACCTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGca                      <  2:1119556/150‑1 (MQ=255)
                                                                                                                                 cAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGcac                     >  2:4852591/1‑150 (MQ=255)
                                                                                                                                 cAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGcac                     >  1:1559181/1‑150 (MQ=255)
                                                                                                                                  aaTTGCTGGTTCAGTTGCTAGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGcaca                    <  1:3576303/150‑1 (MQ=255)
                                                                                                                                   aTTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACAt                   <  2:2497788/150‑1 (MQ=255)
                                                                                                                                    gtGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACCtt                  <  2:3555940/149‑1 (MQ=255)
                                                                                                                                     tGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTg                 <  1:4382935/150‑1 (MQ=255)
                                                                                                                                        tGGTTCAGTTGCTCGGCAAGCGTTGCCAGGATGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCACGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCg              >  1:3416974/1‑150 (MQ=255)
                                                                                                                                              aGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACACCGACGCCGCTTGCCAGCACATTGCCGTACAGt        >  2:596069/1‑150 (MQ=255)
                                                                                                                                                ttGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTaa      >  2:2001133/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCC  >  CP009273/2497326‑2497623

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: