Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,501,665 T→A 5.1% intergenic (+22/+30) PP_3987 → / ← PP_5651 site‑specific recombinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,501,6650TA 5.1% 1005.2 / 27.9 292intergenic (+22/+30)PP_3987/PP_5651site‑specific recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (119/158);  new base A (9/6);  total (128/164)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 302 total aligned reads displayed.

TACGCCCCAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAGTTG  >  NC_002947/4501516‑4501814
                                                                                                                                                     |                                                                                                                                                     
     cccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGAAATTCGAAGCGGGGTACGCTCCCCg                                                                                                                                                  >  1:4188490/1‑150 (MQ=255)
     cccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCAAACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCg                                                                                                                                                  >  2:237567/1‑150 (MQ=255)
      ccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGc                                                                                                                                                 >  1:57644/1‑150 (MQ=255)
        aaGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCtt                                                                                                                                               >  2:5189350/1‑150 (MQ=255)
        aaGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCtt                                                                                                                                               >  2:4102296/1‑150 (MQ=255)
         aGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCTAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCttt                                                                                                                                              <  2:5457444/150‑1 (MQ=255)
         aGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCttt                                                                                                                                              <  1:1893407/150‑1 (MQ=255)
         aGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGAGCGTACCCCGCttc                                                                                                                                              <  1:5400784/150‑2 (MQ=255)
           gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTcc                                                                                                                                                     >  1:2347071/1‑141 (MQ=255)
           gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                            <  1:4708358/150‑1 (MQ=255)
           gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                            <  1:6390196/150‑1 (MQ=255)
           gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                            <  2:3372530/150‑1 (MQ=255)
              gcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGTTAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGca                                                                                                                                         <  2:1760310/150‑1 (MQ=255)
              gcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGca                                                                                                                                         >  1:5237321/1‑150 (MQ=255)
               cAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGcac                                                                                                                                        >  1:6132726/1‑150 (MQ=255)
                   cGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAAtt                                                                                                                                    <  1:6756239/150‑1 (MQ=255)
                     aaaTCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTca                                                                                                                                  <  1:5066816/150‑1 (MQ=255)
                      aaTCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTcac                                                                                                                                 <  1:5192115/150‑1 (MQ=255)
                        tCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACaa                                                                                                                               >  1:5459874/1‑150 (MQ=255)
                           gCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTg                                                                                                                            >  2:3007040/1‑150 (MQ=255)
                             gATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGTTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGAt                                                                                                                          <  2:5328983/150‑1 (MQ=255)
                              aTTTGGGTCACGCGAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGca                                                                                                                                         <  2:363228/134‑1 (MQ=255)
                                  gggTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTAtt                                                                                                                     <  2:3542368/150‑1 (MQ=255)
                                  gggTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGATATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTAtt                                                                                                                     >  1:5560344/1‑150 (MQ=255)
                                   ggTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTg                                                                                                                    >  2:4892286/1‑150 (MQ=255)
                                    gTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGa                                                                                                                   >  1:1675267/1‑150 (MQ=255)
                                    gTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGCACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGa                                                                                                                   <  2:1937231/150‑1 (MQ=255)
                                     tCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAg                                                                                                                  >  2:3136532/1‑150 (MQ=255)
                                       aCGCAAGCCGCCTCTTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCt                                                                                                                >  2:2520632/1‑150 (MQ=255)
                                       aCGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCt                                                                                                                >  2:1049422/1‑150 (MQ=255)
                                        cGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCta                                                                                                               >  1:2295782/1‑150 (MQ=255)
                                         gCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCtat                                                                                                              >  2:4647853/1‑150 (MQ=255)
                                           aaGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATAc                                                                                                            <  1:6635780/150‑1 (MQ=255)
                                           aaGCCGCATATTGGCTCACACTGACAAGCTCATAACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATAc                                                                                                            <  1:6862949/150‑1 (MQ=255)
                                             gccgccTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGc                                                                                                          <  1:2509675/150‑1 (MQ=255)
                                                   tATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTa                                                                                                    <  1:4583019/150‑1 (MQ=255)
                                                   tATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTa                                                                                                    <  2:1259697/150‑1 (MQ=255)
                                                     ttGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATc                                                                                                  >  1:4086019/1‑150 (MQ=255)
                                                       gggTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCaa                                                                                                >  1:1447152/1‑150 (MQ=255)
                                                        ggTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAAt                                                                                               >  1:5271315/1‑150 (MQ=255)
                                                            acacTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTa                                                                                           >  1:4146970/1‑150 (MQ=255)
                                                             cacTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAg                                                                                          >  1:4808890/1‑150 (MQ=255)
                                                               cTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCt                                                                                        >  2:3892141/1‑150 (MQ=255)
                                                                tGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCtt                                                                                       >  1:2911475/1‑150 (MQ=255)
                                                                tGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCtt                                                                                       >  1:1177859/1‑150 (MQ=255)
                                                                  aCAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTtgc                                                                                     >  1:4485028/1‑150 (MQ=255)
                                                                        gcATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGctct                                                                               >  2:783981/1‑150 (MQ=255)
                                                                          aTCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGc                                                                             <  2:5798409/150‑1 (MQ=255)
                                                                          aTCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGc                                                                             >  2:1394861/1‑150 (MQ=255)
                                                                           tCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGcc                                                                            >  2:6448141/1‑150 (MQ=255)
                                                                           tCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGcc                                                                            <  2:1601896/150‑1 (MQ=255)
                                                                           tCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGcc                                                                            <  1:85207/150‑1 (MQ=255)
                                                                           tCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGcc                                                                            >  2:4259697/1‑150 (MQ=255)
                                                                             acacaAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCAc                                                                          >  2:4783160/1‑150 (MQ=255)
                                                                                  aaCGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAat                                                                     <  2:6327256/150‑1 (MQ=255)
                                                                                      gTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAatcatc                                                                 <  2:1127925/150‑1 (MQ=255)
                                                                                      gTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAatcatc                                                                 <  2:688016/150‑1 (MQ=255)
                                                                                       tatataGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAAtcatca                                                                >  1:3701811/1‑150 (MQ=255)
                                                                                         tataGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATa                                                              <  1:3283291/150‑1 (MQ=255)
                                                                                           taGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAgg                                                            <  1:2059460/150‑1 (MQ=255)
                                                                                            aGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGt                                                           >  2:3274737/1‑150 (MQ=255)
                                                                                            aGACGAGTAGGTGAAATCGTACAGCCAACACGAGAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGt                                                           >  2:590908/1‑150 (MQ=255)
                                                                                               cGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTc                                                        <  1:3383664/150‑1 (MQ=255)
                                                                                                 aGTAGGTGAAATCTTACAGCGAACACGTTAGGTAATTCGTAGCGGGGTACGCTGCCCGCTTTGAGCTCTATTCACTACTGATCTATTGAGCTATACGCTTTGTATCAATCTTAGCTTGCTGATCTGCCACCCTATCATCATAGGTTTAga                                                      <  2:1344123/150‑1 (MQ=255)
                                                                                                 aGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCga                                                      <  2:4113322/150‑1 (MQ=255)
                                                                                                  gTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCgag                                                     >  2:2842572/1‑150 (MQ=255)
                                                                                                     ggTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCgagaga                                                  >  1:4911552/1‑150 (MQ=255)
                                                                                                     ggTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCgagaga                                                  <  1:4674063/150‑1 (MQ=255)
                                                                                                      gTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCgatggaa                                                 >  2:3756732/1‑145 (MQ=255)
                                                                                                      gTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGaa                                                 <  1:3864297/150‑1 (MQ=255)
                                                                                                      gTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGaa                                                 >  2:4032374/1‑150 (MQ=255)
                                                                                                       tGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAAt                                                >  1:2770310/1‑150 (MQ=255)
                                                                                                         aaaTCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATAc                                              <  1:1814092/150‑1 (MQ=255)
                                                                                                           aTCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCt                                            <  2:2132224/150‑1 (MQ=255)
                                                                                                            tCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTc                                           >  1:3822619/1‑150 (MQ=255)
                                                                                                             cGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCg                                          <  1:6206831/150‑1 (MQ=255)
                                                                                                               tACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCa                                        >  1:3681528/1‑150 (MQ=255)
                                                                                                                  aGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGc                                     >  1:4248890/1‑150 (MQ=255)
                                                                                                                   gCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGcc                                    >  2:5625374/1‑150 (MQ=255)
                                                                                                                      aaCACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGa                                 <  2:6042675/150‑1 (MQ=255)
                                                                                                                      aaCACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGa                                 <  1:3091783/150‑1 (MQ=255)
                                                                                                                       acacGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAg                                <  1:257728/150‑1 (MQ=255)
                                                                                                                       acacGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAg                                <  1:4011616/150‑1 (MQ=255)
                                                                                                                        cacGAGAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGt                               <  1:4370361/150‑1 (MQ=255)
                                                                                                                                gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                       >  1:2650463/1‑150 (MQ=255)
                                                                                                                                gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                       <  2:714174/150‑1 (MQ=255)
                                                                                                                                gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAACCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                       >  1:2400414/1‑150 (MQ=255)
                                                                                                                                  tgctcaaaGCGGGGAGCGTACCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTgg                     >  1:3986513/7‑150 (MQ=255)
                                                                                                                                   aTTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGa                    <  2:4345613/150‑1 (MQ=255)
                                                                                                                                     tcaaaGCGGGGAGCGTACCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGAcc                  >  2:5126828/4‑150 (MQ=255)
                                                                                                                                     tcaaaGCGGGGAGCGTACCCCGCTTCGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGAcc                  >  2:6281170/4‑150 (MQ=255)
                                                                                                                                        aaGCGGGGAGCGTACCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTg               >  2:3973691/1‑150 (MQ=255)
                                                                                                                                         aGCGGGGTACGCTCCCCGCTTTGAGCACACTTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGc              <  2:228091/150‑1 (MQ=255)
                                                                                                                                         aGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCGTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGc              <  2:4188490/150‑1 (MQ=255)
                                                                                                                                           cGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCa            <  1:3007040/150‑1 (MQ=255)
                                                                                                                                             gggTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCaga          >  2:5536703/1‑150 (MQ=255)
                                                                                                                                              ggTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCagag         >  1:3101627/1‑150 (MQ=255)
                                                                                                                                                   gCTCCCAGCTTTGATCAATATTCACAACTGATCTATTGAGCTATTAGCTTTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAGt    <  1:5906798/150‑1 (MQ=255)
                                                                                                                                                     tCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAGTTg  >  2:1031051/1‑150 (MQ=255)
                                                                                                                                                     tCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAGTTg  >  2:6565606/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACGCCCCAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCTAGTTG  >  NC_002947/4501516‑4501814

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