Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,884,678 T→C 9.8% intergenic (+31/+17) 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,6780TC 9.8% 502.5 / 31.8 174intergenic (+31/+17)serA/PP_5156bifunctional D‑3‑phosphoglycerate dehydrogenase/alpha‑ketoglutarate reductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (77/80);  new base C (7/10);  total (84/90)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 182 total aligned reads displayed.

TGCAGACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCCTTG  >  NC_002947/5884529‑5884827
                                                                                                                                                     |                                                                                                                                                     
 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)
         aCGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGGCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTc                                                                                                                                              >  2:3841979/1‑150 (MQ=255)
           gagaAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTccc                                                                                                                                            <  1:1957259/150‑1 (MQ=255)
           gagaAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTccc                                                                                                                                            >  1:1125566/1‑150 (MQ=255)
            agaAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCt                                                                                                                                           <  2:5278668/150‑1 (MQ=255)
              aaaGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCttt                                                                                                                                         <  2:2372118/150‑1 (MQ=255)
              aaaGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttt                                                                                                                                         <  2:1980155/150‑1 (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:2220508/150‑1 (MQ=255)
                     gTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTAtt                                                                                                                                  >  2:72647/1‑150 (MQ=255)
                     gTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCACTTTTTTTAtt                                                                                                                                  <  2:2220432/150‑1 (MQ=255)
                     gTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCtttt                                                                                                                                        >  1:4817265/1‑144 (MQ=255)
                     gTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCtttt                                                                                                                                        <  2:4817265/144‑1 (MQ=255)
                      ttACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTAttt                                                                                                                                 <  2:2429029/150‑1 (MQ=255)
                                   gATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTg                                                                                                                    >  1:4271255/1‑150 (MQ=255)
                                   gATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTg                                                                                                                    >  2:2515918/1‑150 (MQ=255)
                                    aTGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttt                                                                                                                                         <  1:2199972/128‑1 (MQ=255)
                                    aTGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCACCAGGGCCTCCCttt                                                                                                                                         >  2:2199972/1‑128 (MQ=255)
                                           cGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATcttctt                                                                                                            <  1:2579766/150‑1 (MQ=255)
                                           cGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATcttctt                                                                                                            <  2:1986156/150‑1 (MQ=255)
                                             ccGAGTACTCGGACCTGGCGCAGGAAAAACAGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCATGTTCTAAGGTAGCTGGCGGTATGAAAAAGTGAGGCCATGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTtcg                                                                                                          >  2:5008251/1‑150 (MQ=255)
                                                   aCTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTtcgactcga                                                                                                    <  1:4537357/150‑1 (MQ=255)
                                                     tCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACg                                                                                                  >  1:4453055/1‑150 (MQ=255)
                                                     tCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACg                                                                                                  >  1:1871317/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)
                                                               cgcAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGa                                                                                        <  2:1663761/150‑1 (MQ=255)
                                                               cgcAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGa                                                                                        >  1:1738036/1‑150 (MQ=255)
                                                                gcAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTGCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGaa                                                                                       >  1:4852566/1‑150 (MQ=255)
                                                                   ggAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATgg                                                                                    <  2:4898951/150‑1 (MQ=255)
                                                                   ggAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATgg                                                                                    >  1:2982216/1‑150 (MQ=255)
                                                                    gAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGATGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGc                                                                                   <  2:22448/150‑1 (MQ=255)
                                                                     aaaaaCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGca                                                                                  <  1:908650/150‑1 (MQ=255)
                                                                      aaaaCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGcac                                                                                 >  1:2271230/1‑150 (MQ=255)
                                                                      aaaaCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGcac                                                                                 >  2:2139996/1‑150 (MQ=255)
                                                                       aaaCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGcaca                                                                                <  2:2751367/150‑1 (MQ=255)
                                                                       aaaCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGcaca                                                                                <  2:3343282/150‑1 (MQ=255)
                                                                       aaaCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGcaca                                                                                <  1:909838/150‑1 (MQ=255)
                                                                         aCTGGAAGAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAATGGGAGGCGGTGGTGGCGTCGCTTTTTTTATTTCAGTTCAAGCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATg                                                                              >  2:1536715/1‑150 (MQ=255)
                                                                         aCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATg                                                                              <  1:287305/150‑1 (MQ=255)
                                                                         aCTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATg                                                                              <  1:287204/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)
                                                                                          gCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGaa                                                             >  1:3575004/1‑150 (MQ=255)
                                                                                          gCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGaa                                                             <  1:3444318/150‑1 (MQ=255)
                                                                                            aCCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATc                                                           <  2:378282/150‑1 (MQ=255)
                                                                                             ccATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAACAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCa                                                          <  2:2703420/150‑1 (MQ=255)
                                                                                              cATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCTGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCGCTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAg                                                         <  2:2935115/150‑1 (MQ=255)
                                                                                              cATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAg                                                         >  2:2287071/1‑150 (MQ=255)
                                                                                               ctCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGt                                                        <  2:3747384/149‑1 (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                                                                                                                  <  2:228014/74‑1 (MQ=255)
                                                                                                                 ttCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGAt                                                                                                                  >  1:228014/1‑74 (MQ=255)
                                                                                                                    tAAGGTCGCTGGCGGTATGAAAAATGGAGGACCTTGTGCCCTCACTTTATTTATTTCACTTCAAACGTCATATTCTTCGACACGAAGCTAGGCTTGAATGGCACAAGATACTGGAAACCCAGAATCATTTGAAGGGAGAGCTAGCCCggg                                   >  2:2671418/1‑150 (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)
                                                                                                                            cTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACGGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGt                           >  2:1732578/1‑150 (MQ=255)
                                                                                                                               gCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAAt                        >  2:2444663/1‑150 (MQ=255)
                                                                                                                               gCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAAt                        >  1:107807/1‑150 (MQ=255)
                                                                                                                                  gTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTGCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCt                     <  2:4852566/150‑1 (MQ=255)
                                                                                                                                    aTGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTgg                   <  2:1945972/150‑1 (MQ=255)
                                                                                                                                      aaaaaaGGGAGGCCACCAGGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGtt                 >  2:2809587/2‑150 (MQ=255)
                                                                                                                                       aaacaGGGATGCCCTGGTGGCCTCACTTTTTTTATTTCGATTCAACCGTGATCTTCTTCGACTCGACGATAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGttt                <  2:3745186/150‑1 (MQ=255)
                                                                                                                                       aaaaaGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCGCTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGttt                >  1:3745186/1‑150 (MQ=255)
                                                                                                                                         aaaGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCg              <  1:3277889/150‑1 (MQ=255)
                                                                                                                                         aaaGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCg              <  2:5029071/150‑1 (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:1262482/1‑150 (MQ=255)
                                                                                                                                                   ccTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCCt    >  2:1513865/1‑150 (MQ=255)
                                                                                                                                                    cTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCCtt   <  1:2139996/150‑1 (MQ=255)
                                                                                                                                                    cTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCCtt   <  1:1465923/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCAGACCAACGAGAAAGTCGGTTACGTGGTAATCGATGTTGACGCCGAGTACTCGGACCTGGCGCAGGAAAAACTGCAACAGGTCAAAGGCACCATCCGTTCGCGCGTCCTGTTCTAAGGTCGCTGGCGGTATGAAAAAGGGAGGCCCTGGTGGCCTCCCTTTTTTTATTTCACTTCAACCGTGATCTTCTTCGACTCGACGCTAGGCTTGAATGGCACATGATACTGGTCACCCAGAATCAGTTGCAGGGTGTGCTTGCCCGGGGTCAGGGTAATGCTGGTTTCGGTCTGTGCCTTG  >  NC_002947/5884529‑5884827

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