Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,143,605 G→T 32.0% intergenic (‑579/+19) ppc ← / ← argE phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,143,6050GT32.0% 480.7 / 263.5 290intergenic (‑579/+19)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (101/96);  new base T (53/40);  total (154/136)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 293 total aligned reads displayed.

CAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGC  >  CP009273/4143460‑4143754
                                                                                                                                                 |                                                                                                                                                     
cAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCtt                                                                                                                                                   <  1:3803497/150‑1 (MQ=39)
   tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                                >  2:287801/1‑150 (MQ=39)
     ccGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCg                                                                                                                                              >  1:1802683/1‑150 (MQ=39)
       ggATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                                            >  2:3125811/1‑150 (MQ=21)
        gATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                                           >  1:4264196/1‑150 (MQ=39)
         aTGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCt                                                                                                                                          >  2:2096194/1‑150 (MQ=255)
            cggcgCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTCTCCGGCCTACg                                                                                                                                       >  2:2058705/1‑150 (MQ=255)
                 cgAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCTCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTa                                                                                                                                  <  1:4358710/150‑1 (MQ=255)
                 cgAGCGCCTTATCCGGCCTACGGTTGGGCATCGGTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTa                                                                                                                                  >  2:1302297/1‑150 (MQ=255)
                  gAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTaa                                                                                                                                 <  2:2348203/150‑1 (MQ=255)
                  gAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTaa                                                                                                                                 >  2:3359197/1‑150 (MQ=255)
                   aGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATATGATGCGCAAGCATCGCAGCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAAt                                                                                                                                <  1:3125112/150‑1 (MQ=255)
                     cgcCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGc                                                                                                                              >  1:3377740/1‑150 (MQ=255)
                      gcCTTATCCGGCCTAGGGTTGGGCATCGTTTGCGTCACTGTCGGTCGGATAAGATGCGCATGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACCGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGcc                                                                                                                             <  1:1347712/150‑1 (MQ=14)
                      gcCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGcc                                                                                                                             >  1:2608367/1‑150 (MQ=255)
                        cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTAAAGCCAg                                                                                                                           >  2:30004/1‑150 (MQ=255)
                           aTCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCac                                                                                                                        >  2:4289221/1‑149 (MQ=255)
                            tCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCaaa                                                                                                                       >  2:3406286/1‑150 (MQ=255)
                              cGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCaaaaa                                                                                                                     >  1:3175700/1‑150 (MQ=255)
                                  cTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGt                                                                                                                 <  2:1680759/150‑1 (MQ=255)
                                   tACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTg                                                                                                                >  1:2934478/1‑150 (MQ=255)
                                    aatgtggggCATCTTTTGAGTCTCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGa                                                                                                               <  2:181775/144‑1 (MQ=255)
                                     cGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGaa                                                                                                              >  2:32369/1‑150 (MQ=255)
                                      ggTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAt                                                                                                             <  1:495626/150‑1 (MQ=255)
                                      ggTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAt                                                                                                             >  1:2516346/1‑150 (MQ=255)
                                        ttGGGCTTCGTATGTGTAGCTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTa                                                                                                           <  2:4517288/150‑1 (MQ=17)
                                        ttGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGGATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTa                                                                                                           <  1:4303988/150‑1 (MQ=255)
                                         tGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAtt                                                                                                            >  1:2365252/1‑148 (MQ=255)
                                         tGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAtt                                                                                                            <  2:2365252/148‑1 (MQ=255)
                                          gggCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTAcc                                                                                                         >  1:3992136/1‑150 (MQ=255)
                                           ggCATCGTTTGAGTCTCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCt                                                                                                        <  1:773902/150‑1 (MQ=255)
                                                  tttGAGTCACTGTCGGTCGGATAAGATGCGCACGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTAt                                                                                                 <  1:358843/150‑1 (MQ=255)
                                                  tttGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTAt                                                                                                 <  1:691678/150‑1 (MQ=255)
                                                      aGTCAGTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGt                                                                                             <  1:1644712/150‑1 (MQ=255)
                                                      aGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGt                                                                                             >  2:4332418/1‑150 (MQ=255)
                                                        tCACTGTCGGTCGGATAAGATGCGCAAGGATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTc                                                                                           <  1:4404075/150‑1 (MQ=255)
                                                        tCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTc                                                                                           <  1:3114384/150‑1 (MQ=255)
                                                        acagggTAGGTCGTTTACTATGCGCAAGCATCGATTCCGACTTTATGGCGGCACTGGAGTTTGGCAACCGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTc                                                                                           <  1:4010730/145‑1 (MQ=21)
                                                         cACTGTCGGTCGTATTAGTTGCTCAAGTATCGAATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTcg                                                                                          <  1:101999/150‑1 (MQ=18)
                                                          aCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTcgc                                                                                         >  1:1340000/1‑150 (MQ=255)
                                                          aCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTcgc                                                                                         >  1:1557544/1‑150 (MQ=255)
                                                              tAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAAt                                                                                                                                >  1:2296278/1‑107 (MQ=255)
                                                              tAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAAt                                                                                                                                <  2:2296278/107‑1 (MQ=255)
                                                                  tCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGc                                                                                 >  2:3589964/1‑150 (MQ=255)
                                                                    ggATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCtt                                                                               >  1:1302222/1‑150 (MQ=255)
                                                                     gATAAGATGCGCAAGCATCGCATCCGACATTCGTGCGGCACTGGCGTTTGGCAACCGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTg                                                                              <  1:619460/150‑1 (MQ=255)
                                                                      aTAAGATGCGCAAGCATCGCATCCGACATTAGTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGa                                                                             <  1:3852667/150‑1 (MQ=255)
                                                                        aaGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATa                                                                           >  1:4662763/1‑150 (MQ=255)
                                                                          gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATaaa                                                                         <  2:126666/150‑1 (MQ=255)
                                                                          gATGCGCAAGAATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATaaa                                                                         <  2:2150952/150‑1 (MQ=255)
                                                                           aTGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAAt                                                                                                                                >  2:93944/1‑94 (MQ=255)
                                                                           aTGCGCAAGCATCGCATCCGACATGATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAAc                                                                        <  1:2658419/150‑1 (MQ=255)
                                                                            tccgcCAGCATCGCATCCGACATTCTTGCGGCACTGGAGTGTCGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAAcc                                                                       <  1:2844291/148‑1 (MQ=255)
                                                                             gcgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCg                                                                      >  1:1123233/1‑150 (MQ=255)
                                                                              cgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt                                                                     >  1:2242563/1‑150 (MQ=255)
                                                                               gcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgtg                                                                    >  1:208857/1‑150 (MQ=255)
                                                                               gcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgtg                                                                    <  2:4264196/150‑1 (MQ=255)
                                                                               gcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGGTCGCGGGTGGGCTTGATAAACCgtg                                                                    >  2:4110075/1‑150 (MQ=255)
                                                                                cAAGCATCGCATCCGACATTATTGCGGCACTGGAGGTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAAGGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgtgt                                                                   <  1:982015/150‑1 (MQ=255)
                                                                                 aaGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGtt                                                                  >  2:3254350/1‑150 (MQ=255)
                                                                                     atcgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCa                                                              >  1:1462534/1‑150 (MQ=255)
                                                                                        gcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGat                                                           <  1:548950/150‑1 (MQ=255)
                                                                                          atcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                                                                 <  1:1721261/62‑1 (MQ=35)
                                                                                           tcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATAtt                                                        >  2:3080647/1‑150 (MQ=255)
                                                                                           tcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATAtt                                                        <  1:4752093/150‑1 (MQ=255)
                                                                                           tcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATAtt                                                        <  1:3006567/150‑1 (MQ=255)
                                                                                            ccGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtc                                                       >  2:3523902/1‑150 (MQ=255)
                                                                                               aCATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtcatc                                                    >  2:738955/1‑150 (MQ=255)
                                                                                                cATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtcatca                                                   <  1:2370599/150‑1 (MQ=255)
                                                                                                 attattGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAg                                                  >  2:241269/1‑150 (MQ=255)
                                                                                                   tattGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGt                                                >  2:569437/1‑150 (MQ=255)
                                                                                                   tattGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGt                                                >  2:2959171/1‑150 (MQ=255)
                                                                                                    attGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGtt                                               >  1:1084959/1‑150 (MQ=255)
                                                                                                     ttGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTtg                                              <  1:3260624/150‑1 (MQ=255)
                                                                                                       gCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTtgat                                            >  1:2666279/1‑150 (MQ=255)
                                                                                                        cGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTtgatg                                           >  1:4479307/1‑150 (MQ=255)
                                                                                                         ggCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTtgatga                                          <  2:1618679/150‑1 (MQ=255)
                                                                                                          gCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAg                                         <  1:2796700/150‑1 (MQ=255)
                                                                                                            aCTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGcc                                       >  2:3306077/1‑150 (MQ=255)
                                                                                                            aCTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGcc                                       >  2:1067977/1‑150 (MQ=255)
                                                                                                                    ttGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAAt                               <  1:3870369/150‑1 (MQ=255)
                                                                                                                     tGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAAtt                              >  1:1820799/1‑150 (MQ=255)
                                                                                                                        cAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAg                           <  2:3197440/150‑1 (MQ=255)
                                                                                                                         aaCAGCGCCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGc                          <  2:3992136/150‑1 (MQ=255)
                                                                                                                            aGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAg                       <  2:4639756/150‑1 (MQ=255)
                                                                                                                            aGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAg                       >  2:3814728/1‑150 (MQ=255)
                                                                                                                             gcgcCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAgg                      <  1:1117026/150‑1 (MQ=255)
                                                                                                                             gcgcCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAgg                      <  1:2052887/150‑1 (MQ=255)
                                                                                                                              cgcCGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTAGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGc                     >  1:4540686/1‑150 (MQ=255)
                                                                                                                               gcCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGcc                    <  2:1568956/150‑1 (MQ=255)
                                                                                                                                 cGGATGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGcccc                  >  2:2477763/1‑150 (MQ=255)
                                                                                                                                   gATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCaa                >  2:2623988/1‑150 (MQ=255)
                                                                                                                                    aTGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAac               >  1:3000083/1‑150 (MQ=255)
                                                                                                                                    aTGCGGCGCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAac               >  1:4161974/1‑150 (MQ=255)
                                                                                                                                     tGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAaca              <  2:1090573/150‑1 (MQ=255)
                                                                                                                                        ggcgCGAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCa           <  1:1982937/150‑1 (MQ=255)
                                                                                                                                          cgcgAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGc         <  1:1174231/150‑1 (MQ=255)
                                                                                                                                            cgAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGt       <  2:1137497/150‑1 (MQ=255)
                                                                                                                                             gAGCTCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTc      <  2:1463671/150‑1 (MQ=255)
                                                                                                                                               gcgcCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCgg    <  1:1940643/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGC  >  CP009273/4143460‑4143754

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