Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,713 C→A 26.0% 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,7130CA26.0% 441.4 / 33.5 243intergenic (‑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 (62/95);  new base A (48/22);  total (126/117)
Fisher's exact test for biased strand distribution p-value = 5.53e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

TTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCG  >  CP009273/2170565‑2170861
                                                                                                                                                    |                                                                                                                                                    
ttGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTc                                                                                                                                                      <  2:674089/149‑1 (MQ=255)
 tGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTCa                                                                                                                                                     >  1:2976389/1‑149 (MQ=255)
    cAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaaa                                                                                                                                                  <  1:133867/149‑1 (MQ=255)
     aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGAGTTGaaaaa                                                                                                                                                 >  1:944690/1‑148 (MQ=255)
       cGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAatat                                                                                                                                               <  1:583678/149‑1 (MQ=255)
         gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTGGATTTAAAAATATaa                                                                                                                                             >  1:1973753/1‑149 (MQ=255)
          tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                            >  2:3192350/1‑149 (MQ=255)
           ttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATg                                                                                                                                           <  1:4534649/149‑1 (MQ=255)
              tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                        <  2:293900/149‑1 (MQ=255)
              tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                        <  2:1905419/149‑1 (MQ=255)
              tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTATATTTAAAAAAAAAATTaaa                                                                                                                                        >  2:109527/1‑149 (MQ=255)
              tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAATATaaataaa                                                                                                                                        >  2:1711531/1‑144 (MQ=255)
                 caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTTGATTTAAAAAAATAATTAAaata                                                                                                                                     >  1:1043749/1‑149 (MQ=255)
                   ccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                                   <  2:781887/149‑1 (MQ=255)
                    cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTAAAATAttt                                                                                                                                  >  2:1250580/1‑149 (MQ=255)
                      cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGt                                                                                                                                <  2:1330623/149‑1 (MQ=255)
                         gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtttt                                                                                                                             <  1:1019995/149‑1 (MQ=255)
                           aTCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTa                                                                                                                           >  2:4277941/1‑149 (MQ=255)
                           aTCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTCAAAATATAATGAAAATATTTGtttttt                                                                                                                           >  2:2727391/1‑148 (MQ=255)
                           aTCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAATAAAATTAAAATATTTTTTTTTa                                                                                                                           >  1:1940112/1‑149 (MQ=255)
                           aTCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTa                                                                                                                           <  1:3215766/149‑1 (MQ=255)
                           aTCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTATATTTAAAAAAATAAATAAAAAATttttttttt                                                                                                                           >  1:16629/1‑148 (MQ=255)
                            tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCAGAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaa                                                                                                                          <  2:4478348/149‑1 (MQ=255)
                             cGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaaa                                                                                                                         <  2:1349278/149‑1 (MQ=255)
                              gTTTCGAGAGTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAat                                                                                                                        <  2:3661378/149‑1 (MQ=255)
                               tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAAAATAATTAAAATAATTTTTTTTaaaaa                                                                                                                       >  2:4239872/1‑147 (MQ=255)
                                 tCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGAATTAAAAAAATAAATAAATTATTTTTTTTTTAaaata                                                                                                                     >  1:3536648/1‑148 (MQ=255)
                                     gtttGTGTATATTGAATGCCGGAACCGCATAACCTCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTAGATACCACGTACATTTTCAGATCCTGTCGTTTGTTTTCGATTTCAAAATATATTGAAATTATTTGTTTTTATATTTCgaga                                                                                                                 <  2:2419702/147‑1 (MQ=255)
                                     gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                                 <  2:1060409/149‑1 (MQ=255)
                                     gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                                 <  1:2499862/149‑1 (MQ=255)
                                     gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAATAATTAAAATATTTTTTTTTTAATATAgaga                                                                                                                 >  1:2933950/1‑149 (MQ=255)
                                      aTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                                <  1:1699283/149‑1 (MQ=255)
                                        tgtgAATATTGAATGCCGGAACCGCATATCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATaa                                                                                                              <  2:3335681/149‑1 (MQ=255)
                                        tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTGTTAGATTTAAAAAAATAATGAAATTATTTGTTTTTATATATCgagctca                                                                                                              >  1:3729132/1‑145 (MQ=255)
                                         gtgAATTTTGATTGCCGGAACCGCATATCCGCCGCGCTGTGAGTTGTTCAGGATCTGCTTTGTCGATACCACGTACTTTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATTTCGAGATatc                                                                                                             <  1:3673453/149‑3 (MQ=255)
                                           gAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGa                                                                                                           <  1:2242559/149‑1 (MQ=255)
                                            aaTAGTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGAt                                                                                                          <  2:2497703/149‑1 (MQ=255)
                                             agatTGAATGCCGGATCCGCATAACCGCCGCTCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATc                                                                                                         <  2:4430435/147‑1 (MQ=255)
                                              tatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATca                                                                                                        <  2:2523449/149‑1 (MQ=255)
                                               atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAAATATTTTTTTTTAAATATCGAGATAACGATcac                                                                                                       >  2:3961936/1‑149 (MQ=255)
                                               atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTTTTTTTAAATATAGAGATAACGATcac                                                                                                       >  2:2545603/1‑149 (MQ=255)
                                                ttGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTTAATTATTTTTTTTTAAATATAGCGATAAAGATaaca                                                                                                      >  2:562131/1‑149 (MQ=255)
                                                    aTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                                  <  2:3236800/149‑1 (MQ=255)
                                                     tGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACAATCACAAAAAc                                                                                                 <  2:960173/149‑1 (MQ=255)
                                                       ccGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGa                                                                                               <  2:1180091/149‑1 (MQ=255)
                                                         ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTATATTTAAAAATATAAATTAAATATTTTTTTTTAAAAATAGAGATAAAGATCACaaaaaacaca                                                                                             >  1:1145551/1‑144 (MQ=255)
                                                         ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTAAAAAAATAATTAAATTATTTTTTTTTAAATATATAGATAAATATAACaaaaaaaaca                                                                                             >  1:743565/1‑144 (MQ=255)
                                                          gAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTGTTGGATTTGAAAAAATAAATAAAATATTTTTTTTTTTTTATATAGATGAAGATGACAAAAAAAACaa                                                                                            >  1:803882/1‑149 (MQ=255)
                                                           aaCCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTGAAAATATAATGTAAATATTTTTTTTTAAAAATAGTGATTAAGATAACCAAAACGACaaa                                                                                           >  1:4285796/1‑148 (MQ=255)
                                                            acctaaTAACCGCCGCGGTGTGTTTTGTTCTGCTTGTGATTTGTCTATACCACTAACTTTTTAATTTCCTGTGTTTTTTTTTGGTTTTAAAAATATAATTAATTTTTTTTTTTTTAATTATGTATATAACGATCACAAAAACGACAata                                                                                          <  2:4106002/144‑1 (MQ=255)
                                                             ccGCATAACCGCCGCGCTGTGGGTTGTTCAGCATCTGATTTGTCTATACCACTTACTTTTTCATATCCTGTGGTTTGTTTTCTATTTCAAAATATAATGAATTTTTTTGTTTTTAATTTTCGAGATAACGATCAAAAACACGACAatat                                                                                         <  2:92515/149‑1 (MQ=255)
                                                                 aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTAAAAATATAATTAAAATATTTTTTTTTTAATATAGAGGTAACGATAACAAAAAACACAAAATGaaa                                                                                     >  1:2644220/1‑149 (MQ=255)
                                                                  tAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAATGAAAATATTTTTTTTTAAAAATCGAGATAACGAACAAAAAAAACACAAAAAGaaaa                                                                                    >  2:4100214/1‑149 (MQ=255)
                                                                  tAACCGCCGCGCTGTGCGTTGTTCAGAATCTGCTTTGTCGATACCACGTAAATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACTATCACAAAAACGACAATATGaaaa                                                                                    <  2:3569567/149‑1 (MQ=255)
                                                                     ccgccgCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACAGTTTCATATCCTGTCGTGTGTTTTGGATTTAAAAAAAAAAAGAAATTATTTTTTTTTTAATATAGAGATAACGAACACAAAAAAGACAAAAATAAAatta                                                                                 >  1:1540594/1‑149 (MQ=255)
                                                                      cgccggGGTGTGGGTTGTTCAGCATCTGCTTTGTCGATCCCCCGTAAATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATTAATTTTTTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAattat                                                                                <  1:2677480/149‑1 (MQ=255)
                                                                      cgccgcGCTGTGCTTTTTTCAGAATCTGATTTGTCGATACAACTTAATTTTTAATATCGTGTCTTTTTTTTTCTTTTTCAAATTATAATTAATTTATTTTTTTTTAAATATCTAGATAACGAACAAAAACACGACAATATGAAAattat                                                                                <  1:1492758/149‑1 (MQ=255)
                                                                       gccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAattatt                                                                               <  2:1143767/149‑1 (MQ=255)
                                                                       gccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATTAAAATATTTTTTTTTAAATATAGAGATAACCATAACAAAAACCACAAAATGAAAaatatt                                                                               >  2:1701763/1‑149 (MQ=255)
                                                                            gcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGTTTTCTAAAAATAATGAAAATATTTTTTTTTTAAAATAGAGTTAAAGATAACAAAAAAACCAAAATTGAAATTATTCtagg                                                                          >  2:805991/1‑149 (MQ=255)
                                                                             cTGTGGTTTGTCCAGCATCGGCTTTGTCTAAACCACTTACATTTTCATAACCTTTCTTTTTTTTTCTATTAAAAAATATAATCAAATTTTTTTTTTTTAAATATCGAGATAACTATCACAAAAACGACAATATGAAAATTATTCgagga                                                                         <  1:4718189/149‑1 (MQ=255)
                                                                                tgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTGAAAATATAATTTAATTATTTTTTTTTTAATATAGATATAAAGATCACAAAAAAGACAAAATTGAAATTTTTAGAGGGGTg                                                                      >  1:2690813/1‑149 (MQ=255)
                                                                                        tCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAAAATAATGAAATTATTTTTTTTTAAATATAGAGATAAACATAACAAAAACGACAAAATGAAAATTATTAGAGGAGTGGAAGGCaa                                                              >  2:4264202/1‑149 (MQ=255)
                                                                                         cAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTCAAAATATAATTAAATTATTTTTTTTTTAATATATAGATAAAGATAAAAAAAAAGACAAAATTAAAAATATTAGAGGAGTGTAAGTCaaa                                                             >  2:2166140/1‑149 (MQ=255)
                                                                                          aGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaaaa                                                            <  2:627677/149‑1 (MQ=255)
                                                                                          aGCATCTGCTTTGTCGATACCACGTACAGTTTCATATCCTGTCGTGTTTTTTAGATTTAAAAAAATAATGTAAATATTTTTTTTTTAAAAAAGAGATAAAAATAACAAAAACAAAAAAATTAAAAATATTGTAGGTGTGGAAGGCcaaa                                                            >  1:2270486/1‑149 (MQ=255)
                                                                                             aTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTTGATTTAAAAAAAAAATTAAAATATTTTTTTTTTAATATAGAGTTAAAGATCACAAAAAAGAAAATATTAAAAATATTATATGAGTGAAAGGCaaaaaaa                                                         >  1:1080140/1‑149 (MQ=255)
                                                                                             aTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTGAAAAAAAAATTAAAAAATTTTTTTTTTAATATAGCGAAAACCATAACAAAAACAAAAAAAATAAAAAAATTAGAGGAGGGAAAAGCacaaaaa                                                         >  1:4651389/1‑149 (MQ=255)
                                                                                               cTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTAAAAATATAATGAAAATATTTTTTTTTAAATAACGAGATAAAGATAAAAAAAACGACAATATTAAAATTATTAGTGGGGTGAGAGGGaaaaaaaag                                                       >  2:3521238/1‑147 (MQ=255)
                                                                                                 gCTTTGTCGATACCACGTACAGTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAAAATAATGAAATTAATTGTTTTTAAATATCGAGATAACCATCACAAAAACCACCATATGAAAATTATTCGAGGAGGGAAATGCCAAAAAAAGGc                                                     >  2:1973671/1‑149 (MQ=255)
                                                                                                   tttGTCTATACCACGTACATTTTCATATCCTGTGGTTTTTTTTCGATTTAAAAATATAATTAAATTTTTTTTTTTTAAATATCTATATAACGATAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCt                                                   <  2:1821303/149‑1 (MQ=255)
                                                                                                   tttGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTTGTTTTAAAAATATAATTTAATTATTTTTTTTTTAATATAGATTTAAACATAAAAAAAAAGAAAATATTTAAATTATTATGGGAGTGGAATGCaaaaaaaaaggct                                                   >  2:2080768/1‑143 (MQ=255)
                                                                                                     tGTCGAAACCACGTACATTTCAATATCCTGTCGTTTTTTTTCTATTTAAAAATATAATCAAATTATTTTTTTTTAAATATCTAGATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTcc                                                 <  2:32808/149‑1 (MQ=255)
                                                                                                        cGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAAAAAAATGAAATTATTTTTTTTTTAATATAGAGAAAAAGATCAAAAAAAACACAAAATTAAAATTTTTAGAGGTGTGAAAGGCAAAAAAAAGGCCACCCGa                                              >  2:1421696/1‑149 (MQ=255)
                                                                                                             ccACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATAGAGATAAAGATAACCAAAACGACAAAATTTAAATTATTGGAGGAGTGGTATGCCAAAAAAATGCCTCCCGCTAggg                                         >  2:2999377/1‑149 (MQ=255)
                                                                                                              cACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTGTAAATATAATTGAAATATTTGTTTTTTAATATAGAGTTAACGATAGCAAAAAAGACCCTATGAAAATTTTTGGCGGAGTTGTAGGCCACAAAACGGCCGCCCGAGAGGGa                                        >  2:3787194/1‑149 (MQ=255)
                                                                                                                   aCATTTTCATATCCTGTCGTTTTTTTTTGATTTTAAAATATAATGAAAATATTTTTTTTTAAATATAGAGATAAAGATCACAAAAAACACACTATTAAAATTTTTAGAGGGGTGTAATGCAAAAAAAATGCCGCCCCATATGGCAGccc                                   >  2:1316994/1‑148 (MQ=255)
                                                                                                                     aTTTTAATATCCTGTCTTTTTTTTTATATTAAAAAATATATTAAAATTATTTTTTTTTAAATATCTATATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                 <  1:3485941/149‑1 (MQ=255)
                                                                                                                      ttttCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                <  2:2648273/149‑1 (MQ=255)
                                                                                                                        ttCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTTTTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCa                              <  1:394402/149‑1 (MQ=255)
                                                                                                                                   tCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                   >  2:2991614/1‑149 (MQ=255)
                                                                                                                                    cGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATACCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtgt                  <  2:1587178/149‑1 (MQ=255)
                                                                                                                                     ttttgtttTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                 <  2:532666/148‑1 (MQ=255)
                                                                                                                                     gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTTTTCGAGGAGTGGAAGGCAAAATAACGGCCTCCCGAGCGGGGCGCCGTTGCAAAGTGCGCGGGtt                 >  2:4509539/1‑149 (MQ=255)
                                                                                                                                     gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                 >  1:2300068/1‑149 (MQ=255)
                                                                                                                                     gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                 >  1:2299960/1‑149 (MQ=255)
                                                                                                                                      tttgtttTCTATTAAAAAATATATTTAAATTATTTTTTTTTAAATATCGAGATAACGATCACAACAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt                <  1:330093/149‑1 (MQ=255)
                                                                                                                                       ttgtttTCGTTTTAAAAATTTAATGAAATTTTTTTTTTTTAATTATCGTGATACGGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt               <  1:4167352/149‑1 (MQ=255)
                                                                                                                                       ttgtttTCGATTTCAAAATATAATTAAATTTTTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt               <  2:3303786/149‑1 (MQ=255)
                                                                                                                                       ttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt               >  1:3997325/1‑149 (MQ=255)
                                                                                                                                            ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc          >  1:3188622/1‑149 (MQ=255)
                                                                                                                                              cGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGACAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCTAAGTGCGCGTGTTTTTATGCCg        >  2:3394787/1‑149 (MQ=255)
                                                                                                                                              cGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCg        >  2:1108775/1‑149 (MQ=255)
                                                                                                                                              cGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCg        >  1:3821961/1‑149 (MQ=255)
                                                                                                                                              cGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCg        >  1:2937885/1‑149 (MQ=255)
                                                                                                                                               gATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAACGTGCGCGTGTTTTTATGCCgg       >  2:4571753/1‑149 (MQ=255)
                                                                                                                                               gATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCgg       >  1:3391221/1‑149 (MQ=255)
                                                                                                                                                 tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGACGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGag     >  2:2763908/1‑148 (MQ=255)
                                                                                                                                                   tCAAAATATAATGAAATTATTTGTTTTTAAATATCGATATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGc   <  1:3448720/149‑1 (MQ=255)
                                                                                                                                                    cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg  <  1:201921/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCG  >  CP009273/2170565‑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: