Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,454 C→G 5.6% 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,4540CG 5.6% 297.5 / 6.0 89intergenic (+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 (37/46);  new base G (2/3);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTCCGC  >  NC_000913/256309‑256601
                                                                                                                                                 |                                                                                                                                                   
cATCTTCGCAAAACCGGCTGGTCGTCCGCTGGGTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGc                                                                                                                                                  <  1:1384146/149‑1 (MQ=255)
   attCGAAAAACAGTCTGGTCGTCCTCTGGTGGATGACTATGTTGTTGATATTCCGGAAGATACCTTGATTGACCAGCCGTGGGATATGGGCTTCGTATTCGTCCCGCCAATCTCCGGTGGCTAATCTTTTCAACGCCTGGCACTGCCgg                                                                                                                                               <  2:995892/148‑1 (MQ=255)
    ttCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCggg                                                                                                                                              >  2:83624/1‑149 (MQ=255)
    ttCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCggg                                                                                                                                              >  2:880216/1‑149 (MQ=255)
     tCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGACAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGACGGGc                                                                                                                                             >  1:133600/1‑149 (MQ=255)
      cGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCg                                                                                                                                            <  1:1386024/149‑1 (MQ=255)
       gCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgt                                                                                                                                           <  1:751810/149‑1 (MQ=255)
       gCAAAACCGCCTCGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGAGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgt                                                                                                                                           >  2:895068/1‑149 (MQ=255)
        cAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                                                                          <  2:401665/149‑1 (MQ=255)
        cAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                                                                          >  1:1425731/1‑149 (MQ=255)
        cAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgtt                                                                                                                                          <  2:946337/149‑1 (MQ=255)
         aaaaCCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgttg                                                                                                                                         <  2:216628/149‑1 (MQ=255)
           aaCCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgttgtt                                                                                                                                       >  1:1560620/1‑149 (MQ=255)
             ccGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttgaaaa                                                                                                                                     <  1:177339/149‑5 (MQ=255)
                     tcgtcCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACt                                                                                                                             >  2:154144/1‑149 (MQ=255)
                           gCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGc                                                                                                                       <  2:1100055/149‑1 (MQ=255)
                                 tgatgaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTGTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTa                                                                                                                 <  1:224101/149‑1 (MQ=255)
                                    tgaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACaa                                                                                                              <  1:1327648/149‑1 (MQ=255)
                                      aCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATa                                                                                                            <  1:769706/149‑1 (MQ=255)
                                      aCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATa                                                                                                            <  1:619430/149‑1 (MQ=255)
                                      aCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATa                                                                                                            <  2:647089/149‑1 (MQ=255)
                                         aTGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCAATGCCGGGCGTTGTTCTTTTAAACTTCAGGCGGGTTACAATAGtt                                                                                                         >  1:138416/1‑149 (MQ=255)
                                           gttgttGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTc                                                                                                       <  2:1498547/149‑1 (MQ=255)
                                            ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTcc                                                                                                      >  1:385514/1‑149 (MQ=255)
                                            ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTcc                                                                                                      >  1:441518/1‑149 (MQ=255)
                                             tgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                                                                     <  1:66839/149‑1 (MQ=255)
                                             tgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCAGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa                                                                                                     >  2:1226911/1‑149 (MQ=255)
                                              gttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCag                                                                                                    <  2:42729/149‑1 (MQ=255)
                                                 gATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaa                                                                                                 <  1:154144/149‑1 (MQ=255)
                                                 gATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaa                                                                                                 <  2:498276/149‑1 (MQ=255)
                                                  atatCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                                                                                <  1:625191/149‑1 (MQ=255)
                                                  atatCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                                                                                <  2:1592530/149‑1 (MQ=255)
                                                       ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAtta                                                                                           >  2:775141/1‑148 (MQ=255)
                                                       ccGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTc                                                                                           >  2:200608/1‑149 (MQ=255)
                                                         gCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                                                                         >  1:1392600/1‑149 (MQ=255)
                                                            aGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAg                                                                                      <  1:238748/149‑1 (MQ=255)
                                                             gATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAgg                                                                                     <  2:1266240/149‑1 (MQ=255)
                                                             gATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAgg                                                                                     <  1:1474320/149‑1 (MQ=255)
                                                              aTACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCTATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGc                                                                                    <  1:200608/149‑1 (MQ=255)
                                                              aTACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGc                                                                                    >  2:1060929/1‑149 (MQ=255)
                                                                     caTTGAACAGCCGGGGGATATGGGCGGCGACTTCGTCCCGCCAATCTCCCGCCGCTAAGCTTTTCTACGCCTCTCCCTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATcc                                                                             <  1:99542/148‑1 (MQ=255)
                                                                       ttGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCTCTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCAt                                                                           <  1:1426702/149‑1 (MQ=255)
                                                                       ttGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCAt                                                                           <  2:596231/149‑1 (MQ=255)
                                                                       ttGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCAt                                                                           <  2:1548838/149‑1 (MQ=255)
                                                                       ttGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCAt                                                                           <  2:1019054/149‑1 (MQ=255)
                                                                       ttGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACGTCCGGCGGGTTTCAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCAt                                                                           >  2:1198526/1‑149 (MQ=255)
                                                                          aaCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCTGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGac                                                                        <  1:605205/149‑1 (MQ=255)
                                                                                 gTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCa                                                                 >  1:1098102/1‑149 (MQ=255)
                                                                                  tGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCaa                                                                <  2:1560620/149‑1 (MQ=255)
                                                                                    ggATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAc                                                              <  1:1023701/149‑1 (MQ=255)
                                                                                     gATATGGGCGTCGTAGTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAAcc                                                             >  1:152964/1‑149 (MQ=255)
                                                                                      atatGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCt                                                            >  1:1152433/1‑149 (MQ=255)
                                                                                      atatGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCt                                                            >  1:27688/1‑149 (MQ=255)
                                                                                          gggCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTTACTTCAGGCGGGGTACAATAGTTTCCAGTAACTATTCTGGAGGCTGCATCCATGACACCGGCCAAccggcgc                                                        >  2:1332757/1‑144 (MQ=255)
                                                                                          gggCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGc                                                        >  2:866154/1‑149 (MQ=255)
                                                                                              gtcgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTTACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGaaa                                                    >  1:513808/1‑149 (MQ=255)
                                                                                                cgtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAAcc                                                  <  2:1392600/149‑1 (MQ=255)
                                                                                                  tATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCt                                                >  1:1249177/1‑149 (MQ=255)
                                                                                                     tCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGtt                                             <  1:83624/149‑1 (MQ=255)
                                                                                                     tCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGtt                                             <  1:409526/149‑1 (MQ=255)
                                                                                                     tCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGtt                                             <  1:1757/149‑1 (MQ=255)
                                                                                                           cGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAACCCTGAGCGCAGCCCTGTTCACAcc                                       >  2:1370096/1‑149 (MQ=255)
                                                                                                            gCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAAccc                                      >  1:851572/1‑149 (MQ=255)
                                                                                                               aaTCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGc                                   <  1:1060929/149‑1 (MQ=255)
                                                                                                               aaTCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGc                                   <  1:880216/149‑1 (MQ=255)
                                                                                                                 tctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGcat                                 <  2:1425731/149‑3 (MQ=255)
                                                                                                                 tctcCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCtt                                 >  2:1421220/1‑149 (MQ=255)
                                                                                                                     cGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTaaa                             <  2:1098102/149‑1 (MQ=255)
                                                                                                                     cGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTaaa                             <  1:395218/149‑1 (MQ=255)
                                                                                                                     cGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTATTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATGCATGACAAAGGCAAACCTGATCGAAACCCTGTTAAAACCCCGCTTTaaa                             >  2:628493/1‑149 (MQ=255)
                                                                                                                         cGCTACTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTGAGTCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAATCCTGAGGGAAACCCTGTTCAAACCCCGCTTTAAAGATc                         <  2:624780/149‑1 (MQ=255)
                                                                                                                           cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCt                       >  2:929424/1‑149 (MQ=255)
                                                                                                                           cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCt                       <  2:385994/149‑1 (MQ=255)
                                                                                                                            tAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGCATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAACCCCTGTTCAAACCCCGCTTTAAACATCCTg                      >  1:352399/1‑149 (MQ=255)
                                                                                                                              aTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGaa                    >  2:183843/1‑149 (MQ=255)
                                                                                                                              aTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGAAAACCTGAGCGAAACCCTGTTCAAACACCGCTTTAAACATCCTGaa                    >  1:109315/1‑149 (MQ=255)
                                                                                                                               aaagaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGaaa                   >  2:129020/7‑149 (MQ=255)
                                                                                                                                cTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAAc                  <  2:1371421/149‑1 (MQ=255)
                                                                                                                                     cAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGa             <  1:701147/149‑1 (MQ=255)
                                                                                                                                      aaCGCCTGGCACTGCCGGGCGTTGTTCTTGTTTACTTCAGGCGGGTTACAATAGTTTCCAGGAAGTTTTCTGGAGGCTGCGGCCATGCCCCGGGCAAACCTGAGCGAACCCCGGTTCAAACCCCGCGTTAAACAACCAGAGACCTCGAc            >  1:1406375/1‑149 (MQ=255)
                                                                                                                                       aCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACg           >  1:507473/1‑149 (MQ=255)
                                                                                                                                        cGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGc          <  2:1249177/149‑1 (MQ=255)
                                                                                                                                        cGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGc          >  1:657922/1‑149 (MQ=255)
                                                                                                                                         gCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTTCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCt         >  1:397303/1‑149 (MQ=255)
                                                                                                                                         gCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGCGCGAAACCCTGTTCAACCCCCGCTATAAACATCCTGAACCCTCGCCGCt         >  1:780913/1‑149 (MQ=255)
                                                                                                                                          ccTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTa        <  2:881592/149‑1 (MQ=255)
                                                                                                                                             ggCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTc     >  2:293974/1‑149 (MQ=255)
                                                                                                                                                aCTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTccgc  <  2:385514/149‑1 (MQ=255)
                                                                                                                                                aCTACCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTccgc  >  1:1436611/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CATCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAGTCCGC  >  NC_000913/256309‑256601

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