Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,708 G→T 7.5% intergenic (‑25/+283) 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,7080GT 7.5% 632.0 / 5.3 222intergenic (‑25/+283)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base G (89/105);  new base T (5/23);  total (94/128)
Fisher's exact test for biased strand distribution p-value = 6.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

TGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGG  >  CP009273/2170559‑2170856
                                                                                                                                                     |                                                                                                                                                    
 gCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGa                                                                                                                                                     <  1:1736427/150‑1 (MQ=255)
   aGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAtt                                                                                                                                                   <  1:4472993/150‑1 (MQ=255)
   aGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAtt                                                                                                                                                   >  2:803131/1‑150 (MQ=255)
   aGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAtt                                                                                                                                                   <  2:1051975/150‑1 (MQ=255)
   aGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGATGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAtt                                                                                                                                                   <  2:5111008/150‑1 (MQ=255)
      ttGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCa                                                                                                                                                <  1:1837872/150‑1 (MQ=255)
         gCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaaa                                                                                                                                             <  2:1620602/150‑1 (MQ=255)
         gCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaaa                                                                                                                                             <  2:1802390/150‑1 (MQ=255)
         gCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaaa                                                                                                                                             <  1:2508406/150‑1 (MQ=255)
          cAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAat                                                                                                                                            <  1:5316019/150‑1 (MQ=255)
           aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                           <  1:4609404/150‑1 (MQ=255)
           aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGTTTTGAAAata                                                                                                                                           >  1:1965427/1‑150 (MQ=255)
              ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                        <  2:1877147/150‑1 (MQ=255)
              ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                        <  2:4609871/150‑1 (MQ=255)
              ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATaa                                                                                                                                        >  1:2256059/1‑150 (MQ=255)
               gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                       <  2:171184/150‑1 (MQ=255)
                tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATg                                                                                                                                      <  2:4396549/150‑1 (MQ=255)
                   cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGGTTTTAAAAaataatttaa                                                                                                                                   >  1:95160/1‑146 (MQ=255)
                    tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                  <  2:311005/150‑1 (MQ=255)
                    tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                  <  2:293103/150‑1 (MQ=255)
                    tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                  <  1:5114587/150‑1 (MQ=255)
                     accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATAAAATGaaaat                                                                                                                                 >  2:4694890/1‑148 (MQ=255)
                        accacTTGCATCGTTTCGAGGTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAAATAAattatt                                                                                                                              >  2:549064/1‑150 (MQ=255)
                        accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTAAattatt                                                                                                                              >  2:1252666/1‑150 (MQ=255)
                         ccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTAttt                                                                                                                             <  2:2553854/150‑1 (MQ=255)
                          cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                            <  1:1181620/150‑1 (MQ=255)
                          cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTGAAAATATAATGAAATTATTTg                                                                                                                            >  2:2807685/1‑150 (MQ=255)
                            cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtt                                                                                                                          <  2:2641319/150‑1 (MQ=255)
                            cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTAAAAAAATAATTTAAATAtttttt                                                                                                                          >  1:1369848/1‑150 (MQ=255)
                             ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttt                                                                                                                         <  2:4309435/150‑1 (MQ=255)
                             ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttt                                                                                                                         <  1:2916438/150‑1 (MQ=255)
                              tGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtttt                                                                                                                        <  1:2083450/150‑1 (MQ=255)
                              tGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtttt                                                                                                                        <  2:1670192/150‑1 (MQ=255)
                              tGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtttt                                                                                                                        <  1:3436784/150‑1 (MQ=255)
                              tGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtttt                                                                                                                        <  1:2083405/150‑1 (MQ=255)
                                cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTAAAAAAATAAATAAAATATttttttttt                                                                                                                      >  1:1855397/1‑149 (MQ=255)
                                 aTCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAAAGAAAAAATttttttttta                                                                                                                     >  2:5227507/1‑148 (MQ=255)
                                  tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaaa                                                                                                                    <  1:1999607/150‑1 (MQ=255)
                                  tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaaa                                                                                                                    <  1:1904134/150‑1 (MQ=255)
                                    gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAata                                                                                                                  <  2:3869312/150‑1 (MQ=255)
                                    gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGTTTTCAAAATATAATTTAATTATTTTTTTTTAAatt                                                                                                                  >  1:417703/1‑149 (MQ=255)
                                    gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTATAGTTAAAAAAAAAATTAAAATATTTTTTTTTaaaaa                                                                                                                  >  1:5375816/1‑148 (MQ=255)
                                    gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGGTTTGGATTTCAAAAAAAAAATAAAAAATTTTTTTTTTAata                                                                                                                  >  1:1472005/1‑150 (MQ=255)
                                     tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                 <  2:2698901/150‑1 (MQ=255)
                                     tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                 <  2:3540870/150‑1 (MQ=255)
                                      ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATc                                                                                                                <  1:3510800/150‑1 (MQ=255)
                                        cGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTGTTTGGATTTAAAAAAAAAATTAAAATATTTTTTTTTTAATATAga                                                                                                              >  1:254174/1‑150 (MQ=255)
                                         gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgag                                                                                                             >  2:3903056/1‑150 (MQ=255)
                                         gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgag                                                                                                             <  1:3694327/150‑1 (MQ=255)
                                          agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                            >  2:5534083/1‑150 (MQ=255)
                                          agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGTATTTAAAAAAAAAAATAAAAAAATTTTTTTTAAATAAAgaga                                                                                                            >  1:2626910/1‑150 (MQ=255)
                                           gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATGAAATTATTTTTTTTTTAatatatatta                                                                                                           >  1:2623540/1‑144 (MQ=255)
                                               gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTAAAAAAATAATGAAATTATTTTTTTTTTAAAATATAGATaaac                                                                                                       >  1:2269375/1‑148 (MQ=255)
                                                 gAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTAAAAAAATAAATTAATTATTTTTTTTTTAATAAATAGATaaacat                                                                                                     >  1:4094468/1‑146 (MQ=255)
                                                  aaTATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATACTGAAATTATTTTTTTTTTAATATAGAGATAACGATc                                                                                                    >  2:5263705/1‑150 (MQ=255)
                                                     atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAAAATAATTAAAATATTTTTTTTTAAATATAGAGATAACGATaaca                                                                                                 >  2:3161111/1‑150 (MQ=255)
                                                         aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGGGTTTTAGAATTCAAAAAAAAAAGAAAAAATATTTTTTTAAATAAAGAGATAACCATAACaaaaa                                                                                             >  1:2888698/1‑150 (MQ=255)
                                                             ccGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTAAAAATATAAATAAAATATTTTTTTTTATAAATAGAGAGAAAAATAACAAAAAAGAc                                                                                         >  1:3410764/1‑150 (MQ=255)
                                                             ccGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTAAAAATATAAAGAAAATATTTTTTTTTTAAAATAGAGATAAAGAACACaaaaaacac                                                                                         >  2:1553060/1‑146 (MQ=255)
                                                               ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTGTTTAGATTTAAAAAAATAAATAAAATATTTTTTTTTTAATATAGAGTTAAAGATAACAAAAAAGACaa                                                                                       >  1:1777988/1‑150 (MQ=255)
                                                                    cGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTAAAAAAATAAATAAAATAATTTTTTTTTAATATAGAGATAAAGATAACAAAAAcaacaaaaata                                                                                  >  1:1619132/1‑145 (MQ=255)
                                                                       aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACTTAAATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGaaaa                                                                               <  2:2239353/150‑1 (MQ=255)
                                                                       aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAAAGAAAATATTTGTTTTTAAATAAAGAGAAAACGAAAACAAAAACCACAATATGaaaa                                                                               >  2:4022110/1‑150 (MQ=255)
                                                                         acccgccgCGGTGTGCTTTTTTCAGCATCTGCTTTGTCTATACCACTTAAATTTTAATATCCTTTCTTTTTTTTTCTTTTTCAAAATATAATTAATTTTTTTTTTTTTAAATATCGAGATAACTATCAAAAAAACGACAATATGAAAatt                                                                             <  1:3002752/148‑1 (MQ=255)
                                                                         aaCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAAAATAAATAAAATATTTTTTTTTAAATATAGAGATAAAGATAAAAAAAAAGACAATATTAAAatt                                                                             >  2:2850246/1‑150 (MQ=255)
                                                                          aCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTGAAAAAAAAATTAAAATATTTTTTTTTAAAAAAATATATAAAGAAAACAAAAAAAACAAAATTAAAaata                                                                            >  2:1970058/1‑150 (MQ=255)
                                                                            cgccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTAAAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAattatt                                                                          <  2:3823816/150‑1 (MQ=255)
                                                                             gccgcGCTGTGAGTTGTTCAGCATGTGCTTTGTCTATACCACTTACATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATTAAATTTTTTTTTTTTAAATATCGAGATAACGATCAAAAAAACGACAATATGAAAATTATTc                                                                         <  2:4407338/150‑1 (MQ=255)
                                                                              ccgcGCTGTGCGTGGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATTGCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCg                                                                        <  2:5180664/150‑1 (MQ=255)
                                                                               cgcgCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAATAATTAAATAATTTTTTTTTAAATAAAGAGAGAAACAAAACAAAAAAAACAAAAAGAAAaataatagc                                                                       >  2:948574/1‑145 (MQ=255)
                                                                                 cgcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAAATATTTTTTTTTAAATATAGAGATAAACATCACAAAAAAGACAAAAAGAAAAATAATCgagg                                                                     >  2:192431/1‑150 (MQ=255)
                                                                                   cTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGAATTAAAAAAAAAATGTAAATATTTTTTTTTTAAAAAAGATAAAACCATAACAAAAACAAAAAAAATAAAAATATTCgatgag                                                                   >  2:17348/1‑150 (MQ=255)
                                                                                   cTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAATATAATTAAAATATTTGTTTTTAAATATAGAGATAACGATAAAAAAAAAGAAAAAATGAAAATTAATCgagggg                                                                   >  2:2439156/1‑148 (MQ=255)
                                                                                      ggCGTTGTTCAGCATCTGCTTTGCCGATACCACGTACATTTTCAAATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGa                                                                <  2:4533391/149‑1 (MQ=255)
                                                                                       gCGTTGTTCCGCACCGGATTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATACCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGaa                                                               <  2:3515935/150‑1 (MQ=255)
                                                                                           tgttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAAAATAATGAAATTATTTTTTTTTAAATATAGAGATAACGATCACAAAAAAGACAAAATGAAAATTATTCGAGGAGTGAAAGGc                                                           >  2:4891340/1‑150 (MQ=255)
                                                                                            gttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATAGAGATAACGATCACAAAAAAGACAATATTAAAATTATTAGAGGAGTGAAAGGCa                                                          >  2:866031/1‑150 (MQ=255)
                                                                                                 gCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaaaaaa                                                     <  2:4900168/150‑1 (MQ=255)
                                                                                                     cTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTTAATATATAGATAAAGATCACAAAAAAGACCATATTTAAATTATTAGAGGTGTTAAAGGCaaaaaaaagt                                                 >  2:763989/1‑147 (MQ=255)
                                                                                                           ggTCTATACAACGAACATTTTCTTATCCTGTCGTTTTTTTTCTTTTACAAAATATAATTAATTTATTTGTTTTTAAATATCGAGATAACTATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTccc                                           <  2:2003220/149‑1 (MQ=255)
                                                                                                                aTACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTAAAAAAATAATTAAATTATTTTTTTTTAATTATAGAGATAACGATCACAAAAACCACAAAATGAAAATTATTCGAGGAGGGAAAGGCCAAAAAAAAGCCTccccctag                                      >  2:2136445/1‑145 (MQ=255)
                                                                                                                  aCCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAggg                                    <  2:3609305/150‑1 (MQ=255)
                                                                                                                     ccctaaaTTTTAATATCCTGTGTTTTTTTTCCTATTACAAAATATAATTAATTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAg                                 <  2:956429/144‑1 (MQ=255)
                                                                                                                       ttaatttttAATATCATGTCGTTTTTTTTGTATTACAAAATATATTGAATTTATTTTTTTTTAAATATCTATATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                               <  1:4443811/145‑1 (MQ=255)
                                                                                                                       gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                               <  2:4534478/150‑1 (MQ=255)
                                                                                                                          aaTTTTCAAACCCTGTCGTTTTTTTACGAATAAAAAATATAATTAAATTATTTTTTTTTAAATATCGATATAACGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                            <  1:4522585/149‑1 (MQ=255)
                                                                                                                             tttCATATCCGGTCGTTTGTTTTCTATTTCAAAATATAATGAATTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCa                         <  2:4944734/150‑1 (MQ=255)
                                                                                                                              ttCATATCCTGTCGTTTGTTTTTGATTTCAAAATAAAAATAAAATATTTTTTTTTAAATAAAGAGATAACGATAAAAAAAAACACAAAATTAAAAATATTTGGGGTGTTAAAGGGCAAAAAAAGGCCGCCCCAAAAGGGAGACCTAGCaa                        >  2:5196285/1‑150 (MQ=255)
                                                                                                                                     ccTGTCGTTTGTTTTCGATTTCAAAATATAATGAAAATATTTGTTTTTAAATATCGAGATAACGATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGTAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTATCAAAGTgcgc                 >  2:1370904/1‑150 (MQ=255)
                                                                                                                                          cGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt            >  2:2293306/1‑150 (MQ=255)
                                                                                                                                           gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt           >  2:3669361/1‑150 (MQ=255)
                                                                                                                                            tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt          >  1:207584/1‑150 (MQ=255)
                                                                                                                                             ttgtttTCGTTTTCAAAATATAATGAAATTTTTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt         <  1:5534083/150‑1 (MQ=255)
                                                                                                                                             ttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt         >  1:1507961/1‑150 (MQ=255)
                                                                                                                                              tgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAATCGGCCTCCCGATCGGGAAGCCGTAGCAACGTGCGCGTGTTTTTa        >  2:1474376/1‑150 (MQ=255)
                                                                                                                                              tgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa        >  2:5481987/1‑150 (MQ=255)
                                                                                                                                                ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg      <  2:4455534/150‑1 (MQ=255)
                                                                                                                                                  ttCTATTTAAAAATATAATTAAATTATTTTTTTTTAAATATCGATATAACGAACACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc    <  1:5132615/150‑1 (MQ=255)
                                                                                                                                                  ttCTATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc    <  2:3587000/150‑1 (MQ=255)
                                                                                                                                                    cGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCgg  >  1:657973/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGG  >  CP009273/2170559‑2170856

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: