Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,557,574 T→A 10.0% 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,557,5740TA10.0% 400.8 / 11.1 157intergenic (‑80/+94)PP_4887/PP_4888hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (45/89);  new base A (15/4);  total (64/93)
Fisher's exact test for biased strand distribution p-value = 2.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

GCCGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCA  >  NC_002947/5557425‑5557723
                                                                                                                                                     |                                                                                                                                                     
gCCGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAt                                                                                                                                                       >  2:1184119/1‑150 (MQ=255)
   gTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGc                                                                                                                                                    <  1:1549298/150‑1 (MQ=255)
      cGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGc                                                                                                                                                 <  1:1336994/150‑1 (MQ=255)
           tCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTTATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCATGTGGCGCGAAAGGGCCGCAAAGCcccccccg                                                                                                                                            >  2:3525402/1‑148 (MQ=255)
           tCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCTATGCGGcccttt                                                                                                                                            <  2:1478997/150‑4 (MQ=255)
           tCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCTATGCGGcccttt                                                                                                                                            <  2:362136/150‑4 (MQ=255)
            cACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCTATGCGGccctttc                                                                                                                                           <  1:775024/150‑5 (MQ=255)
             aCTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCa                                                                                                                                          <  2:432600/150‑1 (MQ=255)
              cagtgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTATTGTCATGTGTCGCCGACAGCAAAGCAACACCGCTCAGGTGAGGGAGCGCCCAAATTGCGCGAAAAAGCCCCAAATAGGCCCCGGaaa                                                                                                                                         >  2:1611798/3‑150 (MQ=255)
                 tgCAGGAAATTCCGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAAttt                                                                                                                                      <  2:1637517/150‑1 (MQ=255)
                 tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAAttt                                                                                                                                      <  2:2896562/150‑1 (MQ=255)
                 tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGGGGGGGGGAAGGGGCGCCAAGCGGGCCCCGcccaat                                                                                                                                      >  1:1663694/1‑145 (MQ=255)
                   cAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTTTTTGGCGCAAGGGGCGCCGAGCGGCCCCCGccaaatac                                                                                                                                    >  1:2659957/1‑143 (MQ=255)
                   cAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCCAATTcc                                                                                                                                    >  2:2859256/1‑150 (MQ=255)
                   cAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTcc                                                                                                                                    <  1:1872810/150‑1 (MQ=255)
                     ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCcgc                                                                                                                                  <  2:905396/150‑1 (MQ=255)
                      gAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGAAAGGGCCGCAGAGCGGCCCCCGCCAATTAcgcg                                                                                                                                 >  1:2236102/1‑150 (MQ=255)
                         aTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCgtgt                                                                                                                              <  1:3298312/150‑1 (MQ=255)
                          ttCAGGTCGGCCTGGGCAGCTATTCGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTa                                                                                                                             <  1:3879462/150‑1 (MQ=255)
                            cAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGAAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATg                                                                                                                           <  1:3770951/150‑1 (MQ=255)
                             aGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGc                                                                                                                          <  1:2335009/150‑1 (MQ=255)
                                           aGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCCGCAATTTCCGCGTGTATGCAACAAACCAGGGgc                                                                                                            >  2:3976034/1‑150 (MQ=255)
                                           aGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGGCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGgc                                                                                                            >  2:3976026/1‑150 (MQ=255)
                                              gATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGgctgc                                                                                                         <  2:3278887/150‑1 (MQ=255)
                                                 tAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTAt                                                                                                      <  1:1296869/150‑1 (MQ=255)
                                                        gCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGccc                                                                                               <  2:2893080/150‑1 (MQ=255)
                                                        gCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCCAAAGGGCCCCAAAGCGGCCCCCGCAAATTCCCCGGGGATGCAACAAAACAGGGGCTGCGATGCAGccc                                                                                               >  2:3814495/1‑150 (MQ=37)
                                                        gCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCCCGAAAGGGCCGCAAAGCGGCCCCCGCAAATTACGCGTGTATGCAACAAACCAGGGGCTGCTAAGCAcccc                                                                                               >  2:1035490/1‑150 (MQ=38)
                                                          ggACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCATAGCAGcccca                                                                                             <  2:4343555/150‑3 (MQ=255)
                                                          ggACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCATAGCAGcccca                                                                                             <  2:1523875/150‑3 (MQ=255)
                                                          ggACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAAATTCCGCGTGTATGCAACAAACCTGGGGCTGCGATGCAGCCCtt                                                                                             >  2:2195063/1‑150 (MQ=255)
                                                            aCATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTc                                                                                           <  1:3943351/150‑1 (MQ=255)
                                                            aCATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTAGGGCTGCATAGCAGccccagg                                                                                           <  2:2367743/150‑5 (MQ=38)
                                                                gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGa                                                                                       <  1:673981/150‑1 (MQ=255)
                                                                   cAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAAAGCGGCCCCGGCCAATTCCGCGGGGATGCAACAAACCAGGGGCTGCTATGCAGCCCGTTCGCGacac                                                                                    >  2:4415851/1‑150 (MQ=38)
                                                                     aCTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGAAGGGCCGCCAAGCGGCCCCGGCCCTTTGCGCGTGTGTGGAACAAAACAGGGGCTGCTATGCAGCCCGTTCGCGACACaa                                                                                  >  2:3516889/1‑150 (MQ=37)
                                                                     aCTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGGAATGGCCCCCAAGCGGCCCCCGCCAATTACGCCGGTATTCAACAAAACAGGGGCTGCAATTCACCCCCATcgccgcccaa                                                                                  >  1:1833831/1‑143 (MQ=255)
                                                                      cTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAg                                                                                 <  2:4397008/150‑1 (MQ=255)
                                                                       tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCCGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAgg                                                                                <  1:1355831/150‑1 (MQ=255)
                                                                       tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGc                                                                                                                          >  1:1083559/1‑108 (MQ=255)
                                                                       tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGc                                                                                                                          <  2:1083559/108‑1 (MQ=255)
                                                                       tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAgg                                                                                <  2:2339915/150‑1 (MQ=255)
                                                                        ccTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCGCATGTCGGCCCCCGCCAATTCCGCGGGGATGCCCCAAAACGGGGGCTGCGATGCAGCCCGTTAGCGACACAAGGc                                                                               >  2:1364693/1‑150 (MQ=255)
                                                                         cTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGAATGGGCGCATAACGGCCCCCGCCATTTACGCGTGTATGGAACAAAACTGGGGCTGCTATGCAGCCCTTTCGCGCCAAAAGGcc                                                                              >  2:93224/1‑150 (MQ=255)
                                                                         cTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGc                                                                                                                          >  1:4002368/1‑106 (MQ=255)
                                                                          ttttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCCATTTCCGCGTGTTTGCAACAAACCAGGGGCTGCTATGCAGCCCTTTCGCGACACCAGGccc                                                                             >  2:4194732/1‑150 (MQ=255)
                                                                            ttGTCTTGTGGCGCCGACAGCATAGAAACACCGCTCTGGTGTGGGAGGGGCGTTGGTTCGCGAAGGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTGCGCGACACAAGGCCcac                                                                           <  1:851604/150‑1 (MQ=255)
                                                                               tCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACaaa                                                                        <  2:3581363/150‑1 (MQ=255)
                                                                               tCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGGTTCGCGAAAAGGCCGCCAAGCGGCCCCCGCCAATTACCCGGGTGTGCAACAAACCGGGGGCGGGTATGCCGCCCGTTAGCGTCACAAGGCCCACaaa                                                                        >  2:3830434/1‑150 (MQ=255)
                                                                                cATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCCACAAACCAGGGGCTGCTATGCAGCCCATTCGCGACACCACGCCCCCAAAg                                                                       >  2:2518523/1‑150 (MQ=38)
                                                                                cATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAg                                                                       <  2:3079412/150‑1 (MQ=255)
                                                                                 aTGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGa                                                                                       <  1:1928420/133‑1 (MQ=255)
                                                                                 aTGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGa                                                                                       >  2:1928420/1‑133 (MQ=255)
                                                                                   gtgtCGCCGACTGAATAGAAACACCGCTCTGGTGTGGTTGCGGCCTTGTGTCGCGAACGGGCCGATTAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAc                                                                    <  1:2134443/150‑1 (MQ=37)
                                                                                   gtgtCGCCGACAGATTAGCAACACCGCTCAGGTGTGGGTGGGGCCTTGTGTGGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTTTGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAc                                                                    <  1:2134504/150‑1 (MQ=38)
                                                                                    tgtgGCCTACTGAATAGCAACCCGGCTCTGGTGTGGGTGGGGGTTTGTGTCGCTAATGGGCCGCTTAGCGCCCCGGGCAATTTCCGCGTTTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                                                   <  1:3516889/150‑1 (MQ=25)
                                                                                    tgtCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTc                                                                                           >  1:3729330/1‑126 (MQ=255)
                                                                                    tgtCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTc                                                                                           <  2:3729330/126‑1 (MQ=255)
                                                                                     gtCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGCAAGGGGCCCGTATCCGCCCCGGCCAATTGCGCGTGTATGCCACAAACCAGGGGCAGCGTAGGAGCCCGTTGGCGGCACAACGCCCACCAAGGAccc                                                                  >  1:2958836/1‑149 (MQ=255)
                                                                                         cctaCAGCATAGCAACACCGCTCTGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCa                                                              <  2:199511/147‑1 (MQ=255)
                                                                                              aGCAGAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCtt                                                         <  2:552801/150‑1 (MQ=255)
                                                                                                cATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGAAATGGCCGCCTAGCGGCCCCCGCCATTTGCGCGGGTATGCCACAAAACGGGGGCTGCTGTGCAGCCCGTTAGCGACCCAAGGCCCCCAAAGGAACCCCCAGGCGTcc                                                       >  2:965927/1‑149 (MQ=255)
                                                                                                cATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCa                                                       <  2:3053236/150‑1 (MQ=255)
                                                                                                cATAGCAACACCGCTCAGGTGGGGGAGCGGCCTTGTGTCGCGTAAGGGCCGCATAGCGGCCCCGGCAAGTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCa                                                       <  2:425764/150‑1 (MQ=255)
                                                                                                cATAGAAACCCCGCTCTGGTGTGGGTGCGGCCTTGTGTCGCGAACGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAAAAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCa                                                       <  1:2758512/150‑1 (MQ=37)
                                                                                                 aTAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                                      <  2:265719/150‑1 (MQ=255)
                                                                                                  tAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAATGGCCGCATAGCGGCCCCGGCAAATTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTAGCGACACAAGGCCCACAAAGGACCGCGCCGGCGTCAgg                                                     >  2:4013025/1‑150 (MQ=38)
                                                                                                   aGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCCGCAAATTACGCCAGTTTGCAACAAACCTGGGGCTGCGATGCAGCCCGTTCGCGACACAAGGCCCACCAAGGACCCCGCAGGCGTCAGGc                                                    >  2:1457624/1‑150 (MQ=37)
                                                                                                     cAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                                               >  1:894435/1‑89 (MQ=255)
                                                                                                      aaCACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCCCCTAGCGGCCCCCGCAAATTCCGCGGGTTTGCAACAAACCAGGGGCGGCGCTGCAGCCCGTTCGCGCCACCAGGCCCCCCAAGGACCCCGCCGGCGTCATGCCCg                                                 >  2:627395/1‑150 (MQ=255)
                                                                                                       acacCGCTCTGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                                                <  1:2563235/150‑1 (MQ=255)
                                                                                                        cacCGCTCAGGTGTGGGAGCGGCCTTTTGTGGCGCAAAGGCCGCAAAGCGGCCCCCGCAAATTACGCGGGTTTTGCACAAAACAGGGGGTGCTTATCAGCCCCTTCGCGACAAAAAGCCCACAAAGGGCCCCCCAGGCGTGAGGccccgc                                               >  2:1247423/1‑147 (MQ=255)
                                                                                                        cacCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCCCCAAGCCGCCCCCGCCAATTACCCCTGTATGCAACAAAACAGGGGCTGGGATTCAGCCCATTCGCGACACAAGGCCCCCAAACGACCGCCCCGGCGTACGGCCCGGa                                               >  2:258739/1‑150 (MQ=255)
                                                                                                         acCGCTCAGGTGTGGGAGCGGCCTTTTGTGGCGAAAAGGCCGCCAAGCGGCCCCCGCAAATTCCCCGGTTATGGAAAAAAACAGGGGCTGGGAATCAGCCCATTAGCGACACAACGCCCACAAAGGACCCCGCAGGCGTCACGCCCCGaa                                              >  2:3032067/1‑150 (MQ=255)
                                                                                                           cGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTTTTGAACAAAACGGGGGCTGGGATGCAGCCCGTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCCGAACt                                            >  2:2536022/1‑150 (MQ=37)
                                                                                                             ctctGGTGTGGGAGCGGCCTTGGGTCGCGAACGGGCCGCATAGCGCCCCCGGCAAGTCCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTgg                                          <  1:1272164/150‑1 (MQ=38)
                                                                                                             cTCAGGTGTGGGCGCGGCCTTGGGCCGCGAACGGGCCGCATCGCGGCCCCGGCAATTCCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTgg                                          <  1:918034/150‑1 (MQ=38)
                                                                                                             cTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTgg                                          <  1:2932268/150‑1 (MQ=255)
                                                                                                                  gtgtGGGAGCGGCCTTGTGTCGCTATAGGGCCGCATAGCGGCGCCGGCAATTTCCGCGTGTATGCATCAAACCTGGGGCTTCTTTGCAGGCCTTTCGCGACACAAGGCCCTCAAAGGACCGCGCAGGCCTCAGGCCGGGAACTGGCCCAg                                     <  2:3448986/150‑1 (MQ=38)
                                                                                                                     tgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCt                                  <  2:78757/150‑1 (MQ=255)
                                                                                                                      gggAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg                                 <  2:4130425/150‑1 (MQ=255)
                                                                                                                             gCCTTGGGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATGTCCGCGTGTATGCAACAAACCGGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCa                          <  1:3359366/150‑1 (MQ=255)
                                                                                                                                ttGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCt                       >  1:4500140/1‑150 (MQ=255)
                                                                                                                                   tgtCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAg                    >  2:3384194/1‑150 (MQ=255)
                                                                                                                                     tCGCGACTGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgcc                  <  1:2031366/150‑1 (MQ=255)
                                                                                                                                     tCGCGAATGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGcc                                                                                                <  1:1003944/72‑1 (MQ=255)
                                                                                                                                     tCGCGAACGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgcc                  <  1:276644/150‑1 (MQ=255)
                                                                                                                                     tCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgcc                  >  1:2735798/1‑150 (MQ=255)
                                                                                                                                       gcgAAGGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc                >  2:3657017/1‑150 (MQ=255)
                                                                                                                                       gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                                                <  2:818523/118‑1 (MQ=255)
                                                                                                                                       gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc                <  2:2316290/150‑1 (MQ=255)
                                                                                                                                       gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc                <  2:1764831/150‑1 (MQ=255)
                                                                                                                                         gAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCa              <  2:1576887/150‑1 (MQ=255)
                                                                                                                                              ggCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGac                                                                                      <  2:2164665/73‑1 (MQ=255)
                                                                                                                                              ggCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGac                                                                                      >  1:2164665/1‑73 (MQ=255)
                                                                                                                                                ccGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcatc       >  1:4096242/1‑150 (MQ=255)
                                                                                                                                                   cATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAg    <  1:423575/150‑1 (MQ=255)
                                                                                                                                                    aTAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGc   <  2:1150109/150‑1 (MQ=255)
                                                                                                                                                     tAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCa  <  2:3456534/150‑1 (MQ=255)
                                                                                                                                                     tAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCGGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCa  <  2:585478/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCA  >  NC_002947/5557425‑5557723

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: