Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,454 C→G 17.3% 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,4540CG17.3% 432.2 / 76.2 166intergenic (+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 (69/68);  new base G (17/12);  total (86/80)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

GTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAA  >  NC_000913/256355‑256550
                                                                                                   |                                                                                                
gttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCAc                                                                                                  >  2:1008308/1‑100 (MQ=255)
        cccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCGCTGCCGGGc                                                                                          >  1:2830517/1‑100 (MQ=255)
         ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGcn                                                                                         <  1:3219623/100‑2 (MQ=255)
         ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCg                                                                                         >  1:1966652/1‑100 (MQ=255)
         ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCg                                                                                         <  1:2887215/100‑1 (MQ=255)
          cGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgt                                                                                        >  2:138707/1‑100 (MQ=255)
          cGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgt                                                                                        <  1:2404444/100‑1 (MQ=255)
           gCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                       >  1:2255061/1‑100 (MQ=255)
           gCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                       >  2:3349237/1‑100 (MQ=255)
           gCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                       >  2:2719150/1‑100 (MQ=255)
            cAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttg                                                                                      <  1:311224/100‑1 (MQ=255)
            cAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAAGCTTTTCAACGCCCGGCAGTGCCAGGCgttg                                                                                      <  2:426723/100‑1 (MQ=255)
            cAAGATACCTGGATTGAACAGCCGTGGGATATGGCCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttg                                                                                      <  2:1625135/100‑1 (MQ=255)
             aaGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttga                                                                                     <  1:1460858/100‑2 (MQ=39)
             aaGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttga                                                                                     <  2:874217/100‑2 (MQ=39)
              aGATACCTTGATTGAACAGCCGTGGGATATGCGCGTCGTATTCGTCCCGCCAATATCCGGTCGCTAATCTTTTCAACGCCTGGCCCTGCCGGGCgttgtt                                                                                    <  1:720035/100‑1 (MQ=255)
              aGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttgaa                                                                                    <  1:1779806/100‑3 (MQ=38)
              aGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTACTCCCGCCAATCTCCGGTCGCTAATATTTTCAACGCCTGGCACTGCCGGGCgttg                                                                                      >  2:1090503/1‑98 (MQ=255)
                aTACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCt                                                                                  >  1:2735711/1‑100 (MQ=255)
                aTACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttgaaaa                                                                                  <  1:737227/100‑5 (MQ=38)
                aTACCTCGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCt                                                                                  <  2:1641804/100‑1 (MQ=255)
                  aCCTGGATTGAACAGCCGTCGGATATGGGGGTCGTATTCGTCCCGCCGAGCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                                 <  2:3537169/99‑1 (MQ=255)
                   ccTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtttt                                                                               >  2:3008975/1‑100 (MQ=255)
                    cTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGTTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttttt                                                                              <  2:2663863/100‑1 (MQ=255)
                        aTTTAACAGCCGTCGGATATGGGTGTCGTATTCGTCCCGCCAATCTCCTGTCGCTAATCTTTTCAACCTCTTGCTCTGCCTGCCGTTGTTCTTTTTaa                                                                            <  2:1722868/98‑1 (MQ=255)
                          tGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTc                                                                        <  1:2940816/100‑1 (MQ=255)
                           gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                       <  1:3637862/100‑1 (MQ=255)
                           gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                       >  1:1539785/1‑100 (MQ=255)
                           gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                       >  1:3030119/1‑100 (MQ=255)
                             aCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAgg                                                                     >  1:1486913/1‑100 (MQ=255)
                               aGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAgg                                                                     >  2:3090269/1‑98 (MQ=255)
                                gCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTaagatcggaaga                                                                  >  2:905083/1‑90 (MQ=255)
                                gCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCgg                                                                  >  2:1199913/1‑100 (MQ=255)
                                 ccGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCggg                                                                 <  1:1948500/100‑1 (MQ=255)
                                  cGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGt                                                                <  1:2869114/100‑1 (MQ=255)
                                   gTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGtt                                                               >  2:645769/1‑100 (MQ=255)
                                    tGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTa                                                              <  2:2412860/100‑1 (MQ=255)
                                     gggATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTAc                                                             <  1:2978359/100‑1 (MQ=255)
                                     gggATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTAc                                                             <  1:3334508/100‑1 (MQ=255)
                                      ggATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACa                                                            <  1:1733424/100‑1 (MQ=255)
                                       gATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACaa                                                           <  1:1345778/100‑1 (MQ=255)
                                       gATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACaa                                                           <  2:3484918/100‑1 (MQ=255)
                                        atatGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAAt                                                          <  1:272369/100‑1 (MQ=255)
                                           tGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGt                                                       >  2:2368691/1‑100 (MQ=255)
                                             ggCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGttt                                                     <  1:1008315/100‑1 (MQ=255)
                                              gCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTc                                                    <  1:119708/100‑1 (MQ=255)
                                               cgtcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTcc                                                   <  1:2276798/100‑1 (MQ=255)
                                                gtcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTcc                                                   >  2:1789609/1‑99 (MQ=255)
                                                gtcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                  >  2:1547805/1‑100 (MQ=255)
                                                  cgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagt                                                >  2:3272727/1‑100 (MQ=255)
                                                   gtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagta                                               >  1:451186/1‑100 (MQ=255)
                                                    tATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaa                                              >  2:2062129/1‑100 (MQ=255)
                                                     aTTCGTCCCGCCTATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                             <  1:1723713/100‑1 (MQ=255)
                                                     aTTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                             >  1:772361/1‑100 (MQ=255)
                                                        cGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAt                                          <  1:596867/100‑1 (MQ=255)
                                                        cGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAt                                          >  1:3521988/1‑100 (MQ=255)
                                                        cGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAt                                          >  2:81049/1‑100 (MQ=255)
                                                          tCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTc                                        <  1:3063133/100‑1 (MQ=255)
                                                            ccGCCTATCTCCGGTCGCTAATCTTTTCACCTTCTGGCACTGCCGGTCGTTGTTCTTTTTAATTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCt                                       <  1:347700/99‑1 (MQ=255)
                                                            ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                      >  2:3266704/1‑100 (MQ=255)
                                                            ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                      <  1:81049/100‑1 (MQ=255)
                                                                  aTCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGcn                                <  1:567534/100‑2 (MQ=255)
                                                                   tctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTg                               >  2:1326954/1‑100 (MQ=255)
                                                                   tctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTg                               <  2:1836375/100‑1 (MQ=255)
                                                                   tctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTg                               <  1:2698590/100‑1 (MQ=255)
                                                                    ctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGc                              <  1:1258692/100‑1 (MQ=255)
                                                                     tcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCa                             <  1:3453909/100‑1 (MQ=255)
                                                                        gttCGCTAATCTTTGCAACGCCGGACACTTCCTGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCGGCATcc                          <  2:2314446/98‑1 (MQ=255)
                                                                             cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                     >  1:1120246/1‑100 (MQ=255)
                                                                             cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                     >  2:811494/1‑100 (MQ=255)
                                                                             cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                     >  1:1611606/1‑100 (MQ=255)
                                                                             cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                     <  1:487150/100‑1 (MQ=255)
                                                                                aaaagaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAGCTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGAcaca                  >  2:1985571/8‑100 (MQ=37)
                                                                                 tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAg                 >  2:1156482/1‑100 (MQ=255)
                                                                                   ttttCAACGCCTGGCGCTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGc               <  2:2830497/100‑1 (MQ=255)
                                                                                   ttttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAgg                >  2:1473032/1‑99 (MQ=255)
                                                                                   ttttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGc               >  2:3323271/1‑100 (MQ=255)
                                                                                    tttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCa              <  1:1031085/100‑1 (MQ=255)
                                                                                    gaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCa              >  1:316392/4‑100 (MQ=38)
                                                                                    gaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCa              >  2:2756705/4‑100 (MQ=38)
                                                                                      tCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCaaa            <  1:3640385/100‑1 (MQ=255)
                                                                                      tCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCaaa            >  2:2494627/1‑100 (MQ=255)
                                                                                      acAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCaaa            >  1:992410/2‑100 (MQ=39)
                                                                                      acAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCaaa            >  2:2344007/2‑100 (MQ=39)
                                                                                      acAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCaaa            >  1:1690048/2‑100 (MQ=39)
                                                                                      acAACACCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGTCTACAATAGTTTCCCGTAAGTATTCTGGAGTCTGCATCCATGACACAGGc               >  2:638881/2‑97 (MQ=255)
                                                                                       cAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCCGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCCAAc           >  1:3382404/1‑100 (MQ=255)
                                                                                       cAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAc           >  2:2711012/1‑100 (MQ=255)
                                                                                       cAACGCCTGGCACTGCCGGGAGTTGTTCTTTTTAACTTCAGGCGGATTCAAAAATTTTCCAATAAGTAGTATGCAGGATGCCTCCATGaccacgacagac           >  2:3620235/1‑90 (MQ=255)
                                                                                        aaCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAcc          >  2:2037633/1‑100 (MQ=255)
                                                                                        aaCGCCCGGCAGTGCCTGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTACGTATTCTGGAGGCTGCATCCATGACACACGCAAAcc          >  2:3219474/1‑100 (MQ=255)
                                                                                        aaCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAcc          >  1:2990369/1‑100 (MQ=255)
                                                                                         gcGCCCGGCAGTGCCAGGCGAGGTTCTTTTTAACTTCAAGCGGGTTACAATAGTTCCTAGTAAGTATTCTGGAGGCTCCATCCAAGACACAGGCAAAcc          >  2:755764/2‑99 (MQ=255)
                                                                                         aCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCt         >  2:1845059/1‑100 (MQ=255)
                                                                                            ccTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAg      <  2:479154/100‑1 (MQ=255)
                                                                                             cTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGc     <  1:1849908/100‑1 (MQ=255)
                                                                                             cTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGc     <  2:25554/100‑1 (MQ=255)
                                                                                              tGGCACTGCCGGGCGTTGTTCGTTTTAACCTCAGGCAGGGTACAGCAGGTGCGAGTGAGTATTCTGGCGGCTGCCTCCCGCGCACAGGCAGAACAGAGCg    >  2:1298176/1‑100 (MQ=255)
                                                                                               ggCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGa   >  1:1920463/1‑100 (MQ=255)
                                                                                                gCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGaa  >  1:2086946/1‑100 (MQ=255)
                                                                                                   |                                                                                                
GTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAA  >  NC_000913/256355‑256550

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: