Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,306,540 T→C 6.4% S20G (AGC→GGC)  yjcS ← metallo‑beta‑lactamase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,306,5400TC 6.4% 301.3 / 8.5 94S20G (AGC→GGC) yjcSmetallo‑beta‑lactamase superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (46/42);  new base C (3/3);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

GGGGCGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCG  >  NC_000913/4306392‑4306688
                                                                                                                                                    |                                                                                                                                                    
ggggCGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCt                                                                                                                                                      <  4:354298/149‑1 (MQ=255)
     gATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcgg                                                                                                                                                 <  2:5248/149‑1 (MQ=255)
      aTAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcggc                                                                                                                                                >  3:361714/1‑149 (MQ=255)
      aTAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCagcggc                                                                                                                                                >  2:117763/1‑149 (MQ=255)
        aaaTCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCat                                                                                                                                              <  2:583318/149‑2 (MQ=255)
         aaTCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGt                                                                                                                                             <  1:415047/149‑1 (MQ=255)
          aTCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCatga                                                                                                                                            >  3:681115/1‑146 (MQ=255)
           tCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTa                                                                                                                                           >  2:520991/1‑149 (MQ=255)
           tCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCatgat                                                                                                                                           <  4:681115/149‑5 (MQ=255)
           tCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCatgat                                                                                                                                           <  1:117763/149‑5 (MQ=255)
            ccACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTaa                                                                                                                                          <  4:650109/149‑1 (MQ=255)
            ccACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTaa                                                                                                                                          <  4:121024/149‑1 (MQ=255)
             cACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAAc                                                                                                                                         <  4:361714/149‑1 (MQ=255)
             cACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAAc                                                                                                                                         <  2:840706/149‑1 (MQ=255)
                        tCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAAt                                                                                                                              <  3:236715/149‑1 (MQ=255)
                         cGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATc                                                                                                                             <  2:275734/149‑1 (MQ=255)
                          gTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATcc                                                                                                                            >  4:840534/1‑149 (MQ=255)
                           acAAAATCCTGGCGATCGCTAAAAGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTGTCGCTTCCTCTTTTTCGTTATCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCt                                                                                                                           >  2:323554/2‑149 (MQ=255)
                             aaaaTCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGc                                                                                                                         <  3:515829/149‑1 (MQ=255)
                             aaaaTCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGc                                                                                                                         >  1:126755/1‑149 (MQ=255)
                              aaaTCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCt                                                                                                                        <  3:167000/149‑1 (MQ=255)
                                 tCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCAGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCaa                                                                                                                     <  1:858838/149‑1 (MQ=255)
                                    tGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGACTAAGCTGTTGGGTATACTGAGTGGCGGCTTTCGCATCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGCCAGTTAACGCAATAACTATCCTACTCAAACg                                                                                                                  >  3:571928/1‑149 (MQ=255)
                                       cGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGaa                                                                                                               <  3:911985/149‑1 (MQ=255)
                                        gATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAAt                                                                                                              >  2:871476/1‑149 (MQ=255)
                                         aTCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATa                                                                                                             >  4:542598/1‑149 (MQ=255)
                                              tAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCga                                                                                                        >  4:822405/1‑149 (MQ=255)
                                              tAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCga                                                                                                        <  1:332747/149‑1 (MQ=255)
                                              tAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCCCTTCTCCTTTCGCGTTACCGGTATTTACAATCATGCCGCTGGCGGCAGTTAACGAAATAACAATCCTCCTCAAACGGAATAACCga                                                                                                        >  1:573780/1‑149 (MQ=255)
                                                aaCGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCgaga                                                                                                      <  3:612362/149‑1 (MQ=255)
                                                aaCGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCgaga                                                                                                      <  2:126755/149‑1 (MQ=255)
                                                aaCGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCgaga                                                                                                      <  1:434291/149‑1 (MQ=255)
                                                aaCGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCgaga                                                                                                      <  1:520991/149‑1 (MQ=255)
                                                 aCGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGagag                                                                                                     <  1:828185/149‑1 (MQ=255)
                                                    gTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTa                                                                                                  >  2:766767/1‑149 (MQ=255)
                                                    gTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTa                                                                                                  >  2:82235/1‑149 (MQ=255)
                                                     tAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCgaga                                                                                                      <  3:349210/144‑1 (MQ=255)
                                                     tAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCgaga                                                                                                      >  4:349210/1‑144 (MQ=255)
                                                     tAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTAt                                                                                                 <  3:822405/149‑1 (MQ=255)
                                                      aaTGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTAtt                                                                                                <  1:706888/149‑1 (MQ=255)
                                                               gCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAt                                                                                       >  3:621769/1‑149 (MQ=255)
                                                                 gTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAtct                                                                                     >  1:19370/1‑149 (MQ=255)
                                                                      tCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatt                                                                                >  1:887409/1‑149 (MQ=255)
                                                                      tCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatt                                                                                >  1:783114/1‑149 (MQ=255)
                                                                      tCTGATTAACCTGTTGGGTATATTGAGTGGCGGCGTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCATACGGAATAACCGAGAGTTATTCATATAAATCTCcatt                                                                                <  1:766767/149‑1 (MQ=255)
                                                                       cTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcattt                                                                               >  1:818942/1‑149 (MQ=255)
                                                                         gATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcattta                                                                              >  4:369111/1‑148 (MQ=255)
                                                                         gATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcattta                                                                              <  3:369111/148‑1 (MQ=255)
                                                                          aTTAACCTGATGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTGTTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGGTATTCATATAAATCTCcatttatc                                                                            >  4:28808/1‑149 (MQ=255)
                                                                           ttAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatca                                                                           <  3:353089/149‑1 (MQ=255)
                                                                                 tGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcattta                                                                       >  2:345216/1‑147 (MQ=255)
                                                                                 tGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcattta                                                                       <  1:345216/147‑1 (MQ=255)
                                                                                  gTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGt                                                                    <  4:336914/149‑1 (MQ=255)
                                                                                   ttGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTc                                                                   >  2:63808/1‑149 (MQ=255)
                                                                                     gggTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGt                                                                 >  3:902981/1‑149 (MQ=255)
                                                                                          tatTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAata                                                            >  1:314940/1‑149 (MQ=255)
                                                                                          tatTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAata                                                            >  3:876499/1‑149 (MQ=255)
                                                                                             tGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAAc                                                                                                                   <  3:408901/91‑1 (MQ=255)
                                                                                             tGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAAc                                                                                                                   >  4:408901/1‑91 (MQ=255)
                                                                                             tGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGt                                                         <  3:455993/149‑1 (MQ=255)
                                                                                             tGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGt                                                         >  2:120765/1‑149 (MQ=255)
                                                                                               aGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTaa                                                       <  4:895989/149‑1 (MQ=255)
                                                                                                gTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAAt                                                      >  2:800959/1‑149 (MQ=255)
                                                                                                   gcggcTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCatat                                                                                           <  3:114316/109‑1 (MQ=255)
                                                                                                   gcggcTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCatat                                                                                           >  4:114316/1‑109 (MQ=255)
                                                                                                    cggcTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCAt                                                  <  4:712233/149‑1 (MQ=255)
                                                                                                             cTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATAc                                         <  3:753174/149‑1 (MQ=255)
                                                                                                             cTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATAc                                         <  1:179012/149‑1 (MQ=255)
                                                                                                             cTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATAc                                         >  1:704045/1‑149 (MQ=255)
                                                                                                              ttCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACg                                        >  3:853646/1‑149 (MQ=255)
                                                                                                                  tctTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCAt                                    >  3:127286/1‑149 (MQ=255)
                                                                                                                          gTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCt                            >  3:825781/1‑149 (MQ=255)
                                                                                                                          gTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCt                            >  4:681462/1‑149 (MQ=255)
                                                                                                                           ttAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTGTCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGTTg                           >  2:746152/1‑149 (MQ=255)
                                                                                                                           ttAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTg                           >  3:531320/1‑149 (MQ=255)
                                                                                                                             aGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGc                         <  4:621769/149‑1 (MQ=255)
                                                                                                                              gCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAt                                                                                       <  4:894168/86‑1 (MQ=255)
                                                                                                                              gCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAt                                                                                       >  3:894168/1‑86 (MQ=255)
                                                                                                                                 gTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTa                     >  2:121405/1‑149 (MQ=255)
                                                                                                                                  tATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTAt                    >  1:736236/1‑149 (MQ=255)
                                                                                                                                   aTTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAACCCTGCACAAACGGAATAACCGAGAGTTATTCATATCAATCTCCATTTATCATTTAAGTCGTAAATATGTAATCCATGCAAACTATGTCATTTAAAGCTGGCGTTATc                   <  4:575666/149‑1 (MQ=255)
                                                                                                                                    tttACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCa                  >  1:323524/1‑149 (MQ=255)
                                                                                                                                      tACCATCATGCCGCTGGCGTCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTCAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATa                >  3:761905/1‑149 (MQ=255)
                                                                                                                                         catcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAAtt             >  1:384839/1‑149 (MQ=255)
                                                                                                                                          atcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTg            >  1:22866/1‑149 (MQ=255)
                                                                                                                                           tcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTgg           <  2:783114/149‑1 (MQ=255)
                                                                                                                                           tcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTgg           <  2:736236/149‑1 (MQ=255)
                                                                                                                                           tcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTgg           >  2:807456/1‑149 (MQ=255)
                                                                                                                                           tcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTgg           <  4:129761/149‑1 (MQ=255)
                                                                                                                                           tcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTgg           <  1:323554/149‑1 (MQ=255)
                                                                                                                                           tcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTgg           <  3:681462/149‑1 (MQ=255)
                                                                                                                                           aactGCCGCCAGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTgg           >  3:262781/4‑149 (MQ=255)
                                                                                                                                             atGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGCGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGa         >  4:637671/1‑149 (MQ=255)
                                                                                                                                              tGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGaa        >  2:533501/1‑149 (MQ=255)
                                                                                                                                                   cTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGc   <  4:792129/149‑1 (MQ=255)
                                                                                                                                                    tGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCg  <  4:853646/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGGCGATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATTGCG  >  NC_000913/4306392‑4306688

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: