Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,449 T→C 18.8% intergenic (+14/‑78) gpt → / → frsA xanthine phosphoribosyltransferase; xanthine‑guanine phosphoribosyltransferase/fermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913256,4490TC18.8% 400.6 / 77.5 184intergenic (+14/‑78)gpt/frsAxanthine phosphoribosyltransferase; xanthine‑guanine phosphoribosyltransferase/fermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity
Reads supporting (aligned to +/- strand):  ref base T (87/62);  new base C (22/13);  total (109/75)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

ATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCG  >  NC_000913/256350‑256548
                                                                                                   |                                                                                                   
aTGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCt                                                                                                     >  2:3722817/1‑100 (MQ=255)
   ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGc                                                                                                  >  2:2019955/1‑100 (MQ=255)
    tgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCa                                                                                                 >  2:2217566/1‑100 (MQ=255)
    tgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCa                                                                                                 >  2:4099514/1‑100 (MQ=255)
     gttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCAc                                                                                                <  1:1700462/100‑1 (MQ=255)
     gttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCAc                                                                                                >  2:1220063/1‑100 (MQ=255)
           atCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCgg                                                                                          >  2:4140507/1‑100 (MQ=255)
           atCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCgg                                                                                          <  2:4073760/100‑1 (MQ=255)
            tCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCggg                                                                                         <  2:579570/100‑1 (MQ=255)
            tCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCggg                                                                                         <  2:1834255/100‑1 (MQ=255)
             cccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGc                                                                                        >  1:2246304/1‑100 (MQ=255)
             cccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGc                                                                                        <  1:999522/100‑1 (MQ=255)
              ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCg                                                                                       >  1:2419779/1‑100 (MQ=255)
              ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCg                                                                                       >  1:4140428/1‑100 (MQ=255)
               cGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgt                                                                                      <  1:3226108/100‑1 (MQ=255)
                gCAAGATACCTGGATTGAACAGCCGTGTGATATGGGCGTCGTATTCGTCCCGCCAATCTCCCGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgtt                                                                                     <  1:10084/100‑1 (MQ=255)
                 cAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttg                                                                                    <  1:1587025/100‑1 (MQ=255)
                  aaGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTGATCTTTTCAACGACCGGCAGTGCCAGGCgttg                                                                                    <  1:3980413/99‑1 (MQ=255)
                  aaGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttga                                                                                   <  1:2743328/100‑2 (MQ=39)
                   aGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgttg                                                                                    <  2:3852097/98‑1 (MQ=255)
                    gATACCTGGATTGACCCGCCGTGGGGGATGGGCGTCGTATTCGTCCCGCGAGGCTCCGGTCGCTAATCTTTGCAACGCCTGGCACTGCCGGGCgttg                                                                                    <  2:3933256/97‑1 (MQ=255)
                    gATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGGCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttga                                                                                   <  2:2855203/98‑2 (MQ=38)
                      tACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                               >  2:2608895/1‑100 (MQ=255)
                      tACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                               <  2:3161303/100‑1 (MQ=255)
                      tACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                               >  1:1351031/1‑100 (MQ=255)
                      tACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttgagatc                                                                               >  1:522881/1‑95 (MQ=38)
                        ccTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTAGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtttt                                                                             >  2:2686771/1‑100 (MQ=255)
                         cTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttttt                                                                            <  2:2318105/100‑1 (MQ=255)
                         cTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttttt                                                                            <  1:1120116/100‑1 (MQ=255)
                         cTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttttt                                                                            >  1:1789742/1‑100 (MQ=255)
                         cTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttttt                                                                            <  1:2648289/100‑1 (MQ=255)
                         cTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTCGTTCttttt                                                                            >  2:3541544/1‑100 (MQ=255)
                          tGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTACCGCCAATCTCCGGACGGTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTa                                                                           >  2:3417394/1‑100 (MQ=255)
                           ggATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTaa                                                                          <  1:3865161/100‑1 (MQ=255)
                             aTTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACt                                                                        <  1:802433/100‑1 (MQ=255)
                               tGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTc                                                                      >  1:3503816/1‑100 (MQ=255)
                               tGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTc                                                                      <  2:3299737/100‑1 (MQ=255)
                                gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                     <  1:2104731/100‑1 (MQ=255)
                                gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                     >  2:2719968/1‑100 (MQ=255)
                                gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                     >  1:1748509/1‑100 (MQ=255)
                                 aaCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAg                                                                    >  2:136735/1‑100 (MQ=255)
                                     gCCGTTGGATATGGGCGTCGTATTCGTCGCGCTAAGCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGc                                                                  <  2:721723/98‑1 (MQ=255)
                                     gCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCgg                                                                >  2:506732/1‑100 (MQ=255)
                                       cGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGt                                                              <  1:2800379/100‑1 (MQ=255)
                                       cGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGt                                                              >  1:796824/1‑100 (MQ=255)
                                            gATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCATGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTCCaa                                                         >  1:1678941/1‑100 (MQ=255)
                                                tGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGt                                                     >  2:2543664/1‑100 (MQ=255)
                                                tGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGt                                                     >  2:2239077/1‑100 (MQ=255)
                                                 gggCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGtt                                                    <  2:118755/100‑1 (MQ=255)
                                                    cgtcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTcc                                                 <  1:480534/100‑1 (MQ=255)
                                                    cgtcgtATTCGTCCCGCCAATCACCGGTCGCTAATCTTGTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGTTTAAACTAGTTTcc                                                 >  2:3977693/1‑100 (MQ=255)
                                                      tcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCag                                               <  1:1256895/100‑1 (MQ=255)
                                                      tcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCag                                               <  1:2686778/100‑1 (MQ=255)
                                                       cgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagt                                              >  2:713647/1‑100 (MQ=255)
                                                        gtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagt                                              <  1:631114/99‑1 (MQ=255)
                                                          aTTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                           <  2:3590830/100‑1 (MQ=255)
                                                          aTTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                           >  1:2064895/1‑100 (MQ=255)
                                                          aTTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                           <  1:4099543/100‑1 (MQ=255)
                                                           ttCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaagt                                          <  1:4054187/100‑1 (MQ=255)
                                                             cGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAt                                        >  2:3413267/1‑100 (MQ=255)
                                                              gTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAtt                                       >  1:1912086/1‑100 (MQ=255)
                                                              gTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAtt                                       >  1:3350299/1‑100 (MQ=255)
                                                              gTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAtt                                       >  2:2485850/1‑100 (MQ=255)
                                                                nccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCt                                     >  2:2218594/2‑100 (MQ=255)
                                                                 ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                    >  2:3099437/1‑100 (MQ=255)
                                                                 ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                    >  2:3020228/1‑100 (MQ=255)
                                                                  cGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTgg                                   >  1:3698585/1‑100 (MQ=255)
                                                                   nccAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGa                                  >  1:2247375/2‑100 (MQ=255)
                                                                   gCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGa                                  >  2:2381057/1‑100 (MQ=255)
                                                                       aTCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCt                              <  2:1711789/100‑1 (MQ=255)
                                                                         ctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGc                            <  2:3819846/100‑1 (MQ=255)
                                                                          tcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCa                           >  1:553367/1‑100 (MQ=255)
                                                                             ggTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATcc                        <  1:3884385/100‑1 (MQ=255)
                                                                                cGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATg                     >  1:62282/1‑100 (MQ=255)
                                                                                  cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                   >  1:607299/1‑100 (MQ=255)
                                                                                  cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                   <  1:2668613/100‑1 (MQ=255)
                                                                                  cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                   <  1:3945152/100‑1 (MQ=255)
                                                                                   taaaaagaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGaca                  >  1:1286882/10‑100 (MQ=37)
                                                                                     aaaagaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGAcaca                >  2:2838220/8‑100 (MQ=38)
                                                                                     aTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGAcaca                >  1:893741/1‑100 (MQ=255)
                                                                                      tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTCACTTCAGGCGGGTTACAATACATTCCAGAAAGTATTCTGGAGGCAGCCTCCATGAcaca                >  1:396413/1‑99 (MQ=255)
                                                                                      tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAg               >  2:1623922/1‑100 (MQ=255)
                                                                                       cTTTTGAGCGGCTGGCAGTGCCGGGCGTGGTTCTTTTTAACTTCGGGGGGGTTACAATAGTTTCCCGGAGGTATTCTGGAGGCTGCGTCCCTGAcaca                <  2:2194172/98‑1 (MQ=255)
                                                                                       cTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAgg              <  1:577545/100‑1 (MQ=255)
                                                                                       cTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAgg              >  2:35626/1‑100 (MQ=255)
                                                                                       aagaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTCCGTTAGTTTCCCGTAGGTATTCTGGAGGTTGCATACAtatta                  >  2:2606122/6‑92 (MQ=16)
                                                                                        ttttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCagatcggaaga             >  1:47359/1‑90 (MQ=255)
                                                                                        ttttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGc             <  1:3686609/100‑1 (MQ=255)
                                                                                        agaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAg               >  2:3656224/5‑98 (MQ=38)
                                                                                          aacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGCCGGGTTACAATAGTTTCCAGTAAGTATTCTAGAGGCTGCATCCATGACACAGCCaa           >  2:2790243/3‑100 (MQ=38)
                                                                                           acAACGCCCGGCAGTGCCAGGCGTTGTTGTTTTTAACTTCAGGCGGGTTACAACAGTTCCCCGTAAGTCGTCTGGAGGCTGCATCCATGACACCGGCaaa          >  2:1352662/2‑100 (MQ=25)
                                                                                           acAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCaaa          >  1:364106/2‑100 (MQ=39)
                                                                                            cAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTCCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAc         >  1:2012146/1‑100 (MQ=255)
                                                                                            cAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAc         >  2:1683639/1‑100 (MQ=255)
                                                                                             aaCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAGCTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCAGGAGGCTGCATCCATGACACAGGCAAAcc        >  2:2523880/1‑100 (MQ=255)
                                                                                             aaCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGCGGCTGCATCCATGACACAGGCAAAcc        >  1:1217707/1‑100 (MQ=255)
                                                                                              aCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCt       >  2:3503507/1‑100 (MQ=255)
                                                                                              aCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCTGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCt       >  2:3456892/1‑100 (MQ=255)
                                                                                              aCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCt       >  1:958714/1‑100 (MQ=255)
                                                                                               cGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGCGGCTGCATCCATTACACAGGCAAAcat       >  1:302663/1‑97 (MQ=255)
                                                                                                   |                                                                                                   
ATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCG  >  NC_000913/256350‑256548

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