Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,459 G→A 5.5% intergenic (+24/‑68) 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,4590GA 5.5% 805.3 / 22.0 239intergenic (+24/‑68)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 G (90/134);  new base A (6/7);  total (98/141)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 245 total aligned reads displayed.

ATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTCCGCCGCTTTA  >  NC_000913/256310‑256608
                                                                                                                                                     |                                                                                                                                                     
aTCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCg                                                                                                                                                       >  1:3354403/1‑150 (MQ=255)
 tCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCgg                                                                                                                                                      <  2:896727/150‑1 (MQ=255)
  cTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCggg                                                                                                                                                     <  2:1380235/150‑1 (MQ=255)
   ttCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGc                                                                                                                                                    >  2:3600149/1‑150 (MQ=255)
   ttCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGc                                                                                                                                                    >  2:3430893/1‑150 (MQ=255)
   ttCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGc                                                                                                                                                    <  1:3313476/150‑1 (MQ=255)
   ttCGCAAAACCGGCTGATCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTACTCGTTTCAACGCCTGGCACTGCCGGGc                                                                                                                                                    >  2:3599951/1‑150 (MQ=255)
    tCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCg                                                                                                                                                   <  2:1072387/150‑1 (MQ=255)
      gCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                                                                                 >  2:2715962/1‑150 (MQ=255)
      gCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                                                                                 <  1:1722355/150‑1 (MQ=255)
      gCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgtt                                                                                                                                                 <  1:2469520/150‑1 (MQ=255)
         aaaCCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgttgtt                                                                                                                                              >  2:546761/1‑150 (MQ=255)
            ccGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTcgt                                                                                                                                           >  1:385030/1‑148 (MQ=255)
            ccGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                                                                                                                           <  1:1638142/150‑1 (MQ=255)
             cGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttt                                                                                                                                          <  1:3098410/150‑1 (MQ=255)
              ggCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtttt                                                                                                                                         >  1:2173468/1‑150 (MQ=255)
              ggCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtttt                                                                                                                                         <  1:1073199/150‑1 (MQ=255)
              ggCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAAACTTTTCAACGCCTGGCACTGCCGGGCGGTGTTAtttt                                                                                                                                         >  1:1118494/1‑150 (MQ=255)
                   gtcgtcCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACt                                                                                                                                    >  2:3409700/1‑150 (MQ=255)
                    tcgtcCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACtt                                                                                                                                   >  2:2725842/1‑150 (MQ=255)
                      gtcCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                                                                                 <  2:1325457/150‑1 (MQ=255)
                          gCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCg                                                                                                                             >  1:1615460/1‑150 (MQ=255)
                           cTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCTTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCgg                                                                                                                            <  2:2035651/150‑1 (MQ=255)
                           cTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCgg                                                                                                                            <  1:1957171/150‑1 (MQ=255)
                            tGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCggg                                                                                                                           <  1:2270010/150‑1 (MQ=255)
                            tGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCggg                                                                                                                           <  2:2763905/150‑1 (MQ=255)
                              gTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATATTTTCAACGCCAGGCACAGCCGGGCGGTGTTGTTTTTAACTTAAGGCGGGtt                                                                                                                         >  2:1783632/1‑150 (MQ=255)
                                 gatgaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACa                                                                                                                      >  2:3338130/1‑150 (MQ=255)
                                  atgaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACaa                                                                                                                     <  2:1857392/150‑1 (MQ=255)
                                    gaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATa                                                                                                                   <  1:805788/150‑1 (MQ=255)
                                    gaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATa                                                                                                                   <  1:3633857/150‑1 (MQ=255)
                                    gaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGGTGTTCTTTTTAACTTCAGGCGGGTTACAATa                                                                                                                   >  2:3252070/1‑150 (MQ=255)
                                         tgttgtTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTc                                                                                                              <  2:3549219/150‑1 (MQ=255)
                                           ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                                                                            <  2:3488313/150‑1 (MQ=255)
                                           ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                                                                            <  1:2879476/150‑1 (MQ=255)
                                           ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                                                                            <  1:2816650/150‑1 (MQ=255)
                                              ttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGGTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagta                                                                                                         >  2:3906826/1‑150 (MQ=255)
                                              ttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTAAACGCCAGGCACTGCCGGGCGTTGTTAATTTTAAAATAATGCGGGTTACAATAGTTTACagta                                                                                                         >  2:833970/1‑150 (MQ=255)
                                                gATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                                                                                       <  2:363272/150‑1 (MQ=255)
                                                gATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGAGATGGGCGTCGTATTCGTCCCGCCACTCTCCGGTCGCGAAGCGGTTCAACGCCTGGCACGGCCGGGCGTTGTTCTTTTTAACGTCCGGCGGGTTCCCAGAGTTTCCagtaag                                                                                                       >  2:3983348/1‑150 (MQ=255)
                                                 atatCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTTATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCATGCGGGTTACAATAGTTTGaagtaagt                                                                                                      >  1:3935787/1‑150 (MQ=255)
                                                 atatCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCAGGCCCAGCCGGGCGTTGTTGTTTTTAACATCAAGCCGGTTACAAAAGTTTCaagtaagt                                                                                                      >  1:424440/1‑150 (MQ=255)
                                                      ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCt                                                                                                 <  1:445622/150‑1 (MQ=255)
                                                      ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCt                                                                                                 <  2:2595509/150‑1 (MQ=255)
                                                       cGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTGATCTTTTAAACGCCTGGCACAGCCGGGCGTTGTTGTTTTTTAATTAAGGCGGGTGAAAATAGTTTGCAGTAAGTAttagt                                                                                                >  1:1126731/1‑147 (MQ=255)
                                                         cAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGa                                                                                              <  2:3845538/150‑1 (MQ=255)
                                                           aGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAAATTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAgg                                                                                            >  1:3760729/1‑150 (MQ=255)
                                                                 cTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGGAAGTATTCTGGAGGCTGCAt                                                                                      >  2:305493/1‑150 (MQ=255)
                                                                  tGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCATGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATc                                                                                     >  2:2217927/1‑150 (MQ=255)
                                                                     ttagaccAGCCGTGGGATTTGGGCGTCGACTTCGTCCCGACAATCTCCGGTCGCTAATCATTTCAAAGCCGGGCACTGCCGGGAGTTGTTCATTAGAACTTCATGCGTGTTTAACTAGTTTCCAGTAAGTATTCTGGTGTCTGCATCGAt                                                                                  <  2:2299777/144‑1 (MQ=255)
                                                                     aTTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCag                                                                                  <  1:1525888/150‑2 (MQ=255)
                                                                     aTTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCAt                                                                                  <  1:984702/150‑1 (MQ=255)
                                                                     aTTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCAt                                                                                  <  2:2349251/150‑1 (MQ=255)
                                                                         aaCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCACTCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGGAAGTATTCTGGCGGCTGCATCCATGaca                                                                              >  2:2605639/1‑150 (MQ=255)
                                                                           cAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGAcaca                                                                            <  2:2474104/150‑1 (MQ=255)
                                                                            aGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAg                                                                           <  2:2856389/150‑1 (MQ=255)
                                                                            aGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGGTGTTCTTTTTAAATTCAGTCGGGTTACAATAGTTTACCGTAAGTATTCTGGGGGCTGCATCCAAGACACAg                                                                           >  1:2018888/1‑150 (MQ=255)
                                                                             gCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAgg                                                                          <  2:1591078/150‑1 (MQ=255)
                                                                             gCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAgg                                                                          <  2:1591189/150‑1 (MQ=255)
                                                                              ccGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGc                                                                         <  2:3021689/150‑1 (MQ=255)
                                                                              ccGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGc                                                                         <  2:1565082/150‑1 (MQ=255)
                                                                              ccGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGc                                                                         <  1:2649698/150‑1 (MQ=255)
                                                                                      tatGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGa                                                                 <  2:3002518/150‑1 (MQ=255)
                                                                                              tcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAAcc                                                         <  2:1970228/150‑1 (MQ=255)
                                                                                              tcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAAcc                                                         <  1:233455/150‑1 (MQ=255)
                                                                                                gtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCt                                                       <  2:3130621/150‑1 (MQ=255)
                                                                                                gtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCt                                                       >  2:687548/1‑150 (MQ=255)
                                                                                                   ttCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGCGGCTGCATCCAAGACACAGGCAAACCTGAGCGAAACCCTGtt                                                    >  2:514555/1‑150 (MQ=255)
                                                                                                   ttCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAACCCTGAGCGAAACCCTGtt                                                    <  2:1201142/150‑1 (MQ=255)
                                                                                                   ttCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGtt                                                    <  2:777746/150‑1 (MQ=255)
                                                                                                   ttCGTCCCGCCAATCTCCGGTCGCTAATCGTTTCAACGCCTGGCACTGCCGGGCGGTGTTCTTTTTAACTTCAGGCGGGTTACAATAGGTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTgat                                                    >  2:654824/1‑148 (MQ=255)
                                                                                                       tCCCGCCAATCTCCGGGCGCTAAGCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAACTGAGCGAAACCCTGTTCaaa                                                >  2:1634837/1‑150 (MQ=255)
                                                                                                           gCCAATCTCCGGTCGCTAATCTTTTCAAAGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAAcccc                                            >  2:938545/1‑150 (MQ=255)
                                                                                                            ccAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCg                                           <  1:865665/150‑1 (MQ=255)
                                                                                                                      gTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAATCCCTGTTCAAACCCCGCTTTAAACAt                                 <  2:343538/150‑1 (MQ=255)
                                                                                                                        cGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATcc                               <  1:1243621/150‑1 (MQ=255)
                                                                                                                         gCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCt                              <  2:2120759/150‑1 (MQ=255)
                                                                                                                         gCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCt                              >  2:2564796/1‑150 (MQ=255)
                                                                                                                         gCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCt                              <  2:3117365/150‑1 (MQ=255)
                                                                                                                          cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTg                             <  1:1380281/150‑1 (MQ=255)
                                                                                                                           taaaattttCAACGACTGGCACTGCCGGGCGTTGTTCTTTATAACTTCAGTCGGGTTACAATAGTTTCCCGTAAGTATTATGGAGGCTGCATACATGACACAGGAAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGa                            >  2:2292717/6‑150 (MQ=255)
                                                                                                                           tAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGGTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGa                            >  2:1574648/1‑150 (MQ=255)
                                                                                                                           tAATCTTTTCAACGCCTGGCACTGCCGGGCGTGGTTATTTTTAAAATAAGGCGGGTTACAATATTTTACAGAAAGTATTCTTGAGGCTGCATCCAAGACACAGGCAAACCAGAGCGAAAACCAGTTCAAACCCCCCGTTAAACaaacaga                            >  1:2912987/1‑144 (MQ=255)
                                                                                                                           tAATCTTCTCAACGCCTGGCACTGCCGGGCGTTGTTCTATTGAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAACCCCTGTTCAAACCCCGCTTTAAACATCCTGa                            <  1:3989149/150‑1 (MQ=255)
                                                                                                                            aaTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGaa                           <  1:687410/150‑1 (MQ=255)
                                                                                                                                 tttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTAATGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCAGTTCAAACACCGCTTTAAACATACTGaaaacac                      >  1:2250791/1‑146 (MQ=255)
                                                                                                                                  ttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGGTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCg                     >  2:3700785/1‑150 (MQ=255)
                                                                                                                                  ttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAAAATAAGGCGGGTTACAATATTTTCCAGTAAGTATTCTGGAGGCTGCATCCAAGACACAGGCAAACCTGAGCGAAACCCAGTTCAAACCCCCCATTAAACATCCAGAAACCACg                     >  1:2113026/1‑150 (MQ=255)
                                                                                                                                  ttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAAAATAAGGCGGGTTACAATATTTTACAGTAAGTATTATGGAGGCTGCATCCATGACACAGGCAAACCAGAGCGAAAACCAGTTCAAACCCCCCGTTAAAAATACAGAAACCTCg                     >  1:2298082/1‑150 (MQ=255)
                                                                                                                                    cAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGAc                   <  1:3409700/150‑1 (MQ=255)
                                                                                                                                     aaCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCCGTACGTCTTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCGTTAAACATCCTGAAACCTCGACg                  >  2:3188323/1‑150 (MQ=255)
                                                                                                                                      aCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGc                 >  2:2821316/1‑150 (MQ=255)
                                                                                                                                      aCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGc                 >  1:604577/1‑150 (MQ=255)
                                                                                                                                        gCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTa               <  2:1118576/150‑1 (MQ=255)
                                                                                                                                           tGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTc            >  2:1161155/1‑150 (MQ=255)
                                                                                                                                            ggCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTcc           <  1:1290367/150‑1 (MQ=255)
                                                                                                                                               aCTGCCGGGCGTTGTTATTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCCGAAACCTCGACGCTAGTccgcc        <  2:1216147/150‑1 (MQ=255)
                                                                                                                                                 tGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATTGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTACACATCCTGAAACCTCGACGCTAGTCcgccgc      <  1:746460/150‑1 (MQ=255)
                                                                                                                                                 tGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTCcgccgc      <  2:2174738/150‑1 (MQ=255)
                                                                                                                                                     gggCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTCCGCCGCTTTa  <  1:833970/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTCCGCCGCTTTA  >  NC_000913/256310‑256608

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