Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,646,916 A→G 6.7% intergenic (‑46/+20) PP_3211 ← / ← PP_3212 ABC transporter ATP‑binding protein/Rieske 2Fe‑2S family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,646,9160AG 6.7% 793.7 / 28.2 251intergenic (‑46/+20)PP_3211/PP_3212ABC transporter ATP‑binding protein/Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base A (119/115);  new base G (8/9);  total (127/124)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

TAGGTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACACAGG  >  NC_002947/3646767‑3647063
                                                                                                                                                     |                                                                                                                                                   
tAGGTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCa                                                                                                                                                     >  1:8982913/1‑150 (MQ=255)
   gTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTATATGCCCCGGCGCACCg                                                                                                                                                  <  1:5338484/150‑1 (MQ=255)
     tttttCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCggg                                                                                                                                                <  2:7966225/150‑1 (MQ=255)
      ttttCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCgggg                                                                                                                                               >  1:6830944/1‑150 (MQ=255)
        ttCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGTGCGCCGGGGca                                                                                                                                             <  2:7623405/150‑2 (MQ=255)
        ttCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGTGCGCCGGGGca                                                                                                                                             <  1:7357842/150‑2 (MQ=255)
        gtCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGTGCGCCGGGGca                                                                                                                                             <  2:5662759/149‑2 (MQ=255)
         tCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGTGCGCCGGGGcat                                                                                                                                            <  1:2002646/150‑3 (MQ=255)
         tCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGt                                                                                                                                            <  1:2630429/150‑1 (MQ=255)
          ccACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGtt                                                                                                                                           <  2:7337027/150‑1 (MQ=255)
          ccACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTACGCCCCGGTGCGCCGGGGcata                                                                                                                                           <  1:3295841/150‑4 (MQ=255)
           cACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTg                                                                                                                                          >  2:7855574/1‑150 (MQ=255)
           cACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTg                                                                                                                                          >  2:476508/1‑150 (MQ=255)
                  ttGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAgt                                                                                                                                   <  1:5481126/150‑1 (MQ=255)
                  ttGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAgt                                                                                                                                   <  1:8434086/150‑1 (MQ=255)
                     gCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTc                                                                                                                                <  2:307323/150‑1 (MQ=255)
                     gCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTc                                                                                                                                >  2:2941958/1‑150 (MQ=255)
                      cGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCa                                                                                                                               <  1:2204332/150‑1 (MQ=255)
                       gAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAg                                                                                                                              >  1:6875982/1‑150 (MQ=255)
                       gAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAg                                                                                                                              >  2:2675465/1‑150 (MQ=255)
                          cAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTgg                                                                                                                           >  1:8895189/1‑150 (MQ=255)
                          cAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTgg                                                                                                                           >  2:6576497/1‑150 (MQ=255)
                           aaaGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGc                                                                                                                          >  2:3435973/1‑150 (MQ=255)
                              ggggTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGa                                                                                                                       >  1:8221095/1‑150 (MQ=255)
                              ggggTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGa                                                                                                                       >  2:6362275/1‑150 (MQ=255)
                              ggggTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGa                                                                                                                       <  1:771661/150‑1 (MQ=255)
                                  tgtCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGaa                                                                                                                   >  1:2277660/1‑150 (MQ=255)
                                    tCCACACCGCGGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCGTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGGTGGGCACGTGTCAGTGGCGGAGGAaca                                                                                                                 >  2:7649272/1‑150 (MQ=255)
                                    tCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAaca                                                                                                                 >  2:5215205/1‑150 (MQ=255)
                                     ccACACCGATCGTGTTCGGCAGGGCCTGCGGCTCTACGTCGGCCTGGGATAGGGTCTCCATCAACATAGCGCTCTCCCCGTAGGTCATCAGCGGTGGTGTATGCACAGGCGCACCGGGGCGTTCGGCAAGTGTCAGTGGCTGAGGAacac                                                                                                                >  1:1173278/1‑150 (MQ=255)
                                      cacaCCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACt                                                                                                               >  1:4225889/1‑150 (MQ=255)
                                       acacCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCACGTGTCAGTGGCTGAGGAACACtt                                                                                                              >  2:114240/1‑150 (MQ=255)
                                       acacCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACtt                                                                                                              >  1:3698488/1‑150 (MQ=255)
                                        cacCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTc                                                                                                             <  1:6793903/150‑1 (MQ=255)
                                         acCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGa                                                                                                                       <  2:5193924/139‑1 (MQ=255)
                                         acCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGa                                                                                                                       >  1:5193924/1‑139 (MQ=255)
                                         acCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCACGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCg                                                                                                            >  1:6656588/1‑150 (MQ=255)
                                             aTGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGctct                                                                                                        <  1:873745/150‑1 (MQ=255)
                                             aTGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGctct                                                                                                        >  1:583592/1‑150 (MQ=255)
                                                  gtacGGCAGGGCCTGCGGCTGGACGTCGGCATTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCa                                                                                                   >  2:1775122/4‑150 (MQ=255)
                                                   ttCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCAt                                                                                                  >  1:1382587/1‑150 (MQ=255)
                                                       gCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGc                                                                                              >  2:4122936/1‑150 (MQ=255)
                                                       gCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGc                                                                                              <  2:666408/150‑1 (MQ=255)
                                                          gggCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTc                                                                                           <  2:6371284/150‑1 (MQ=255)
                                                             ccTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAg                                                                                        <  1:6240293/150‑1 (MQ=255)
                                                                 cGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGa                                                                                    <  1:7570695/150‑1 (MQ=255)
                                                                 cGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGa                                                                                    <  2:4225889/150‑1 (MQ=255)
                                                                    cTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGAccc                                                                                 >  2:6189839/1‑150 (MQ=255)
                                                                     tGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCa                                                                                >  1:7552941/1‑150 (MQ=255)
                                                                     tGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCGCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCTAGTGTCAGTGTCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGTCCCa                                                                                <  1:2494707/150‑1 (MQ=255)
                                                                           tCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGtt                                                                          <  1:6500104/150‑1 (MQ=255)
                                                                               cgttGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCg                                                                      >  1:1100984/3‑150 (MQ=255)
                                                                                cTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGa                                                                                                                       >  2:2085920/1‑100 (MQ=255)
                                                                                   ggATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACgg                                                                  >  1:1231678/1‑150 (MQ=255)
                                                                                   ggATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACgg                                                                  <  1:8190981/150‑1 (MQ=255)
                                                                                     aTAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGAt                                                                <  1:5739029/150‑1 (MQ=255)
                                                                                          gTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATaacca                                                           >  2:1436095/1‑150 (MQ=255)
                                                                                           tCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATaaccaa                                                          >  1:4550199/1‑150 (MQ=255)
                                                                                              gCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACAATTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAAccc                                                       <  2:744119/150‑1 (MQ=255)
                                                                                               cATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCg                                                      >  2:4501624/1‑150 (MQ=255)
                                                                                                aTCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCgg                                                     >  2:2558716/1‑150 (MQ=255)
                                                                                                aTCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCgg                                                     <  1:1514529/150‑1 (MQ=255)
                                                                                                 tCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCggg                                                    >  1:725255/1‑150 (MQ=255)
                                                                                                 tCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTATATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCggg                                                    >  2:8565757/1‑150 (MQ=255)
                                                                                                     cATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTg                                                <  1:7259304/150‑1 (MQ=255)
                                                                                                     cATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTAAGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGCGGACCACTTCGCTCTCGCCAGGGGCTTCAAGCAGACCCAGGAGTTGCCGACGGATAAACACCCCGGGCTTg                                                >  1:4116301/1‑150 (MQ=255)
                                                                                                       tAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTc                                              <  1:3985330/150‑1 (MQ=255)
                                                                                                          cgcTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGa                                           <  2:5876474/150‑1 (MQ=255)
                                                                                                            ctctCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGCCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATg                                         >  1:2673589/1‑150 (MQ=255)
                                                                                                             tctcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATgg                                        >  2:3794698/1‑150 (MQ=255)
                                                                                                               tcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCa                                      <  2:8221095/150‑1 (MQ=255)
                                                                                                               tcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCa                                      >  2:9056825/1‑150 (MQ=255)
                                                                                                                 cccGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtg                                    >  2:486944/1‑150 (MQ=255)
                                                                                                                  ccGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAATCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtgt                                   >  2:5499160/1‑150 (MQ=255)
                                                                                                                  ccGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtgt                                   >  1:3901953/1‑150 (MQ=255)
                                                                                                                   cGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtgtg                                  <  1:9067764/150‑1 (MQ=255)
                                                                                                                   cGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtgtg                                  <  1:5595391/150‑1 (MQ=255)
                                                                                                                     tGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCa                                <  2:3556286/150‑1 (MQ=255)
                                                                                                                       ggTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGt                              <  1:1693186/150‑1 (MQ=255)
                                                                                                                       ggTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGt                              <  2:2875721/150‑1 (MQ=255)
                                                                                                                        gTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGtt                             <  2:8385377/150‑1 (MQ=255)
                                                                                                                         tCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGa                                                                     >  2:2563070/1‑109 (MQ=255)
                                                                                                                         tCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTc                            >  1:6996581/1‑150 (MQ=255)
                                                                                                                                 ggtggtGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTgc                    >  2:8838813/1‑150 (MQ=255)
                                                                                                                                  gtggtgTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTgcg                   >  1:3327480/1‑150 (MQ=255)
                                                                                                                                   tggtgTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCagt                                                                              >  1:8918958/1‑90 (MQ=255)
                                                                                                                                     gtgtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTgcggcg                >  1:1865848/1‑150 (MQ=255)
                                                                                                                                      tgtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGcggcgg               <  1:4366376/150‑1 (MQ=255)
                                                                                                                                      tgtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGcggcgg               <  2:2277660/150‑1 (MQ=255)
                                                                                                                                       gtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCggg              >  1:8624055/1‑150 (MQ=255)
                                                                                                                                       gtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCggg              >  1:5272551/1‑150 (MQ=255)
                                                                                                                                       caacgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCggg              >  1:4839420/5‑150 (MQ=255)
                                                                                                                                       caacgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCggg              >  1:2993139/5‑150 (MQ=255)
                                                                                                                                       caacgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCggg              >  1:3189255/5‑150 (MQ=255)
                                                                                                                                        aacgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGgt             >  2:3424247/4‑150 (MQ=255)
                                                                                                                                            ccccGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGc                                                                        >  2:7961396/1‑87 (MQ=255)
                                                                                                                                            ccccGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGc                                                                        <  1:7961396/87‑1 (MQ=255)
                                                                                                                                            ccccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCg         >  2:6346185/1‑150 (MQ=255)
                                                                                                                                             cccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGa        >  2:3305314/1‑150 (MQ=255)
                                                                                                                                                   gcACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACACAgg  >  1:5703254/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TAGGTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACACAGG  >  NC_002947/3646767‑3647063

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: