Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,080,491 T→C 7.3% M350V (ATG→GTG)  lysC ← aspartokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,080,4910TC 7.3% 194.3 / 4.6 68M350V (ATG→GTG) lysCaspartokinase
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base C (2/3);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

TATGCAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACGGTGAAGGTGAAGTCGGTGGTGTTA  >  NC_002947/5080344‑5080636
                                                                                                                                                   |                                                                                                                                                 
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                                                                                                               <  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                                         <  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                   >  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)
                                                                                                                                                gCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACGGTGAAGGTGAAGTCGGTGGTGTTa  <  5:102937/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TATGCAGCGCGCGTACGGCCAGCTCCAGGTACTTCTCTTCGAGCACCACCGAAACCTTGATCTCGGAAGTGGAGATCATCTGGATGTTGATGCTCTCCTTGGCCAGGGCTTCGAACATGCAGCTGGCAACACCGGCGTGCGAACGCATGCCGACGCCGACGATCGAGACCTTGGCGATCTTGGTGTCGCCGATCACTTCACGGGCACCGATTTCGCGCGCGGTGTTTTCCAGCACGCTCTGCGCCTTCTCGTACTCATTGCGGTGCACGGTGAAGGTGAAGTCGGTGGTGTTA  >  NC_002947/5080344‑5080636

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