Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,707 C→A 16.1% intergenic (‑24/+284) 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,7070CA16.1% 416.3 / 7.6 199intergenic (‑24/+284)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base C (31/107);  new base A (27/4);  total (88/111)
Fisher's exact test for biased strand distribution p-value = 2.18e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

ATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCG  >  CP009273/2170558‑2170855
                                                                                                                                                      |                                                                                                                                                     
aTGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTc                                                                                                                                                       <  1:3931481/150‑1 (MQ=255)
 tGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCg                                                                                                                                                      <  1:68567/150‑1 (MQ=255)
 tGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCg                                                                                                                                                      <  1:2495455/150‑1 (MQ=255)
  gCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGa                                                                                                                                                     <  1:587036/150‑1 (MQ=255)
    aGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAtt                                                                                                                                                   <  2:2034091/150‑1 (MQ=255)
    aGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAtt                                                                                                                                                   <  1:1425126/150‑1 (MQ=255)
        tGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaa                                                                                                                                               <  1:3121849/150‑1 (MQ=255)
         ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTaaaa                                                                                                                                              >  1:4300432/1‑150 (MQ=255)
          gCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaaa                                                                                                                                             <  1:3815552/150‑1 (MQ=255)
          gCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaaa                                                                                                                                             <  2:2160898/150‑1 (MQ=255)
           cAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAat                                                                                                                                            <  1:2077791/150‑1 (MQ=255)
           cAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAat                                                                                                                                            <  1:3429503/150‑1 (MQ=255)
            aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                           <  2:70003/150‑1 (MQ=255)
             gCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAatat                                                                                                                                          <  1:3779200/150‑1 (MQ=255)
             gCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAatat                                                                                                                                          <  2:1418018/150‑1 (MQ=255)
              cGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAtata                                                                                                                                         <  2:3208805/150‑1 (MQ=255)
               ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                        <  1:528340/150‑1 (MQ=255)
                gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                       <  2:4275174/150‑1 (MQ=255)
                 tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATg                                                                                                                                      <  2:494736/150‑1 (MQ=255)
                 tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACAGTTTCATATCCTGTCGTGTGGTTTTGATTTCAAAATATAATg                                                                                                                                      >  1:892431/1‑150 (MQ=255)
                   tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                    <  1:4232332/150‑1 (MQ=255)
                    cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                   <  1:392124/150‑1 (MQ=255)
                      accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATTAAatt                                                                                                                                 >  2:1697004/1‑150 (MQ=255)
                         accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                              <  1:3021804/150‑1 (MQ=255)
                           cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                            <  2:1909909/150‑1 (MQ=255)
                           cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                            <  1:1463304/150‑1 (MQ=255)
                           cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                            <  2:239962/150‑1 (MQ=255)
                             cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTTGAATTTCAAAAATTATTAAAAtttttttt                                                                                                                           >  1:3720255/1‑148 (MQ=255)
                              ttGCATCGTTTAGAGATTGTTAATATTGAATGCCTGAACCTCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttt                                                                                                                         <  2:1162512/150‑1 (MQ=255)
                              ttGAATCGTTTCGAGATTGTGAATATTGAATTCCGTAACCGCATAACCGCCGAGCTGTTCGTTGTTCAGCATCTGCTTTTTCGATACCACGTTCATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATTAAATTATTTGttt                                                                                                                         <  1:2263480/150‑1 (MQ=255)
                               tGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTATATTTAAAAAAATAATAAAAAAAtttttttt                                                                                                                        >  1:3319944/1‑150 (MQ=255)
                                gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTCAAAAAATAAATAAATTATtttttttt                                                                                                                       >  2:1119086/1‑150 (MQ=255)
                                     gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATGAAATTATTTGTTTTTAAata                                                                                                                  >  2:1772703/1‑150 (MQ=255)
                                     gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCTATTTCAAAATATAATGAAATTTTTTTTTTTTTAata                                                                                                                  >  1:1724997/1‑150 (MQ=255)
                                      tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                 <  2:1706476/150‑1 (MQ=255)
                                       ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATc                                                                                                                <  1:2308078/150‑1 (MQ=255)
                                       gtCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATc                                                                                                                <  1:1556197/149‑1 (MQ=255)
                                         cGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTAAAAATATAATTTAATTAATTTTTTTTAAatataca                                                                                                              >  2:3075640/1‑147 (MQ=255)
                                          gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGCTTTCAAAATATAATTAAATTATTTTTTTTTAAATTTagag                                                                                                             >  2:2705708/1‑150 (MQ=255)
                                           agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                           <  2:3317504/150‑1 (MQ=255)
                                           agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAAAATAATGAAATTATTTTTTTTTAAATATAgaga                                                                                                           >  2:4198817/1‑150 (MQ=255)
                                              ttGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGAATTGAAAAAAAAAATAAATTATTTTTTTTTAAATATAGAGaaaa                                                                                                        >  2:3399254/1‑150 (MQ=255)
                                               tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAc                                                                                                       <  2:2486686/150‑1 (MQ=255)
                                               tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAc                                                                                                       <  2:3720598/150‑1 (MQ=255)
                                                gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTTGATTTAAAAAAATAAATTAATTAATTTTTTTTAATTATAGAGATaaag                                                                                                      >  1:3135745/1‑148 (MQ=255)
                                                gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTAAAATATTTTTTTTTAAATATAGAGATaaag                                                                                                      >  2:1403288/1‑148 (MQ=255)
                                                gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACg                                                                                                      <  1:4165932/150‑1 (MQ=255)
                                                 tgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGa                                                                                                     <  1:3405943/150‑1 (MQ=255)
                                                  gAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGTTTTGAAAAAATAAATAAAATATTTTTTTTTTTATATAGTGGTAACGAt                                                                                                    >  1:3716891/1‑150 (MQ=255)
                                                   aaTATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATc                                                                                                   <  2:3172288/150‑1 (MQ=255)
                                                   aaTATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATc                                                                                                   <  2:2250553/150‑1 (MQ=255)
                                                    atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATca                                                                                                  <  2:3475352/150‑1 (MQ=255)
                                                        tGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAAATATTTGTTTTTAAATATAGAGAGAACGATCACaaa                                                                                              >  2:2821439/1‑150 (MQ=255)
                                                            tGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACg                                                                                          <  2:3236417/150‑1 (MQ=255)
                                                            tGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAATATAATTTAATTATTTTTTTTTTAATATAGAGATAAATATaaaaaaaaag                                                                                          >  2:639737/1‑148 (MQ=255)
                                                             gCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGa                                                                                         <  2:2101935/150‑1 (MQ=255)
                                                                 gAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGGTTTAAAAAAATAAATAAAATTTTTTTTTTTTAATATATAGATAAAGATaaaaaaaaacaccaa                                                                                     >  1:3320505/1‑143 (MQ=255)
                                                                    ccGCATCACCGCCGGGCTGTGGGTTTTTCTGCTTGTGCTTTGTCTATACAACTTAAATTTTAATATCCTGTCTTTTGTTTTCTTTTTAAAAATATAATTAATTTTTTTTTTTTTAAATATCGATTTAACTATCACAAAAACGACAATATg                                                                                  <  2:4105296/150‑1 (MQ=255)
                                                                      gCATAACCGCGGCGCTGTGCTTTGTTCAGCTTCTGGTTTGTCTATACCACTTACTTTTTCATATCCTGTCGTTTTTTTTCGTTTTCAAAATATATTGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGaa                                                                                <  2:3872800/150‑1 (MQ=255)
                                                                      gCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACAGTTTCATATCCTGTCGTGTGTTTTCGATTTAAAAATATAATGAAAATATTTTTTTTTAAATATATAGTATAAGATACCAAAAACAACAAAATGaa                                                                                >  2:2607395/1‑150 (MQ=255)
                                                                         tacccgcggCGCGGTGCTTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTGGTTTGTTTTCTATTAAAAAATATAATTAATTTTTTTTTTTTTAAATATCTTGATAACGATAAAAAACACGACAATATGAAAat                                                                             <  1:4208005/147‑1 (MQ=255)
                                                                         tAACCGCCGGGGTGTGCGTTGTTAAGCTTCTGGTTTGTCTATACCACGTAAATTTTCATATCATGTCGTTTGTTTTCTATTTAAAATTATAATTAAATTTTTTTTTTTTAAATTTCTATATAACGATCACAAAAACGACAATATGAAAat                                                                             <  2:4041539/150‑1 (MQ=255)
                                                                         tAACCGCCGCGCTGTGGGTTTTTCAGCATTTGCTTTGTCGATACCCCTTAAATTTTCATATCCTGTGTTTTGTTTTCGATTAAAAAATATATTTAATTTTTTTTTTTTTAAATATCTATATAACGATCACAAAAACGACAATATGAAAat                                                                             <  2:2181803/150‑1 (MQ=255)
                                                                          aaCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGGTTTAAAAAAAAAATTTAATTATTTTTTTTTAAATATAGAGGTTACGAACACAAAAAAGACAAAATGagaaat                                                                            >  1:3572728/1‑148 (MQ=255)
                                                                          aaCCGCCGCGCTGTGCGTTGTCCAGCATCTGCTTTGTCTATACCACGTACTTTTAAATATCCTTTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAatt                                                                            <  2:4335389/150‑1 (MQ=255)
                                                                           aCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGGTTTGGATTTAAAAATATAATTAAAATATTTTTTTTTTAAAATAGAGATGACGATCACAAAAACAACAAAATGAAAaata                                                                           >  1:2089058/1‑150 (MQ=255)
                                                                              gccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTc                                                                        <  2:1423091/150‑1 (MQ=255)
                                                                                    cTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTAAAAATATAATTAAAATATTTTTTTTTTAAAATAGCGATAAAGATAAAAAAAAACAAAAAATTAAAAATATTAgaggag                                                                  >  1:2512430/1‑150 (MQ=255)
                                                                                    cTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGGGTTTTTTTTTTTTAAATATAATTTAAATATTTTTTTTTTTATATTAAGATAAAGATCAAAAAAAAAACAAAAATAAAaatattattggtg                                                                  >  1:2944850/1‑143 (MQ=255)
                                                                                     tgtgCGTTGTACACCATCTGCTTTTTCCACACCACTAACATTTTCATACCCTGTCTTTTTTTTTCTATTAAAAAATATAATTAATTTATTTTTTTTTAAATATCGAGAAAACTAACACAAACACGACAATATGAAAATTATTCGAGGAGt                                                                 <  1:1222719/150‑1 (MQ=255)
                                                                                      gtgCGTTGTTCAGCATCTGCTTGGTCGATACCACGTACATTTTCATATCCGGTCGTTTGTTTTCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTg                                                                <  2:1351616/150‑1 (MQ=255)
                                                                                       tgCGTTGTTCAGCATCTGCTTGGTCGATACCACGTACATTTTCATATCCTGGCGTTTTTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGa                                                               <  2:1010578/150‑1 (MQ=255)
                                                                                             gttCAGCATCTGCTTTGTCTATACCCCGTAAATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAATTTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCa                                                         <  2:3556273/150‑1 (MQ=255)
                                                                                             gttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTAAAAAAAAAAATAAAATATTTTTTTTTAAATATCGAGATAACGATAACAAAAAAGAAAAAAAGAAAAATATTAGAGGAGTTGAAGGCa                                                         >  2:1692722/1‑150 (MQ=255)
                                                                                              ttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTAAAAAAAAAATTAAAAAAATTTTTTTTTAATATAGAGAAAAACAAAACAAAAAAAACAAAAAGAAAATTAATAGAGGACTGTAAGGCaa                                                        >  1:1366525/1‑150 (MQ=255)
                                                                                               tCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAATATAATGAAATTATTTTTTTTTAAATATAGAGATAACGAAAACAAAAAAGACAAAATTAAAATTATTAGTGGGTTGAAAAGCaaa                                                       >  2:1595012/1‑150 (MQ=255)
                                                                                                  tcATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCTTTTGTTTTCTATTACAAAATATAATTAAATTATTTTTTTTTAAATATCGAGATAACGATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaaaaaa                                                    <  2:2997888/149‑1 (MQ=255)
                                                                                                               cGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTAAATTATTTTTTTTTTAATATCGAGATAAAGATAACAAAAAAGACAATATTAAAATTATTGGAGGAGTGTAAGGGAAAAAAAAGGCCTCCCGaa                                       >  2:3568260/1‑149 (MQ=255)
                                                                                                               cGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTCAAAATATAATGAAAATATTTTTTTTTAAATATCGAGTTAACGATCACAAAAACGACACTATGAAAATTATTAGAGGGGTGAAAGGCCAAAAAACCGCCACCCGag                                       >  2:1933523/1‑149 (MQ=255)
                                                                                                                   aCCACGTACATTTTCATATCCTGTCGTGTGTTTTTGATTTGAAAATATAATTAAAATATTTTTTTTTTAATTTAGAGTTTACGATAACCAAAACGACAATATTAAAATTTTTAGGGGGGGGTAATGGAAAAAAAAGGCCTGCCGCTTggg                                   >  2:4163751/1‑150 (MQ=255)
                                                                                                                   aCCACGTACATTTTCATATCCTGTCGTGTGTTTTGGTTTTGAAAAAATAATTTAAATATTTTTTTTTTAATATGGAGATTAACATGACAAAAAAGACCAAATTAAAATTATTTGGGGGGGGTAATGCCAAAAAAAGGCCGCCCGGTTggg                                   >  2:571123/1‑150 (MQ=255)
                                                                                                                     cACTTACTTTTTCTTATCTTGTCTTTTTTTTTCTTTTTAAAAATATAATTAATTTTTTTTTTTTTAAATATCTTGATACCGAACAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGaa                                 <  1:1015078/150‑1 (MQ=255)
                                                                                                                     cACGTAAATTTTCATATCCTGTCTTTTTTTTTATATTTAAAAATATAATTAAATTATTTTTTTTTAAATATCTATATAACTAACAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGaa                                 <  1:320681/150‑1 (MQ=255)
                                                                                                                      aCGTACATTTTCATATCCTGTCGTGTGTTTTAGCTTTAAAAAAATAATGAAATTATTTGTTTTTTAATATCGAGATTACGATAACAAAAACGACAATATGACAATTATTGGAGGAGTGGAAGGCACAAAAAACGCCACCCGCTAGGGGAg                                >  2:2767416/1‑150 (MQ=255)
                                                                                                                       cGTACATTTTCATATCCTGTCGTGTTGTTTTGATTTAAAAATAAAATGAAAATATTTTTTTTTTAAAATAGAGATAACGATAAAAAAAAACACCAAAATAAAATTATTAGGGGTGTGGAAGGCAAAAAAAACGCCACCCCAGAGGGGAGc                               >  2:780183/1‑150 (MQ=255)
                                                                                                                        gTACATTTTAAGACCCTGTCTTTTGTTTTCGATTCAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                              <  2:2212713/150‑1 (MQ=255)
                                                                                                                           cATTTTCATATCCTGTCGTGTGTTTTGGATTTACAAAAAAAATTAAAATATTTTTTTTTTAATTTAGAGATAACGATCAAAAAAAATACAATATTAAAATTATTAGTGGGGTTTAAAGGAAAAAAAAAGCCGCCCCCTAGGGGCGCCgga                           >  2:1187951/1‑148 (MQ=255)
                                                                                                                            tttttCATATCCTGGGGTTTTTTTGCGATTTCAAAATATAATTAATTTTTTTTTTTTTAAATTTCGATATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                          <  1:126072/149‑1 (MQ=255)
                                                                                                                               ttCATACCCTTTCTTTTTTTTTCTTTTTAAAAATATAATTAAATTTTTTTTTTTTAAATATCTATATAACTAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaa                       <  1:272193/150‑1 (MQ=255)
                                                                                                                                         gTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGAAAATATGAAAATTATTCGAGTTGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg             >  2:3891515/1‑150 (MQ=255)
                                                                                                                                           tgtttgtttTCTATTTAAAAATATAATTAAATTATTTGTTTTTAAATATCGATATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt           <  2:4038696/149‑1 (MQ=255)
                                                                                                                                           cGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt           >  1:3505448/1‑150 (MQ=255)
                                                                                                                                           cGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt           >  1:2076723/1‑150 (MQ=255)
                                                                                                                                            gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt          >  1:945345/1‑150 (MQ=255)
                                                                                                                                             tttgtttTCGATATCAAAAAATAATGAAAATATTTGTTTTTAAATAACGAGATAACGATCACAAAAACGACAATATGAAAATAATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCCATATGGAAGCCGTAGCAAAGAGCGAGTGttat         >  1:3239958/1‑148 (MQ=255)
                                                                                                                                               ttgtttTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt        <  2:3629053/150‑1 (MQ=255)
                                                                                                                                                tgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa       >  2:3130310/1‑150 (MQ=255)
                                                                                                                                                 gtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt      >  2:1259770/1‑150 (MQ=255)
                                                                                                                                                    ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAGATTATTGGTGGAGTGTAAGGCAATAAAACGGCGTCGGGATGGGGTAGCCGTAGCTAGGTGCGGGGGtttgtttggg   >  1:2194285/1‑143 (MQ=255)
                                                                                                                                                    ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAACGTGCGCGTGTTTTTATGcc   >  1:1704703/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCG  >  CP009273/2170558‑2170855

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: