Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,080,486 A→T 8.8% R351R (CGT→CGA)  lysC ← aspartokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,080,4860AT 8.8% 186.0 / 8.3 68R351R (CGT→CGA) lysCaspartokinase
Reads supporting (aligned to +/- strand):  ref base A (30/32);  new base T (3/3);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

GGTATGCAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACGGTGAAGGTGAAG  >  NC_002947/5080342‑5080624
                                                                                                                                                |                                                                                                                                          
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                                                                                                                             <  3:13679/147‑7 (MQ=255)
             cgTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCATGCGTTCGCACgccggtgttg                                                                                                                             >  4:13679/1‑141 (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                                                                                                   <  3:13925/142‑1 (MQ=255)
                                            gCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGgt                                                                                                   >  4:13925/1‑142 (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                             <  3:171454/149‑1 (MQ=255)
                                                                                                           gggCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTAc                             <  1:298869/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       >  4:85347/1‑149 (MQ=255)
                                                                                                                                 aaCACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACggtgaagg       >  8:183433/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)
                                                                                                                                                |                                                                                                                                          
GGTATGCAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACGGTGAAGGTGAAG  >  NC_002947/5080342‑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: ‑