Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,594,150 T→C 8.1% intergenic (‑97/+48) 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,1500TC 8.1% 193.5 / 7.3 74intergenic (‑97/+48)PP_4921/thiCNCS1 family transporter/phosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base T (35/33);  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 = 1.76e-01

GCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCA  >  NC_002947/5594002‑5594297
                                                                                                                                                     |                                                                                                                                                    
gctgctGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGt                                                                                                                                                      >  4:78150/1‑149 (MQ=255)
gctgctGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGt                                                                                                                                                      <  4:65010/149‑1 (MQ=255)
     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                                                                      <  2:337308‑M2/103‑65 (MQ=255)
                                                                                                                             aaaGCCGGTGCCATGCACCGCGGCGCCTGTTCGCGGGCAcgcccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgtttttt                                                                      >  1:337308‑M2/1‑39 (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                               <  2:242969/129‑1 (MQ=255)
                                                                                                                                          tGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGa                               >  1:242969/1‑129 (MQ=255)
                                                                                                                                          tGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCACTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCg           <  7:17933/149‑1 (MQ=255)
                                                                                                                                            gCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAAcc                                                             <  5:192659/98‑1 (MQ=255)
                                                                                                                                            gCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAAcc                                                             >  6:192659/1‑98 (MQ=255)
                                                                                                                                              aCCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATc        >  2:111768/1‑149 (MQ=255)
                                                                                                                                                  cattgcttgctTCGCGGCCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTgg    <  2:354857/146‑1 (MQ=255)
                                                                                                                                                    gTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCa  <  5:40897/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCA  >  NC_002947/5594002‑5594297

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: