Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,884,677 C→A 10.1% intergenic (+30/+18) serA → / ← PP_5156 bifunctional D‑3‑phosphoglycerate dehydrogenase/alpha‑ketoglutarate reductase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,884,6770CA10.1% 540.7 / 38.4 175intergenic (+30/+18)serA/PP_5156bifunctional D‑3‑phosphoglycerate dehydrogenase/alpha‑ketoglutarate reductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (75/81);  new base A (8/10);  total (84/91)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 182 total aligned reads displayed.

CTGCAGACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCCTT  >  NC_002947/5884528‑5884826
                                                                                                                                                     |                                                                                                                                                     
cTGCAGACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGccc                                                                                                                                                       <  2:3107920/150‑1 (MQ=255)
 tGCAGACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCt                                                                                                                                                      >  2:1727915/1‑150 (MQ=255)
  gCAGACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGGAAAAGGGAGGCCCtg                                                                                                                                                     <  1:2332490/150‑1 (MQ=255)
  gCAGACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCtg                                                                                                                                                     <  2:3621624/150‑1 (MQ=255)
  gCAGACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCtg                                                                                                                                                     <  2:1205354/150‑1 (MQ=255)
   cAGACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCtgg                                                                                                                                                    >  2:1326496/1‑150 (MQ=255)
          aCGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAATGGGAGGCCCTGGTGGCCTc                                                                                                                                             >  1:2478491/1‑150 (MQ=255)
          aCGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTc                                                                                                                                             >  2:2423882/1‑150 (MQ=255)
          aCGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGGCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTc                                                                                                                                             >  2:3841979/1‑150 (MQ=255)
            gagaAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTccc                                                                                                                                           >  1:1125566/1‑150 (MQ=255)
            gagaAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTccc                                                                                                                                           <  1:1957259/150‑1 (MQ=255)
             agaAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCt                                                                                                                                          <  2:5278668/150‑1 (MQ=255)
               aaaGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttt                                                                                                                                        <  2:1980262/150‑1 (MQ=255)
               aaaGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttt                                                                                                                                        <  2:1980155/150‑1 (MQ=255)
                 aGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCttttt                                                                                                                                      >  2:2195557/1‑150 (MQ=255)
                 aGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCttttt                                                                                                                                      >  1:3883554/1‑150 (MQ=255)
                 aGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttttt                                                                                                                                      <  2:1074572/150‑1 (MQ=255)
                 aGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttttt                                                                                                                                      <  1:3093066/150‑1 (MQ=255)
                    cGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTa                                                                                                                                   <  1:1714112/150‑1 (MQ=255)
                    cGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTa                                                                                                                                   <  1:4141667/150‑1 (MQ=255)
                      gTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTAtt                                                                                                                                 >  2:72647/1‑150 (MQ=255)
                      gTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTAtt                                                                                                                                 <  2:2220508/150‑1 (MQ=255)
                      gTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCtttt                                                                                                                                       >  1:4817265/1‑144 (MQ=255)
                      gTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCtttt                                                                                                                                       <  2:4817265/144‑1 (MQ=255)
                           gTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACt                                                                                                                            >  2:2918922/1‑150 (MQ=255)
                             ggTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAACCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTc                                                                                                                          <  2:2236849/150‑1 (MQ=255)
                             ggTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTc                                                                                                                          >  2:2579766/1‑150 (MQ=255)
                              gTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttttt                                                                                                                                      <  2:221978/137‑1 (MQ=255)
                                    gATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTg                                                                                                                   >  1:4271255/1‑150 (MQ=255)
                                     aTGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttt                                                                                                                                        <  1:2199972/128‑1 (MQ=255)
                                            cGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATcttctt                                                                                                           <  1:2579766/150‑1 (MQ=255)
                                              ccGAGTACTCGGACCTGGCGCAGGAAAAACAGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCATGTTCTAAGGTAGCTGGCGGTATGAAAAAGTGAGGCCATGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTtcg                                                                                                         >  2:5008251/1‑150 (MQ=255)
                                                    aCTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTtcgactcga                                                                                                   <  1:4537357/150‑1 (MQ=255)
                                                      tCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAATGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACg                                                                                                 >  2:2365644/1‑150 (MQ=255)
                                                      tCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGCTCTTCTTCGACTCGACg                                                                                                 >  2:1986323/1‑150 (MQ=255)
                                                      tCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACg                                                                                                 >  1:1871317/1‑150 (MQ=255)
                                                      tCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACg                                                                                                 >  1:4453055/1‑150 (MQ=255)
                                                      tCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCACTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACg                                                                                                 >  1:1299762/1‑150 (MQ=255)
                                                        ggACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCt                                                                                               <  1:3890168/150‑1 (MQ=255)
                                                         gACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTa                                                                                              <  1:2423882/150‑1 (MQ=255)
                                                           ccTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAgg                                                                                            <  1:1646230/150‑1 (MQ=255)
                                                             tGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCt                                                                                          >  2:2700550/1‑150 (MQ=255)
                                                                cgcAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGa                                                                                       >  1:1738036/1‑150 (MQ=255)
                                                                 gcAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTGCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGaa                                                                                      >  1:4852566/1‑150 (MQ=255)
                                                                    ggAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATgg                                                                                   >  1:2982216/1‑150 (MQ=255)
                                                                    ggAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATgg                                                                                   <  2:4898951/150‑1 (MQ=255)
                                                                      aaaaaCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGca                                                                                 <  1:908650/150‑1 (MQ=255)
                                                                       aaaaCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGcac                                                                                >  1:2271230/1‑150 (MQ=255)
                                                                        aaaCTGCAACAGGTCAAAGGCACCCTCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGcaca                                                                               <  2:2751360/150‑1 (MQ=255)
                                                                        aaaCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGcaca                                                                               <  1:909838/150‑1 (MQ=255)
                                                                        aaaCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGcaca                                                                               <  2:3343282/150‑1 (MQ=255)
                                                                          aCTGGAAGAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAATGGGAGGCGGTGGTGGCGTCGCTTTTTTTATTTCAGTTCAAGCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATg                                                                             >  2:1536715/1‑150 (MQ=255)
                                                                          aCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATg                                                                             <  1:287204/150‑1 (MQ=255)
                                                                          aCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATg                                                                             <  1:287305/150‑1 (MQ=255)
                                                                           cTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGa                                                                            <  1:615348/150‑1 (MQ=255)
                                                                            tGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGAt                                                                           >  1:2106070/1‑150 (MQ=255)
                                                                              cAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATAc                                                                         <  1:4010708/150‑1 (MQ=255)
                                                                               aaCAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCt                                                                                                                                              <  1:4390594/80‑1 (MQ=255)
                                                                               aaCAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCt                                                                                                                                              >  2:4390594/1‑80 (MQ=255)
                                                                                     tCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCAc                                                                  <  1:2511699/150‑1 (MQ=255)
                                                                                     tCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttt                                                                                                                                        <  2:2771570/80‑1 (MQ=37)
                                                                                     tCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttt                                                                                                                                        >  1:2771570/1‑80 (MQ=37)
                                                                                      cAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCAcc                                                                 <  2:624059/150‑1 (MQ=255)
                                                                                       aaaGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGAGACTGGTCAccc                                                                >  2:2776949/1‑150 (MQ=255)
                                                                                           gCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGaa                                                            <  1:3444318/150‑1 (MQ=255)
                                                                                           gCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGaa                                                            >  1:3575004/1‑150 (MQ=255)
                                                                                             aCCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATc                                                          <  2:378282/150‑1 (MQ=255)
                                                                                               cATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAg                                                        >  2:2287043/1‑150 (MQ=255)
                                                                                                ctCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGt                                                       <  2:3747384/149‑1 (MQ=255)
                                                                                                aTCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGt                                                       >  2:2477739/1‑150 (MQ=255)
                                                                                                 tCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGtt                                                      <  1:3787848/150‑1 (MQ=255)
                                                                                                         cgcgTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGgtg                                              >  2:3970861/1‑150 (MQ=255)
                                                                                                            gTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGc                                           <  1:1727915/150‑1 (MQ=255)
                                                                                                               cTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTg                                        <  1:334883/150‑1 (MQ=255)
                                                                                                                tGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTg                                                                       <  1:4859611/118‑1 (MQ=255)
                                                                                                                tGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGc                                       >  2:4697579/1‑150 (MQ=255)
                                                                                                                  ttCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGAt                                                                                                                 >  1:228014/1‑74 (MQ=255)
                                                                                                                         gTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGc                                           <  2:130414/137‑1 (MQ=255)
                                                                                                                         gTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCa                              >  1:414876/1‑150 (MQ=255)
                                                                                                                          tCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAg                             <  1:3341182/150‑1 (MQ=255)
                                                                                                                           cGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAgg                            <  1:2014277/150‑1 (MQ=255)
                                                                                                                            gCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAggg                           >  1:4088177/1‑150 (MQ=255)
                                                                                                                                gCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAAt                       >  1:107807/1‑150 (MQ=255)
                                                                                                                                gCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAAt                       >  2:2444663/1‑150 (MQ=255)
                                                                                                                                   gTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTGCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCt                    <  2:4852566/150‑1 (MQ=255)
                                                                                                                                        aaacaGGGATGCCCTGGTGGCCTCACTTTTTTTATTTCGATTCAACCGTGATCTTCTTCGACTCGACGATAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGttt               <  2:3745186/150‑1 (MQ=255)
                                                                                                                                        aaaaaGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCGCTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGttt               >  1:3745186/1‑150 (MQ=255)
                                                                                                                                          aaaGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCg             >  2:2249956/1‑150 (MQ=255)
                                                                                                                                          aaaGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCg             <  2:5029071/150‑1 (MQ=255)
                                                                                                                                          aaaGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCg             <  1:3277889/150‑1 (MQ=255)
                                                                                                                                          aaaGGGAGGCCACCAGGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCg             >  2:1894012/1‑150 (MQ=255)
                                                                                                                                              ggAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCt         <  1:5327684/150‑1 (MQ=255)
                                                                                                                                                aGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCtgt       <  1:5198366/150‑1 (MQ=255)
                                                                                                                                                aGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCtgt       <  2:5165071/150‑1 (MQ=255)
                                                                                                                                                  gCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGc     <  1:3970861/150‑1 (MQ=255)
                                                                                                                                                  gCCCTGGTGGCCTCCCTTATTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGc     <  2:5326907/150‑1 (MQ=255)
                                                                                                                                                    ccTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCCt   >  2:1513865/1‑150 (MQ=255)
                                                                                                                                                    ccTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCCt   >  2:1262482/1‑150 (MQ=255)
                                                                                                                                                     cTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCCtt  <  1:1465923/150‑1 (MQ=255)
                                                                                                                                                     cTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCCtt  <  1:2139996/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCAGACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCCTT  >  NC_002947/5884528‑5884826

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: