Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,375,484 A→T 5.8% L160H (CTT→CAT) ‡ menC ← O‑succinylbenzoyl‑CoA synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,375,4840AT 5.8% 293.0 / 5.4 83L160H (CTT→CAT) ‡menCO‑succinylbenzoyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base A (32/46);  new base T (3/2);  total (35/48)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGCCAT  >  NC_000913/2375336‑2375631
                                                                                                                                                    |                                                                                                                                                   
ggcaatgCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCaa                                                                                                                                                     <  3:194428/149‑1 (MQ=255)
  caatgCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTNAAGa                                                                                                                                                   <  3:751727/149‑1 (MQ=255)
        cGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCaaa                                                                                                                                             <  4:514896/149‑1 (MQ=255)
         ggTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCaaa                                                                                                                                             <  1:845055/148‑1 (MQ=255)
         ggTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCaaa                                                                                                                                             <  2:502838/148‑1 (MQ=255)
         ggTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAAt                                                                                                                                            <  3:296940/149‑1 (MQ=255)
          gTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATg                                                                                                                                           >  4:235621/1‑149 (MQ=255)
                 gggCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATcc                                                                                                                                    <  4:646195/149‑1 (MQ=255)
                 gggCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATcc                                                                                                                                    <  3:201357/149‑1 (MQ=255)
                    cAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGa                                                                                                                                 <  4:596249/149‑1 (MQ=255)
                    cAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGa                                                                                                                                 >  3:459688/1‑149 (MQ=255)
                    cAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGa                                                                                                                                 >  3:731087/1‑149 (MQ=255)
                     aaaCGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGaa                                                                                                                                >  3:341626/1‑149 (MQ=255)
                     aaaCGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGa                                                                                                                                 <  2:491879/148‑1 (MQ=255)
                       aCGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAAtt                                                                                                                              <  3:410662/149‑1 (MQ=255)
                          cTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAAtt                                                                                                                              <  3:596610/146‑1 (MQ=255)
                          cTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGcc                                                                                                                           <  3:757782/149‑1 (MQ=255)
                          cTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGc                                                                                                                            <  2:863521/148‑1 (MQ=255)
                          cTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTCGCGTCAAGACGCAAATGCAGATCCGGAAtt                                                                                                                              >  4:596610/1‑146 (MQ=255)
                            cgcgAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCt                                                                                                                          >  1:759083/1‑148 (MQ=255)
                            cgcgAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCt                                                                                                                          >  2:559781/1‑148 (MQ=255)
                            cgcgAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCt                                                                                                                          <  1:906632/148‑1 (MQ=255)
                            cgcgAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCt                                                                                                                          >  1:334319/1‑148 (MQ=255)
                             gcgAATCATCGCGGGTTTTGAACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTcc                                                                                                                        <  4:684401/149‑1 (MQ=255)
                               gAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAg                                                                                                                      <  3:324969/149‑1 (MQ=255)
                                aaTCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAg                                                                                                                      >  2:376967/1‑148 (MQ=255)
                                   catcGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCaa                                                                                                                   <  1:366627/148‑1 (MQ=255)
                                        cgGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGAtt                                                                                                             <  3:369741/149‑1 (MQ=255)
                                           gTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTca                                                                                                           <  1:14098/148‑1 (MQ=255)
                                                gCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCaccacc                                                                                                      <  2:437654/148‑1 (MQ=255)
                                                         ttCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTcgc                                                                                            >  4:363121/1‑149 (MQ=255)
                                                             agaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCAc                                                                                        <  4:73195/149‑1 (MQ=255)
                                                                aaCGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGc                                                                                     <  4:4321/149‑1 (MQ=255)
                                                                   gcgATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCtt                                                                                   <  1:304409/148‑1 (MQ=255)
                                                                             cgcgATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTcaccac                                                                                                       <  1:606760/118‑1 (MQ=255)
                                                                             cgcgATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTcaccac                                                                                                       >  2:606760/1‑118 (MQ=255)
                                                                              gcgATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATccc                                                                        <  1:170101/148‑1 (MQ=255)
                                                                               cgATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCAc                                                                      <  4:459688/149‑1 (MQ=255)
                                                                                gATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGa                                                                                                                                                   >  4:695661/1‑71 (MQ=255)
                                                                                gATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGa                                                                                                                                                   <  3:695661/71‑1 (MQ=255)
                                                                                 aTAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACt                                                                     >  1:571878/1‑148 (MQ=255)
                                                                                     tCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAc                                                                >  4:455236/1‑149 (MQ=255)
                                                                                     tCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTAGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGa                                                                 >  2:744835/1‑148 (MQ=255)
                                                                                               cGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACt                                                      <  3:143684/149‑1 (MQ=255)
                                                                                                 tATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACttt                                                    >  3:570278/1‑149 (MQ=255)
                                                                                                   tttGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTtc                                                   <  2:334319/148‑1 (MQ=255)
                                                                                                   tttGGCAAACTGCTGACCTTTCAGCGGTGTACAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCAt                                                                                                    >  4:884886/1‑99 (MQ=255)
                                                                                                   tttGGCAAACTGCTGACCTTTCAGCGGTGTACAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCAt                                                                                                    <  3:884886/99‑1 (MQ=255)
                                                                                                    ttGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTTCAATCCCACTTTGACCTTCGCCACTTtct                                                  <  2:280243/148‑1 (MQ=255)
                                                                                                    ttGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATc                                                                          >  1:315746/1‑124 (MQ=255)
                                                                                                    ttGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATc                                                                          <  2:315746/124‑1 (MQ=255)
                                                                                                      ggCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCg                                                <  1:392930/148‑1 (MQ=255)
                                                                                                       gCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGcc                                              <  4:624780/149‑1 (MQ=255)
                                                                                                       gCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGcc                                              >  4:819776/1‑149 (MQ=255)
                                                                                                          aaCTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTgg                                           >  4:532687/1‑149 (MQ=255)
                                                                                                          aaCTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTg                                            >  1:308029/1‑148 (MQ=255)
                                                                                                                 gACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAc                                                                >  4:756765/1‑121 (MQ=255)
                                                                                                                 gACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAc                                                                <  3:756765/121‑1 (MQ=255)
                                                                                                                 gACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCt                                    <  3:363121/149‑1 (MQ=255)
                                                                                                                    cTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCa                                 >  4:662626/1‑149 (MQ=255)
                                                                                                                           cGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGa                          <  4:731087/149‑1 (MQ=255)
                                                                                                                            ggTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGa                          <  2:172491/148‑1 (MQ=255)
                                                                                                                            ggTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGa                          <  1:744835/148‑1 (MQ=255)
                                                                                                                            ggTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGa                          <  2:170592/148‑1 (MQ=255)
                                                                                                                            ggTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAg                         <  3:235621/149‑1 (MQ=255)
                                                                                                                            ggTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAg                         <  3:662626/149‑1 (MQ=255)
                                                                                                                            ggTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAg                         <  4:341626/149‑1 (MQ=255)
                                                                                                                               gtCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTg                                            <  4:693476/127‑1 (MQ=255)
                                                                                                                               gtCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTg                                            >  3:693476/1‑127 (MQ=255)
                                                                                                                               gtCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGAt                      >  3:425613/1‑149 (MQ=255)
                                                                                                                                 ccAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATc                     <  1:559781/148‑1 (MQ=255)
                                                                                                                                   aGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGa                  >  4:478319/1‑149 (MQ=255)
                                                                                                                                    ggCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGAt                 >  3:743546/1‑149 (MQ=255)
                                                                                                                                     gcgcGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGAt                 >  1:831840/1‑148 (MQ=255)
                                                                                                                                     gcgcGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATc                >  4:126830/1‑149 (MQ=255)
                                                                                                                                      cgcgATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATc                >  2:628212/1‑148 (MQ=255)
                                                                                                                                         gATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTc             <  2:104189/148‑1 (MQ=255)
                                                                                                                                           tttGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCg           >  1:93379/1‑148 (MQ=255)
                                                                                                                                           tttGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCg           >  2:629850/1‑148 (MQ=255)
                                                                                                                                           tttGCGTCTTGACGCAAATACAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCgg          >  4:517350/1‑149 (MQ=255)
                                                                                                                                            ttGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCgg          >  1:573070/1‑148 (MQ=255)
                                                                                                                                              gCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATc       <  3:387566/149‑1 (MQ=255)
                                                                                                                                                 tCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGcc    <  4:570278/149‑1 (MQ=255)
                                                                                                                                                  cAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGCCa   >  3:336360/1‑149 (MQ=255)
                                                                                                                                                   aaGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGccgt  <  4:650301/149‑3 (MQ=255)
                                                                                                                                                   aaGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGCCa   <  2:275665/148‑1 (MQ=255)
                                                                                                                                                    aGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGCCAt  >  1:730827/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGCCAT  >  NC_000913/2375336‑2375631

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑