Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,596,083 T→A 5.9% intergenic (‑57/+88) 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,596,0830TA 5.9% 341.6 / 8.4 112intergenic (‑57/+88)PP_4921/thiCNCS1 family transporter/phosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base T (52/50);  new base A (5/2);  total (57/55)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 121 total aligned reads displayed.

GGAAAACCACAGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGC  >  NC_002947/5595934‑5596232
                                                                                                                                                     |                                                                                                                                                     
ggaagtgcgtgagtacgctgccacgatcgaggccgtggatgtgacggtcgagcaaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCtt                                                                                                                                                       >  2:80220‑M1/132‑150 (MQ=255)
ggaagtgcgtgagtacgctgccaagatcgaggccgtggatgtgacggtcgagcaaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCtt                                                                                                                                                       >  2:80450‑M1/132‑150 (MQ=255)
ggaagtgcgtgagtacgctgccaagatcgaggccgtggatgtgacggtcgagcaaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCtt                                                                                                                                                       >  2:80210‑M1/132‑150 (MQ=255)
ggAAAACCACAGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTT                                                                                                                                                       <  2:3713492‑M1/150‑1 (MQ=255)
ggAAAACCACAGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTT                                                                                                                                                       <  2:1123920‑M1/150‑1 (MQ=255)
 gaagtgcgtgagtacgctgccaagatcgaggccgtggatgtgacggtcgagcaaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTg                                                                                                                                                      >  1:4194‑M1/131‑150 (MQ=255)
  aagtgcgtgagtacgctgccaagatcgaggccgtggatgtgacggtcgagcaaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGc                                                                                                                                                     >  1:4517827‑M1/130‑150 (MQ=255)
  aaaaCCACAGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGC                                                                                                                                                     >  2:651578‑M1/1‑150 (MQ=255)
  aaaaCCACAGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGC                                                                                                                                                     >  2:4202663‑M1/1‑150 (MQ=255)
      gcgtgagtacgctgccaagatcgaggccgtggatgtgacggtcgagcaaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCcgc                                                                                                                                                 <  1:4901732‑M1/25‑1 (MQ=255)
      gcgtgagtacgctgccaagatcgaggccgtggatgtgacggtcgagcaaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCAAGCCcgc                                                                                                                                                 >  1:1460538‑M1/126‑150 (MQ=255)
      gcgtgagtacgctgccaagatcgaggccgtggatgtgacggtcgagcaaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCAAGCCcgc                                                                                                                                                 >  1:41007‑M1/126‑150 (MQ=255)
       cacaGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGGTTGCCCGCG                                                                                                                                                >  1:2255902‑M1/1‑150 (MQ=255)
        gtgagtacgctgccaagatcgaggccgtggatgtgacggtcgagcaaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCAAGcccgctc                                                                                                                                               >  2:837486‑M1/124‑148 (MQ=255)
           gTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGcacc                                                                                                                                            >  2:2712130‑M1/1‑146 (MQ=255)
           gTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCAAGCCCGCtccca                                                                                                                                            <  1:1055578‑M1/150‑6 (MQ=255)
               aaaCAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCgg                                                                                                                                        <  2:4285426/150‑1 (MQ=255)
                aaCAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGc                                                                                                                                       >  2:4185235/1‑150 (MQ=255)
                aaCAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCGTGCCCGCGAAcccccgc                                                                                                                                       >  2:1424753‑M1/1‑143 (MQ=255)
                   aGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCaaa                                                                                                                                    <  1:1818101/150‑1 (MQ=255)
                      tCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGcc                                                                                                                                 <  1:2990978/150‑1 (MQ=255)
                      tCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGcc                                                                                                                                 <  1:2990969/150‑1 (MQ=255)
                         cgAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCTTTGCCGGt                                                                                                                              <  1:502342/150‑1 (MQ=255)
                         cgAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCTTTGCCGGt                                                                                                                              <  1:502423/150‑1 (MQ=255)
                         cgAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCAGGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGt                                                                                                                              <  2:2774020/150‑1 (MQ=255)
                         cgAGCACCGAAAATGCGCTGCTGGATGGGGATTGGAGGGGCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGt                                                                                                                              <  1:637255/150‑1 (MQ=255)
                            gCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGcc                                                                                                                           <  2:3862509/150‑1 (MQ=255)
                            gCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGcc                                                                                                                           <  2:1281583/150‑1 (MQ=255)
                                cGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGc                                                                                                                       >  2:1949713/1‑150 (MQ=255)
                                 gAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGGTTTGCCGCCAAAACCGGCAAAGCCGGGGGCATGCa                                                                                                                      >  1:989227/1‑150 (MQ=255)
                                 gAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCTTTGCCGGTgttcgcgg                                                                                                                      <  2:399504/150‑8 (MQ=39)
                                    aaTGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCg                                                                                                                   <  1:3942652/150‑1 (MQ=255)
                                        cgcTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttg                                                                                                               <  2:1726381‑M1/150‑32 (MQ=255)
                                         gctgctGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtg                                                                                                              <  2:2204863/150‑1 (MQ=255)
                                         gctgctGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtg                                                                                                              <  2:4591385/150‑1 (MQ=255)
                                         gctgctGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtg                                                                                                              <  2:583266/150‑1 (MQ=255)
                                         gcggcgggATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCCTGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtg                                                                                                              <  2:583237/144‑1 (MQ=255)
                                            gctgGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtgctt                                                                                                           <  1:1614641/150‑1 (MQ=255)
                                             ggtcgagcaaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCAAgcccgctcc                                                                                                                                              >  1:1312853‑M1/87‑111 (MQ=255)
                                             ctgGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggt                                                                                                          <  2:1188533‑M1/150‑37 (MQ=255)
                                               ggATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCt                                                                                                        >  2:130584‑M2/1‑149 (MQ=255)
                                                   ggggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggcttt                                                                                                                <  1:4546763‑M1/138‑31 (MQ=255)
                                                   ggggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggcttt                                                                                                                >  2:4546763‑M1/1‑108 (MQ=255)
                                                   ggggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTACttcgc                                                                                                    >  2:2304632‑M2/1‑145 (MQ=255)
                                                   ggggATTGGATGGTCGGGAGAGAACGGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgc                                                                                                    >  2:4323630‑M2/1‑145 (MQ=255)
                                                    gggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcg                                                                                                   <  2:3013376‑M1/150‑44 (MQ=255)
                                                    gggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcg                                                                                                   >  1:659008‑M2/1‑144 (MQ=255)
                                                    gggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcg                                                                                                   >  1:3224768‑M2/1‑144 (MQ=255)
                                                    caaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTcgcg                                                                                                   >  1:1025031‑M1/80‑150 (MQ=255)
                                                     ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgg                                                                                                  <  2:1839723‑M1/150‑45 (MQ=255)
                                                     ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgg                                                                                                  <  2:3707504‑M1/150‑45 (MQ=255)
                                                     ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgg                                                                                                  <  2:2333247‑M1/150‑45 (MQ=255)
                                                     ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgg                                                                                                  <  2:2333402‑M1/150‑45 (MQ=255)
                                                     ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgg                                                                                                  <  2:3707601‑M1/150‑45 (MQ=255)
                                                     ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgg                                                                                                  >  1:1627349‑M2/1‑143 (MQ=255)
                                                     ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcga                                                                                                  >  2:4222300‑M2/1‑143 (MQ=255)
                                                     ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCACGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgg                                                                                                  <  2:3878095‑M1/150‑45 (MQ=255)
                                                     aaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                  >  2:3856190‑M1/79‑150 (MQ=255)
                                                     aaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                  >  1:4170545‑M1/79‑150 (MQ=255)
                                                     aaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                  >  2:1008468‑M1/79‑150 (MQ=255)
                                                     aaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                  >  2:2790285‑M1/79‑150 (MQ=255)
                                                     aaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                  >  2:1765056‑M1/79‑150 (MQ=255)
                                                     aaggcatgcgtgagcaggccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                  >  2:4931257‑M1/79‑150 (MQ=255)
                                                      gATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAACCGGTGCCATGCACAGCAGTTCTTGCttcgcggg                                                                                                 <  2:742020‑M2/150‑9 (MQ=255)
                                                         tGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGTTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaacggctttgccggtgttcgcgggcaa                                                                                              <  2:495267‑M1/150‑49 (MQ=255)
                                                         tGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgggcaa                                                                                              <  2:495451‑M1/150‑49 (MQ=255)
                                                         tGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCAGACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaa                                                                                              >  1:36112‑M2/1‑139 (MQ=255)
                                                                            tggctggGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCCCAGTGCGTGCTTCGCGGGCAAGCCCGCGCCCACAGGGTaa                                                                           >  2:1419449/1‑150 (MQ=255)
                                                                                      aTGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCCGCAAAGCCCGTGCCCAGCACCGCAGTGCGTGCTTCGCGGGCAAGCCCGCACCCACACGGTAAGGTGCTGCgg                                                                 >  2:3490275/1‑150 (MQ=37)
                                                                                      aTGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCGTGCCCGCGCACACCGGCAAAGCCCGTGCCATGCACCGCCGGGCTGGCTTAGCGGGCCAGCCCGCTCCCCCAGGGGAAGGTGCTGCgg                                                                 >  1:2849019/1‑150 (MQ=25)
                                                                                        gtcatggtcatgTCCGCAGGAGTGATGTTAAGTGGCAGACACCACCCTGTGGGTGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGGGCCATGCACCGCCGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGtt                                                               <  2:4435303/150‑1 (MQ=37)
                                                                                          catggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTa                                                             <  2:2962334/150‑1 (MQ=255)
                                                                                          catggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTa                                                             <  2:5011642/150‑1 (MQ=255)
                                                                                           atggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAACCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAg                                                            >  2:4901732/1‑150 (MQ=255)
                                                                                             ggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCCCGAACCCCGGCAAACCCGGTGCCCAGCACCGCAGTGCGTGCGTCGCGGGCCAGCCCGCGCCCACAGGGGAAGGGGCGGCGGGGTAGAc                                                          >  2:1025527/1‑150 (MQ=37)
                                                                                               tcatgTCCGCTGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGTGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCCCCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACtt                                                        <  2:2688508/150‑1 (MQ=255)
                                                                                                catgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGGGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACttt                                                       <  2:3138220/150‑1 (MQ=255)
                                                                                                 atgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTtg                                                      >  2:135937/1‑150 (MQ=255)
                                                                                                   gTCGGCAGGATTGATGTTAAGTGGCATACACCACCCTGTGGGCGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTtgtg                                                    <  2:3873626/150‑1 (MQ=255)
                                                                                                      cGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTa                                                 <  2:3814213/150‑1 (MQ=255)
                                                                                                                aTGTTAAGTGGCATACACCACCCTGTGGGAGCGGGGTTGGCCGCGCACACCGGCAAAACCCGGGCCATGCACCGCCGGGCGTGCGTCGCGGGCCAGCCCGCTCCCCCCCGGTAAGGGGCTGCGGGTTTGACGTTGTGGTACAgcgggcgg                                       >  1:1137022/1‑148 (MQ=255)
                                                                                                                aTGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCAAGCCCGCGAAAACCCGCAAAGCCCGTGCCCAGCACCCCACTGCTTGCGTCGCCGGCAACCCCCCTCCCACCGGGTGAGGGGCGGCGGTTTAGACGTTGTGGGACAgcgggcag                                       >  2:1404515/1‑148 (MQ=255)
                                                                                                                     aaGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTc                                  >  2:408378/1‑150 (MQ=255)
                                                                                                                     aaGCGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCCCAGCCGGTGCAATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTc                                  <  2:2633629/150‑1 (MQ=255)
                                                                                                                       ggggCATACACCCCCCTGGGGGAGCGGGCTGGCCCGCGACCACCGGCAACGCCGGTGCCATGCCCCGCAGGGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCt                                <  2:3208169/148‑1 (MQ=37)
                                                                                                                       gTGGCATACACCCCCCTGTGGGTGGGGGGTTGCCCGCGACCACCGGCAAAGCCGGTGCCATGCACGGCATTGGTTGCTTCGCGGGCACGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCt                                <  1:4134613/150‑1 (MQ=25)
                                                                                                                          gaacacacCACCGTGGGGGCGGGGGCTGGCCCGCCACCCCCGCCAACGCGGGGGCCATCCCCGGCGGTGCTTGCTTCGCGGCCACCCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGc                             <  1:135937/146‑1 (MQ=255)
                                                                                                                          gCATACACCACCCTGTGGTCGCGGGCTTGCCCGCGAACACCGGAAAAGCCGGTGCCATGCACCGCCGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGc                             <  1:742020/150‑1 (MQ=255)
                                                                                                                          gCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCCGTGCGTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGc                             >  2:4944561/1‑150 (MQ=255)
                                                                                                                          gCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggca                                                                                               >  2:4266178‑M2/1‑74 (MQ=255)
                                                                                                                            aTACACCCCCCTGGGGGCGCGGGCTGGCCCGCCACCACCGGCAACGCCGGGGCCATCCACCGACGGGCGGGCTCCGCGGCAACCCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCgg                           <  1:3021377/150‑1 (MQ=255)
                                                                                                                                     CCTGTGGGAGCGGGCTTTCCCCCCAAAAACGGCAAAACCCGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCACCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCAGGCTGCCTTCCTGGCGGAAccgcacg                  <  2:1583552‑M1/150‑4 (MQ=255)
                                                                                                                                          GGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCTTGCACCGCTGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTg             >  2:3608031‑M1/1‑150 (MQ=255)
                                                                                                                                            GAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCt           <  2:4923357‑M1/150‑1 (MQ=255)
                                                                                                                                               CGGGCTTGCCCGCGAACACCGGCAACGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCAc        >  1:4185235‑M1/1‑150 (MQ=255)
                                                                                                                                                    ttGCCCGCGAACACCGGCAAAGCCGTTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaacatgacatcacccagccagttctctcccgaccatcca   <  2:495267‑M2/150‑103 (MQ=255)
                                                                                                                                                    ttGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatcca   <  2:495451‑M2/150‑103 (MQ=255)
                                                                                                                                                    ttGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatcca   <  1:4176813‑M2/150‑103 (MQ=255)
                                                                                                                                                     tGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgggcaagcccgctcccacagggtggtgtatgc                                                                    <  1:4266178‑M1/84‑75 (MQ=255)
                                                                                                                                                     TGCCCGCGAACACCGGCAACGCGGGTGCCATGCACCGCTGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGc  >  2:2255902‑M1/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGAAAACCACAGTGAAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGC  >  NC_002947/5595934‑5596232

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: