Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,080,480 G→A 8.5% H353H (CAC→CAT)  lysC ← aspartokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,080,4800GA 8.5% 193.0 / 6.7 71H353H (CAC→CAT) lysCaspartokinase
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base A (3/3);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

GATCGAACGCGGTATGCAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACGGTGAAGGTGAAG  >  NC_002947/5080332‑5080624
                                                                                                                                                    |                                                                                                                                                
gATCGAACGCGGTATGCAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCg                                                                                                                                                  >  7:87948/1‑149 (MQ=255)
 aTCGAACGCGGTATGCAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGt                                                                                                                                                 <  6:94104/149‑1 (MQ=255)
    gAACGCGGTATGCAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGc                                                                                                                                               <  2:301238/148‑1 (MQ=255)
    gAACGCGGTATGCAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGc                                                                                                                                               <  2:164815/148‑1 (MQ=255)
          ggTATGCAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCATGCGTTCGCa                                                                                                                                        >  6:167032/1‑149 (MQ=255)
            tATGCAGCGCGCGTACGACCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATg                                                                                                                                      <  4:57854/149‑1 (MQ=255)
                cAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCATGCGTTCGCACgccg                                                                                                                                   <  3:273238/148‑1 (MQ=255)
                cAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCATGCGTTCGCACgccg                                                                                                                                   <  5:167032/148‑1 (MQ=255)
                   cgcgcgTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATgccgacgc                                                                                                                               <  1:178875/149‑1 (MQ=255)
                       cgTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCATGCGTTCGCACgccggtgttg                                                                                                                             >  4:13679/1‑141 (MQ=255)
                       cgTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCATGCGTTCGCACgccggtgttg                                                                                                                             <  3:13679/147‑7 (MQ=255)
                        gTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATgccgacgccgac                                                                                                                           >  2:153254/1‑148 (MQ=255)
                            ggCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATc                                                                                                                       >  4:171454/1‑148 (MQ=255)
                            ggCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATCCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATc                                                                                                                       >  7:200545/1‑148 (MQ=255)
                             gCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCga                                                                                                                     <  8:290636/149‑1 (MQ=255)
                              ccAGCTACAGCTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATAATCTGGAAGTTGATCATCTCCTTGGCCAGGCCTTAGAACATGCAGCTGGCAACACCGGCGTGAGAAAGCATGCCGACGCCGACCATc                                                                                                                       >  3:3675/1‑146 (MQ=255)
                                aGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGAc                                                                                                                  >  7:176347/1‑149 (MQ=255)
                                      aGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGc                                                                                                            >  7:176495/1‑149 (MQ=255)
                                      aGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGc                                                                                                            >  7:24941/1‑149 (MQ=255)
                                          aCTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATc                                                                                                        >  3:206881/1‑149 (MQ=255)
                                              ctctTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGGTGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCtt                                                                                                      <  4:112688/147‑1 (MQ=255)
                                              ctctTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCtt                                                                                                      <  2:138636/147‑1 (MQ=255)
                                              ctctTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTg                                                                                                     <  4:36517/148‑1 (MQ=255)
                                                    gAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGc                                                                                              <  8:176347/149‑1 (MQ=255)
                                                      gCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGgt                                                                                                   >  4:13925/1‑142 (MQ=255)
                                                      gCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGgt                                                                                                   <  3:13925/142‑1 (MQ=255)
                                                          caccGAAACCTTGATCTCGGAAGTGGAGGTCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCa                                                                                        >  5:287949/1‑149 (MQ=255)
                                                          caccGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCa                                                                                        >  7:9386/1‑149 (MQ=255)
                                                             cGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACtt                                                                                     >  7:249557/1‑149 (MQ=255)
                                                                aaCCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCAc                                                                                  <  8:9386/149‑1 (MQ=255)
                                                                  ggttGTTCTCGGAAGTGGAGATCATCTGGATGTTGAGGGTCTCCTTGTCCCGGGCTTGTTACTTGGTGATGTCAAAACGGGGGTGGGAACTAATGCCGACGCCGACGATCGAGACCTTGTCGATCTTGGTGTCGCCGATCACTTGACgt                                                                                <  4:3675/147‑2 (MQ=255)
                                                                             gAATTGGAGATCATCTGGATGTTGATGCTCTCCTTTGCCAGGGCTTCTAACATGCAGCTTGCAACACCGGCTTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTc                                                                     <  8:200545/149‑1 (MQ=255)
                                                                                   gagaTCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTcgcgcg                                                                <  8:87948/148‑1 (MQ=255)
                                                                                    agaTCATCTGTATGTTGATGCTCTCCTTGGCCAGGGGTTCGAACATGCAGCTGGCAACAGCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTcgcgcgcg                                                              <  1:130415/149‑1 (MQ=255)
                                                                                      atcatcTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTTGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGgt                                                            <  4:27168/149‑1 (MQ=255)
                                                                                          tcTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGttt                                                        >  1:110078/1‑149 (MQ=255)
                                                                                          tcTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGttt                                                        >  2:241699/1‑149 (MQ=255)
                                                                                           cTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGtttt                                                       >  1:216195/1‑149 (MQ=255)
                                                                                              gATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTc                                                      <  3:254736/147‑1 (MQ=255)
                                                                                              gATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTc                                                      <  7:18715/147‑1 (MQ=255)
                                                                                                 gTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCa                                                 >  6:102937/1‑149 (MQ=255)
                                                                                                      tGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGctc                                            >  8:264759/1‑149 (MQ=255)
                                                                                                        ctctCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGctct                                           <  8:24941/148‑1 (MQ=255)
                                                                                                          ctcCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTc                                                                     <  3:227221/120‑1 (MQ=255)
                                                                                                          ctcCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTc                                                                     >  4:227221/1‑120 (MQ=255)
                                                                                                          ctcCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTgcg                                        >  1:46649/1‑149 (MQ=255)
                                                                                                            ccTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGcc                                      >  5:152828/1‑149 (MQ=255)
                                                                                                             cTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTgcgc                                       <  4:206881/147‑1 (MQ=255)
                                                                                                             cTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTgc                                         <  6:85129/145‑1 (MQ=255)
                                                                                                              ttGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCtt                                    >  7:129612/1‑149 (MQ=255)
                                                                                                                ggCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTtc                                   <  8:176495/148‑1 (MQ=255)
                                                                                                                 gCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCTCCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTtct                                  <  8:249557/148‑1 (MQ=255)
                                                                                                                   cAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCg                                >  4:35998/1‑148 (MQ=255)
                                                                                                                    aGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTa                              >  3:163939/1‑149 (MQ=255)
                                                                                                                     gggCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTAc                             <  1:298869/149‑1 (MQ=255)
                                                                                                                     gggCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTAc                             <  3:171454/149‑1 (MQ=255)
                                                                                                                      ggCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTAc                             >  2:308361/1‑148 (MQ=255)
                                                                                                                        cTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCATTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCa                          <  1:241699/149‑1 (MQ=255)
                                                                                                                           cGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCAtt                        <  2:216195/148‑1 (MQ=255)
                                                                                                                             aaCATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCg                     <  1:308361/149‑1 (MQ=255)
                                                                                                                                  gCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCa                >  8:120598/1‑149 (MQ=255)
                                                                                                                                        ggCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACggtga          <  6:152828/149‑1 (MQ=255)
                                                                                                                                         gCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACggtg           >  3:253673/1‑147 (MQ=255)
                                                                                                                                          cAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACggtgaag        <  6:287949/149‑1 (MQ=255)
                                                                                                                                           aaCACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACggtgaagg       >  8:183433/1‑149 (MQ=255)
                                                                                                                                           aaCACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACggtgaagg       >  4:85347/1‑149 (MQ=255)
                                                                                                                                            gcgtcGGCATGCGTTCGCACGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACggtgaaggt      >  4:109414/5‑149 (MQ=255)
                                                                                                                                            acacCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACggtgaaggt      <  2:46649/149‑1 (MQ=255)
                                                                                                                                             cacCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACggtgaaggtg     <  7:264759/149‑1 (MQ=255)
                                                                                                                                              acCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACggtgaaggtg     >  5:260785/1‑148 (MQ=255)
                                                                                                                                              acCGGCGTGCGAACGCATGCCGACGCAGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGAACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACggtgaaggtg     >  6:269290/1‑148 (MQ=255)
                                                                                                                                                cGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACGgtgaaggtgaag  <  1:153254/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GATCGAACGCGGTATGCAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACGGTGAAGGTGAAG  >  NC_002947/5080332‑5080624

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: ‑