Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,713 C→A 24.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,7130CA24.9% 187.9 / 14.5 97intergenic (‑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 (21/44);  new base A (19/8);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 1.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCG  >  CP009273/2170564‑2170861
                                                                                                                                                     |                                                                                                                                                       
gTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCATCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTc                                                                                                                                                         <  2:1644262/150‑1 (MQ=255)
gTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTc                                                                                                                                                         <  2:5004811/150‑1 (MQ=255)
   ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaa                                                                                                                                                      <  2:5081861/150‑1 (MQ=255)
   ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTaaaa                                                                                                                                                      >  2:228042/1‑150 (MQ=255)
   ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTaaaa                                                                                                                                                      >  2:395322/1‑150 (MQ=255)
    gCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaaa                                                                                                                                                     <  1:3659101/150‑1 (MQ=255)
     cAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAat                                                                                                                                                    <  1:2510569/150‑1 (MQ=255)
     cAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAat                                                                                                                                                    <  1:691824/150‑1 (MQ=255)
        cGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAtata                                                                                                                                                 <  1:1362027/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                               <  1:723019/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                               <  1:722820/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                               <  2:4487195/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                               <  2:2669681/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                               <  2:4486145/150‑1 (MQ=255)
           tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAATATAATg                                                                                                                                              >  2:1987495/1‑150 (MQ=255)
           tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTCAAAATATAATg                                                                                                                                              >  2:4616500/1‑150 (MQ=255)
            ttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATaaata                                                                                                                                             >  2:3480247/1‑147 (MQ=255)
             tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                            <  1:4818774/150‑1 (MQ=255)
             tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTATATTTAAAAAAAAAATTaa                                                                                                                                            >  1:2522353/1‑150 (MQ=255)
              cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                           <  2:4463116/150‑1 (MQ=255)
               tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                          <  2:3694056/150‑1 (MQ=255)
               tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                          <  1:4762431/150‑1 (MQ=255)
                accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAatt                                                                                                                                         <  1:3083802/150‑1 (MQ=255)
                    ccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTAttt                                                                                                                                     <  1:3927872/150‑1 (MQ=255)
                     cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                                    <  1:603958/150‑1 (MQ=255)
                        ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttt                                                                                                                                 <  1:4739741/150‑1 (MQ=255)
                        ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttt                                                                                                                                 <  1:4739655/150‑1 (MQ=255)
                            aTCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaa                                                                                                                             <  2:22190/150‑1 (MQ=255)
                            aTCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAATAAATAAATTATtttttttttt                                                                                                                             >  2:3504050/1‑148 (MQ=255)
                             tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaaa                                                                                                                            <  1:3203460/150‑1 (MQ=255)
                              cGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAat                                                                                                                           <  1:974359/150‑1 (MQ=255)
                                tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                         <  1:345687/150‑1 (MQ=255)
                                  tCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCg                                                                                                                       <  2:1509140/150‑1 (MQ=255)
                                  tCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCg                                                                                                                       <  2:416786/150‑1 (MQ=255)
                                    gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgag                                                                                                                     <  2:485418/150‑1 (MQ=255)
                                      gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                                   <  1:2771144/150‑1 (MQ=255)
                                       aTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATAAAATGAAAATATTTTTTTTTAAAAATCCCGATa                                                                                                                  >  2:3870585/1‑150 (MQ=255)
                                        ttGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATaa                                                                                                                 <  2:4187871/150‑1 (MQ=255)
                                         tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAc                                                                                                                <  2:1299867/150‑1 (MQ=255)
                                           tgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGTTTTAAAAATATAATGTAATTTTTTTTTTTTTAATATCTAGAttaaga                                                                                                              >  2:3273150/1‑145 (MQ=255)
                                           tgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGTAATTATTTTTTTTTTAATATAGAGATAACGa                                                                                                              >  2:214503/1‑150 (MQ=255)
                                           tgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTAAATTAATTTTTTTTAAAAATAGAGATAAcca                                                                                                              >  2:307510/1‑148 (MQ=255)
                                           tgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGGTTTAGATTTAAAAATATAAAGAAAATATTTTTTTTTAAATATACAGATAAAGa                                                                                                              >  2:5080309/1‑150 (MQ=255)
                                           tgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTAAAAATATAAAGAAATTATTTTTTTTTTAATATAGGGTTAACGa                                                                                                              >  2:213609/1‑150 (MQ=255)
                                            gAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATTAAAATATTTGTTTTTAAATATCGAGATAACGAt                                                                                                             >  2:767628/1‑150 (MQ=255)
                                            gAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAAAGAAATTATTTTTTTTTAAATATAGAGATAACGAt                                                                                                             >  2:767565/1‑150 (MQ=255)
                                              atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTGTTTTGTTTTGAAAAAATAATTTAATTATTTTTTTTTTAATTTGTAGATTACgcata                                                                                                           >  1:1012205/1‑146 (MQ=255)
                                              atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTTTTTTTAAATATAGGGATAAAGATca                                                                                                           >  2:2566909/1‑150 (MQ=255)
                                              atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGGTTTTTGTTTGAATATATAATTTAAATATTTTTTTTTTTATTAAGAGATAAAGATca                                                                                                           >  1:1012937/1‑150 (MQ=255)
                                               tatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcac                                                                                                          <  1:1108668/150‑1 (MQ=255)
                                               tatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcac                                                                                                          <  1:547573/150‑1 (MQ=255)
                                                atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAGTGAAATTATTTGTTTTTAAATATCGAGATAACGATcaca                                                                                                        <  2:377956/150‑1 (MQ=255)
                                                atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcaca                                                                                                        <  2:377881/150‑1 (MQ=255)
                                                          ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTCAAAAAAAAAAGAAAATATTTTTTTTTTAATATAGAGATTACGATAACAAAAACGACac                                                                                              >  2:2271383/1‑149 (MQ=255)
                                                                gaaTAACCGCCGCGCTGTGCGTTTTTCAGCATCTTCTTTGTCGATACCAGTTAAATTTTCATATCCTGTCTTTTGTTTTCTATTTAAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCAAAAAAACGACAATATGaa                                                                                        <  2:4464957/148‑1 (MQ=255)
                                                                 cATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTAAAAATATAGATAAAGATAAAAAAAACGACAAAATTaaa                                                                                       >  2:1693891/1‑150 (MQ=255)
                                                                      ccgccgCGGTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGAAAATTTACATATCCTTTCGTTTTTTTTCGATTTAAAAATATAATTAAATTATTTTTTTTAAAATATCTAGATAACTATAAAAAAAACGACAATATGAAAattat                                                                                  <  1:3480247/150‑1 (MQ=255)
                                                                          cgcgGTGTGCTTTGTTCAGAATCTGCTTTGTCGATACCACGTACATTTTAATATCCTGTCGTTTTTTTTCGTTTTAAAAATATATTTAAATTATTTTTTTTTAAATATCGAGATACCGATCACAAAAACGACAATATGAAAATTATTCga                                                                              <  2:2776971/150‑1 (MQ=255)
                                                                          cgcgCTGTGCGTTGTTCAGAATCGGCTTTGTCGAAACCACGTACATTTACATATCCTGTCGTTTTTTTACTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCga                                                                              <  2:2777110/150‑1 (MQ=255)
                                                                              cTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTAAAAATATAATTTAAATATTTTTTTTTAAATATAGAGGTTAAGATAAAAAAAAAGAAAAAATTAAAATTATTAgagggg                                                                          >  1:2167970/1‑148 (MQ=255)
                                                                                gtgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTGAAAAAAAAATTAAAATAATTTTTTTTAAATATAGAGATAAACAAAACAAAAAAGAAAAAAATAAAAATATTAgaggaggg                                                                        >  1:3057885/1‑148 (MQ=255)
                                                                                gtgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAAAAAAGAAAATATTTTTTTTTAAATAACGATATAAACATCACAAAAAAGACAAAATTAAAAATATTCGTGGtgtg                                                                        >  1:3056662/1‑150 (MQ=255)
                                                                                 tgCGTTGTTCAGCATCTGCTTTGTCTATACCACGTAAATTTTCATATCCTGTCTTTTTTTTTCTATTTCAAAATAAAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGa                                                                      <  2:1266266/150‑1 (MQ=255)
                                                                                    gttgttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTTAATTATTTGTTTTTAAATATAGATATAAAGATCAAAAAAAAAACAAAATTAAAATTATTAGAGGAGTGAAAg                                                                  >  2:1934109/1‑150 (MQ=255)
                                                                                              aTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTTGTTTTAAAAAAATAATTAAAATATTTTTTTTTTTATTTAGAGATTAAGATAAAAAAAAAGACAAAAAGAAAAATATTATAGGAGTTAAATTGCaaaaaaa                                                         >  2:1021149/1‑150 (MQ=255)
                                                                                               tCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAAAAAAAAAAAAATATTTTTTTTTAAAAAAAGAGATAAACATAACAAAAAAAAAAAAAAGTAAAATATTAGATGAGTGTAATGACaaaaaaaa                                                        >  2:498811/1‑149 (MQ=255)
                                                                                               tCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTAAAAAAAAAATTAAATTATTTTTTTTTAAAAATAGTGAAAAAGAAAAAAAAAAAAAAAAAATTTAAAATATTAGAGGATTGTAAGGGaaaaaaaaa                                                       >  1:2103648/1‑148 (MQ=255)
                                                                                               tCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTGTGTTTTTAAAAAAAAATGAAATTATTAGTTTTTATATATAGATATAAACAAGACAAAAAACAAAAAAATTAAAATATACGAGTATTGAAATCaaaaaaaaaa                                                       >  1:2103640/1‑148 (MQ=255)
                                                                                                 tgatttGTCGATACAACTTACATTTTAATATCCTTTCGTTTTTTTTCTATTTAAAAATATATTTAAATTATTTTTTTTTAAATATCTATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGc                                                     <  2:2522353/147‑1 (MQ=255)
                                                                                                         cGATACCACTTATATTTTCTTATCATGTCGTTTTTTTTCTTTTTAAAAATATAATTAAATTATTTGTTTTTAATTATCGATATAACGATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGAt                                             <  2:4467665/150‑1 (MQ=255)
                                                                                                          gATACCACGTACATTTTCATATCCTGTCGTTTTTTTTGGATTTCAAAAAATAATGAAAATATTTTTTTTTAAATATAGAGATAAAGAACACAAAAACGACAAAATGAAAATTATTAGAGGCGTGTAAGGCAAAAAAAAGGCCTCCCGCTa                                            >  2:5223672/1‑150 (MQ=255)
                                                                                                          gATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTTAAATATTTTTTTTTAAATATAGAGATAACGATAACAAAAAACACAATATTAAAAATATTCGAGGAGTGAAAGGCCAAAAAAAGGCCACCCGATa                                            >  2:5224471/1‑150 (MQ=255)
                                                                                                          gATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAAAAAAATGAAATTATTTTTTTTTTAATATAGAGATAAAGATAACAAAAAAGAAAAAATTAAAATTATTAGAGGAGTGTAATGCAAAAAAAACGCCAccccaga                                            >  2:5224608/1‑146 (MQ=255)
                                                                                                           aTACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATAGATATAAAGATCACAAAAAAGAAAATATGAAAATTATTAGTGGATTGAAAGGCAAAAAAAAGGGCTCCCGCTat                                           >  2:2010191/1‑149 (MQ=255)
                                                                                                           aTACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAAAAAAATGAAATTATTTTTTTTTATATATCGAGATAAAGATAACAAAAACGACAATATGAAAATTATTAGAGGTGTGAAAGGGCAAAAAACTGCCTACCGGTAg                                           >  2:3495883/1‑150 (MQ=255)
                                                                                                            tACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTGAAAATATAATTTAATTATTTTTTTTTTTATATCGGGTTTAAGATAACAAAAAACACAAAATGGGCAGTATTGGGGGGGTGGATGGGAAAAAAAAGGCCGCCCgcgagg                                          >  2:4057998/1‑145 (MQ=255)
                                                                                                            tACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAATAATGTAAATATTTTTTTTTTAATATCGAGATGACGATCACAAAAAAGACAAAATGAAAAATATTAGAGGAGGGACATGCGAAAAAACGGCCTCCCGAGAgg                                          >  2:4057948/1‑150 (MQ=255)
                                                                                                                 cttaaaTTTTAATATCCTGTGTTTTTTTTTCTATTAAAAAATATATTGAAATTATTTTTTTTTAAATATCTATATACCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGc                                     <  1:2431808/145‑1 (MQ=255)
                                                                                                                  ttaatttttCATATCCTGTCTTTTGTTTTATATTTAAAAATATAATTAATTTTTTTTTTTTTAAATATCTATATAACTATAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                                    <  1:4616500/145‑1 (MQ=255)
                                                                                                                  ttaatttttCATATCCTGTCGTTTTTTTTCTTTTTCAAAATATAATGAATTTATTTTTTTTTAAATTTCTAGATAACGATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                                    <  1:1987495/145‑1 (MQ=255)
                                                                                                                        tttCATATCCTGTCGTGTTTTTTAGATTTAAAAATATAAAGAAATTATTTTTTTTTAAAAATAGAGCTAACGATAACAAAAAAGACAAAAATAAAAATATTAGAGGAGTGAAAGGCAAAAAAAAGGCCAACCCATATGGAAGCCGTAGCa                              >  2:5166828/1‑150 (MQ=255)
                                                                                                                        ttaaatatCCTGTCGTTTGTTTTCTATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCa                              <  2:2979368/146‑1 (MQ=255)
                                                                                                                           cATATCCTGTCGTTTGTTTTCTATTACAAAATATAATTAAATTATTTGTTTTTAAATATCTATATAACGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAg                           <  1:2554921/150‑1 (MQ=255)
                                                                                                                             tttCCTGTCGTTTTTTTTCGATTACAAAATAAAATGAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                         <  2:1144225/148‑1 (MQ=255)
                                                                                                                              atCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGATATACCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgc                        <  1:2596408/150‑1 (MQ=255)
                                                                                                                                     cGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                 <  2:2455835/150‑1 (MQ=255)
                                                                                                                                        ttgtttTCGATTTCAAAATATATTTAAATTATTTTTTTTTAAATATCGATATAACGATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt              <  1:395322/150‑1 (MQ=255)
                                                                                                                                        ttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt              >  1:2959834/1‑150 (MQ=255)
                                                                                                                                        ttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt              >  1:3971814/1‑150 (MQ=255)
                                                                                                                                         tgtttTCTATAACAAAATATATTCAATTTATTTTTTTTTAAATATCGATAAAACTAAAAAAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa             <  1:4422524/150‑1 (MQ=255)
                                                                                                                                          gtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt            >  1:3341732/1‑150 (MQ=255)
                                                                                                                                             ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc         >  1:1847877/1‑150 (MQ=255)
                                                                                                                                             ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc         >  1:1847854/1‑150 (MQ=255)
                                                                                                                                             ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATGCGAGGAGTGAAAGGCAAAATAACGGCCTCCCGATAGGGAAGCCGTAGCACAGTGCGCGTGTGGTTATGcc         >  1:448815/1‑150 (MQ=255)
                                                                                                                                                 aTTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGCGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt     <  2:4713653/150‑1 (MQ=255)
                                                                                                                                                 aTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt     <  1:2575872/150‑1 (MQ=255)
                                                                                                                                                  tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg    >  1:4631949/1‑150 (MQ=255)
                                                                                                                                                   ttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGc   >  1:811197/1‑150 (MQ=255)
                                                                                                                                                    tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGTGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg  >  2:71171/1‑150 (MQ=255)
                                                                                                                                                    tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGGTTTTATGCCGGATGCg  >  2:73531/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCG  >  CP009273/2170564‑2170861

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: