Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,449 T→C 7.9% 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,4490TC 7.9% 291.4 / 10.6 101intergenic (+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 (45/48);  new base C (3/5);  total (48/53)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

TTTGTCACCATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTC  >  NC_000913/256301‑256598
                                                                                                                                                    |                                                                                                                                                     
tttGTCACCATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTg                                                                                                                                                      <  2:221492/150‑1 (MQ=255)
   gTCACCATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCa                                                                                                                                                   <  2:1056406/150‑1 (MQ=255)
    tCACCATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTATTCTTTTCAACGCCTGGCAc                                                                                                                                                  >  1:913801/1‑150 (MQ=255)
    tCACCATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCAc                                                                                                                                                  >  2:1643433/1‑150 (MQ=255)
     cACCATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACt                                                                                                                                                 >  2:713781/1‑150 (MQ=255)
        cATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGcc                                                                                                                                              <  1:66442/150‑1 (MQ=255)
        cATCTTCGCAAAACCGGCTGGGCGTCCGCTGTTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGcc                                                                                                                                              <  1:1103608/150‑1 (MQ=255)
          tCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGTGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCgg                                                                                                                                            >  2:1111831/1‑150 (MQ=255)
           cTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGTCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCggg                                                                                                                                           >  1:1533603/1‑150 (MQ=255)
            ttCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCGAATCTCCGGTCGCTATTCTTTTCAACGCCTGGCACTGCCgcgc                                                                                                                                          >  2:1711243/1‑150 (MQ=255)
              cGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGTCGTAGTATTCGTCACGCCAATCTCAGGTCGCTAATCTTTTCAACGCCTGGCACTGCCgggatt                                                                                                                                        >  1:36321/1‑147 (MQ=255)
               gCAAAACCGGCTGTTCGTCCGCGGGTTTATGACTATGTTGTTGATATCCGGCAAGATACATGGATTGAACAGCCGTTGGATATGGTCGTCGTATTCGTCCCGCCAATCTCCGGTCGGTAATCTTTTCAACGCCCGGCAGTGCCAGGCgtt                                                                                                                                       <  2:1044990/150‑1 (MQ=255)
               gCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGGCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                                                                       >  1:1664906/1‑150 (MQ=255)
                cAAATCCGGCTGGTCGTCCGCTGGTTGAGGACTAAGATGTTCATATCCCGCAATATACCTCGATTGAACAGCCGTGGGGTATGGGCGACGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgttg                                                                                                                                      <  2:388800/150‑1 (MQ=255)
                cAAAACCGTCTGGTCGGCCGCTGGTTGATGACTATGTTGTTGAAATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttg                                                                                                                                      <  1:1376724/150‑1 (MQ=255)
                cAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTCATATCCCGCAAGATACCTGGATTGAACAGCCGAGGGATAAGGGCGTCGTATTAGTCCCGCCAATCACCTCTCGCAAATCTTTTCAAAGCCTGGAACTGCCGGGCgttt                                                                                                                                      >  1:313651/1‑149 (MQ=255)
                   aaCCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTc                                                                                                                                   <  1:1657766/150‑1 (MQ=255)
                   aaCCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCCCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTc                                                                                                                                   <  2:379452/150‑1 (MQ=255)
                    aCCGGCTGGTCGTCCGCTTGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGGTATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttgaaaa                                                                                                                                  <  2:542626/150‑5 (MQ=255)
                     ccGGCGGGGCGTCCGCTGGTTGATGACTATGTTGTTGCTATCCCGCCAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttgaaaag                                                                                                                                 <  1:1804776/150‑6 (MQ=255)
                       ggCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGGCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttgaaaagat                                                                                                                               <  1:1392291/150‑8 (MQ=255)
                        gCTGGTCGTCCGCTGGTTGACGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttttt                                                                                                                              >  2:1134840/1‑150 (MQ=255)
                           gttcgtcCGCTGTTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGTCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAAc                                                                                                                           <  2:574391/148‑1 (MQ=255)
                           ggTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCAGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGATCGCTAATCTTTTCAACGCCTGGCACTGCCGCGCGTTGTTCTTTTTAAc                                                                                                                           >  1:265955/1‑150 (MQ=255)
                                tcCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAg                                                                                                                      <  1:1543513/150‑1 (MQ=255)
                                      ggTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGt                                                                                                                >  1:1424573/1‑150 (MQ=255)
                                      ggTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGt                                                                                                                <  2:29819/150‑1 (MQ=255)
                                      ggTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGt                                                                                                                >  2:17905/1‑150 (MQ=255)
                                         tgatgaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTAc                                                                                                             <  2:1416621/150‑1 (MQ=255)
                                         tgatgaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAAGCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTAc                                                                                                             <  1:1217624/150‑1 (MQ=255)
                                          gatgaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAAGTTCAGGCGGGTTACa                                                                                                            >  2:334608/1‑150 (MQ=255)
                                           atgaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACaa                                                                                                           <  1:630941/150‑1 (MQ=255)
                                           atgaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACaa                                                                                                           >  2:196073/1‑150 (MQ=255)
                                            tgaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGTGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAAt                                                                                                          <  2:1602283/150‑1 (MQ=255)
                                            tgaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAAt                                                                                                          <  2:639756/150‑1 (MQ=255)
                                              aCTCTGGTGTTGATATCCCTCACGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAg                                                                                                        <  1:325319/150‑1 (MQ=255)
                                              aCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGTGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAg                                                                                                        <  1:187891/150‑1 (MQ=255)
                                              aCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTTATCTTTTCAACGCCTGGCACTGCCCGGCGTTGGTGTTTTTAAAATCAGGCGCGGTTCAATaa                                                                                                        >  1:1495023/1‑149 (MQ=255)
                                              aCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAg                                                                                                        <  1:1453005/150‑1 (MQ=255)
                                                 aTGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGttt                                                                                                     <  2:1511672/150‑1 (MQ=255)
                                                   gttgttGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTcc                                                                                                   <  1:1535728/150‑1 (MQ=255)
                                                      gttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagt                                                                                                >  1:1645389/1‑150 (MQ=255)
                                                        tgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaa                                                                                              >  2:486561/1‑150 (MQ=255)
                                                        tgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTATTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaa                                                                                              >  1:1361041/1‑150 (MQ=255)
                                                         gATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGGCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                                                                             <  1:104653/150‑1 (MQ=255)
                                                           tatCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaagta                                                                                           <  2:1332715/150‑1 (MQ=255)
                                                             tCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAtt                                                                                         >  1:658971/1‑150 (MQ=255)
                                                              cccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTc                                                                                        <  2:184772/150‑1 (MQ=255)
                                                               ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTTATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTTACATCAGGGGGGTTACAATAGTTTCCAGTAAGTATTCt                                                                                       >  1:678379/1‑150 (MQ=255)
                                                               ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCt                                                                                       >  2:1196549/1‑150 (MQ=255)
                                                                cGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                                                                      >  2:1456751/1‑150 (MQ=255)
                                                                    aGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAgg                                                                                  >  1:1236968/1‑150 (MQ=255)
                                                                    aGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGGCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTATCTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAgg                                                                                  <  2:1312166/150‑1 (MQ=255)
                                                                     gATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGc                                                                                 <  1:1543316/150‑1 (MQ=255)
                                                                       tACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTg                                                                               <  2:966849/150‑1 (MQ=255)
                                                                        aCCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGc                                                                              <  1:1035789/150‑1 (MQ=255)
                                                                         ccTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCa                                                                             >  1:508702/1‑150 (MQ=255)
                                                                           tGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATc                                                                           >  2:381470/1‑150 (MQ=255)
                                                                               ttGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATg                                                                       <  2:1794400/150‑1 (MQ=255)
                                                                                 gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                                                                     <  1:751827/150‑1 (MQ=255)
                                                                                   aCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGacac                                                                   <  1:462440/150‑1 (MQ=255)
                                                                                       ccGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGc                                                               <  2:1606239/150‑1 (MQ=255)
                                                                                       ccGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGc                                                               >  1:1748352/1‑150 (MQ=255)
                                                                                           gggATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATATTTTCAACGCCTGGCACTGCCGGGCGTTGTTGTTTTTAACTTCAGGCGGGTTTCAATAGTTTCCAGTAAGTATTGTGGAGGCTGCATACATGACACAGGCAAAc                                                           >  1:1538583/1‑150 (MQ=255)
                                                                                            ggATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAcc                                                          >  2:1727343/1‑150 (MQ=255)
                                                                                              atatGGGCGTCGTATTCGTCCCGCAAATCTCCGGTCGCTAATCTTTTCATCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTg                                                        <  2:36321/150‑1 (MQ=255)
                                                                                                 tGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGc                                                     <  1:1206443/150‑1 (MQ=255)
                                                                                                 tGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGGTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGc                                                     >  1:1629358/1‑150 (MQ=255)
                                                                                                      gtcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAAc                                                <  2:584288/150‑1 (MQ=255)
                                                                                                        cgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAAccc                                              >  1:1188738/1‑150 (MQ=255)
                                                                                                          tATTCGTCCCGCCAATCTCCGGTCTCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTg                                            <  1:1648354/150‑1 (MQ=255)
                                                                                                          tATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTg                                            >  1:1645519/1‑150 (MQ=255)
                                                                                                            ttCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGtt                                          <  2:1664906/150‑1 (MQ=255)
                                                                                                              cGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCa                                        <  1:1430681/150‑1 (MQ=255)
                                                                                                              cGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCa                                        >  2:783855/1‑150 (MQ=255)
                                                                                                               gTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCAGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCaa                                       >  1:317169/1‑150 (MQ=255)
                                                                                                                 cccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAAc                                     >  1:599028/1‑150 (MQ=255)
                                                                                                                 cccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATACTATCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAAc                                     <  1:284340/150‑1 (MQ=255)
                                                                                                                  ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAAcc                                    <  1:196073/150‑1 (MQ=255)
                                                                                                                   cGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGGTGTTCTTTTTTACTTCAGGCGGGGTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCAGAGCGAAACCCTGTTCAAAccc                                   >  1:1778557/1‑150 (MQ=255)
                                                                                                                   agCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGTGGATGCATCCATGAAGCAGGCAAACCTGAGCGAAACCCTGTTCAAAccc                                   >  2:60962/2‑150 (MQ=255)
                                                                                                                      cAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGc                                >  1:934508/1‑150 (MQ=255)
                                                                                                                           tcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTaa                           <  2:209244/150‑1 (MQ=255)
                                                                                                                           tcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTaa                           >  2:965400/1‑150 (MQ=255)
                                                                                                                            ccGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTTAAATAGTTTCCAATAAGTATTAAGGTGGCTGGATCCAAGCCACAGGCAAACCGTAGCGAAAACCAGGTCAACACCCCCTGTaaa                          >  1:1590377/1‑150 (MQ=255)
                                                                                                                             cGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTTAGCGAAACCCTGTTCAAACCCCGCTTTAAAc                         <  1:1196549/150‑1 (MQ=255)
                                                                                                                                 cGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATcc                     >  1:452456/1‑150 (MQ=255)
                                                                                                                                  gCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCt                    <  1:822468/150‑1 (MQ=255)
                                                                                                                                   cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTg                   <  1:713781/150‑1 (MQ=255)
                                                                                                                                   cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTAAAACCCCGCTTTAAACATCCTg                   >  1:1594329/1‑150 (MQ=255)
                                                                                                                                   cTAAGCTTTTCATCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTg                   <  2:265955/150‑1 (MQ=255)
                                                                                                                                     aaTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGaa                 <  1:381470/150‑1 (MQ=255)
                                                                                                                                     aaTCTTTACAACGCCTGGCACTGCCGGGCGTGGTTCATTTTAACTTCAGGCGGGTTTCAATAGTTTCCAGTAAGTATTCTGCTGCCTCAAACCATGACAAAGGCAAACCTTAGCGAAACCCTGTTCAAACCACGCTTTAAAAATCCTGaa                 >  2:1508672/1‑150 (MQ=255)
                                                                                                                                         ttttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCGGCATCCATGACACAGGCAAGCCTGAGCGAGACCCTGTTCACCCCCCGCTTTAACCATCCTGTAACCt             >  2:1308683/1‑150 (MQ=255)
                                                                                                                                          gaacAACGCCCGGCAGTGCCAGGCGTTGTACTGTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTTCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTAAAACCCCGCGTTACACATCCTGAACCCta            >  1:1758079/4‑149 (MQ=255)
                                                                                                                                           ttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCg           <  1:1111831/150‑1 (MQ=255)
                                                                                                                                            acAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGa          >  1:436602/2‑150 (MQ=255)
                                                                                                                                              aaCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTATCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGACAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGAcc        >  1:640939/1‑149 (MQ=255)
                                                                                                                                              aaCGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGAcc        >  2:434347/1‑149 (MQ=255)
                                                                                                                                                    tGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTc  >  1:1632970/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTGTCACCATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTC  >  NC_000913/256301‑256598

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