Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,498,634 A→G 5.6% intergenic (+18/+34) 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,498,6340AG 5.6% 644.2 / 16.8 197intergenic (+18/+34)PP_3987/PP_5651site‑specific recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (92/93);  new base G (5/6);  total (97/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 199 total aligned reads displayed.

GACATACGCCCCAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCT  >  NC_002947/4498485‑4498782
                                                                                                                                                     |                                                                                                                                                    
gACATACGCCCCAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTa                                                                                                                                                      >  1:4436300/1‑150 (MQ=255)
       gCCCCAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTccc                                                                                                                                               >  1:63255/1‑150 (MQ=255)
        ccccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTcccc                                                                                                                                              >  1:652323/1‑150 (MQ=255)
         cccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCg                                                                                                                                             >  1:3372497/1‑150 (MQ=255)
         cccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCg                                                                                                                                             >  1:3626053/1‑150 (MQ=255)
         cccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCg                                                                                                                                             <  2:1396796/150‑1 (MQ=255)
          ccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGc                                                                                                                                            >  2:3388201/1‑150 (MQ=255)
          ccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGc                                                                                                                                            >  1:3879165/1‑150 (MQ=255)
           cAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCt                                                                                                                                           >  1:430787/1‑150 (MQ=255)
            aaGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCtt                                                                                                                                          >  1:2676481/1‑150 (MQ=255)
             aGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGAGCGTACCCCGCttc                                                                                                                                         <  1:3484358/150‑2 (MQ=255)
             aGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCAAAGCGGGGAGCGTACCCCGCttc                                                                                                                                         <  1:1116863/150‑2 (MQ=255)
               gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACCCGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                       <  1:2386513/150‑1 (MQ=255)
               gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                       <  2:2153471/150‑1 (MQ=255)
               gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                       <  1:4391116/150‑1 (MQ=255)
               gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                       <  1:1739282/150‑1 (MQ=255)
               gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                       <  1:1840171/150‑1 (MQ=255)
               gcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGa                                                                                                                                       <  1:313486/150‑1 (MQ=255)
                cggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAg                                                                                                                                      <  2:2843481/150‑1 (MQ=255)
                   cAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGcac                                                                                                                                   >  2:3300926/1‑150 (MQ=255)
                   cAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGcac                                                                                                                                   >  1:2741674/1‑150 (MQ=255)
                    aaGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCAAAGCGGGGAGCGTACCCCGCTTCGaattac                                                                                                                                  <  2:2291406/150‑6 (MQ=255)
                      gCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAAt                                                                                                                                <  1:3023430/150‑1 (MQ=255)
                       cGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAAtt                                                                                                                               <  1:1099455/150‑1 (MQ=255)
                        gAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTc                                                                                                                              >  1:3508560/1‑150 (MQ=255)
                              ggCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACt                                                                                                                        >  1:821526/1‑150 (MQ=255)
                              ggCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACt                                                                                                                        >  1:3918534/1‑150 (MQ=255)
                                 gATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGAt                                                                                                                     <  1:3407836/150‑1 (MQ=255)
                                      gggTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTAtt                                                                                                                <  2:1776651/150‑1 (MQ=255)
                                      gggTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTAtt                                                                                                                <  2:2178872/150‑1 (MQ=255)
                                        gTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGa                                                                                                              <  1:1237946/150‑1 (MQ=255)
                                         tCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAg                                                                                                             >  2:4209892/1‑150 (MQ=255)
                                         tCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAg                                                                                                             >  2:690706/1‑150 (MQ=255)
                                         tCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAg                                                                                                             >  1:4159239/1‑150 (MQ=255)
                                            cGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCta                                                                                                          <  1:1193544/150‑1 (MQ=255)
                                               aaGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATAc                                                                                                       <  1:767538/150‑1 (MQ=255)
                                                aGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACg                                                                                                      >  1:4060652/1‑150 (MQ=255)
                                                    gccTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCAtt                                                                                                  >  1:2335442/1‑150 (MQ=255)
                                                            ggTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAAt                                                                                          <  2:3727243/150‑1 (MQ=255)
                                                            ggTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAAt                                                                                          >  2:3544190/1‑150 (MQ=255)
                                                                acacTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTa                                                                                      <  1:1791629/150‑1 (MQ=255)
                                                                  acTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGc                                                                                    >  2:2430822/1‑150 (MQ=255)
                                                                  acTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGc                                                                                    >  2:4463141/1‑150 (MQ=255)
                                                                    tGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCtt                                                                                  >  2:2809173/1‑150 (MQ=255)
                                                                       cAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTtgct                                                                               <  2:1996793/150‑1 (MQ=255)
                                                                            gcATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGctct                                                                          <  2:976563/150‑1 (MQ=255)
                                                                             cATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTg                                                                         >  2:2635332/1‑150 (MQ=255)
                                                                               tCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACGCGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGcc                                                                       >  2:2094611/1‑150 (MQ=255)
                                                                               tCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCAGcc                                                                       >  1:250128/1‑150 (MQ=255)
                                                                                 acacaAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCAc                                                                     <  1:1515025/150‑1 (MQ=255)
                                                                                   acaAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCAccc                                                                   >  2:3425404/1‑150 (MQ=255)
                                                                                   acaAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCAccc                                                                   >  1:1378175/1‑150 (MQ=255)
                                                                                     aaaCGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCaa                                                                 >  2:2764567/1‑150 (MQ=255)
                                                                                       aCGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCAGAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAatc                                                               <  1:1322276/150‑1 (MQ=255)
                                                                                        cGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAatca                                                              <  2:1413103/150‑1 (MQ=255)
                                                                                            atataGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATcatcat                                                          <  2:1154425/150‑1 (MQ=255)
                                                                                             tataGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATa                                                         <  2:1579179/150‑1 (MQ=255)
                                                                                              ataGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAg                                                        <  1:2831311/150‑1 (MQ=255)
                                                                                              ataGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTGTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAg                                                        >  1:228140/1‑150 (MQ=255)
                                                                                                  aCGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGttt                                                    >  1:980950/1‑150 (MQ=255)
                                                                                                  aCGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGttt                                                    <  2:2397660/150‑1 (MQ=255)
                                                                                                    gAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCg                                                  >  2:4090314/1‑150 (MQ=255)
                                                                                                     aGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCga                                                 <  1:3525027/150‑1 (MQ=255)
                                                                                                     aGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCga                                                 <  2:1261204/150‑1 (MQ=255)
                                                                                                       tAGGTGAAATAGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCTCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCgaga                                               <  2:151767/150‑1 (MQ=255)
                                                                                                        aGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCgagag                                              <  1:1346190/150‑1 (MQ=255)
                                                                                                          gTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGaa                                            >  1:4410342/1‑150 (MQ=255)
                                                                                                            gAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATa                                          >  1:2517596/1‑150 (MQ=255)
                                                                                                                tCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTc                                      <  2:652323/150‑1 (MQ=255)
                                                                                                                tCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTc                                      <  1:1763530/150‑1 (MQ=255)
                                                                                                                 cGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCg                                     <  1:658254/150‑1 (MQ=255)
                                                                                                                  gTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGc                                    >  2:1930140/1‑150 (MQ=255)
                                                                                                                      aGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTGGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGc                                <  1:1972745/150‑1 (MQ=255)
                                                                                                                      aGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGc                                <  1:1675659/150‑1 (MQ=255)
                                                                                                                        ccAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGcca                              <  2:3190467/150‑2 (MQ=255)
                                                                                                                        ccAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCt                              <  1:512491/150‑1 (MQ=255)
                                                                                                                         cAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTg                             >  1:1186594/1‑150 (MQ=255)
                                                                                                                         cAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTg                             >  2:1002659/1‑150 (MQ=255)
                                                                                                                          aaCACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGa                            >  2:1201660/1‑150 (MQ=255)
                                                                                                                          aaCACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGa                            <  1:1062517/150‑1 (MQ=255)
                                                                                                                           acacGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAg                           <  2:3626053/150‑1 (MQ=255)
                                                                                                                               gATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGc                       >  2:76911/1‑150 (MQ=255)
                                                                                                                                   ggTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCtt                   >  1:4006975/1‑150 (MQ=255)
                                                                                                                                   ggTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCtt                   >  1:3207006/1‑150 (MQ=255)
                                                                                                                                   ggTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCtt                   <  1:3657238/150‑1 (MQ=255)
                                                                                                                                    gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATTCGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                  <  2:2711953/150‑1 (MQ=255)
                                                                                                                                    gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                  >  2:3040639/1‑150 (MQ=255)
                                                                                                                                    gTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCttt                  >  2:1800390/1‑150 (MQ=255)
                                                                                                                                     tAATTCGAAGCGGGGTACGCTCCCCTCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTg                 <  2:3870282/150‑1 (MQ=255)
                                                                                                                                        ttCGAACCGGGCTACGCTCCCCGCTTTGAGCACAATACACAACTCATCTATTGAGCTATACGCAATGTATCAATCTTAGCTTGCAGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGAc              <  2:3191628/150‑1 (MQ=255)
                                                                                                                                           gAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCtt           >  2:64390/1‑150 (MQ=255)
                                                                                                                                           gAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAACCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCtt           >  1:2072423/1‑150 (MQ=255)
                                                                                                                                            aaGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTg          >  1:1428449/1‑150 (MQ=255)
                                                                                                                                            aaGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTg          <  2:4435923/150‑1 (MQ=255)
                                                                                                                                            aaGCGGGGAGCGTACCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTg          >  1:3787489/1‑150 (MQ=255)
                                                                                                                                             aGCGGGGTTCGCTCCCCGCTTTGAGCACATTTCACTACTGATCTTTTGATCTATACGCATTGTATCAATATTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTTGCTCTTTGGACGTTGc         <  1:734851/150‑1 (MQ=255)
                                                                                                                                             aGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGc         <  1:2817699/150‑1 (MQ=255)
                                                                                                                                               cGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCa       <  2:3372497/150‑1 (MQ=255)
                                                                                                                                                 gggTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCaga     <  1:307618/150‑1 (MQ=255)
                                                                                                                                                    tACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCt  <  1:896194/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACATACGCCCCAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTGCCAGAGCT  >  NC_002947/4498485‑4498782

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: