Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,497,475 C→A 5.0% 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 5.0% 1132.2 / 19.5 297E164* (GAG→TAG) fryAputative PTS enzyme, Hpr component/enzyme I component/enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base C (165/117);  new base A (8/7);  total (173/124)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 297 total aligned reads displayed.

ATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCC  >  CP009273/2497327‑2497623
                                                                                                                                                    |                                                                                                                                                    
  gatCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTAGGc                                                                                                                                                   >  1:984950/1‑150 (MQ=255)
     cgccgcCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAg                                                                                                                                                >  1:1379221/1‑150 (MQ=255)
     cgccgcCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAg                                                                                                                                                <  2:3532488/150‑1 (MQ=255)
      gccgccCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGc                                                                                                                                               <  1:4876233/150‑1 (MQ=255)
      gccgccCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGc                                                                                                                                               >  2:2634737/1‑150 (MQ=255)
          cccccAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTg                                                                                                                                           <  2:2569066/150‑1 (MQ=255)
            cccAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGcc                                                                                                                                         >  2:3748780/1‑150 (MQ=255)
              cAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAg                                                                                                                                       <  1:3027888/150‑1 (MQ=255)
              cAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAg                                                                                                                                       >  2:3467615/1‑150 (MQ=255)
                gCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGGTCAGCTGCTCGGCAAGCGTTGCCAGGc                                                                                                                                     >  2:574604/1‑150 (MQ=255)
                  cccTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCtg                                                                                                                                   >  1:2885563/1‑150 (MQ=255)
                    cTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCtgtg                                                                                                                                 >  1:4103687/1‑150 (MQ=255)
                        tGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTcc                                                                                                                             <  1:5054667/150‑1 (MQ=255)
                           tGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGc                                                                                                                          >  1:4027165/1‑150 (MQ=255)
                            gTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCACTTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGcc                                                                                                                         >  2:671285/1‑150 (MQ=255)
                             ttCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCg                                                                                                                        <  2:4472611/150‑1 (MQ=255)
                             ttCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCg                                                                                                                        >  2:5251088/1‑150 (MQ=255)
                                tGTCATCNGNCNACGGATAGTGCCTGCAAATTCATCATCCTGACTCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCCCGCAATTGCGGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGNt                                                                                                                     <  1:1136501/150‑1 (MQ=255)
                                 gTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTg                                                                                                                    <  2:4059570/150‑1 (MQ=255)
                                    atcaGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGaa                                                                                                                 >  1:4291425/1‑150 (MQ=255)
                                     tcaGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAAt                                                                                                                >  1:3238409/1‑150 (MQ=255)
                                     tcaGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAAt                                                                                                                >  2:3422705/1‑150 (MQ=255)
                                          cGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTgc                                                                                                           <  1:3165315/150‑1 (MQ=255)
                                           gACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTgcg                                                                                                          >  1:274590/1‑150 (MQ=255)
                                             cGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCa                                                                                                        >  1:1690574/1‑150 (MQ=255)
                                              ggATATTGCCTGCAAATTCATCATCCTGAATCAGCGACACATGGGCGCTGTGGATAGTTTTGCTTTCGCCGTCACGCTCACGCACTTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCAAGCCGGGTGGAATCTTGCGCAc                                                                                                       >  1:4050216/1‑150 (MQ=255)
                                               gATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACt                                                                                                      >  2:5089152/1‑150 (MQ=255)
                                               gATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACt                                                                                                      <  2:414664/150‑1 (MQ=255)
                                               gATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTAGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACt                                                                                                      <  2:2342117/150‑1 (MQ=255)
                                               gATATTGCCTGCAAATTCATCATCATGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACt                                                                                                      >  2:1667295/1‑150 (MQ=255)
                                               gATAGTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCGGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACt                                                                                                      <  2:4819991/150‑1 (MQ=255)
                                                atatTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTg                                                                                                     <  1:1860969/150‑1 (MQ=255)
                                                atatTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTg                                                                                                     >  1:2047681/1‑150 (MQ=255)
                                                 tatTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTgg                                                                                                    >  1:1327948/1‑150 (MQ=255)
                                                  atTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGc                                                                                                   >  2:1312477/1‑150 (MQ=255)
                                                    tGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCgg                                                                                                 <  2:3370482/150‑1 (MQ=255)
                                                      ccTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCgggg                                                                                               <  1:519912/150‑1 (MQ=255)
                                                       cTGCAAATTCATCATCCTGAATCGGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGa                                                                                              >  2:1539971/1‑150 (MQ=255)
                                                        tGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGAt                                                                                             <  1:2557266/150‑1 (MQ=255)
                                                         gCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGAtt                                                                                            >  2:5321063/1‑150 (MQ=255)
                                                         gCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGAtt                                                                                            <  2:110150/150‑1 (MQ=255)
                                                          cAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTg                                                                                           >  2:3153372/1‑150 (MQ=255)
                                                           aaaTTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGc                                                                                          >  2:4880649/1‑150 (MQ=255)
                                                             aTTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGACTCTTGCGCACTGGCGGGGATTGccc                                                                                        >  2:2912963/1‑150 (MQ=255)
                                                               tcatcaTCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGa                                                                                      >  2:1710420/1‑150 (MQ=255)
                                                                catcatCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGAt                                                                                     >  2:1181646/1‑150 (MQ=255)
                                                                catcatCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGAt                                                                                     >  2:1930228/1‑150 (MQ=255)
                                                                catcatCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGAt                                                                                     >  2:2763959/1‑150 (MQ=255)
                                                                 atcatcCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATa                                                                                    >  1:84197/1‑150 (MQ=255)
                                                                 atcatcCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAAGTGCTGGTTCAGTTGCTAGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATa                                                                                    >  2:3610262/1‑150 (MQ=255)
                                                                    atcCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACt                                                                                 >  2:1813052/1‑150 (MQ=255)
                                                                    atcCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACt                                                                                 >  2:755938/1‑150 (MQ=255)
                                                                       cTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTTCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTc                                                                              <  1:4685008/150‑1 (MQ=255)
                                                                       cTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTc                                                                              >  1:892909/1‑150 (MQ=255)
                                                                       cTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTc                                                                              <  1:4414440/150‑1 (MQ=255)
                                                                         gAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGa                                                                            >  1:2955706/1‑150 (MQ=255)
                                                                           aTCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTAGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAgg                                                                          <  2:4694543/150‑1 (MQ=255)
                                                                               gCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGt                                                                      <  1:2641340/150‑1 (MQ=255)
                                                                                cGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTc                                                                     <  2:2704536/150‑1 (MQ=255)
                                                                                 gACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCg                                                                    >  2:2095741/1‑150 (MQ=255)
                                                                                    aaaTGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCtt                                                                 >  1:2172578/1‑150 (MQ=255)
                                                                                           cgcTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACa                                                          >  2:4094472/1‑150 (MQ=255)
                                                                                           cgcTGAGGAGAGTTTTGCGTTCGCCGTCACGCTCACGCTCTTGCTGGTTCAGTTGCTCGGCTAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACa                                                          <  1:4969300/150‑1 (MQ=255)
                                                                                            gcTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAg                                                         >  1:5329153/1‑150 (MQ=255)
                                                                                                aGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTc                                                     >  2:960398/1‑150 (MQ=255)
                                                                                                 ggATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCa                                                    <  2:4347406/150‑1 (MQ=255)
                                                                                                  gATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAg                                                   >  2:3229131/1‑150 (MQ=255)
                                                                                                  gATAGTTTTGCTTTCGCCGTCACGCTCACGCAAGTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAg                                                   >  2:2232296/1‑150 (MQ=255)
                                                                                                     aGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCAGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGt                                                >  2:3932281/1‑150 (MQ=255)
                                                                                                        tttGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTAcc                                             >  2:3288216/1‑150 (MQ=255)
                                                                                                        tttGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTAcc                                             >  2:1633703/1‑150 (MQ=255)
                                                                                                         ttGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTAccc                                            >  2:1415425/1‑150 (MQ=255)
                                                                                                         ttGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTAccc                                            >  1:4234938/1‑150 (MQ=255)
                                                                                                          tGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCa                                           <  2:3735716/150‑1 (MQ=255)
                                                                                                           gCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCac                                          <  1:3467615/150‑1 (MQ=255)
                                                                                                            cTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacg                                         >  1:3369138/1‑150 (MQ=255)
                                                                                                                 gCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgccgac                                    <  1:3614206/150‑1 (MQ=255)
                                                                                                                 gCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgccgac                                    <  1:2610864/150‑1 (MQ=255)
                                                                                                                  ccGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgccgacg                                   <  2:4971359/150‑1 (MQ=255)
                                                                                                                   cGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgccgacgc                                  <  1:3748780/150‑1 (MQ=255)
                                                                                                                    gTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCacgccgacgcc                                 >  1:2916415/1‑150 (MQ=255)
                                                                                                                       acgctcacgcAATCGCTGGTTCAGTGGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCTCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCt                              <  1:164657/150‑1 (MQ=255)
                                                                                                                           tcacgcAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGcc                          >  2:5180113/1‑150 (MQ=255)
                                                                                                                            cacgcAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCa                         <  1:1964307/150‑1 (MQ=255)
                                                                                                                              cgcAATTGCTGGTTCAGTTGCTCGGCAAGCGTAGCCAGGCTGTGCTCCAGCCGCGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGACGAGGCTGACGCTTTGTAACAGGGTCAGCGAACCCACGCCTACGCCGCTAGCCAGc                       >  2:4498565/1‑150 (MQ=255)
                                                                                                                                cAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGcac                     >  1:2382004/1‑150 (MQ=255)
                                                                                                                                 aaTTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGcaca                    <  2:4285927/150‑1 (MQ=255)
                                                                                                                                  aTTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACAt                   >  2:4955387/1‑150 (MQ=255)
                                                                                                                                    tGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTg                 <  2:869223/150‑1 (MQ=255)
                                                                                                                                      cTGGTTCAGTTGCTCGGCTAGCGGTGCCAGGCTGTGCTCCAGCCGGGTGGACTCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGGACCCACGCCGACGCCGCTTGCCAGCACATTGcc               <  1:103252/150‑1 (MQ=255)
                                                                                                                                       tGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCg              <  2:4583979/150‑1 (MQ=255)
                                                                                                                                       tGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCg              >  1:3968797/1‑150 (MQ=255)
                                                                                                                                       tGGTTCAGTTGCTAGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCg              <  2:3140755/150‑1 (MQ=255)
                                                                                                                                        ggTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGt             >  2:3033217/1‑150 (MQ=255)
                                                                                                                                             aGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGt        >  1:1721169/1‑150 (MQ=255)
                                                                                                                                             aGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGt        <  1:786196/150‑1 (MQ=255)
                                                                                                                                              gTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTa       >  2:2066195/1‑150 (MQ=255)
                                                                                                                                              gTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTa       >  1:1451998/1‑150 (MQ=255)
                                                                                                                                                 gCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAAt    <  1:4777685/150‑1 (MQ=255)
                                                                                                                                                  cTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATc   <  1:5250879/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTTCAGTTGCTCGGCAAGCGTTGCCAGGCTGTGCTCCAGCCGGGTGGAATCTTGCGCACTGGCGGGGATTGCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCC  >  CP009273/2497327‑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: