Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,234 C→A 15.7% 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
*NC_0009132,177,2340CA15.7% 502.7 / 16.6 199intergenic (‑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 (76/81);  new base A (29/10);  total (108/91)
Fisher's exact test for biased strand distribution p-value = 4.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

TTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGG  >  NC_000913/2177086‑2177383
                                                                                                                                                    |                                                                                                                                                       
ttGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCa                                                                                                                                                        <  1:666967/150‑1 (MQ=255)
  ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCGGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaa                                                                                                                                                      <  2:5019415/150‑1 (MQ=255)
  ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaa                                                                                                                                                      <  1:4011214/150‑1 (MQ=255)
  ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaa                                                                                                                                                      <  2:2170905/150‑1 (MQ=255)
  ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaa                                                                                                                                                      <  1:2509399/150‑1 (MQ=255)
   gCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCGTTGTCGATACCACGTACAGTTTCATATCCTGTCGTGTGTTTTCGATTTCaaaa                                                                                                                                                     >  2:2825694/1‑150 (MQ=255)
    cAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAat                                                                                                                                                    <  1:1850207/150‑1 (MQ=255)
     aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                                   <  1:2655987/150‑1 (MQ=255)
     aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                                   <  2:2269144/150‑1 (MQ=255)
      gCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAatat                                                                                                                                                  >  1:4678859/1‑150 (MQ=255)
         gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                               <  2:3264775/150‑1 (MQ=255)
          tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATg                                                                                                                                              >  1:5004973/1‑150 (MQ=255)
           ttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGa                                                                                                                                             <  2:1636568/150‑1 (MQ=255)
            tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                            >  1:200670/1‑150 (MQ=255)
             cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                           >  1:2138663/1‑150 (MQ=255)
               accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAatt                                                                                                                                         >  2:3838017/1‑150 (MQ=255)
               accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTaaaat                                                                                                                                         >  2:3966706/1‑148 (MQ=255)
               accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCAGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGAATTCAAAAAATAATGAAata                                                                                                                                         >  1:3927291/1‑149 (MQ=255)
                ccaccaCTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAatta                                                                                                                                        >  2:3431114/1‑150 (MQ=255)
                 caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattat                                                                                                                                       <  2:2702743/150‑1 (MQ=255)
                  accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                                      <  1:2688574/150‑1 (MQ=255)
                   ccacTTGCATCGTTTCGATATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTAttt                                                                                                                                     <  1:2716603/150‑1 (MQ=255)
                   ccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTAttt                                                                                                                                     <  1:4814413/150‑1 (MQ=255)
                   ccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTAttt                                                                                                                                     <  1:4692515/150‑1 (MQ=255)
                   ccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTAttt                                                                                                                                     >  2:3303132/1‑150 (MQ=255)
                         gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATGAAATTATtttttttt                                                                                                                               >  2:2518947/1‑150 (MQ=255)
                           aTCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaa                                                                                                                             <  2:4485751/150‑1 (MQ=255)
                           aTCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaa                                                                                                                             <  2:1346378/150‑1 (MQ=255)
                            tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATGAAAATATTTTTTTTTaaa                                                                                                                            >  1:3102051/1‑150 (MQ=255)
                               tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTGAAAAAATAATGAAAATATTTTTTTTTAAatat                                                                                                                         >  1:3073674/1‑150 (MQ=255)
                               tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                         <  1:967374/150‑1 (MQ=255)
                                ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                                >  2:166672/1‑126 (MQ=255)
                                ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                                <  1:166672/126‑1 (MQ=255)
                                ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATc                                                                                                                        <  1:4881436/150‑1 (MQ=255)
                                 tCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCg                                                                                                                       <  2:2605483/150‑1 (MQ=255)
                                  cGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTCAAAATATAATGAAAATATTTTTTTTTAAatataaa                                                                                                                      >  1:4376936/1‑147 (MQ=255)
                                   gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTTTTTTTAAATATagag                                                                                                                     >  1:749473/1‑150 (MQ=255)
                                   gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAATAATTAAAATATTTGTTTTTAAATATagag                                                                                                                     >  2:4330812/1‑150 (MQ=255)
                                    agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                                    <  1:2852940/150‑1 (MQ=255)
                                     gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                                   <  2:1885040/150‑1 (MQ=255)
                                      aTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATGAAATTAATTTTTTTTAAATATCGAGATa                                                                                                                  >  2:4989840/1‑150 (MQ=255)
                                         gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACg                                                                                                               >  2:5032304/1‑150 (MQ=255)
                                         gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAAATATTTTTTTTTAAATATCGAGATaaag                                                                                                               >  2:4808295/1‑148 (MQ=255)
                                         gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAAAAAAAAGAAATTATTTTTTTTTTAATATCGTGATAACg                                                                                                               >  1:3569250/1‑150 (MQ=255)
                                           gAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGAt                                                                                                             >  2:171976/1‑150 (MQ=255)
                                           gAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTGTAAATATCGAGATAACGAt                                                                                                             >  2:3852764/1‑150 (MQ=255)
                                           gAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTGGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAGTTAAAAATCAAATTAAATTATTTGTTGGTAAAGAGCGATATAAAGAt                                                                                                             >  2:2774289/1‑150 (MQ=255)
                                            aaaaaagAAAGCCGGAACCGCAAAACCCCCGCGCTCTGCGTACTACACCAACAGCATTCACGATACCACCTACAAATTCATATCCTCTCGATAGATTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATc                                                                                                            <  2:119351/144‑1 (MQ=255)
                                             atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATca                                                                                                           <  1:3551394/150‑1 (MQ=255)
                                             atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATca                                                                                                           <  1:4848921/150‑1 (MQ=255)
                                               atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCGTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATAAAATTTAATTATTTTTTTTTAAATATAGAGATAAAGATcaca                                                                                                         >  1:1854227/1‑150 (MQ=255)
                                                ttGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGTTTTAAAAATAAAATTAAAATATTTTTTTTTAAATTTAGAGATAAAAATAACaa                                                                                                        >  1:4311695/1‑150 (MQ=255)
                                                ttGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATGAAATTATTTTTTTTTAAATATAGAGATAAAGATCACaa                                                                                                        >  1:3797583/1‑150 (MQ=255)
                                                 tGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAATATAATGAAAATATTTTTTTTTAAATATCGAGATAAAGATCACaaa                                                                                                       >  1:3216262/1‑150 (MQ=255)
                                                   aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                                     >  2:2047957/1‑150 (MQ=255)
                                                   aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTGAAAATATTATGAAATTATTTGTTTTTTTATATGGAGTGTACGATGACcaaaa                                                                                                     >  1:3906481/1‑150 (MQ=255)
                                                   aaTGCCCCAACCGCATAACCGCCGCGCTGTGCGATGTACACCATCTGCTTTGTCGATACCACCTACATTTTCATATCCAGTCCTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGATTTTAAATATCGAGATAACGATCACaaaaa                                                                                                     <  2:1370989/150‑1 (MQ=255)
                                                          gAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTGAAAATATAAATAAATTATTTTTTTTTTAATATAGAGATAAACATAAAAAAAACGaacaa                                                                                              >  1:1864981/1‑146 (MQ=255)
                                                          gAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATTAAAATATTTTTTTTTAAATATCGAGATAACGATCACAAAAAAGACaaa                                                                                             >  2:917028/1‑149 (MQ=255)
                                                              cGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTTAATATAGAGATAAAGATCACAAAAAAGACAatatta                                                                                         >  1:5000469/1‑148 (MQ=255)
                                                               gCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGaa                                                                                        >  2:3026559/1‑150 (MQ=255)
                                                               gCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAAAGATAACAAAAACGACAATATGaa                                                                                        >  2:4569362/1‑150 (MQ=255)
                                                                cATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATAGAGATAAAGATCACAAAAAAGACAATATGaaa                                                                                       >  2:4013724/1‑150 (MQ=255)
                                                                           cgcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTAATCgagg                                                                            >  2:3783878/1‑150 (MQ=255)
                                                                            gcTGGGCTTTGTTCCGCATCTGCTTTGTCGAACCCACGTACATTTTCATATCCTGGCGTTTGTTTTCGATTAAAAAATATAATCAAATTTTTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCgagga                                                                           <  1:3148438/150‑1 (MQ=255)
                                                                               gggCGTTGTTCAGCTTCTGTTTTGTCTCCACCACTTAATTTTTCATATCCTGTCGTTTTTTTTCTATTTCAAAATATAATTAAATTATTTTTTTTTAAATATCGAGACAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTg                                                                        <  1:5016527/148‑1 (MQ=255)
                                                                                    ttgttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAgg                                                                  <  2:4687479/150‑1 (MQ=255)
                                                                                        tCAGCATCTGCTTTGCCTAAACCACGAACATTTTCATATCCTGTCTTTTTTTTTCTATTAAAAAATATAATCAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaaa                                                              <  1:3650431/150‑1 (MQ=255)
                                                                                         cAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTCAAAATATAATGAAAATATTTTTTTTTAAAAATATAGATAACGAAAACAAAAAAGACAATATTTAAATTATTAGGGGGGTGAAAGGCaaaa                                                             >  2:4368702/1‑150 (MQ=255)
                                                                                           gCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATATAAAAAAAAAAATAAAATAATTTTTTTTTAATAAAGAGAAAAAGCACACAAAAAAGACACAAATGAAAAAATTAGAGGGGTGAAAGGCaaaaaa                                                           >  1:2418707/1‑150 (MQ=255)
                                                                                             aTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTTTTTTTTTTTAAATATCTAGATAACTATAACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAc                                                         <  1:289119/150‑1 (MQ=255)
                                                                                              tCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAATATAATTACATTATTTTTTTTTTAATAACGTGATAAACCTCCCAAAAAAGACAAGATGTAAAAAATTCTGGGCGTGTGAGGCaaaaaaaag                                                        >  1:2744141/1‑148 (MQ=255)
                                                                                                  cTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATAGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGCGGGGTGAAAGGCAAAAAAACGGCCt                                                    >  2:1469956/1‑150 (MQ=255)
                                                                                                   tttGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAAAATAATGAAAATATTTTTTTTTAAATATAGAGATAACGATCACAAAAACCACAAAATGAAAAATATTCGAGGAGTGAAAGGCAAAAAAAAGGCCTc                                                   >  2:2316027/1‑150 (MQ=255)
                                                                                                    ttGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTTGATTTAAAAAAAAAATTAAAATATTTTTTTTTAAATATATATATAAAGATCACAAAAAACACAAACTGAAAAATATTAGGGGAGTGTAAGGGCaaaaaaatgccacc                                                   >  1:4540375/1‑143 (MQ=255)
                                                                                                         gaaaCCACTTACTTTTAAATATCTTTTGTTTTTTTTTCTATTAAAAAATATAATCAATTTTTTTTTTTTTAAATATCTATATAACTAACAAAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGAAAAAAAACGGCCTCCCGATa                                             <  1:1631713/147‑1 (MQ=255)
                                                                                                         gATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCACAAAAACGGCCTCCCGATa                                             >  2:2831533/1‑150 (MQ=255)
                                                                                                                   aCATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAAGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAAGGCCTCCCGATAGGGAAGCCGt                                   >  2:3671459/1‑150 (MQ=255)
                                                                                                                    cATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                  <  2:1085791/150‑1 (MQ=255)
                                                                                                                     aTTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATAGAGATAACGATCACAAAAACGACAATATTTAAATTATTGGAGGGGTGTAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                 >  2:1207369/1‑150 (MQ=255)
                                                                                                                       gtacATAGCCTGCCTTTTGTTTTCTATTTAAAAATATACTGAAATTATTTTTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGttttt                               <  1:4237583/147‑5 (MQ=255)
                                                                                                                         tCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaaa                             <  2:598668/150‑1 (MQ=255)
                                                                                                                         tCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATAGAGATAACGATCACAAAAACGACAATATGAAAAATATTCGTGGAGTGTAAGGCCAAAAAAAGGCCACCCGAGAGGGAATCCGTAGCaaa                             >  2:1530408/1‑150 (MQ=255)
                                                                                                                         acATATCCTGTCGTTTGTTTACGATTTCAAAATATAATGAAATTTTTTGTTTTTAAATATCTAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaaa                             <  2:1405625/149‑1 (MQ=255)
                                                                                                                               ccTTTCTTTTGTTTTCGATTACAAAATATAATGAAATTTTTTTTTTTTAAATATCGATATAACGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                       <  1:2028909/150‑1 (MQ=255)
                                                                                                                               ccTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGATTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                       <  2:3068701/150‑1 (MQ=255)
                                                                                                                               ccTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                       <  2:3216262/150‑1 (MQ=255)
                                                                                                                               ccTGTCGTTTGTTTTCGATTTATAAATATAATGTAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                       <  2:2905653/150‑1 (MQ=255)
                                                                                                                                     gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt                 <  2:2996427/150‑1 (MQ=255)
                                                                                                                                     gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt                 <  2:4149813/150‑1 (MQ=255)
                                                                                                                                      tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                >  2:3871130/1‑150 (MQ=255)
                                                                                                                                      tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                <  1:3296408/150‑1 (MQ=255)
                                                                                                                                          ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg            >  1:2451846/1‑150 (MQ=255)
                                                                                                                                           tttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGttttttt                             <  2:4345742/132‑7 (MQ=255)
                                                                                                                                               gATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa       >  1:1183214/1‑150 (MQ=255)
                                                                                                                                               gATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa       >  2:659402/1‑150 (MQ=255)
                                                                                                                                                aTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt      >  1:1493521/1‑150 (MQ=255)
                                                                                                                                                  ttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGc    >  2:1299984/1‑150 (MQ=255)
                                                                                                                                                  ttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAAGACACTATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACTGCAACCCGAGAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAGGc    >  1:3744315/1‑150 (MQ=255)
                                                                                                                                                    cAAAATATAATCAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg  <  1:4219212/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
TTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGG  >  NC_000913/2177086‑2177383

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: