Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,713 C→A 27.9% intergenic (‑30/+278) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,170,7130CA27.9% 284.0 / 34.5 194intergenic (‑30/+278)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base C (49/75);  new base A (45/16);  total (103/91)
Fisher's exact test for biased strand distribution p-value = 1.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

GTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATG  >  CP009273/2170564‑2170859
                                                                                                                                                     |                                                                                                                                                   
gTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTc                                                                                                                                                     <  2:3469104/150‑1 (MQ=255)
gTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTc                                                                                                                                                     <  1:3022080/150‑1 (MQ=255)
 ttGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCa                                                                                                                                                    <  1:4640559/150‑1 (MQ=255)
 ttGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCa                                                                                                                                                    <  2:3095657/150‑1 (MQ=255)
 ttGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCa                                                                                                                                                    <  2:3095680/150‑1 (MQ=255)
  tGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTaaa                                                                                                                                                   >  2:649982/1‑150 (MQ=255)
      aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                               <  1:1497834/150‑1 (MQ=255)
      aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGTTTTAaaaaaa                                                                                                                                               >  1:1483802/1‑148 (MQ=255)
         ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAAAATaa                                                                                                                                            >  2:1327290/1‑150 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                           <  1:773376/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                           <  2:4013450/150‑1 (MQ=255)
              cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                       <  2:5075758/150‑1 (MQ=255)
              cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGTTTTAAAAATATaaataaa                                                                                                                                       >  1:2629973/1‑145 (MQ=255)
               tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATGaaag                                                                                                                                      >  2:3525824/1‑149 (MQ=255)
                   accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                                  <  1:3915799/150‑1 (MQ=255)
                   accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTAAaatatt                                                                                                                                  >  1:4206090/1‑150 (MQ=255)
                     cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                                <  2:908908/150‑1 (MQ=255)
                        ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAAATATtttttt                                                                                                                             >  2:1884503/1‑150 (MQ=255)
                          gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttttt                                                                                                                           <  1:3756246/150‑1 (MQ=255)
                            aTCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTaa                                                                                                                         >  2:1684028/1‑150 (MQ=255)
                              cGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTCAAAATATAATTAAAATATTTTTTTTTAAat                                                                                                                       >  2:1914318/1‑150 (MQ=255)
                                tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                     <  1:3637281/150‑1 (MQ=255)
                                 ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAAAATAAATAAAATATTTTTTTTTAAatata                                                                                                                    >  2:4059725/1‑149 (MQ=255)
                                   cGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATAga                                                                                                                  >  2:3704782/1‑150 (MQ=255)
                                   cGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATAga                                                                                                                  >  2:3704731/1‑150 (MQ=255)
                                     agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                                <  2:1131649/150‑1 (MQ=255)
                                     agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTTTTTTTAAATATCgaga                                                                                                                >  2:2530264/1‑150 (MQ=255)
                                     agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATAgaga                                                                                                                >  2:1501398/1‑150 (MQ=255)
                                     agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTTAATTATTTTTTTTTTAATATctata                                                                                                                >  2:938581/1‑146 (MQ=255)
                                       aTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATa                                                                                                              <  2:4837478/150‑1 (MQ=255)
                                       aTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAAAATAATTAAATTATTTTTTTTTTAAATTAGAGATa                                                                                                              >  2:1610438/1‑150 (MQ=255)
                                       aTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTAAATTATTTTTTTTTTTATATAGatatt                                                                                                              >  2:4205422/1‑149 (MQ=255)
                                           tgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGa                                                                                                         >  2:2843075/1‑150 (MQ=255)
                                           tgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTGAAAATATAATGAAATTATTTGTTTTTATATATCGAGGTAAcca                                                                                                          >  2:122587/1‑148 (MQ=255)
                                                     aTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAc                                                                                               <  2:2190953/150‑1 (MQ=255)
                                                          ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTGAAAATATAATGAAATTATTTGTTTTTTAATATTGCGGTTACGATACAAAAAAAGACaa                                                                                          >  1:3343937/1‑150 (MQ=255)
                                                          ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAAAATAATTAAATTATTTGTTTTTTAATTTAGAGATAACGATCACAAAAACggaaa                                                                                          >  1:2517088/1‑146 (MQ=255)
                                                          ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTGAAAATATAAATAAAATATTTTTTTTTTAATATAGAGAGAAAGATAAAAAAAAACACaa                                                                                          >  1:683401/1‑150 (MQ=255)
                                                           gAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTGGATTTAAAAAAATAATGAAATTATTTTTTTTTAAATATAGAGATAAAGATAAAAAAAAcaaaaaa                                                                                         >  1:791947/1‑144 (MQ=255)
                                                             aCCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATAACAAAAACGACAatat                                                                                       >  2:417903/1‑150 (MQ=255)
                                                              ccGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGAAAAAATg                                                                                      >  2:3746844/1‑150 (MQ=255)
                                                               cGCATAACCGCCGCGCTGTGCTTTGTTCAGCATCTGATTTGTCGATACCACTTAAATTTTCATATCCTGTCTTTTGTTTTCTATTTAAAAATATAATTAAATTTTTTTTTTTTAAATTTCGAGATAACTAACACAAAAACGACAATATGa                                                                                     <  1:304200/150‑1 (MQ=255)
                                                               cGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATACAAAGAAATTATTTTTTTTTTAATATAGAGATAACCATCACAAAAACGACAaaatta                                                                                     >  2:548011/1‑148 (MQ=255)
                                                               cGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAATTAAATTATTTTTTTTTAAATATAGAGATAAAGATCAAAAAAACCACAaaatta                                                                                     >  2:4993685/1‑148 (MQ=255)
                                                                gaaTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTACAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGaa                                                                                    <  2:4629619/148‑1 (MQ=255)
                                                                gCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATAGAGATAACGATCAAAAAAACGACAATATGaa                                                                                    >  2:3389506/1‑150 (MQ=255)
                                                                  aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGTATTTAAAAATATAATTAAATTATTTTTTTTTTAATTTAGAGATAACGATAACAAAaaagaaaaaaagaaaa                                                                                  >  1:4913081/1‑150 (MQ=255)
                                                                   tAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAATATAATTAAATTAATTTTTTTTAAAAATAGAGATAAACAAAAAAAAAAAGACAATATGtaaaa                                                                                 >  1:2552838/1‑149 (MQ=255)
                                                                   tAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTAAATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAat                                                                                 <  2:660630/150‑1 (MQ=255)
                                                                    acccgccgCGCTGTGCGTGGTTCAGCATCTGCTTTGTCTATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAatt                                                                                <  2:4547788/148‑1 (MQ=255)
                                                                    aaCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAATATAATTAAATTATTTTTTTTTAAATATAGAGATTAAGATCACAAAAACCAACATATGaaaaat                                                                                >  2:4068789/1‑148 (MQ=255)
                                                                    aaCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAatt                                                                                <  2:2577624/150‑1 (MQ=255)
                                                                    aaCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTTGTTTTAAAAAAAAAAATTAAATATTTTTTTTTTAATATAGAGAAAAAGAAACAAAAAAACACAATAAGAAAatt                                                                                >  1:483612/1‑150 (MQ=255)
                                                                    aaCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAATATAAAAATAAATTATTTTTTTTTTAATATAGAGAGTGAGATAAAAAAAAAAACAAAATGAAAata                                                                                >  1:623399/1‑149 (MQ=255)
                                                                       cgccgcGCTGTGGGTTGTTCATCATCTGCTTTGTCTATACCACGTACATTTTCATAGCCGGTCGTTTGTTTTCGATTTAAAAATATAATGAATTTTTTTGTTTTTAAATTTCGATAAACCGATCAAAAAAACGACAATATGAAAattatt                                                                             <  2:29411/150‑1 (MQ=255)
                                                                        gccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTAATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTc                                                                            <  2:2430455/150‑1 (MQ=255)
                                                                         ccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTAAAAAAAAAATTAAAATATTTTTTTTTAAAAATAGAGATAACCATCAAAAAAAAGACAAAAAGAAAaatattag                                                                           >  1:4997265/1‑148 (MQ=255)
                                                                           gcgcTTGGCTTTGTTCAGCTTCTGCTTTGTCGATACCACGTACTTTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATATTGAAATTATTTGTTTTTAAATATCGTGATAACGATCACAAAAACGACAATATGAAAATTATTCgag                                                                         <  2:5005400/150‑1 (MQ=255)
                                                                                 tgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATAGAGATAACGATCACAAAAACCACAATATGAAAATTATTAGAGGAGTGa                                                                   >  2:100118/1‑150 (MQ=255)
                                                                                  gCGTTTTTCACCATCTGCTTTGTCGATACCACTTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGaa                                                                  <  2:1483802/150‑1 (MQ=255)
                                                                                  gCGTTGTTCAGCATCTGCTTTGTCTATCCACCGTACATTTTCATAGCCTGTCTTTTGTTTTCGATTACAAAATATAATGAAATTATTTTTTTTTAAATATCGATATAACTAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGaa                                                                  <  1:4100747/150‑1 (MQ=255)
                                                                                   ctttttacAGAATCTGTTTTGTCTATCCCCCTTAATTTTTAATACCCTGTCTTTTTTTTTCTTTTACAAAATATAATTAAATTATTTTTTTTTAAATATCTATATAACTATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGaaa                                                                 <  1:932968/143‑1 (MQ=255)
                                                                                   cGTTGTTCAGCATCTGCTGTGTCTATACCCCTTAAATTTTCATATCCGGTCGTTTGTTTCCGATTTCAAAATATAATGAAATTTTTTGTTTTTAAATATCGAGAAAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGaaa                                                                 <  2:1796984/150‑1 (MQ=255)
                                                                                     ttgttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTAAAAATATAATGAAATTATTTGTTTTTTAATATAGAGATTAAGATAACAAAAACGACAATATTAAAATTATTAGATGAGTGAAAgg                                                               >  2:4430569/1‑150 (MQ=255)
                                                                                     ttgttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAAATATTTGTTTTTAAATATAGAGATAACGATCACAAAAACGACAATATGAAAATTATTAGAGGAGTGAAAgg                                                               >  2:2995201/1‑150 (MQ=255)
                                                                                     ttgttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAAAGAAATTAATTTTTTTTAAAAATCAAGATAAAGAACACAAAAACCACAAAATTAAAAATATTAGAGGAGTGAAAgg                                                               >  2:257112/1‑150 (MQ=255)
                                                                                          catcATGTGCTTTTTCGATACCACGTAATTTTTCATATCCTGTCTTTTGTTTTCTATTTCAAAATATAATTAATTTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaaaa                                                          <  1:4082257/147‑1 (MQ=255)
                                                                                            gCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTTGATTTAAAAAAATAAATAAAATATTTTTTTTTTAATATCGATATAAACATAAAAAAAACCACAAAAAGAAAATAATTAGATGTGTGAAATGCaaaaaa                                                        >  1:3130504/1‑150 (MQ=255)
                                                                                                    tttGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTAAAAATATAAAGAAATTATTTTTTTTTTAAAATAGAGTTAAAGATCACAAAAAAGAAAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAAGGCCTc                                                >  1:4754450/1‑150 (MQ=255)
                                                                                                     ttGTCGCTACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTcc                                               <  2:167417/150‑1 (MQ=255)
                                                                                                      tGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCTATTTAAAAATATAATGAAATTATTTTTTTTAAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTccc                                              <  2:4730242/150‑1 (MQ=255)
                                                                                                      ggTCGATCCCACGTACATTTTCATATCCTGTCTTTTGTTTTCGATTACAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTccc                                              <  2:4450735/149‑1 (MQ=255)
                                                                                                       gTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCg                                             <  2:4336193/150‑1 (MQ=255)
                                                                                                        tCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTATGTTTGAAAAAAAAATTAAAATATTTTTTTTTTAAAATAGAGATTAAGAAAACAAAAAAAACAAAATTAAAAATATTAGGGGAGTGAAAGGCAAAAAAAAAGCCAcccca                                            >  1:1715589/1‑148 (MQ=255)
                                                                                                          gATACCCCGTACATTTACATATCCTGTCGTTTGTTTTCTATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                          <  2:2411473/150‑1 (MQ=255)
                                                                                                          gATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTGAAAATATAATGAAAATATTTTTTTTTATATTTAGAGTTAAAGATCACAAAAAAGACCAAATTTAAATTATTAGAGGAGTGTAAGGGAAAATAAAGGCCAACCGaga                                          >  2:2455146/1‑148 (MQ=255)
                                                                                                          gATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTTAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                          >  2:3435363/1‑150 (MQ=255)
                                                                                                          gATACCACGTACATTTTCATATCCTGTCGTTGGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                          <  2:2411460/150‑1 (MQ=255)
                                                                                                               caattACATTTTAATATCCTTTCGTTTTTTTTCTATTACAAAATATAATCAAATTATTTGTTTTTAAATATCTATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGaa                                     <  1:4196509/146‑1 (MQ=255)
                                                                                                                  ttACTTTTTCATACCCTGTCTTTTGTTTTCTATTACAAAATATAATTAAATTATTTTTTTTAAAATTTCGTGATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                                  <  1:1091433/149‑1 (MQ=255)
                                                                                                                   tACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAATAATGAAATTATTTTTTTTTTAATATAGGGAGAACGATAACAAAAAAGACAAAATGGAAATTATTAGAGGGGGGGAAGGCAAAAAAAAGGCCGCCCCGGAGGGGAGCCg                                 >  2:2326676/1‑150 (MQ=255)
                                                                                                                     cATTTTCATATCCTGTCTTTTGTTTTCTATTTAAAAATATAATTAAATTATTTTTTTTTAAATATCGAGATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                               <  2:923053/150‑1 (MQ=255)
                                                                                                                      aTTTTCATATCCTGTCGTTTGTTTTCGATTACAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                              <  2:2234500/150‑1 (MQ=255)
                                                                                                                        tttCATACCCTGTCTTTTGTTTTCTATTAAAAAATATAATTAAATTTTTTGTTTTAAAATATCGATATAACTATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCa                            <  1:502251/150‑1 (MQ=255)
                                                                                                                         ttCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaa                           <  2:4760798/150‑1 (MQ=255)
                                                                                                                          tCATCTCCTGTCTTTTTTTTTCTATTTAAAAATATAATCAAATTATTTTTATTTAATTATCTATATAACTAACAAAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaaa                          <  1:1811907/150‑1 (MQ=255)
                                                                                                                          tCATATCCTGTCGTTTGTTTTGGTTTTGAAAATATAATGTAATTATTTTTTTTTAAATATGGCAATTAAGATGACAAAAACGAACAAATGTACAATTTTAGAGGGGTGTAAGGCCAAAAAAAGGCCTACCGAGAGGGGAGCCggtgcaca                          >  2:4942165/1‑143 (MQ=255)
                                                                                                                             tatCCTGTCTTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                       <  2:4486088/150‑1 (MQ=255)
                                                                                                                                  tGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgt                  >  2:691569/1‑150 (MQ=255)
                                                                                                                                   gTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                 >  2:1994879/1‑150 (MQ=255)
                                                                                                                                     cGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt               >  1:3664449/1‑150 (MQ=255)
                                                                                                                                     cGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt               >  2:4950602/1‑150 (MQ=255)
                                                                                                                                       tttgtttTCGATTTCAAAATATAATGAATTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt             <  2:4057063/150‑1 (MQ=255)
                                                                                                                                        ttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt            >  2:3284298/1‑150 (MQ=255)
                                                                                                                                           ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg         >  2:1366024/1‑150 (MQ=255)
                                                                                                                                             ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc       >  2:4151258/1‑150 (MQ=255)
                                                                                                                                               cGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCgg     >  1:334106/1‑150 (MQ=255)
                                                                                                                                                 ttttAAAAATATAATTAATTTATTTGTTTTTAAATATCGATATACCGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt   <  1:1914318/149‑1 (MQ=255)
                                                                                                                                                  tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg  <  1:4311935/150‑1 (MQ=255)
                                                                                                                                                  tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg  >  1:2437163/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATG  >  CP009273/2170564‑2170859

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