Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,594,153 T→C 5.6% 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 5.6% 201.6 / 2.1 73intergenic (‑100/+45)PP_4921/thiCNCS1 family transporter/phosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base T (29/38);  new base C (2/2);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

GATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCG  >  NC_002947/5594009‑5594300
                                                                                                                                                |                                                                                                                                                    
gATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCt                                                                                                                                                  <  6:214056‑M2/149‑2 (MQ=255)
  tGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtgcttgcccg                                                                                                                                                >  2:77929/1‑146 (MQ=255)
  tGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttc                                                                                                                                                <  8:247565‑M2/149‑4 (MQ=255)
  tGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttc                                                                                                                                                <  8:252768‑M2/149‑4 (MQ=255)
     ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcg                                                                                                                                             >  3:85033‑M2/1‑143 (MQ=255)
     ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcg                                                                                                                                             >  7:247565‑M2/1‑143 (MQ=255)
     ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcg                                                                                                                                             >  7:252768‑M2/1‑143 (MQ=255)
      aggcatgcgtgagcaggccgagcggttcctccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                                                            >  8:336530‑M1/78‑149 (MQ=255)
      aggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                                                            >  4:355137‑M1/78‑149 (MQ=255)
       ggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCggg                                                                                                                                           >  6:343034‑M1/77‑149 (MQ=255)
        ttGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggc                                                                                                                                          >  5:214056‑M2/1‑140 (MQ=255)
            aTGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGTAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagc                                                                                                                                      >  1:234533‑M2/1‑136 (MQ=255)
                tCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttagcgggcaagcccgc                                                                                                                                  <  1:232543‑M2/149‑18 (MQ=255)
                   ggAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcc                                                                                                                               <  6:149547‑M2/149‑21 (MQ=255)
                   ggAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcc                                                                                                                               <  6:159542‑M2/149‑21 (MQ=255)
                    gagagaACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTTCCCGCGAACACCGGCAAATCCGGTGCCATTCCCCGCAGTGCTTGCttcgcgggcaagcccgctcc                                                                                                                               <  7:54158‑M2/148‑21 (MQ=255)
                    gagagaACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttagcgggcaaccccgctcc                                                                                                                               <  1:326043‑M2/148‑21 (MQ=255)
                        cgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGc                                                                                                                                      <  7:337672‑M1/78‑1 (MQ=255)
                        cgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGc                                                                                                                                      >  8:337672‑M1/60‑137 (MQ=255)
                        caaCTGGCTGGGTTAGGACTTGTCCAGGCCCGCAGGAGGGATGTTAAGGGGCTTACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACCCCGCCAAACCCGGTGCCATGCACCGACGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAg                                                                                                                         <  8:272028/148‑1 (MQ=18)
                         gagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGcc                                                                                                                                     >  4:322295‑M1/59‑137 (MQ=255)
                         gagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGcc                                                                                                                                     <  3:322295‑M1/79‑1 (MQ=255)
                         ccgatcttgggggatgtCTGGGTCATGTCCGCCGGTGTTATTTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCTAACACCGGAAAAGCCGGTGCCATGCACCGCAGTGCTTGCtgcgcgggcaagcccgctcccacag                                                                                                                          >  2:12151‑M2/8‑123 (MQ=255)
                                 gggTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccccag                                                                                                                          <  1:12151‑M2/140‑26 (MQ=255)
                                 gggTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCCAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTGCGCGGGCAAGCCCGCTCCCACAGGGTAAGGtg                                                                                                                <  3:318133/149‑1 (MQ=255)
                                  ggTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCTAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGtg                                                                                                                <  2:254654/148‑1 (MQ=255)
                                      aTGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTTCCCGCGAACAACCGCAAAGCCGGTGCCATGCACCGCAGTGCGTGCTTAGCCGGCAACCCCGCTCCCCCAGGGTAATGTgcggcg                                                                                                           >  1:153349/1‑149 (MQ=37)
                                        gtcatggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCATGCCCCCGAACCCCCGCAAAACCGGTGCCATGCACCACAGTGCATGCTTAGCGGGCAAGCCCGCTCCCACAGGGTAAGGTTCAGCGGt                                                                                                         >  4:134927/1‑149 (MQ=25)
                                         tcatggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGtt                                                                                                        <  1:61833/149‑1 (MQ=255)
                                          catggtcatgGCCGCAGGAGTCATGTTACGGGGCATACACCACCCGGGGGCAGCTGGCTTGCCCGCGAACACCGCCAACGCCGGTGCCCTGCACCGCAGTGCTGGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGttt                                                                                                       <  7:96407/149‑1 (MQ=25)
                                            tggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCCCCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCCCAGGGTAAGGTGCTGCGGTTTAg                                                                                                     >  5:84178/1‑149 (MQ=255)
                                                 atgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCCCCGCAGTGCTTGCTTCGCGGGCAAGCCCGCGCCCACAGGGTAAGGTGCTGCGGGGTAGACttt                                                                                                >  3:275737/1‑149 (MQ=255)
                                                 atgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTTCTTCGCGGGCAAGCCCGCTCCCACAGGGGAAGGTGCTGCGGTTTAGCCtgt                                                                                                >  3:280530/1‑147 (MQ=255)
                                                 atgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACttt                                                                                                >  1:12170/1‑149 (MQ=255)
                                                      cGCCGGAGGGATGTTAAGCGCCATCCACCACCCGGTGGGCGCGGGCTTCCCCGCGAACACCGGCAACCCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTgg                                                                                            <  2:48096/148‑1 (MQ=25)
                                                                 tgtaaggTGCCATACACCACCGTGTGGGAGCGGGCTTCCCGGCTAACACCGGCAAAGCGGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAgctggctg                                                                                <  3:134927/143‑1 (MQ=37)
                                                                   ttAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGGCCGCGAACACCGGCAAAGCCGGTGGCCTGCACCGCCGTGCTTGCGTCGCGGGGCAGCCCGCTCCCCCAGGGGAAGGTGGGGGGGGTTTGACGTTGTGGGACAgctggcgggc                                                                              >  4:271454/1‑145 (MQ=21)
                                                                     aaGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCt                                                                             <  5:149547/148‑1 (MQ=255)
                                                                      cgTGGCATCCCCCACCCTGTGGGAGCGGGCTTGCCCGCTAACACCGCCAAAGCCGGGGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCtt                                                                            <  3:74289/147‑1 (MQ=37)
                                                                           cATACACCACCCTGGGGGAGCGGGCTTGCCCGCTAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGc                                                                      <  2:11531/149‑1 (MQ=255)
                                                                            aTACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttagccggcaagcccgctcccacagggtggtgtatgccacttaacatcacaccggcggacatgaccaagacaacacc                                                                      <  4:356022‑M2/149‑78 (MQ=255)
                                                                                   ACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGc                                                              >  2:12170‑M1/1‑149 (MQ=255)
                                                                                     ccTGTGGGAGCGGGCTTGGCCGCGAACACCGGCAAAGCCGGGGCCATGCCCCGCAGGGGTTGCTTCGCGGGGAAGCCCGGGCGCCCAGGGTGAGGTGGGGGGGTTTTGGCTTTGTGGTAGCGCTGGCGGCCTTCCTGGCGGCGCCGCTc                                                            >  6:81441/1‑149 (MQ=16)
                                                                                      CTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCg                                                           >  5:159542‑M1/1‑149 (MQ=255)
                                                                                      CTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCGCCCACAGGGTAAGGTGGGGCGGTTTAGACTTTGTGGTACAGGTGGCGGCCTTACTGGCGGAACCGCTCg                                                           <  2:75834‑M1/149‑1 (MQ=255)
                                                                                          GGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCt                                                       <  2:8041‑M1/149‑1 (MQ=255)
                                                                                               CGGGCTTGCCCGCTAACACCGGCAAAGCGGGTGCCTTGCACCGCAGTGATTGCTTGTCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCtct                                                  >  2:326043‑M1/1‑148 (MQ=255)
                                                                                                 GGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACg                                                <  3:51633‑M1/149‑1 (MQ=255)
                                                                                                  GCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGc                                               >  8:260431‑M1/1‑149 (MQ=255)
                                                                                                      gCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaa                                            <  6:117689‑M2/149‑104 (MQ=255)
                                                                                                      gCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaa                                            >  3:356022‑M2/1‑46 (MQ=255)
                                                                                                       cccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctc                                                        <  5:47142‑M2/136‑92 (MQ=255)
                                                                                                       cccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctc                                                        >  6:47142‑M2/1‑45 (MQ=255)
                                                                                                        ccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCGGCGCCTGTTTCGCGGGCAcgcccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgttttttgttgcctcgctaacgaatcatgactccct                                            <  7:134574‑M2/147‑94 (MQ=255)
                                                                                                        CCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCt                                         <  2:50424‑M1/149‑1 (MQ=255)
                                                                                                         CGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCtt                                        <  5:27939‑M1/149‑1 (MQ=255)
                                                                                                         CGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCtt                                        <  5:22509‑M1/149‑1 (MQ=255)
                                                                                                                acCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTgg                                                                       >  2:148885/1‑111 (MQ=255)
                                                                                                                acCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTgg                                                                       <  1:148885/111‑1 (MQ=255)
                                                                                                                  cGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCCCAGGGTAAGGTGCTGCGGGTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCGCCCGCATGCCTTGc                                      >  5:331920/1‑142 (MQ=255)
                                                                                                                     cAAAGCCGGTGCCATGCACCGCGGCGCCTGTTTCGCGGGCAcgcccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgttttttgttgcctcgctaacgaatcatgactcccttagaacgctatcaag                             >  8:134574‑M2/1‑41 (MQ=255)
                                                                                                                      aaaGCCGGTGCCATGCACCGCGGCGCCTGTTCGCGGGCAcgcccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgttttttgttgactcgctaacgaataatgaatccctaagaacgctataaagc                            <  8:120605‑M2/148‑110 (MQ=255)
                                                                                                                           cGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCa                      >  7:287853/1‑149 (MQ=255)
                                                                                                                            ggTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCAGCAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTg                                       >  1:269972/1‑132 (MQ=255)
                                                                                                                            ggTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTg                                       <  2:269972/131‑1 (MQ=255)
                                                                                                                               gCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGa                                  >  2:31358/1‑133 (MQ=255)
                                                                                                                               gCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGa                                  <  1:31358/133‑1 (MQ=255)
                                                                                                                                ccATGCACCGCGGCGCCTGTTCGCGGGCAcgcccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgttttttgttgcctcgctaacgaatcatgactcccttagaacgctatcaagcagatctgaaa                  >  7:120605‑M2/1‑29 (MQ=255)
                                                                                                                                  aTGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTc               <  5:81441/149‑1 (MQ=255)
                                                                                                                                   tGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCCCGGCCTc               >  5:371286/1‑148 (MQ=255)
                                                                                                                                          ccgTGCTTGCTTCGCGGGCAAGCCCGCTCCCCCAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTgg       <  4:275737/147‑1 (MQ=255)
                                                                                                                                          cAGTGCTTGCTTCGCGGTCAAGCCCGCTCCCACAGGGTAAGTTGCTGCGGGTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTgg       <  4:280530/149‑1 (MQ=255)
                                                                                                                                           aGTGCTGGCTGCGCGGGCAAGCCCGCTCCCACAGGGTAAGGGGCTGCGGTTTAGACTTTGTTGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGc      <  7:63527/149‑1 (MQ=255)
                                                                                                                                            gTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCa     <  6:371286/149‑1 (MQ=255)
                                                                                                                                             tgcttgctTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAg    <  4:51633/149‑1 (MQ=255)
                                                                                                                                               cttgcttCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGc   <  2:232543/148‑1 (MQ=255)
                                                                                                                                               cttgcttCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCg  >  6:68699/1‑149 (MQ=255)
                                                                                                                                                ttgcttCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCg  >  8:184118/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCG  >  NC_002947/5594009‑5594300

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: