Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,594,153 T→C 8.3% intergenic (‑100/+45) PP_4921 ← / ← thiC NCS1 family transporter/phosphomethylpyrimidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,594,1530TC 8.3% 191.9 / 7.5 74intergenic (‑100/+45)PP_4921/thiCNCS1 family transporter/phosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base C (3/3);  total (38/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 = 3.64e-01

TGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGC  >  NC_002947/5594007‑5594299
                                                                                                                                                   |                                                                                                                                                   
tgGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCGGTGCCATGCACCGCAGtgcttg                                                                                                                                                   >  4:112127/1‑149 (MQ=255)
  gATGGGGATTGGATGGTCGTGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGTCATGCACCGCAGtgcttg                                                                                                                                                   <  6:381744/147‑1 (MQ=255)
        gATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                                                           <  2:298675/149‑1 (MQ=255)
             gATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAg                                                                                                                                      >  2:40412/1‑149 (MQ=255)
                    gggAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTc                                                                                                                               <  3:78150/149‑1 (MQ=255)
                     ggAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGCACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTAGCGGGCAAGCCCGCTcc                                                                                                                              >  2:297174/1‑149 (MQ=255)
                         ccaaCTGGCTGGGGCATGGCATGGTCATGTCCGCAGGAGTGACGTTACGGGGCATACACCACCCGGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCcaca                                                                                                                          <  5:316070/147‑1 (MQ=25)
                            tgtgggtggTTGCTGTAATGGTCATGCCCGCAGTTGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCTAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAgg                                                                                                                        <  8:182270/146‑1 (MQ=21)
                            aCTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGGGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAggg                                                                                                                       <  7:293234/149‑1 (MQ=255)
                              tggctggGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATTCAACGGAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAgggg                                                                                                                      >  8:263684/1‑147 (MQ=255)
                                  tggGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAAACCTGCAAATCCGGTTCCATGAAACGCAGTACTTGCTTAGCCGGCAAGCCCGCTCCCACACGGTAAGGt                                                                                                                 >  7:364049/1‑149 (MQ=37)
                                         tGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGCACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCCCGGGCAAGCCCGCTCCCACAGGGTAAGGTGGTGCgg                                                                                                          >  5:213378/1‑149 (MQ=255)
                                         tGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCACGCTCCCACAGGGTAAGGTGCTGCgg                                                                                                          >  7:187692/1‑149 (MQ=255)
                                         tGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGACAGCGAACACCGGCAAAGCCGGTGCCATGCACCGCACTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCgg                                                                                                          >  7:182629/1‑149 (MQ=255)
                                          gtcatggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGATCCCACAGGGTAAGGTGCTGCGGt                                                                                                         >  4:286016/1‑149 (MQ=255)
                                           gaatgtccatgTCCGACGTAGTCATTTTAAGTGGCATACACCACCCTGGGGGAGCGGGCTGGCCCGCTAACACCGGCAAAGCCGGTGCCATGCACCGCTGTGCTTGCTTCGCGGGAAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGtt                                                                                                        <  8:187692/142‑1 (MQ=21)
                                                      tCCGCAGGAGTGTTGTTAAGGGGCATACACCACCCTGGGGGAGCGGGCTGGCCCGCGAACACCGGCAAAGCCGGTGCCATGCCCCGCAGTGCTTGCTGCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTtgt                                                                                              <  1:344487/148‑1 (MQ=37)
                                                        cGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTgg                                                                                            >  1:260016/1‑148 (MQ=255)
                                                             gAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTTCCCGCGAACACCGGCCAAGCCCGTGCCATGCACCGCAGTGCTTGCTTAGCCGGCAAGCCCGCAACCCCAAGGTAAGGTGCGGCGGTTTTGACGTTGTGGTACAgc                                                                                      >  4:383130/1‑149 (MQ=25)
                                                                tGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCCCTCCCACAgg                                                                                                                        >  1:340273/1‑112 (MQ=255)
                                                                 ttttttAAGTGCCATACACCACCTTGTGGGAGGGGGCTTGCCCGCTAACACGGTCAATGCGGGTGCATTGCACCGCAGTGCTTGCTTCGCGGGAAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAgctgg                                                                                   <  4:145100/144‑1 (MQ=255)
                                                                  aTGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGGAAAGCCGGTGCCATGCACCGCAGTGGTTGCTTCGCGGGCCAGCCCGCTCCCCCAGGGTAAGGTGCTGCGGGTTAGACTTTGTGGTACAgctggct                                                                                 >  1:354857/1‑149 (MQ=37)
                                                                    gTTAAGTGGCATACACCACCCTGTGGGAGCGGGCGTTCCCCCCCCCACCCGCCAAGCCGGTGGCCTGGACCCGCGTGCATGCATAGCGGGCCACCCCGCGCCCACAGGGGAAGGTGCTGCGGGTTGGAAGTTGTGGTTCACCAGGCTGc                                                                               >  8:168933/1‑149 (MQ=255)
                                                                       aaGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCt                                                                             <  2:260016/148‑1 (MQ=255)
                                                                        aGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTc                                                                           >  4:312125/1‑149 (MQ=255)
                                                                         gTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGc                                                                                                                                         <  8:181752/86‑1 (MQ=255)
                                                                          tGGCATACACCACCGTGTGGGTGCGGGCTTGCCCGCTAACACCGGAAAAGCCGGTGCCATGCACCGCAGTGCTGGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCt                                                                         <  7:168933/149‑1 (MQ=37)
                                                                            gCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCCGCAAAGCCCGTGCCCTGCACCGCAGTGCTTTCTTAGCGGGCAAGCCCGCGACCACCGGGTAAGGTGCTGCGGttt                                                                                                       >  6:116274/1‑117 (MQ=25)
                                                                            gCATACACCACCCTGTGGGAGCGGGCTGGCCCGCGAACACCGGCAACGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGttt                                                                                                       <  5:116274/117‑1 (MQ=255)
                                                                             cATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGc                                                                     <  3:112127/149‑1 (MQ=255)
                                                                              aTACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCg                                                                    >  1:56148/1‑149 (MQ=255)
                                                                                 caccacCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGTGGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGaa                                                                 <  8:166724/149‑1 (MQ=255)
                                                                                 caccacCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTGCCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGaa                                                                 <  1:297174/149‑1 (MQ=255)
                                                                                  accaccCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAAc                                                                <  7:263684/149‑1 (MQ=255)
                                                                                    caccCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCg                                                              >  2:91471/1‑149 (MQ=255)
                                                                                     accCTGTGTGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAAGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGc                                                             >  6:182946/1‑149 (MQ=255)
                                                                                            gggAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGcc                                                       >  1:74693/1‑149 (MQ=255)
                                                                                              gtgCGGGCTTGCCCGCTAACACCGGCAAAGCGGGTGCCAGGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTATGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGc                                                    <  5:128721/147‑1 (MQ=255)
                                                                                              gAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGc                                                    <  5:355306/149‑1 (MQ=255)
                                                                                               aGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCt                                                   >  6:280755/1‑149 (MQ=255)
                                                                                                    gCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACg                                               >  6:40897/1‑148 (MQ=255)
                                                                                                        gCCCGCGAACACGGGCAACCCGGGGGCATTCCACCGCTGTGCTGGCTTGGCGGTCATGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGc                                          <  8:160625/149‑1 (MQ=255)
                                                                                                          ccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCGGCGCCTGTTTCGCGGGCAcgcccgctcccaccgactgattgccattaagcgggcctatcgctacagtgacggcccgttttttgttgcctcgctaacgaatcatgactccct                                           <  5:78080‑M2/147‑94 (MQ=255)
                                                                                                          ccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCATTCGTGGCGGAACCGATCGGCCTGATCACGCATGCCt                                        >  1:195689/1‑149 (MQ=255)
                                                                                                           cgcgAACACCGGCAAAGCCGGTGCCATGCACCGCGGCGCCTGTTTCGCGGGCAcgcccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgttttttgttgcctcgctaacgaatcatgactccct                                           <  3:249817‑M2/146‑94 (MQ=255)
                                                                                                           cgcgAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCGGCCTTCCTGTGCGGAACCGCTCGGCCTGCTCACGCATGCCt                                        <  1:40412/149‑1 (MQ=255)
                                                                                                             ctaaCACCGGCAAAGCCGGTGCCATGCCCCGCAGTGCTTGCTTCGCGGGAAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGc                                     <  6:213378/147‑1 (MQ=255)
                                                                                                               aaCACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTc                                   <  5:280755/149‑1 (MQ=255)
                                                                                                                 cacCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCgc                                 >  4:64515/1‑148 (MQ=255)
                                                                                                                  acCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGTGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGAc                                >  8:150741/1‑149 (MQ=255)
                                                                                                                    cGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGc                                          <  5:59676/137‑1 (MQ=255)
                                                                                                                    cGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGc                                          >  6:59676/1‑137 (MQ=255)
                                                                                                                    cGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGAcc                               <  5:252537/148‑1 (MQ=255)
                                                                                                                    cGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGAcc                               >  6:252537/1‑148 (MQ=255)
                                                                                                                    cGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCg                              >  2:306326/1‑149 (MQ=255)
                                                                                                                     ggCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTGTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACAGCTCGGCCTGCTCACGCATGCCTTGCTCGAcagg                             >  3:44277/1‑146 (MQ=255)
                                                                                                                        aaaGCCGGTGCCATGCACCGCGGCGCCTGTTCGCGGGCAcgcccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgtttttt                                                                        >  1:337308‑M2/1‑39 (MQ=255)
                                                                                                                        aaaGCCGGTGCCATGCACCGCGGCGCCTGTTCGCGGGCAcgcccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgtttttt                                                                        <  2:337308‑M2/103‑65 (MQ=255)
                                                                                                                         aaGCCGGTGCCATGCACCGCCGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTcaca                         <  3:286016/149‑1 (MQ=255)
                                                                                                                         aaGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTcaca                         <  8:244242/149‑1 (MQ=255)
                                                                                                                            ccGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATc                       >  4:252048/1‑148 (MQ=255)
                                                                                                                             cGGTGCCATCCACCGCAGTGCGTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCa                     <  5:318569/149‑1 (MQ=255)
                                                                                                                               gTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCACTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACg                   >  8:17933/1‑149 (MQ=255)
                                                                                                                                tGCCATGCACCGCGGCGCCTGTTTCGCGGGCAcgcccgctcccaccgactgattgccattaagcgggcctctcgctacagtgacggcccgttttttgttgcctcgctaacgaatcatgactcccttagaacgctatcaagcagatctga                   >  6:78080‑M2/1‑32 (MQ=255)
                                                                                                                                tGCCATGCACCGCGGCGCCTGTTTCGCGGGCAcgcccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgttttttgttgcctcgctaacgaatcatgactcccttagaacgctatcaagcagatctga                   >  3:136600‑M2/1‑32 (MQ=255)
                                                                                                                                 gCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGc                 <  7:161597/149‑1 (MQ=255)
                                                                                                                                  ccATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGcc                <  1:306326/149‑1 (MQ=255)
                                                                                                                                   cATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGCACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCt               >  4:299736/1‑149 (MQ=255)
                                                                                                                                   cATGCACAGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCt               >  7:257389/1‑149 (MQ=255)
                                                                                                                                     tGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGa                                 >  1:242969/1‑129 (MQ=255)
                                                                                                                                     tGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGa                                 <  2:242969/129‑1 (MQ=255)
                                                                                                                                     tGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCACTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCg             <  7:17933/149‑1 (MQ=255)
                                                                                                                                       gCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAAcc                                                               >  6:192659/1‑98 (MQ=255)
                                                                                                                                       gCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAAcc                                                               <  5:192659/98‑1 (MQ=255)
                                                                                                                                         aCCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATc          >  2:111768/1‑149 (MQ=255)
                                                                                                                                             cattgcttgctTCGCGGCCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTgg      <  2:354857/146‑1 (MQ=255)
                                                                                                                                               gTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCa    <  5:40897/149‑1 (MQ=255)
                                                                                                                                                  cttgcttCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACTCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGc  >  5:292796/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGC  >  NC_002947/5594007‑5594299

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: