Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,560,601 T→A 5.2% intergenic (‑80/+94) PP_4887 ← / ← PP_4888 hypothetical protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,560,6010TA 5.2% 715.9 / 8.8 252intergenic (‑80/+94)PP_4887/PP_4888hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (97/131);  new base A (14/4);  total (116/136)
Fisher's exact test for biased strand distribution p-value = 5.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

CCGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCA  >  NC_002947/5560453‑5560747
                                                                                                                                                    |                                                                                                                                                   
ccGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTTTGTAGCGAAAGGGCCGCATa                                                                                                                                                    >  2:679231/1‑150 (MQ=255)
ccGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTTTGGGGGGAAAGGGGCGCAtt                                                                                                                                                    >  1:132761/1‑149 (MQ=255)
  gTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGc                                                                                                                                                  <  1:6440669/150‑1 (MQ=255)
  gTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGc                                                                                                                                                  <  2:6965233/150‑1 (MQ=255)
  gTACGCTGGCACTGGGCAGGACATTCAGGGCGGCCGGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTGTGTCATGTGTCGCCGACAGCATAGCAACCCCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATCGc                                                                                                                                                  <  2:2606188/150‑1 (MQ=255)
   gaCGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCg                                                                                                                                                 <  2:5007538/149‑1 (MQ=255)
      gCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCTATGCGGcc                                                                                                                                              <  1:6127612/150‑1 (MQ=255)
      gCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGcc                                                                                                                                              >  2:3176857/1‑150 (MQ=255)
        tGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGcccc                                                                                                                                            <  2:5612243/150‑1 (MQ=255)
          tCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCgg                                                                                                                                          <  2:2361320/150‑1 (MQ=255)
          tCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCgg                                                                                                                                          <  2:2849713/150‑1 (MQ=255)
              tgtgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCact                                                                                                                                      >  2:7415139/1‑148 (MQ=255)
                  cAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCCGCAAtttac                                                                                                                                  >  2:7359637/1‑148 (MQ=255)
                      aaaTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCgt                                                                                                                              <  2:784908/150‑1 (MQ=255)
                          tCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGAATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTAt                                                                                                                          <  1:7023109/150‑1 (MQ=255)
                           cAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATg                                                                                                                         <  2:4779286/150‑1 (MQ=255)
                            aGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGc                                                                                                                        <  2:1785724/150‑1 (MQ=255)
                            aGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCATTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGc                                                                                                                        <  1:1617331/150‑1 (MQ=255)
                            aGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGATAAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGc                                                                                                                        <  1:7280099/150‑1 (MQ=255)
                              gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGAAGGGCCGCATAGCGGCCCCCGCCAGTTCCGCGGGTTTGcaa                                                                                                                      >  2:141644/1‑150 (MQ=39)
                              gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGcaa                                                                                                                      <  1:3599036/150‑1 (MQ=255)
                              gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGAATAGCAACACCGCTCCGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGcaa                                                                                                                      <  1:3149804/150‑1 (MQ=255)
                                      gggCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                                              <  1:446857/150‑1 (MQ=255)
                                       ggCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                                             <  1:823750/150‑1 (MQ=255)
                                       ggCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                                             <  2:2037018/150‑1 (MQ=255)
                                       ggCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGAATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                                             <  1:1526613/150‑1 (MQ=255)
                                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCCGCCCCCGCAATTTACCCGTGTATGGAACCAAACAggg                                                                                                            >  2:2727051/1‑150 (MQ=38)
                                          aGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCGCCTAGCGGCCCCGGCAAATTCCGCGGGTATGCAACAAACCGGGGgc                                                                                                          >  2:1647100/1‑150 (MQ=39)
                                             gATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGgctgc                                                                                                       <  2:3295989/150‑1 (MQ=255)
                                             gATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGgctgc                                                                                                       <  2:60518/150‑1 (MQ=255)
                                               ttAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGAAGGGCCCCCTAGCGGCCCCCGCCATTTCCGCGGGGATGGAACAAACCTGGGGCTGCTa                                                                                                     >  2:5209931/1‑150 (MQ=255)
                                                  gTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCCTAGCGGCCCCGGCCCATTACGCGGGGATGCACCAAACCGGGggcggcaaggc                                                                                                  >  2:6469959/1‑145 (MQ=37)
                                                  gTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCcgc                                                                                                                                >  2:1175358/1‑120 (MQ=255)
                                                  gTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCcgc                                                                                                                                <  1:1175358/120‑1 (MQ=255)
                                                   tCGGGCGGACATGAGCAACTCCTTTTGTCTTGTGTCGCCGACAGCATAGCAACACCGCTCTGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCa                                                                                                 <  1:20474/150‑1 (MQ=255)
                                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCCCCTAGCGGCCCCGGCAAATTCCGCGTGTATGCAACAAACCAGGGGCTGCTATGCa                                                                                                 >  2:7230698/1‑150 (MQ=255)
                                                       gCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGGCGCGAAATGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCAGGGGCGGCTATTCAGccc                                                                                             >  2:1209612/1‑150 (MQ=255)
                                                          gACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCttt                                                                                          <  1:1988252/150‑1 (MQ=255)
                                                           aCATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTc                                                                                         <  1:5367340/150‑1 (MQ=255)
                                                             aTGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCCGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTcgc                                                                                       <  1:389243/150‑1 (MQ=255)
                                                             aTGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTcgc                                                                                       <  2:4275700/150‑1 (MQ=255)
                                                             aTGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTcgc                                                                                       <  2:2443127/150‑1 (MQ=255)
                                                             aTGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTcgc                                                                                       <  1:4563362/150‑1 (MQ=255)
                                                              tGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGAAAGGGCCGCATAGCGGCCCCGGCCCATTGCGCGTGTATTCAACAAAACTGGGGCTGCTATGGAGCCCGTTagcg                                                                                      >  2:654593/1‑150 (MQ=37)
                                                               gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAAAGCGGCCCCCGCAAATTCCGCGTGTATGCAACAAAACAGGGGCTGCTATTCAGCCCATTAGCGa                                                                                     >  2:3514229/1‑150 (MQ=38)
                                                                 gCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGaca                                                                                   <  1:7196089/150‑1 (MQ=255)
                                                                 gCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGaca                                                                                   <  2:4564194/150‑1 (MQ=255)
                                                                   aaCTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTAGGGCTGCATAGCAGccccaggt                                                                                        <  1:4676737/143‑6 (MQ=38)
                                                                      tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCCAAAGGGCCGCCTAGCCGCCCCCGCAAATTAACCGGGTATTCAACAAACCAGGGGCTGCTATTCAGCCCATTCGCGACACaaag                                                                              >  2:5593914/1‑148 (MQ=255)
                                                                      tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCTTAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCt                                                                                            <  1:5527694/136‑1 (MQ=255)
                                                                      tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCt                                                                                            >  2:5527694/1‑136 (MQ=255)
                                                                      tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAAATTCCGCGGGTATGCAACAAACCAGGGGCTGCGATGCAGCCCATTCCCGACCCAAgg                                                                              >  2:3779200/1‑150 (MQ=38)
                                                                        cTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAAATTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGcc                                                                            >  2:2990301/1‑150 (MQ=255)
                                                                         ttttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGAATGGCCGCATAGCGGCCCCGGCCAGTTCCGCGGGTATGCAACAAAACGGGGGCTGCTATGGAGCCCATTAGCGACACAAGGccc                                                                           >  2:3904251/1‑150 (MQ=38)
                                                                         ttttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCGCCTATCGGCCCCGGCAAATTGCGCGTGTATGCCACAAACCTGGGGCGGCTATGCAGCCCATTGGCGTCACAAGGccc                                                                           >  2:3132157/1‑150 (MQ=38)
                                                                         ttttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGc                                                                                                  <  1:4777707/127‑1 (MQ=255)
                                                                         ttttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCCCATAGCGGCCCCCGCAAATTACGCCTGTATGCAAAAAAACACGGGCTGCTATGAAGCCCATTAGCGACACAACGccc                                                                           >  2:5914802/1‑150 (MQ=37)
                                                                             gTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAATGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACaa                                                                       <  2:7181917/150‑1 (MQ=255)
                                                                              tCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCGATGCAGCCCGTTCGCCACACAAGGCCCACaaa                                                                      >  2:3651833/1‑150 (MQ=255)
                                                                                   tgtCGCCGAAAGCATAGCAACACCGCTCTGGTGTGGGAGCGGCATTTTGTCGCTAATGGGCCGCATTGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAGGGAcc                                                                 <  1:3640118/150‑1 (MQ=38)
                                                                                     tCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGCAAGGGCCGCCTAGCGGCCCCGGCCAGTTGCGCGGGTTTGCCACAAACCGGGGGCTGCTATGCCGCCCGTTCGCGCCACAACGCCCCCCAAGGGCcgc                                                               >  2:3760578/1‑150 (MQ=255)
                                                                                        ccGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCa                                                            <  2:4804589/150‑1 (MQ=255)
                                                                                         cGACAGCATCGAAACCCGGGTGGGGTGTGGGCGGGGCTTTGTGTCGCGAACGGGCCCCATCGGGCCCCCGGAAATTTCCGGGTGTATGCAAAAAACCGGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAg                                                           <  1:4994366/150‑1 (MQ=255)
                                                                                             aggaTAGCAACACCGCTCAGGTGTGGGAGCGGCATTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCATCCCTTTCGCGACACAAGGCCCACAAAGGACCGAGCAGGCtt                                                       <  2:5179083/147‑1 (MQ=255)
                                                                                             aGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCa                                                                                                 <  1:63893/108‑1 (MQ=255)
                                                                                                aTAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGTTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCGTCAg                                                    >  2:6162257/1‑150 (MQ=255)
                                                                                                 tAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAATGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTTTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAgg                                                   <  1:1106992/150‑1 (MQ=255)
                                                                                                   gCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAACGGCCGCAAAGCGGCCCCGGCAAATTCCGCGTGTATGCAACAAAACGGGGGCGGCTAAGCAGCCCGTTCGCGCCACAACGCCCACAAAGGACCGCGCCGGCGTCAGGcc                                                 >  2:5840422/1‑150 (MQ=25)
                                                                                                    cAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGGAATGGCCGCATAGCGGCCCCCGCCAATTGCGCGTGTATTCAACAAACCAGGGGCTGCTAATCAGCCCGTTCGCGACACAAGGCCCACAAAAGACCCCGCAGGCGTCAGGccc                                                >  2:3732746/1‑150 (MQ=255)
                                                                                                     aaCACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAAAGCGGCCCCCGCAAATTGCGCGGGTATGCAACAAAACCGGGGCTGCGATGCAGCCCGTTCGCGGCACCACGCCCACCAAGGAACGCGCCCGCGTCAGGCCCg                                               >  2:5658222/1‑150 (MQ=21)
                                                                                                        acCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCGCCTAGCGGCCCCCGCCAGTTCCGCGGGTCTGCCACAAACCGGGGGCTGCGATGCCGCCCGGTCGCGCCACAAGGCCCACCCAGGCCCCCGCCGGCGTCAGGCCCggga                                            >  2:2338384/1‑148 (MQ=255)
                                                                                                        acCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGc                                                                                               >  1:2934698/1‑99 (MQ=255)
                                                                                                           gCTCTGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTg                                         <  2:5018898/150‑1 (MQ=255)
                                                                                                            cTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTgg                                        <  2:1502069/150‑1 (MQ=255)
                                                                                                             acAGGTGTGGCAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGCCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGc                                       <  2:2211661/149‑1 (MQ=255)
                                                                                                              ctggTGGGGGCGGGGCCTTGTGTCGCTAACGGGCCGCATCGGGGCCCCGGCAATTTCCGCGTGTATGCAACAACCCTGGGGCGGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGcc                                      <  1:2694402/148‑1 (MQ=255)
                                                                                                                        aGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAAAAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGgcg                           <  1:2381575/150‑1 (MQ=255)
                                                                                                                           ggCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCCGCAATTTCCGCGTGTATGCAACAAACCAGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGgcgcgc                        >  2:2735486/1‑150 (MQ=255)
                                                                                                                            gCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCa                       <  2:4744655/150‑1 (MQ=255)
                                                                                                                             ccTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGACAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAg                      >  2:943257/1‑150 (MQ=255)
                                                                                                                              cTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGaa                                            >  1:969965/1‑128 (MQ=255)
                                                                                                                              cTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGc                     >  2:2865136/1‑150 (MQ=255)
                                                                                                                                 gtgtCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                  <  2:132761/150‑1 (MQ=255)
                                                                                                                                 gtgtCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                  <  2:4111117/150‑1 (MQ=255)
                                                                                                                                 gtgtCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCATTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                  <  2:6287051/150‑1 (MQ=255)
                                                                                                                                   gtCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgc                >  1:6057290/1‑150 (MQ=255)
                                                                                                                                    tCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACAGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgcc               >  2:762061/1‑150 (MQ=255)
                                                                                                                                    gcgcgACTGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgcc               <  1:3373530/149‑1 (MQ=255)
                                                                                                                                      gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc             >  1:504546/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCa           >  2:255390/1‑150 (MQ=255)
                                                                                                                                          aaGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGa         >  1:2951426/1‑150 (MQ=255)
                                                                                                                                             ggCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatca      <  1:3177091/150‑1 (MQ=255)
                                                                                                                                              gCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcat     >  1:5871052/1‑150 (MQ=255)
                                                                                                                                              gCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcat     >  2:2278444/1‑150 (MQ=255)
                                                                                                                                               ccGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGGCTCAGCTGCGCGCCAgcctggg    <  1:4152999/150‑7 (MQ=255)
                                                                                                                                               ccGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcatc    >  1:5234617/1‑150 (MQ=255)
                                                                                                                                               ccGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcatc    <  2:4718515/150‑1 (MQ=255)
                                                                                                                                                cGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCCGATCATcc   >  2:6517715/1‑150 (MQ=255)
                                                                                                                                                cGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATcc   >  2:5338934/1‑150 (MQ=255)
                                                                                                                                                 gCATAGCGGCCCCGGCAATTTCCGCGTGTAGGCAACAAACCTGGGGCTGCTATGCAGCCCTTACGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCa  <  1:1281282/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCA  >  NC_002947/5560453‑5560747

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