Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,594,148 A→G 8.2% intergenic (‑95/+50) 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,1480AG 8.2% 188.2 / 6.9 77intergenic (‑95/+50)PP_4921/thiCNCS1 family transporter/phosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base A (36/31);  new base G (3/3);  total (40/37)
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.90e-01

GCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATC  >  NC_002947/5594002‑5594291
                                                                                                                                                   |                                                                                                                                                
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                                                                >  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                         <  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)
                                                                                                                                                   |                                                                                                                                                
GCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATC  >  NC_002947/5594002‑5594291

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: