Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,454 C→G 19.0% intergenic (+19/‑73) 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,4540CG19.0% 464.6 / 90.8 183intergenic (+19/‑73)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 C (86/61);  new base G (21/14);  total (108/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

GTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAAC  >  NC_000913/256355‑256552
                                                                                                   |                                                                                                  
gttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCAc                                                                                                    >  2:1220063/1‑100 (MQ=255)
gttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCAc                                                                                                    <  1:1700462/100‑1 (MQ=255)
      atCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCgg                                                                                              >  2:4140507/1‑100 (MQ=255)
      atCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCgg                                                                                              <  2:4073760/100‑1 (MQ=255)
      atCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGAGCGTCGTATTCGTCCCGCCAATTACCGGTTGCTAATCTTTTCAACGCCTGGCACTGCCgg                                                                                              >  1:2656695/1‑100 (MQ=255)
       tCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCggg                                                                                             <  2:1834255/100‑1 (MQ=255)
       tCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCggg                                                                                             <  2:579570/100‑1 (MQ=255)
         ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCg                                                                                           >  1:4140428/1‑100 (MQ=255)
          cGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgt                                                                                          >  1:3447083/1‑100 (MQ=255)
           gCAAGATACCTGGATTGAACAGCCGTGTGATATGGGCGTCGTATTCGTCCCGCCAATCTCCCGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgtt                                                                                         <  1:10084/100‑1 (MQ=255)
            gatcttACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttg                                                                                        <  2:2518675/95‑1 (MQ=255)
            cAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttg                                                                                        <  1:1587025/100‑1 (MQ=255)
             aaGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTGATCTTTTCAACGACCGGCAGTGCCAGGCgttg                                                                                        <  1:3980413/99‑1 (MQ=255)
             aaGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttga                                                                                       <  2:28428/100‑2 (MQ=39)
             aaGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttga                                                                                       <  1:2743328/100‑2 (MQ=39)
              aGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgttgt                                                                                       <  2:3846727/99‑1 (MQ=255)
              aGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttga                                                                                       <  2:2528050/99‑2 (MQ=39)
                 tACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                                   >  1:1351031/1‑100 (MQ=255)
                 tACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                                   <  2:3161303/100‑1 (MQ=255)
                 tACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                                   >  2:2608895/1‑100 (MQ=255)
                 tACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttgagatc                                                                                   >  1:2518681/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)
                     tGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTa                                                                               <  1:366071/100‑1 (MQ=255)
                     tGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTa                                                                               <  1:487514/100‑1 (MQ=255)
                     tGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTACCGCCAATCTCCGGACGGTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTa                                                                               >  2:3417394/1‑100 (MQ=255)
                     tGGATTGAACAGCCGTGGGAGATGGGCTTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTa                                                                               <  2:2398964/100‑1 (MQ=255)
                      ggATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTaa                                                                              <  1:3865161/100‑1 (MQ=255)
                      ggATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTaa                                                                              <  1:3959950/100‑1 (MQ=255)
                          tGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTc                                                                          >  1:3503816/1‑100 (MQ=255)
                          tGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTc                                                                          <  2:3299737/100‑1 (MQ=255)
                          tGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGGCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTc                                                                          >  1:3929599/1‑100 (MQ=255)
                           gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                         >  2:2719968/1‑100 (MQ=255)
                           gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                         <  1:2103700/100‑1 (MQ=255)
                           gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                         >  1:1748509/1‑100 (MQ=255)
                            aaCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAg                                                                        >  2:136735/1‑100 (MQ=255)
                              cAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGc                                                                      >  1:37289/1‑100 (MQ=255)
                                gCCGTTGGATATGGGCGTCGTATTCGTCGCGCTAAGCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGc                                                                      <  2:721723/98‑1 (MQ=255)
                                gCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCgg                                                                    >  2:506732/1‑100 (MQ=255)
                                  cGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGt                                                                  >  1:796824/1‑100 (MQ=255)
                                  cGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGt                                                                  <  1:2800379/100‑1 (MQ=255)
                                     gggATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTAc                                                               <  2:3608231/100‑1 (MQ=255)
                                           tGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGt                                                         >  2:2239077/1‑100 (MQ=255)
                                             ggCGTCGTATTCGTCCCGCCAAGCTACGCTCGATAATCTTATCCACGCCCAGCACTGTCGGGCGTTGTTCTTTTTCACTTCAGGCGAGATAc                                                               >  2:3936147/1‑92 (MQ=255)
                                               cgtcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTcc                                                     <  1:480534/100‑1 (MQ=255)
                                                 tcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCag                                                   <  1:1256895/100‑1 (MQ=255)
                                                  cgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagt                                                  >  2:713647/1‑100 (MQ=255)
                                                   gtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagt                                                  <  1:631114/99‑1 (MQ=255)
                                                     aTTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                               <  1:3690886/100‑1 (MQ=255)
                                                     aTTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                               >  1:2064895/1‑100 (MQ=255)
                                                     aTTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                               >  2:3001478/1‑100 (MQ=255)
                                                     aTTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                               <  1:4099543/100‑1 (MQ=255)
                                                      ttCGTCCCGTCAATCTCCGGTCGTTAATCTTTTCAACGCCTGGCCCTGCCGGGCTTTTTTCTTTTTATTTTCCGGCGGGTTACAATAGTTTCCagtaagt                                              <  2:3495353/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                                           >  2:3583024/1‑100 (MQ=255)
                                                           nccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCt                                         >  2:2218594/2‑100 (MQ=255)
                                                            ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                        >  2:391193/1‑100 (MQ=255)
                                                            ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                        >  2:3099437/1‑100 (MQ=255)
                                                             cGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTgg                                       <  1:532930/100‑1 (MQ=255)
                                                              nccAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGa                                      >  1:2247375/2‑100 (MQ=255)
                                                              gCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGa                                      >  2:2381057/1‑100 (MQ=255)
                                                                  aTCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCt                                  >  2:4099140/1‑100 (MQ=255)
                                                                  aTCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCt                                  <  2:1711789/100‑1 (MQ=255)
                                                                   tctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTg                                 >  1:2674313/1‑100 (MQ=255)
                                                                      ccGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCAt                              >  2:2604530/1‑100 (MQ=255)
                                                                       cGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATc                             <  2:277945/100‑1 (MQ=255)
                                                                        ggTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATcc                            >  1:6111/1‑100 (MQ=255)
                                                                           cGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATg                         <  1:331034/100‑1 (MQ=255)
                                                                              taaaaagaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGaca                      >  1:1286882/10‑100 (MQ=37)
                                                                               aaTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGacac                     >  2:19200/1‑100 (MQ=255)
                                                                               aaTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGacac                     >  2:3505702/1‑100 (MQ=255)
                                                                                aaaagaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGAcaca                    >  1:3573217/8‑100 (MQ=37)
                                                                                 tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTCACTTCAGGCGGGTTACAATACATTCCAGAAAGTATTCTGGAGGCAGCCTCCATGAcaca                    >  1:396413/1‑99 (MQ=255)
                                                                                 tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAg                   <  1:218610/100‑1 (MQ=255)
                                                                                 tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAg                   >  2:1623922/1‑100 (MQ=255)
                                                                                  cTTTTGAGCGGCTGGCAGTGCCGGGCGTGGTTCTTTTTAACTTCGGGGGGGTTACAATAGTTTCCCGGAGGTATTCTGGAGGCTGCGTCCCTGAcaca                    <  2:2194172/98‑1 (MQ=255)
                                                                                  cTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAgg                  >  2:35626/1‑100 (MQ=255)
                                                                                  cTTTTCAACGCCTGGCACTGCCGGGAGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAgg                  >  1:561494/1‑100 (MQ=255)
                                                                                  aagaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTGCCAGTCAGTATTCTGGAGGCTGCATCCATGACACAgg                  >  2:2821632/6‑100 (MQ=37)
                                                                                   ttttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGc                 <  1:900080/100‑1 (MQ=255)
                                                                                       cAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAc             >  2:3697523/1‑100 (MQ=255)
                                                                                       cAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTCCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAc             >  1:2012146/1‑100 (MQ=255)
                                                                                       cAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAc             >  2:1683639/1‑100 (MQ=255)
                                                                                        aaCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAcc            >  1:2730342/1‑100 (MQ=255)
                                                                                        aaCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAcc            <  1:1623926/100‑1 (MQ=255)
                                                                                        aaCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGCGGCTGCATCCATGACACAGGCAAAcc            >  1:1217707/1‑100 (MQ=255)
                                                                                        aaCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAcc            >  2:1196721/1‑100 (MQ=255)
                                                                                        aaCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAcc            >  2:1180958/1‑100 (MQ=255)
                                                                                        aaCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAcc            >  1:841840/1‑100 (MQ=255)
                                                                                         aCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCt           >  1:959219/1‑100 (MQ=255)
                                                                                          cGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGCGGCTGCATCCATTACACAGGCAAAcat           >  1:302663/1‑97 (MQ=255)
                                                                                          cGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTg          <  1:601963/100‑1 (MQ=255)
                                                                                              tGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCg      >  1:626251/1‑100 (MQ=255)
                                                                                                 cACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGaaa   >  2:3450105/1‑100 (MQ=255)
                                                                                                 cACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGaaa   >  1:47049/1‑100 (MQ=255)
                                                                                                  aCTGCCGGGCGTTGTTCTTTTTAACTTCAGGCTGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAAc  <  2:1333910/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAAC  >  NC_000913/256355‑256552

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: