Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,713 C→A 12.3% 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,7130CA12.3% 539.2 / 9.2 228intergenic (‑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 (66/117);  new base A (29/11);  total (100/128)
Fisher's exact test for biased strand distribution p-value = 3.43e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

GTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGG  >  CP009273/2170564‑2170862
                                                                                                                                                     |                                                                                                                                                     
 ttGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCa                                                                                                                                                      <  2:2480456/150‑1 (MQ=255)
 ttGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCa                                                                                                                                                      <  2:3987965/150‑1 (MQ=255)
   ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaa                                                                                                                                                    <  1:4072166/150‑1 (MQ=255)
      aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                                 <  2:729375/150‑1 (MQ=255)
      aGCGGGTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                                 <  2:4380947/150‑1 (MQ=255)
       gCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAatat                                                                                                                                                <  1:4333258/150‑1 (MQ=255)
       gCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAatat                                                                                                                                                <  1:606318/150‑1 (MQ=255)
       gCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAatat                                                                                                                                                <  2:4857393/150‑1 (MQ=255)
       gCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAatat                                                                                                                                                <  1:1697473/150‑1 (MQ=255)
        cGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAtata                                                                                                                                               <  2:4670058/150‑1 (MQ=255)
        cGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAtata                                                                                                                                               <  1:2812699/150‑1 (MQ=255)
         ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTAAAAATATaa                                                                                                                                              >  1:3145370/1‑150 (MQ=255)
         ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                              <  2:2742864/150‑1 (MQ=255)
         ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCAGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                              <  1:1057812/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                             <  2:4971518/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                             <  1:4439654/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                             <  1:3277197/150‑1 (MQ=255)
              cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                         <  1:1670924/150‑1 (MQ=255)
              cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                         >  2:3620913/1‑150 (MQ=255)
                accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGAATTCAAAATATAATTAAatt                                                                                                                                       >  2:455289/1‑150 (MQ=255)
                   accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                                    >  2:4273041/1‑150 (MQ=255)
                        ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttt                                                                                                                               <  2:1948656/150‑1 (MQ=255)
                          gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttttt                                                                                                                             <  2:1419899/150‑1 (MQ=255)
                           cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTa                                                                                                                            >  2:1145603/1‑150 (MQ=255)
                             tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaaa                                                                                                                          <  2:3463675/150‑1 (MQ=255)
                               gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAata                                                                                                                        <  2:3167221/150‑1 (MQ=255)
                                tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                       >  2:3729950/1‑150 (MQ=255)
                                    gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAAATATTTTTTTTTAAATATCgag                                                                                                                   >  1:3853659/1‑150 (MQ=255)
                                     agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                                  <  2:4152005/150‑1 (MQ=255)
                                     agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAAAAAAATGAAATTATTTTTTTTTAAATATAgaga                                                                                                                  >  2:1758020/1‑150 (MQ=255)
                                       aTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATa                                                                                                                <  1:1974734/150‑1 (MQ=255)
                                       aTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATa                                                                                                                <  2:674723/150‑1 (MQ=255)
                                         tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAc                                                                                                              <  1:2091960/150‑1 (MQ=255)
                                         tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAc                                                                                                              <  1:4965326/150‑1 (MQ=255)
                                         tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAc                                                                                                              >  2:700444/1‑150 (MQ=255)
                                          gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACg                                                                                                             >  1:422287/1‑150 (MQ=255)
                                           tgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGa                                                                                                            >  2:2054270/1‑150 (MQ=255)
                                              atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATca                                                                                                         <  1:1748780/150‑1 (MQ=255)
                                              atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTGGTTTTCGATTTAAAAAAATAATGAAATTATTTTTTTTTAAATATCGAGATAAAGATca                                                                                                         >  1:2272173/1‑150 (MQ=255)
                                                atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcaca                                                                                                       <  1:1391346/150‑1 (MQ=255)
                                                atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcaca                                                                                                       <  2:1095112/150‑1 (MQ=255)
                                                    aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTTAAATATTTGTTTTTTAATATCGAGTTAAAGATaaaaaaaa                                                                                                   >  1:1040640/1‑150 (MQ=255)
                                                     aTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAc                                                                                                  <  2:5080467/150‑1 (MQ=255)
                                                      tGCCGGAACCTCATAACCGCCGCGCTGTGCGTAGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCACATTAAAAAAAAAAAGAAAAAATTTTTTTTTATATATAAACATAAAGATCAaaaaaaag                                                                                                 >  1:4812894/1‑148 (MQ=255)
                                                          ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACaa                                                                                             >  2:4177291/1‑150 (MQ=255)
                                                            aaCCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAata                                                                                           <  2:3440157/150‑1 (MQ=255)
                                                              ccGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTGGATACCCCGTACATTTTCATATCCTGTCGTTTGTTTTCGTTTTCAAAATATAATGAAATTATTTGTTTTTAAATTTCGAGATAACGATCACAAAAACGACAATATg                                                                                         <  2:3231893/150‑1 (MQ=255)
                                                              ccGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAACACGACAATATg                                                                                         <  2:120176/150‑1 (MQ=255)
                                                               cGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGa                                                                                        <  2:551128/150‑1 (MQ=255)
                                                               cGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTAAATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGa                                                                                        <  2:3104567/150‑1 (MQ=255)
                                                                gCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTGAAAATATAAATAAATTATTTTTTTTTTAATATAGAGATAACGAACACAAAAACCACAATATGaa                                                                                       >  1:3050405/1‑150 (MQ=255)
                                                                gCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGaa                                                                                       <  2:3278713/150‑1 (MQ=255)
                                                                 cATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAATATAATTAAAATATTTTTTTTTAAATATAGAGATAAAGATAAAAAAAAAGAAAATATGaaa                                                                                      >  2:1825404/1‑150 (MQ=255)
                                                                    aaCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAAATATTTGTTTTTAAATATCGAGATAACGATAAAAAAAAAGACAATATTAACatt                                                                                   >  2:4469614/1‑150 (MQ=255)
                                                                      ccgccgCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATAGAGATAACGATCACAAAAACGACAATATGAAAattat                                                                                 >  1:3412587/1‑150 (MQ=255)
                                                                       tgccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAGTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATCTGGAAatgttt                                                                                >  2:2652484/2‑146 (MQ=255)
                                                                           gcgcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAAATATTTGTTTTTAAATATCGAGATAACGATAACAAAAAAGACAAAATTAAAATTATTCgag                                                                            >  2:2574755/1‑150 (MQ=255)
                                                                             gcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCTATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAAGATCACAAAAACGACAATAAGAAAATTATTCgagga                                                                          >  2:2028508/1‑150 (MQ=255)
                                                                             gcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTGAAAATATAATGAAATTATTTTTTTTTTAATATCGAGATAACGATCACCAAAACGACCATATGAAAATTATTGgaggt                                                                          >  2:21488/1‑149 (MQ=255)
                                                                             gcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATAAAATGAAAATATTTTTTTTTAAATATAGAGATAACCATAACAAAAACGACAAAATGAAAATTAATAgagga                                                                          >  1:2877786/1‑150 (MQ=255)
                                                                              cTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCgaggag                                                                         >  2:4627374/1‑150 (MQ=255)
                                                                               tgtgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATCTCCTGTCGTTGGTTTTCGATTGAAAAATATCATGAAATTAATTGTTTGTAAATATCGAGATAACGAACCCAAAAACCACAAAAAGAAAATTATTCgaggact                                                                        >  2:1566041/1‑148 (MQ=255)
                                                                               tgtgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGTTTTCAAAATATAATGAAAATATTTTTTTTTAAATATAGAGATAAAGATAAAAAAAAACACAAAATTAAAATTATTAGAGGGGt                                                                        >  1:4536889/1‑150 (MQ=255)
                                                                                gtgCGTTGTTCAGCATCTGCTTTGTCGATCCCCCGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATCAAATTATTTGTTTTTAAATATCGAGATAACGAACACAAACACGACAATATGAAAATTATTCGAGGAGTg                                                                       <  1:2725043/150‑1 (MQ=255)
                                                                                gtgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTg                                                                       <  2:1438008/150‑1 (MQ=255)
                                                                                 tgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTAGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTTAATATCGAGATAACGATCACAAAAACGACAAAATGAAAATTATTCGAGGAGTGa                                                                      >  1:296654/1‑150 (MQ=255)
                                                                                    gttgttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGtaat                                                                   >  2:1448203/1‑149 (MQ=255)
                                                                                      ttttCAGCATCTGCTTTGTCGATACCACTTACTTTTTCATATCCTGTCTTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGc                                                                 <  1:1267250/148‑1 (MQ=255)
                                                                                       gttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATAACAAAAACGACAATATTAAAATTATTAGAGGAGTGAAAGGCa                                                                >  2:66265/1‑150 (MQ=255)
                                                                                        ttCAGCATCTGCTTTGTCGAAACCCCTTAAATTTTAATAACCTGTCGTTTGTTTACTATTAAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACTATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaa                                                               <  1:3462669/150‑1 (MQ=255)
                                                                                        ttCAGCATCTGATTTGTCGAAACCACGTACATTTTCATATCCTGTGTTTTTTTTTGTATTTCATAATATAATTAATTTATTTGTTTTTAAATATCGTTTTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaa                                                               <  1:2094040/150‑1 (MQ=255)
                                                                                        ttCAGAATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaa                                                               <  1:2546764/150‑1 (MQ=255)
                                                                                         tCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATTAAAATATTTTTTTTTAAATATAGAGATAACGATTACAAAAAAGAAAAAATTAAAATTATTAGAGGAGTGTAATGCaaa                                                              >  1:684242/1‑150 (MQ=255)
                                                                                               tCTGCTTTGTCTATACAAAGAAATTTTTAATATCCTGTCTTTTTTTTTCGATTAAAAAATATAATTAAATTATTTGTTTTTAATTATCGAGATAACTATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACg                                                        <  1:4171069/150‑1 (MQ=255)
                                                                                               tCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACg                                                        <  2:288232/150‑1 (MQ=255)
                                                                                                    tttGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTc                                                   <  2:4747139/150‑1 (MQ=255)
                                                                                                    tttGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTc                                                   <  2:822851/150‑1 (MQ=255)
                                                                                                    tttGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTGAAAAAATAATTTAAATAATTTTTTTTTAATATAGAGATAAACATAACAAAAACCACACAAAGAAAAATTTTAGAGGAGTTGAAGGCaaaaaaaatgccac                                                   >  1:4582212/1‑143 (MQ=255)
                                                                                                     ttGTCGATACCACTTAATTTTTATTTTCCTGTCGTTTGTTTTCGATTAAAAATTATATTTATATTATTTTTTTTTAATTATCTATATAACTATAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTcc                                                  <  2:1343615/150‑1 (MQ=255)
                                                                                                      tGTCTATACCACGTACTTTTTCTTTTCCTTTGGTTTGTTTTGTTTTTAAAAATATATTTAAATTTTTTTTTTTTAAATATCGATTTAAAGATCAAAAAAACGACAATATTAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTccc                                                 <  1:3283836/150‑1 (MQ=255)
                                                                                                       gTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAAAAATGAAAATATTTTTTTTTAAATATAGAGATAAAGATAAAAAAAAAAAAAATATTAAAATTAATAGAGGCGTGTAAGGCCAAAAAAACGCCAcccc                                                >  1:3564177/1‑149 (MQ=255)
                                                                                                         ctaTACAACTTAAATTTTAATATCCTTTCTTTTTTTTTCGATTTCAAAATATATTTAATTTTTTTGTTTTTAAATATCTATAAAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGAt                                              <  1:1791349/148‑1 (MQ=255)
                                                                                                          taTACCACGTAATTTTTAATTTCCTTTCTTTTGTTTTCGATTTCAAATTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                             <  2:1074042/149‑1 (MQ=255)
                                                                                                          gATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                             <  2:1754065/150‑1 (MQ=255)
                                                                                                               cACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGTTAACGATCACAAAAACGACAATATTAAAATTATTAGAGGAGTTAAAGGCAAAAAAAAGGCCTCCCGATAgggta                                        >  2:1963820/1‑148 (MQ=255)
                                                                                                                aCGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAg                                       <  2:4849971/150‑1 (MQ=255)
                                                                                                                   tACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCg                                    <  2:3483400/150‑1 (MQ=255)
                                                                                                                      aTTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                 <  1:1720697/150‑1 (MQ=255)
                                                                                                                           cATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAg                            <  2:503556/150‑1 (MQ=255)
                                                                                                                               tCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATAGAGATAACGATAACAAAAAAGACAAAATGAAAATTATTCGAGGAGTGAAATGCAAAAAAAACGCCTACCCATATGGTAGCCCTAGCAAACTgcg                        >  2:10050/1‑150 (MQ=255)
                                                                                                                                 cTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                      >  2:1731747/1‑150 (MQ=255)
                                                                                                                                   gTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                    <  2:2409970/150‑1 (MQ=255)
                                                                                                                                       tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                >  1:4278790/1‑150 (MQ=255)
                                                                                                                                          gtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt             >  2:881936/1‑150 (MQ=255)
                                                                                                                                          gtttTCGATTTCAAAAGATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt             <  2:2703125/150‑1 (MQ=255)
                                                                                                                                             ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc          <  2:3803798/150‑1 (MQ=255)
                                                                                                                                              tCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTGTTTATGCCg         >  2:1902186/1‑150 (MQ=255)
                                                                                                                                                 aTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt      <  1:4273041/150‑1 (MQ=255)
                                                                                                                                                  tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg     <  1:5115390/150‑1 (MQ=255)
                                                                                                                                                     cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg  >  1:88731/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGG  >  CP009273/2170564‑2170862

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: