Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,949,938 T→G 9.9% Q26P (CAG→CCG)  PP_3483 ← type II secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,949,9380TG 9.9% 396.4 / 2.9 150Q26P (CAG→CCG) PP_3483type II secretion system protein
Reads supporting (aligned to +/- strand):  ref base T (95/35);  new base G (7/12);  total (103/47)
Fisher's exact test for biased strand distribution p-value = 2.86e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 152 total aligned reads displayed.

TCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCTGTG  >  NC_002947/3949815‑3950087
                                                                                                                           |                                                                                                                                                     
tCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACCTCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGAc                                                                                                                             >  2:5996013/1‑150 (MQ=255)
  tgTGGCTTGCAGAGGGGCCGGGGCCTGGCGGTCGGGGGGGTCGGCGGAAATTTTTGCGGGGTTCTTGGGGGGCTGGCCCGCTCCCACCGGTACATCACCGGTCTGTGGGCTGCTTACATCGTGTGGGAGGGGGCATGCCCGCGAAGACGt                                                                                                                           <  2:4935180/149‑1 (MQ=255)
    tGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCTGCATACATCCGGTGGGAGCGGGCATGCCCGCGAAGACGTCg                                                                                                                         >  2:7159067/1‑150 (MQ=255)
      gCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCAGTGGGAGCGGGCATGCCCGCGAAGACCTCGGt                                                                                                                       >  2:1014852/1‑150 (MQ=255)
      gCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCGCCCACAGGTACCTCCCCGGGCGGTGGGCGGCCTACCTCCGGGGGGAGCGGGGCGGGCCGCGCAGACGTCggg                                                                                                                       >  1:2594187/1‑149 (MQ=255)
        ttGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACCTCACAGGTCGGTGGGCTGGCTGAATCCGGTTGGAGCGGGCCTGCCCGCGGAGACGGCGGTGc                                                                                                                     >  1:2935684/1‑150 (MQ=255)
        ttGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGACGACGTCGGTGc                                                                                                                     >  2:7543137/1‑150 (MQ=255)
        ttGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGc                                                                                                                     >  2:5937727/1‑150 (MQ=255)
        ttGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGc                                                                                                                     >  2:910157/1‑150 (MQ=255)
        ttGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCGCCCACAGGTACATCCCAGGTCTGGGGGCTGCCTACATCCTGTGGGAGCGGGCCGGCCCCCGGAGACCTCggggg                                                                                                                     >  1:8693562/1‑147 (MQ=255)
               ggggcctgggcctGGAGGTGGGGGGTGTAGGGGTAATTTTTTGGTGGGTTCTTCGGGGGCTTGCCCGCGCCCACCGGTACCGCACAGGTCTGTGGGCTGAATACATCCTGTGGGCGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCttt                                                                                                              <  2:6995951/149‑1 (MQ=255)
               cGGGCCGGGGCCTGGTGGGCGCGGGGGTCGGCGGAATTTTTGGCTGGGTTCTTCGCGGGCAGGCCCGCGCCCACCGGTACACCACAGGTCTGTGGGCTGCATACACCCGGGGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCttt                                                                                                              <  2:6995826/150‑1 (MQ=255)
                gggCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACCGGTACATCACAGGTCGGTGGGCTGCATACATCCGGTGAGAGCGGGCCAGCCCGCGAAGACGGCGGGGCAGGCTTTg                                                                                                             >  1:4191425/1‑150 (MQ=255)
                 ggcgggggCCAGGAGGGCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGc                                                                                                            <  2:6661216/145‑1 (MQ=255)
                     ggggCCCGGAGGGCGGGGGGGTCGGCGGAATTTTTGGCTGGGTTGTGCGCGGGCAGGCCCGCTCCCACCGGTACATCACCGGTCTGTGGGCTGCATACATCCTGGGGGCGCGGGAATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGAt                                                                                                        <  2:3384271/149‑1 (MQ=255)
                        gCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTcc                                                                                                     >  2:3853252/1‑150 (MQ=255)
                         ccAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTGTGCAGATTCCg                                                                                                    >  2:6373746/1‑150 (MQ=255)
                          cAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCGGTGGGAGCGGGCCTGCCCGCGACGACGTCGGTGCAGGCTTTGCAGATTCggg                                                                                                   >  2:209825/1‑150 (MQ=255)
                          cAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACAACCTGTGGGAGCGGGCATGCCCCCGAAGACCTCGGTGCAGGCTTTGCAGATTCCgg                                                                                                   >  2:209787/1‑150 (MQ=255)
                           aggaggTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGACATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCggg                                                                                                  >  1:4229746/1‑150 (MQ=255)
                            ggaggTCGCGGGGTGTGGGGGAAATTTTTGCTGGGTTCTTGGGGGGCTGGCCCGCTCCCCCCGGTCCATCCCCGGTCTGTGGGCTGCATCCCGCGGGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCgggg                                                                                                 <  2:7218530/150‑1 (MQ=255)
                                   gcgGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCGGTGGGAGCGGGCCTGCCCCCGCAGAAGTCGGTGCAGGCGTTGCAGAATCCGGGGTTGCATc                                                                                          >  1:5915771/1‑150 (MQ=39)
                                              ggtaattttttGGTGGTTTGTGCGGGGCCTGGCCCGCCCCCACGGTTACCCCCCCGGTCGGGGGGCTGCATCCACCCTGGGGGCGGGGGATTGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTg                                                                               <  2:773220/144‑1 (MQ=255)
                                               ggAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTgg                                                                              >  2:4330040/1‑150 (MQ=255)
                                               ggAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCGGCATACATCCTGTGGGAGCGGGCCAGCCCGCGACGACGGCGGTGCAGGCGTTGCAGATGCCGGGGGTGCCTCGAATGCGTCTgg                                                                              >  1:5814524/1‑150 (MQ=38)
                                                gAAATTTTGGGTGTGTTCTTGGCGGTGCTGCCGGCTCCCCCAGTTACATCCCTGGTCTGTGGTCTGCATACCTGCTTTGTGAGCGGGAATGCCCGCGTTGACGTCGGTGCAGGCTTTGCAGATTCCGGGGGTGCATCGAATGCTTCTGGc                                                                             <  2:7790515/150‑1 (MQ=255)
                                                 aaaTTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCa                                                                            >  2:8784307/1‑150 (MQ=255)
                                                                     gcgGACATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAg                                                        <  2:8924791/150‑1 (MQ=255)
                                                                       gggCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAggg                                                      >  1:2165993/1‑150 (MQ=255)
                                                                       gggCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAggg                                                      >  1:5236790/1‑150 (MQ=255)
                                                                       gggCAGGCCCGCTCCCCCCGGTCCCTCCCCGGTCTGTGGGCTGCATCCATCGGGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAggg                                                      <  2:8971578/150‑1 (MQ=255)
                                                                       gggCAGGCCCGCTCCCACCGGTCCACCCCAGGTCTGTGGGCTGCATACAGCCTGTGGGCGGGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAggg                                                      <  2:1211725/150‑1 (MQ=255)
                                                                             gCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGTCCTGCGACTCCGCTCAGGGGCGCGt                                                >  2:8683169/1‑150 (MQ=255)
                                                                              ccggCTCCCACTGGTAAATCACTGGTGTGTGGGCTTCATACATCGTGTGGGTGGGGGCATGCCCGCTAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTc                                               <  2:1163944/147‑1 (MQ=255)
                                                                              cccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTc                                               >  1:4722140/1‑150 (MQ=255)
                                                                              cccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTc                                               >  1:4722144/1‑150 (MQ=255)
                                                                               ccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                                              >  2:102949/1‑150 (MQ=255)
                                                                               ccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                                              >  2:1209677/1‑150 (MQ=255)
                                                                               ccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                                              >  1:2911042/1‑150 (MQ=255)
                                                                               ccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                                              >  2:1988551/1‑150 (MQ=255)
                                                                               ccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                                              >  2:5013018/1‑150 (MQ=255)
                                                                               ccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                                              >  2:5013013/1‑150 (MQ=255)
                                                                               ccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCGCAGGGGCGCGTCg                                              >  1:8038627/1‑150 (MQ=255)
                                                                                cGCTCCCACGGTTACCCCCCGGGGCTGTGGGCTGAATACATCCTGTGGGAGGGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGc                                             <  1:2837858/150‑1 (MQ=255)
                                                                                cGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGc                                             >  1:1132189/1‑150 (MQ=255)
                                                                                cGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGc                                             >  1:4610633/1‑150 (MQ=255)
                                                                                cGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGc                                             >  1:4678620/1‑150 (MQ=255)
                                                                                 gatCCAACAGGTAAAACACTCGTCTGGGGCCTGAATACATCCTTTGGGAGCGGGATTGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCa                                            <  2:2300644/148‑1 (MQ=255)
                                                                                 gCTCCCACGGGGACATCACTGGTCTGTGGTCTGCATACCTCGTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCa                                            <  2:8739279/150‑1 (MQ=255)
                                                                                 gCTCCCACCGGTACCCCCCCGGTCGGTGGGCTGCATACACCCGGGGGGCGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCa                                            <  2:1062905/150‑1 (MQ=255)
                                                                                   gcccACTGGTACATCCCAGGCCTGTGGGCTGAATACATCCTGTGGGAGCGGGAATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAgg                                          <  2:62995/149‑1 (MQ=255)
                                                                                     ccACCGGTACATCTCTGGTCTGTGGGCTGAATCCATCCTGTGGGAGCGGGCATGCCCGGGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggct                                        <  2:2256843/150‑1 (MQ=255)
                                                                                     ccACCGGTACAGCACAGGTCTGTGGTCTGCATCCATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCGGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggct                                        <  2:2935684/150‑1 (MQ=255)
                                                                                     ccACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggct                                        >  2:360907/1‑150 (MQ=255)
                                                                                     ccACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggct                                        >  1:1100897/1‑150 (MQ=255)
                                                                                     ccACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggct                                        >  2:360883/1‑150 (MQ=255)
                                                                                     ccACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggct                                        >  2:7413134/1‑150 (MQ=255)
                                                                                        ccggTACATCCCCGGTCTGTGGGCTGCATACAGCCGGGGGGCGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggctggc                                     <  2:3879770/148‑1 (MQ=255)
                                                                                         cggTACCTCCCCGGTGGGTGGGCTGCATACATCCTGGGGGTGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggctggct                                    <  2:893791/149‑1 (MQ=255)
                                                                                           gTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCg                                  >  2:4435511/1‑150 (MQ=255)
                                                                                           gTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCg                                  >  1:8780673/1‑150 (MQ=255)
                                                                                           gTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCg                                  >  2:5703784/1‑150 (MQ=255)
                                                                                           gTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCg                                  >  1:2714505/1‑150 (MQ=255)
                                                                                            tACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCgg                                 >  2:2895124/1‑150 (MQ=255)
                                                                                             tcATCACAGGTCTGTGGGCTGCTTACATCCTGTGGGCGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGc                                <  2:6041826/149‑1 (MQ=255)
                                                                                               atccccggTCGGTGGGATGAATACATCCTGTGGGCGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCa                              <  1:6765065/144‑1 (MQ=255)
                                                                                               atccccggTCGGTGGCCTGCATACATCCGGTGGGCGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCa                              <  1:6764955/144‑1 (MQ=255)
                                                                                                tCACCGGTCTGTGGGCTGAATACATCCTGTGGGCGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAg                             <  2:3021463/150‑1 (MQ=255)
                                                                                                tCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAg                             >  1:5009820/1‑150 (MQ=255)
                                                                                                tCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAg                             >  2:2304766/1‑150 (MQ=255)
                                                                                                 cacaGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAgg                            >  2:1185108/1‑150 (MQ=255)
                                                                                                 aacaGGTCTGTGTTATGCATACATCCTGTGGGAGCGGGCATGCCCGCGAATAAGTATGTGCATGCTTTGCAGATTCAGGGGTTGCATAGAATGCTTCTTGAATGGCCTGCGACTCCGCTCAGGGGCGCGTCTCAGGCTGTCTCGGCCAgg                            >  2:683112/2‑150 (MQ=255)
                                                                                                  acaGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAggg                           >  1:6240060/1‑150 (MQ=255)
                                                                                                   caGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAgggg                          >  2:7943797/1‑150 (MQ=255)
                                                                                                   caGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCGTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAgggg                          >  2:7944299/1‑150 (MQ=255)
                                                                                                     ggTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTg                        >  1:1461373/1‑150 (MQ=255)
                                                                                                        cTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCgg                     >  2:3359588/1‑150 (MQ=255)
                                                                                                         tgtgGGCTGCATACATCCTGTGGGAGCGGGTATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCggg                    >  2:5597083/1‑150 (MQ=255)
                                                                                                         tgtgGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCAGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCggg                    >  2:4316742/1‑150 (MQ=255)
                                                                                                           tgGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGgcg                  >  1:991734/1‑150 (MQ=255)
                                                                                                            gggCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGgcgc                 >  2:4432136/1‑150 (MQ=255)
                                                                                                             ggcggaaTACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGcgcc                <  2:1708600/144‑1 (MQ=255)
                                                                                                             ggCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGcgcc                >  2:980540/1‑150 (MQ=255)
                                                                                                             ggCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGcgcc                >  1:2418690/1‑150 (MQ=255)
                                                                                                             ggCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGcgcc                >  2:4730548/1‑150 (MQ=255)
                                                                                                                 gCATCCATCCGGTGGGCGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTa            <  2:424475/150‑1 (MQ=255)
                                                                                                                 gCATACCCCCGGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTa            <  2:424474/150‑1 (MQ=255)
                                                                                                                 gCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTa            <  2:7250090/150‑1 (MQ=255)
                                                                                                                   aTACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCGAGGGGTGCGGGCGCGGCTAtt          >  2:7295207/1‑150 (MQ=255)
                                                                                                                   aTACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTAtt          >  2:2414305/1‑150 (MQ=255)
                                                                                                                   aTACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTAtt          >  1:6160198/1‑150 (MQ=255)
                                                                                                                   aTACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTAtt          >  2:9101039/1‑150 (MQ=255)
                                                                                                                   aTACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGGTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTAtt          >  1:7923097/1‑150 (MQ=255)
                                                                                                                      cATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTccc       >  2:7373983/1‑150 (MQ=255)
                                                                                                                      cATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTccc       >  2:7373537/1‑150 (MQ=255)
                                                                                                                      cATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTccc       >  1:7815459/1‑150 (MQ=255)
                                                                                                                       aTCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTcccc      >  2:8855535/1‑150 (MQ=255)
                                                                                                                         ccTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTccc       <  2:8665796/147‑1 (MQ=255)
                                                                                                                           tgtgGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCtgtg  >  1:8142977/1‑150 (MQ=255)
                                                                                                                           tgtgGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCtgtg  >  1:7230212/1‑150 (MQ=255)
                                                                                                                           |                                                                                                                                                     
TCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCTGTG  >  NC_002947/3949815‑3950087

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