Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,985,587 T→A 5.5% intergenic (‑32/+39) PP_4387 ← / ← flgE hypothetical protein/flagellar hook protein FlgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,985,5870TA 5.5% 979.1 / 30.7 293intergenic (‑32/+39)PP_4387/flgEhypothetical protein/flagellar hook protein FlgE
Reads supporting (aligned to +/- strand):  ref base T (129/147);  new base A (11/5);  total (140/153)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

Only 100 of 299 total aligned reads displayed.

ACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAG  >  NC_002947/4985439‑4985736
                                                                                                                                                    |                                                                                                                                                     
  acaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGc                                                                                                                                                    <  2:852634/150‑1 (MQ=255)
  acaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGc                                                                                                                                                    <  2:4759131/150‑1 (MQ=255)
  acaGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGc                                                                                                                                                    <  1:4361507/150‑1 (MQ=255)
    aGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGGAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCgg                                                                                                                                                  >  1:1422994/1‑150 (MQ=255)
    aGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCgg                                                                                                                                                  >  2:3422920/1‑150 (MQ=255)
     ggTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGc                                                                                                                                                 <  2:1984380/150‑1 (MQ=255)
      gTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGACCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGcc                                                                                                                                                <  1:3761930/150‑1 (MQ=255)
       tAGCGGCCGGTGACCAGGCTTTGCCATTGATTGTTGGGGGCGGGGGCTGCGGGGGCAGGGAGATCCAGCTGCAGCTAGGGGATGTAGGGCAGTTACGAGAGGCGGTTCTTTCGGGTGAAGGTTGTGATTATTGGGGCCGGAATGCGGccg                                                                                                                                               <  2:1039326/150‑2 (MQ=255)
           ggCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAAt                                                                                                                                           <  1:1122662/150‑1 (MQ=255)
           ggCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAAt                                                                                                                                           >  1:1896324/1‑150 (MQ=255)
              cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCaagaat                                                                                                                                        <  1:5695063/150‑5 (MQ=255)
              cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTAc                                                                                                                                        <  2:670510/150‑1 (MQ=255)
              cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTAc                                                                                                                                        >  1:6698715/1‑150 (MQ=255)
              cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTAc                                                                                                                                        <  2:2467354/150‑1 (MQ=255)
              cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTAc                                                                                                                                        <  1:4070063/150‑1 (MQ=255)
              cGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTATTGGGGCCGCAAAGCGGCCCCaagaat                                                                                                                                        <  2:4049160/150‑5 (MQ=255)
                 tGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACaaa                                                                                                                                     <  2:4134804/150‑1 (MQ=255)
                 tGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACaaa                                                                                                                                     <  1:4806870/150‑1 (MQ=255)
                 tGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACaaa                                                                                                                                     <  2:5308912/150‑1 (MQ=255)
                    ccAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAacac                                                                                                                                  >  1:6257597/1‑150 (MQ=255)
                     cAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACAcc                                                                                                                                 >  1:4633983/1‑150 (MQ=255)
                        gCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCaaa                                                                                                                              <  2:684294/150‑1 (MQ=255)
                        gCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCaaa                                                                                                                              <  2:159208/150‑1 (MQ=255)
                        gCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCaaa                                                                                                                              <  1:723291/150‑1 (MQ=255)
                        gCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCaaa                                                                                                                              <  1:1621990/150‑1 (MQ=255)
                          tttGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGc                                                                                                                            >  2:5885397/1‑150 (MQ=255)
                                    aTAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatc                                                                                                                  <  1:3331488/150‑1 (MQ=255)
                                       gAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAatcatc                                                                                                               <  2:2641167/150‑1 (MQ=255)
                                        aaCGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAAtcatca                                                                                                              <  1:335051/150‑1 (MQ=255)
                                           ggggcggggGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGt                                                                                                           >  2:2800552/1‑150 (MQ=255)
                                             ggcggggGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAAACAACAGGTCa                                                                                                         >  1:844695/1‑150 (MQ=255)
                                                  gggcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGa                                                                                                    >  1:5366128/1‑150 (MQ=255)
                                                  gggcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGACATCTGa                                                                                                    >  2:4926014/1‑150 (MQ=255)
                                                  gggcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCAACAGGTCATCTGa                                                                                                    >  1:2910217/1‑150 (MQ=255)
                                                   ggcagcgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGaa                                                                                                   <  1:1408747/150‑1 (MQ=255)
                                                        cgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgatg                                                                                              <  2:6154663/150‑1 (MQ=255)
                                                        cgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgatg                                                                                              <  2:5366128/150‑1 (MQ=255)
                                                        cgggggcagGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgatg                                                                                              <  1:5672471/150‑1 (MQ=255)
                                                          ggggcagGCAGAACCAGCACCAGCAAGGTGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGTGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGt                                                                                            <  1:2573381/150‑1 (MQ=255)
                                                            ggcaggcaGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCt                                                                                          >  2:1502147/1‑150 (MQ=255)
                                                             gcaggcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGAAAGCCGCTTCATTCGCGTGAACCTTCACAATCTTCGGGACGCTTTGCCGCCCCAATTACAAACACCAATGCCAAGACACACAGCAGGTAATCTGAAAGATGGTCTg                                                                                         >  2:5780274/1‑150 (MQ=255)
                                                              caggcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGc                                                                                        <  2:164342/150‑1 (MQ=255)
                                                              caggcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGc                                                                                        <  1:3546333/150‑1 (MQ=255)
                                                              caggcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGc                                                                                        <  1:4140916/150‑1 (MQ=255)
                                                               aggcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCa                                                                                       >  2:2949186/1‑150 (MQ=255)
                                                               aggcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCa                                                                                       >  1:3969823/1‑150 (MQ=255)
                                                               aggcagAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAAATACAAACACCAAAGCCAAAACAATCATCAGGTCATCTGAATGATGGTCTGCa                                                                                       >  1:3092081/1‑150 (MQ=255)
                                                                       ccagcaccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTg                                                                               >  2:5218167/1‑150 (MQ=255)
                                                                        cagcaccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGc                                                                              <  2:3379681/150‑1 (MQ=255)
                                                                           caccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTAAAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCttt                                                                           >  1:2447328/1‑150 (MQ=255)
                                                                             ccagcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCa                                                                         <  1:4606145/150‑1 (MQ=255)
                                                                               agcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGt                                                                       <  1:6149763/150‑1 (MQ=255)
                                                                               agcaAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGt                                                                       <  1:5024330/150‑1 (MQ=255)
                                                                                gcaAGGGGAACAAGGCCAGTAACGACAGCAGCTTCATTCTCGTGAACATTCACATTCTTGGGGCCGCTTTGCGGCCCCATTTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGTTGGTCTGCATGATGGTGCTTTCAGTg                                                                      <  2:3667417/150‑1 (MQ=255)
                                                                                 caAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAAAACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTgg                                                                     >  1:1792980/1‑150 (MQ=255)
                                                                                    ggggAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGaaa                                                                  <  1:2066042/150‑1 (MQ=255)
                                                                                    ggggAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGaaa                                                                  >  2:1259598/1‑150 (MQ=255)
                                                                                       gAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATgg                                                               >  2:1759864/1‑150 (MQ=255)
                                                                                       gAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATgg                                                               >  1:6239539/1‑150 (MQ=255)
                                                                                           aaGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCGTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCtt                                                           <  2:190401/150‑1 (MQ=255)
                                                                                              gCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGc                                                        >  2:3947171/1‑150 (MQ=255)
                                                                                              gCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGc                                                        >  2:5942910/1‑150 (MQ=255)
                                                                                              gCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGc                                                        >  1:486154/1‑150 (MQ=255)
                                                                                              gCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTAAAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGc                                                        >  1:2964064/1‑150 (MQ=255)
                                                                                                cAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCAt                                                      >  2:5055229/1‑150 (MQ=255)
                                                                                                       gACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTg                                               <  1:4096378/150‑1 (MQ=255)
                                                                                                        aCAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTgg                                              <  2:1115077/150‑1 (MQ=255)
                                                                                                          aGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACCCCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTa                                            <  2:5484851/150‑1 (MQ=255)
                                                                                                          aGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTa                                            <  2:1484705/150‑1 (MQ=255)
                                                                                                          aGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTa                                            <  2:4610394/150‑1 (MQ=255)
                                                                                                            ccGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGCCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGt                                          <  1:4926014/150‑1 (MQ=255)
                                                                                                            ccGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAAATACAAACACCAAAGCCAATAAAAAAATCAGGTCATCTGAATTATGGTCTGCATGATTGTGCTTTAAGTGGAAATGGTCTTGGCATTGGCCTGGTAGt                                          >  1:4428011/1‑150 (MQ=255)
                                                                                                            ccGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCGTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGTTGGTGCTTTCAGGGGAAAAGGTCTTGGCATTGGCCAGGTAGt                                          >  1:2020730/1‑150 (MQ=255)
                                                                                                             cGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGtt                                         <  1:5494038/150‑1 (MQ=255)
                                                                                                              gCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTg                                        <  2:5163808/150‑1 (MQ=255)
                                                                                                                  cATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACCCCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGctct                                    <  1:1036960/150‑1 (MQ=255)
                                                                                                                  cATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCGGGTAGTTGctct                                    >  1:1512792/1‑150 (MQ=255)
                                                                                                                   aTTCGCGTGAACCTTCACATTCTTGGGGCCCCTTTGCGGCCCCCATTACAAACACCAAATCCAATAAAAACATCAGGTAATAGGTAAGATGGTCTGCATGATGGTGCTTTAAGTGGAAAAGGTCGTGGCAATGGCCTGGTAGGTGCTCTg                                   >  1:1998902/1‑150 (MQ=255)
                                                                                                                     tCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTgcg                                 <  1:6641843/150‑1 (MQ=255)
                                                                                                                          tGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATAAACAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGGCTTGGCCTTGGCCTGGTAGTTGCTCTGCGCCTTg                            >  1:3496791/1‑150 (MQ=255)
                                                                                                                            aaCCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGAt                          <  2:4964077/150‑1 (MQ=255)
                                                                                                                              ccTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCa                        <  2:3956874/150‑1 (MQ=255)
                                                                                                                              ccTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCa                        <  1:5771614/150‑1 (MQ=255)
                                                                                                                               cTTCACATTCTTGGGGCCGCTTTGCGGCCCCAAATACAAAAACCAAAGCCAAAACAATCAACAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCAGTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg                       >  1:6477511/1‑150 (MQ=255)
                                                                                                                               cTTCACAGTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAg                       <  1:6082146/150‑1 (MQ=255)
                                                                                                                                ttCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACACTCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAgg                      <  1:5901433/150‑1 (MQ=255)
                                                                                                                                 tCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGt                     <  2:6616422/150‑1 (MQ=255)
                                                                                                                                 tCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGt                     <  2:6270551/150‑1 (MQ=255)
                                                                                                                                     aTTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAc                 >  2:6417458/1‑150 (MQ=255)
                                                                                                                                     aTTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAc                 >  2:5624526/1‑150 (MQ=255)
                                                                                                                                     aTTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAc                 >  1:234131/1‑150 (MQ=255)
                                                                                                                                     aTTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAc                 <  1:801490/150‑1 (MQ=255)
                                                                                                                                      ttCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTAcc                <  1:6294612/150‑1 (MQ=255)
                                                                                                                                       tCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCa               <  2:3258924/150‑1 (MQ=255)
                                                                                                                                       aattGGGGCCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAAACATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCCTTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCa               >  1:2639036/3‑150 (MQ=255)
                                                                                                                                        ctgggggCCGCTTTGCGGCCCCAATTACAAACCCCAACGCCAATCCAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAg              <  2:456519/147‑1 (MQ=255)
                                                                                                                                        attGGGGCCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAatgatg                                                                                              >  2:218370/2‑70 (MQ=38)
                                                                                                                                         ttGGGGCCGCAAAGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGc             >  1:2283344/1‑150 (MQ=255)
                                                                                                                                              gCCGCTTTGCGGCCCCAATTACAAACCCCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGcc        <  1:5268560/150‑1 (MQ=255)
                                                                                                                                              gCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGcc        >  2:4278520/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGGGCCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAG  >  NC_002947/4985439‑4985736

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: