Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,713 C→A 30.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,7130CA30.3% 171.7 / 19.1 104intergenic (‑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 (16/48);  new base A (23/9);  total (47/57)
Fisher's exact test for biased strand distribution p-value = 1.80e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

TGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGG  >  CP009273/2170566‑2170862
                                                                                                                                                    |                                                                                                                                                      
tGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGTTTTCaa                                                                                                                                                      <  2:1327074/150‑1 (MQ=255)
 ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaa                                                                                                                                                     <  2:3508162/150‑1 (MQ=255)
     gCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAatat                                                                                                                                                 <  2:1770711/150‑1 (MQ=255)
      cGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAtata                                                                                                                                                <  2:664397/150‑1 (MQ=255)
      cGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAtata                                                                                                                                                <  1:3339807/150‑1 (MQ=255)
       ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATaa                                                                                                                                               >  2:3351698/1‑150 (MQ=255)
       ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                               <  2:2538735/150‑1 (MQ=255)
        gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGCTTTCAAAATATAAt                                                                                                                                              >  2:2037505/1‑150 (MQ=255)
        gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTCAAAATATAAt                                                                                                                                              >  2:4422906/1‑150 (MQ=255)
         tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATg                                                                                                                                            <  2:4241735/150‑1 (MQ=255)
         tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATg                                                                                                                                            <  1:4119311/150‑1 (MQ=255)
         tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATg                                                                                                                                            <  1:2797883/150‑1 (MQ=255)
          ttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGa                                                                                                                                           <  1:4276180/150‑1 (MQ=255)
          ttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGa                                                                                                                                           <  1:4276187/150‑1 (MQ=255)
          ttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGGTTTTGATTTAAAAATATAAtta                                                                                                                                            >  1:1665897/1‑148 (MQ=255)
           tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                          <  1:129485/150‑1 (MQ=255)
           tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                          <  1:2285329/150‑1 (MQ=255)
           tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                          <  1:2285113/150‑1 (MQ=255)
           tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                          <  1:129511/150‑1 (MQ=255)
            cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                         <  1:3069425/150‑1 (MQ=255)
             tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGTTTTAAAAAAATaaataaaa                                                                                                                                         >  1:3319112/1‑144 (MQ=255)
              accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTTGGTTTTAAAAAATaaattaaat                                                                                                                                        >  1:2036536/1‑143 (MQ=255)
              accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTCGGTTTGAAAATATAAtttgaat                                                                                                                                        >  1:3909142/1‑144 (MQ=255)
                  ccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTAttt                                                                                                                                   <  1:2901812/150‑1 (MQ=255)
                    acTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGt                                                                                                                                 <  2:2630223/150‑1 (MQ=255)
                    acTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGt                                                                                                                                 <  2:2631416/150‑1 (MQ=255)
                      ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttt                                                                                                                               <  1:2733062/150‑1 (MQ=255)
                       tGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtttt                                                                                                                              <  2:224996/150‑1 (MQ=255)
                       tGCATCGTTTCGAGATTGAGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtttt                                                                                                                              <  2:224816/150‑1 (MQ=255)
                        gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGTAAATATtttttttt                                                                                                                             >  2:3412478/1‑150 (MQ=255)
                        gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTAAAATATtttttttt                                                                                                                             >  2:2191463/1‑150 (MQ=255)
                          aTCGTTTCGAGATTGTTAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaa                                                                                                                           <  2:3556465/150‑1 (MQ=255)
                            cGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTTTATTTAAAAATATTAATAAAtttttttttttttttt                                                                                                                         >  1:4191188/1‑146 (MQ=255)
                             gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAata                                                                                                                        <  2:458690/150‑1 (MQ=255)
                             gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACAGTTTCATATCCTGACGTTTGTTTTAGATATATAAATATAATGTAAATATATTTTTTTAAata                                                                                                                        >  1:2013441/1‑150 (MQ=255)
                                 cGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTAAAAAAAAAAATAAAAATTTTTTTTTTTTatatata                                                                                                                      >  1:605481/1‑149 (MQ=255)
                                 cGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGGTTTAAAAAAAAAAATAAAATATTTTTTTTTtaaaaaata                                                                                                                    >  1:605192/1‑144 (MQ=255)
                                 cGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGTATTTAAAAAAATAATAAAATTATTTTTTTTTTAaaatata                                                                                                                    >  1:605271/1‑147 (MQ=255)
                                   agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                                  <  1:376262/150‑1 (MQ=255)
                                   agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                                  <  1:3252264/150‑1 (MQ=255)
                                      ttGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATGAAATTATTTGTTTTTTAATATAGAGATaa                                                                                                               >  2:2360282/1‑150 (MQ=255)
                                       tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAc                                                                                                              <  2:1445450/150‑1 (MQ=255)
                                       tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAc                                                                                                              <  2:3836316/150‑1 (MQ=255)
                                            atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATca                                                                                                         <  1:3013024/150‑1 (MQ=255)
                                             tatTGACTGCCGGAACCTCATAACCGCCGCTCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACTTAAATTTTCATATCCTGTCTTTTGTTTTAGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcac                                                                                                        <  1:2123355/150‑1 (MQ=255)
                                              atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcaca                                                                                                       <  2:857156/150‑1 (MQ=255)
                                              atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcacc                                                                                                       >  2:3131013/1‑149 (MQ=255)
                                              atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTAAAAATATAATTAAAATATTTTTTTTTTAATATAGAGATAAAGAtaaaa                                                                                                       >  2:3349841/1‑146 (MQ=255)
                                              atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGTTTTCAAAAAATAATTAAAATATTTTTTTTTTTATATAGAGGTTACGATaaca                                                                                                       >  2:3989950/1‑150 (MQ=255)
                                              atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATGTTCATATCCTGTCGTTTGTGTTAGATTTGAAAAAATAAATAAATAAATTGTTTTTAAATATAGAGATGActaaaaaa                                                                                                       >  2:3349832/1‑143 (MQ=255)
                                               ttGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaa                                                                                                      <  2:4040585/150‑1 (MQ=255)
                                               ttGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaa                                                                                                      <  1:2724442/150‑1 (MQ=255)
                                               ttGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAATATAATGAAATTATTTTTTTTTATATTTCGAGATAACGGTCTCaa                                                                                                      >  2:2984818/1‑150 (MQ=255)
                                                 gAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGTTTTTAAAAATAAAAATAAAATATTTTTTTTTTTATATATAGATAAACATCaaaaaa                                                                                                    >  1:1319328/1‑150 (MQ=255)
                                                 gAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTAAAAATATAAAGAAAATATTTTTTTTTAAATATATAGATTACGATCACaaaa                                                                                                    >  1:2997265/1‑150 (MQ=255)
                                                 gAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTCAAAAAAAAATTAAATTATTTTTTTTTTTATATAGAGGGAAAGATCACaaaa                                                                                                    >  1:627522/1‑150 (MQ=255)
                                                   aTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAc                                                                                                  <  2:849488/150‑1 (MQ=255)
                                                   aTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTTGATTTAAAAAAATAATTAAAATATTTTTTTTTAAAAATAGAGATAAAAAACAaaaaaaa                                                                                                  >  2:3676275/1‑149 (MQ=255)
                                                   aTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTTAAAATATAATTAAAATATTTGTTTTTAAATATAGAGATAACCAACACAAAAAc                                                                                                  >  2:3676725/1‑150 (MQ=255)
                                                       cGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTAAAAATATAATGAAAATATTTTTTTTTTAATATAGAGATAACGATCACAAAAACGACa                                                                                              >  2:1729816/1‑150 (MQ=255)
                                                        ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTTAATATCGAGATAACGATCACAAAAACGACaa                                                                                             >  2:2692681/1‑150 (MQ=255)
                                                          aaCCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTGTTCGAGTTCAAAATATAATGAAATTATTTGTTTTTAAATATAGAGATAACGATCACAAAAACGACAata                                                                                           >  2:3699414/1‑150 (MQ=255)
                                                            ccGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATTAAATTATTTTTTTTTTAATATAGAGATAAAGATAACAAAAACGACAATATg                                                                                         >  2:2348259/1‑150 (MQ=255)
                                                            ccGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTGAAAAAAAAATTAAAATATTTTTTTTTTTAAATAGATATTAAGATCACAAAAAAGACaaaaag                                                                                         >  1:653696/1‑146 (MQ=255)
                                                              gaaTAACCGCCGCGCTGTGCTTTGTTCAGCATATGATTTGTCGATACAACTTACATTTTCATATCCTGTCTTTTTTTTTCTATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGATATAACTATCACAAAAACGACAATATGaa                                                                                       <  2:1404571/148‑1 (MQ=255)
                                                                 tAACCGCCGCGCTGTGCTTTGTTCAGCATCTTCTTTGTCTATACCACTTACTTTTTAATATCCTGTCTTTTTTTTTCTATTTAAAAATATAATGAAATTTTTTTTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAat                                                                                    <  1:1160104/150‑1 (MQ=255)
                                                                   cccgccgCGGTGTGCGTTGTTCAGCATCTGCTTTGTCTATACCACTTAAATTTTAATATCCTGTGGTTTGTTTTCTATTTAAAAATATAATGAATTTATTTTTTTTTAAATATCGAGATAACGATCACAAACACGACAATATGAAAatta                                                                                  <  2:442312/149‑1 (MQ=255)
                                                                   cccgccgCGCTGTGCGTTGTTCCGAATCTGCTTTGTCTATACCACTTAAATTTTCATATCCTGTCGTTTGTTTTCTATTACAAAATATAATGAATTTATTTTTTTTTAAATATGGAGATAACTAACACAAACACGACAATATGAAAatta                                                                                  <  2:442359/149‑1 (MQ=255)
                                                                     cgccgcGCTGTGCGTTTTTCAGCTTCTGCTTTGTCGATACCACGTACATTTTAATATCGTGTCTTTTGTTTTCTATTACAAAATTTAATGAAATTATTTGTTTTTAAATATCGATATAACCATCACAAACACGACAATATGAAAattatt                                                                                <  2:2452934/150‑1 (MQ=255)
                                                                           gcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCgagga                                                                          <  2:2715472/150‑1 (MQ=255)
                                                                             tgtgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTTGATTTAAAAAAAAATTTTAATTATTTGTTTTTTAATATCGAGTTAACCATCTCAAAAAAGCCAAAAAGAAAATTATAAGAGGGGt                                                                        >  2:829631/1‑150 (MQ=255)
                                                                                 cGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTTGGTTTAAAAATAAAATTAAAATATTTTTTTTTTAATATATTTATAAAAATAACAAAAAAAACAAAAATAAAATTATTAGTGGAGTGtaa                                                                    >  1:1736570/1‑150 (MQ=255)
                                                                                        cAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTGAAAATAAAATTAAAATATTTTTTTTTAAATATAGAGAAAAACAAAACAAAAAACACCAAATGAAAAAAATTATAGGAGTGGAACGCaaaa                                                             >  1:2696485/1‑150 (MQ=255)
                                                                                           cATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTTAAAAAAAAATTAAAATAATTTTTTTTAAATAAAGAGATAAACAAAACAAAAAACACAAAATTAAAAATATTCTATGAGTGAAAGGCcaaaaaa                                                          >  1:3161467/1‑150 (MQ=255)
                                                                                            ttCTGCTTTGTCTATACCACTTACTTTTTCATATCCTGTCGTTTTTTTTCGATTTCAAATTATATTGAATTTATTTGTTTTTAAATATCGATATAACGAACACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAc                                                         <  2:61746/149‑1 (MQ=255)
                                                                                            aTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCTTTTTCAAAATATAATGAATTTTTTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAc                                                         <  1:771234/150‑1 (MQ=255)
                                                                                            aTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATCAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAc                                                         <  2:4222364/150‑1 (MQ=255)
                                                                                            aTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTTTATTTAAAAAAATAAAAAAAATATTTTTTTTTAAAAATATAGTTAAAGAAAAAAAAAAAAAAAAAATTAAAAATATTCGAGGAGTGAAATGCaaaaaaaa                                                         >  2:1637045/1‑149 (MQ=255)
                                                                                            aTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTAAAAATATAAATTAAATATTTTTTTTTAAAATACGAGATAAACAAAACAAAAAACAAAAAAATAAAAATAATATTGGCGTGAAAGGCcaaaaaaa                                                         >  2:1637026/1‑149 (MQ=255)
                                                                                            aTCTGCTTTGTCGATACCACATACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTTTTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAc                                                         <  2:381923/150‑1 (MQ=255)
                                                                                             tCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGCTTTAAAAATATAAATAAAATATTTTTTTTTTAAAATACTTATAAAAAAAACAAAAAAAACAAAAATAAAAAAATTAGAGGAGTGTAAGGCaaaaaaaac                                                        >  2:459750/1‑148 (MQ=255)
                                                                                              cTGCTTTGTCCAAACAACTAACATTTTAATATCCTTTCTTTTTTTTTCGATTAAAAAATAAAAAGAAATTATTTTTTTTAAAATATCGAGAAAACGATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACgg                                                       <  2:478498/150‑1 (MQ=255)
                                                                                              cTGATTTTTCTATACCAATTAAATTTACATATCCTTTCGTTTTTTTTCTATTTAAAAATATAATTAATTTATTTGTTTTTAAATATCTAGAAAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACgg                                                       <  2:478563/150‑1 (MQ=255)
                                                                                                      tCGAAACCACGTAAATTTTCATATCCTGTCGTTTTTTTTCTATTTCAAAATATAATCAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGa                                               <  2:3842212/150‑1 (MQ=255)
                                                                                                              ccttACTTTTTCATTGCGTGTGGTTTTTTTTCTATTAAAAAATATAATGAATTTTTTTTTTTTTAAATATCGAGATACCGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAg                                       <  2:825411/147‑1 (MQ=255)
                                                                                                                  aCATTTTCATATCCTGTCGTTTGTTTTGGATTTAAAAAAAAAATTAAAATATTTTTTTTTTAATATAGAGATAAAGATAACAAAAACGCAAATATGAAAATTATTCGAGGCGTGAAAGGCCAAAAAAAGGCCCCCCGCTTGGGCAGccct                                   >  2:3038998/1‑148 (MQ=255)
                                                                                                                  aCATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATCAAATTATTTGTTTTTAAATATCGAGATAACGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGt                                   <  2:201144/150‑1 (MQ=255)
                                                                                                                  aCATTTTCATATCCTGTCGTGTTTTTTCGATTTCAAAATATAATGAAAATATTTTTTTTTAAATATAGAGATAACCATCACAAAAAACCAAATATGAAAATTATTCGCGGAGTGTAAGGCAAAAAAACGGCCAACCCGGAGGGCAGCCGt                                   >  2:3038974/1‑150 (MQ=255)
                                                                                                                   atttttAATATCCTGTCGTTTTTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                  <  2:3350215/148‑1 (MQ=255)
                                                                                                                       ttCATATCGTTTTTTTTGTTTTCTTTTTAAAAATATAATTAATTTTTTTTTTTTTAAATATCTATATAACTAACAAAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaa                              <  2:3585454/150‑1 (MQ=255)
                                                                                                                             tCCTGTCGTTTTTTTTAGATTTCAAAAAATAATGAAATTATTTTTTTTTAAATATCGAGAAAACGAAAACAAAAAAGACAAAATTTAAATTATTCGAGGAGTGAAAGGCAAAAAAAAGGCCACCCGATAGGGAAGCCCTAGCAAAATgcg                        >  2:701024/1‑150 (MQ=255)
                                                                                                                                tgtgtttttttttCTATTTAAAAATATATTTAATTTTTTTTTTTTTAATTTTCTATATAACTAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgt                     <  2:1203182/145‑1 (MQ=255)
                                                                                                                                tGTCGTTTTTTTTCTATTTAAAAATATAATTAAATTATTTTTTTTTAAATATCGAGATAACTATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgt                     <  2:3598053/150‑1 (MQ=255)
                                                                                                                                    gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGACGCCGTAGCAAAGTGCGCGTGttt                 >  2:4018226/1‑150 (MQ=255)
                                                                                                                                    gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt                 >  1:4137930/1‑150 (MQ=255)
                                                                                                                                     tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                >  2:4647953/1‑150 (MQ=255)
                                                                                                                                         ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg            >  1:2211302/1‑150 (MQ=255)
                                                                                                                                               aTTTCAAAATTTAATGAAATTTTTTGTTTTTAAATATCGAGATAACGTTCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt      <  2:1732940/150‑1 (MQ=255)
                                                                                                                                                tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg     >  1:4538798/1‑150 (MQ=255)
                                                                                                                                                    cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg  >  2:4243752/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGG  >  CP009273/2170566‑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: