Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,594,156 C→T 5.5% intergenic (‑103/+42) 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,1560CT 5.5% 202.0 / 2.2 74intergenic (‑103/+42)PP_4921/thiCNCS1 family transporter/phosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base C (35/33);  new base T (2/2);  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.90e-01

GATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCGTA  >  NC_002947/5594015‑5594302
                                                                                                                                              |                                                                                                                                                   
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                                  >  6:252537/1‑148 (MQ=255)
                                                                                                            cGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGAcc                                  <  5:252537/148‑1 (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)
                                                                                                                                             gcttCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACTCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCGTa  <  1:91471/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCGTA  >  NC_002947/5594015‑5594302

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: