Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,748,466 G→C 6.8% intergenic (+35/+64) PP_5042 → / ← PP_5043 hypothetical protein/PhoPQ‑activated pathogenicity‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,748,4660GC 6.8% 528.6 / 21.3 156intergenic (+35/+64)PP_5042/PP_5043hypothetical protein/PhoPQ‑activated pathogenicity‑like protein
Reads supporting (aligned to +/- strand):  ref base G (77/68);  new base C (7/4);  total (84/72)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

TAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCTGACCGGG  >  NC_002947/5748317‑5748614
                                                                                                                                                      |                                                                                                                                                    
tAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcgg                                                                                                                                                      >  1:4328116/1‑150 (MQ=255)
tAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcgg                                                                                                                                                      >  1:3929343/1‑150 (MQ=255)
 aaaaGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcggg                                                                                                                                                     <  1:1161013/150‑1 (MQ=255)
 aaaaGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcggg                                                                                                                                                     >  2:3320005/1‑150 (MQ=255)
 aaaaGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcggg                                                                                                                                                     <  2:257356/150‑1 (MQ=255)
  aaaGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcgggc                                                                                                                                                    >  1:1056018/1‑150 (MQ=255)
  aaaGCCGGCTGCGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcgggc                                                                                                                                                    >  2:3757333/1‑150 (MQ=255)
   aaGCCGGCGGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcgggcg                                                                                                                                                   <  2:217378/150‑1 (MQ=255)
    aGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcgggcgg                                                                                                                                                  >  1:4417053/1‑150 (MQ=255)
        ggCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                              <  1:3862459/150‑1 (MQ=255)
           tgtgCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                           <  1:596284/150‑1 (MQ=255)
           tgtgCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                           <  2:2602270/150‑1 (MQ=255)
           tgtgCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                           <  1:3334275/150‑1 (MQ=255)
             tgCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGAtt                                                                                                                                         >  2:1849765/1‑150 (MQ=255)
                aCCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTTGATTCtg                                                                                                                                      >  2:1227831/1‑150 (MQ=255)
                aCCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCtg                                                                                                                                      >  2:2781554/1‑150 (MQ=255)
                aCCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGCTCGAtttat                                                                                                                                      <  1:2532538/150‑4 (MQ=255)
                  cATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCtgtg                                                                                                                                    >  1:2012791/1‑150 (MQ=255)
                     ttATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCg                                                                                                                                 >  2:684205/1‑150 (MQ=255)
                      tATGACACGCGCCGGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCGGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGGGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGa                                                                                                                                <  2:2116876/150‑1 (MQ=255)
                       aTGACACGCGCCCGGTTACCTGCCGCACCCCCCCCCAGGGGGGGCCACGGCGGTGGTAGTGCGCCGAACTGCCCTTGGGGGTTGGGCGTTGAGTCTTGCGCCGGGGTTTAAATCGAGGGCGTGGGGGGGGGCGATCGATTCTTTGGCGAt                                                                                                                               <  2:3875078/150‑1 (MQ=255)
                        tGACACGCGCCCGTATACCTGCCGCAACCCCCCCAAGGTGGGGCCACGGCCGGGGTTTTGCGCATACAACCCCATGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGGGCGGGCGGCGCTCGATTCTGTGGCGATg                                                                                                                              <  1:727621/150‑1 (MQ=255)
                            aaGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGa                                                                                                                          >  2:1326471/1‑150 (MQ=255)
                              gcgcCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGc                                                                                                                        >  2:3578739/1‑150 (MQ=255)
                                   cGGATACCTGCCGCACCCCCCCCCCGGTCGGGCCACGGCCGGGGTATTGCGCCTACACGCCCACGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGAATAAATCGCGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGc                                                                                                                   <  2:1519526/150‑1 (MQ=255)
                                   cGGATACCTGCCGCAACCACCCCACGGTCGGGCCACGGCCGGGGTATTGCGCCTACACGCCCACGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGc                                                                                                                   <  1:3321232/150‑1 (MQ=255)
                                        aCCTGCCGCAACCACCCCACGGTCGGGCCACGGCCGGGGTATTGCGCCTACACGCCCACGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCCAGCTCGGCAtcgc                                                                                                              <  2:1143822/150‑4 (MQ=255)
                                          cTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGc                                                                                                            >  2:3079200/1‑150 (MQ=255)
                                          cTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGc                                                                                                            >  1:2698522/1‑150 (MQ=255)
                                           tGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCa                                                                                                           >  2:3136053/1‑150 (MQ=255)
                                            gCCCCAACCACCCCACGGGCGGGCCCCGGCGGGGTTTTGGCGCCTACACCCCCACGGTGGTGGGGGTTTGTTTCTTGCGCGGGGGAATAAATCGCGCGCCGCGGGGGGGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCAc                                                                                                          <  1:3424612/150‑1 (MQ=255)
                                             ccGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTTGGCATGAAGCACt                                                                                                         >  2:4092374/1‑150 (MQ=255)
                                                cAACCCCCCCAGGGGGGGGCCCGGGCGGGGTTTTGGCGCCTCCCCGCCCAGGGTGGTGGGGGGTTGCGTGTTGCGCGGGGGCATAAATCGCGCGCCGGGGGGGGGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGaa                                                                                                      <  1:922418/150‑1 (MQ=255)
                                                cAACCCCCCCAAGGTGGGGCCACGGCGGGGTTATTGCGCCTACAAGCCCATGGTGGTTGGGCGTTGTTTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGaa                                                                                                      <  1:4089153/150‑1 (MQ=255)
                                                cAACCACCCCCCGGGCGGCCCCCGGCGGGGTTTTTGCGCCTACACCCCCATGTGGGTTGGGCTTTGCGTGTGGCGCGGGGCAATAATCCGCGCGCCGGGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGaa                                                                                                      <  1:169162/150‑1 (MQ=255)
                                                  accaccCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAAcc                                                                                                    >  2:4011072/1‑150 (MQ=255)
                                                   ccaccCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAAccc                                                                                                   >  2:4347733/1‑150 (MQ=255)
                                                     accCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCAc                                                                                                 >  1:4112371/1‑150 (MQ=255)
                                                       cccAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCg                                                                                               >  1:1697897/1‑150 (MQ=255)
                                                        ccAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGc                                                                                              >  1:55917/1‑150 (MQ=255)
                                                         cAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCt                                                                                             >  1:3728219/1‑150 (MQ=255)
                                                         cAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCt                                                                                             >  2:1568487/1‑150 (MQ=255)
                                                         cAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCt                                                                                             >  1:3721817/1‑150 (MQ=255)
                                                           cggccGGGCCACGGCGGGGTTATTGCGCCAACACCCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCGGGGGAATAACTCGCGCGCCGCGGGGGGGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTcc                                                                                           <  1:951428/146‑1 (MQ=255)
                                                                gggCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCAttt                                                                                      <  2:2100708/150‑1 (MQ=255)
                                                                gggCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCAttt                                                                                      >  1:3942106/1‑150 (MQ=255)
                                                                 ggCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTc                                                                                     >  1:3612967/1‑150 (MQ=255)
                                                                  gCCACGGCGGGGTTATTGCGCCTACAAGCCCACGGTGGTTGGGCGTTGCGTCTTGCGCGGGGGAATAAATCGCGCGCCGCGGGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTca                                                                                    <  1:2116034/150‑1 (MQ=255)
                                                                  gCCACGGCCGGGTTTTTGCGCCTACAAGCCCATGGTGGTTGGGCGTTTTTTCTTGCGCCGGGGAATAAATCGCGCGCCGCGGGGGGGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTca                                                                                    <  1:3045042/150‑1 (MQ=255)
                                                                  gCCACGGCCGGGTTTTTGCGCCTACAAGCCCATGGTGGTTGGGCGTTGTGTCTTGCGCCGGGGAATAAATCGCGCGCCGCGGGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTca                                                                                    <  1:2579092/150‑1 (MQ=255)
                                                                            gggTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCgggg                                                                          >  2:4507376/1‑150 (MQ=255)
                                                                             ggTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGa                                                                         >  1:4333367/1‑150 (MQ=255)
                                                                              gTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGACGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGaa                                                                        <  1:2958291/150‑1 (MQ=255)
                                                                               tATTGCGCCTACACGCCCTAGGTGGTTGGGGGTTGTGTATTGCGCGGGGGCTTAAATCGTGCGCCGGGGGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAg                                                                       <  1:1344490/150‑1 (MQ=255)
                                                                                        cTACAAGCCCAAGGTGGTTTGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCAGTTCAGCCAGCGGGGAAGCTCCGGCa                                                               >  2:3489977/1‑150 (MQ=255)
                                                                                           cAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGc                                                            >  1:3910100/1‑150 (MQ=255)
                                                                                              gCCCACGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGAATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCAcc                                                         <  2:2498910/150‑1 (MQ=255)
                                                                                                ccAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCAcccc                                                       >  2:1951273/1‑150 (MQ=255)
                                                                                                ccAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCAcccc                                                       >  1:3636568/1‑150 (MQ=255)
                                                                                                   aGGGGGTGGGGCTTTGAGTCTGGCGCCGGGGCAAAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCt                                                    <  1:2328007/150‑1 (MQ=255)
                                                                                                     ggggtTGGGCGTTGTGTCTTGCGCCGGGGAATAAATCGCGCGCCGCGGGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGa                                                  <  1:2399101/148‑1 (MQ=255)
                                                                                                       ggtTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATc                                                <  2:2012791/150‑1 (MQ=255)
                                                                                                        gtTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATcc                                               >  2:2533771/1‑150 (MQ=255)
                                                                                                              cGTTGTGTCTTGCGCCGGGGCATAAATCGTGCGCCGGGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCtggt                                         <  1:1811341/150‑1 (MQ=255)
                                                                                                              cGTTGAGTCTTGCGCCGGGGAATAATCCGCGCGCGGCGGGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCtggt                                         <  1:4554/150‑1 (MQ=255)
                                                                                                               gTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCtggtg                                        <  2:289151/150‑1 (MQ=255)
                                                                                                                ttGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCtggtgg                                       >  2:1544222/1‑150 (MQ=255)
                                                                                                                  gAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGAt                                     <  2:222596/150‑1 (MQ=255)
                                                                                                                   aGTCTTGCGCCGGGGCATAAATCGCGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATg                                    <  1:3171676/150‑1 (MQ=255)
                                                                                                                     tCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGaa                                                                                                      <  2:3531590/82‑1 (MQ=255)
                                                                                                                     tCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGaa                                                                                                      >  1:3531590/1‑82 (MQ=255)
                                                                                                                     tCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCt                                  >  2:1691735/1‑150 (MQ=255)
                                                                                                                       ttGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCa                                <  2:2486631/150‑1 (MQ=255)
                                                                                                                       ttGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCa                                <  2:3929343/150‑1 (MQ=255)
                                                                                                                       ttGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCa                                >  2:2746550/1‑150 (MQ=255)
                                                                                                                       ttGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCa                                <  2:2517184/150‑1 (MQ=255)
                                                                                                                         gcgcCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCaa                                                              <  1:264040/118‑1 (MQ=255)
                                                                                                                         gcgcCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCaa                                                              >  2:264040/1‑118 (MQ=255)
                                                                                                                          cgcCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGa                             <  1:3590973/150‑1 (MQ=255)
                                                                                                                          cgcCGGGGAATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGa                             <  1:1849765/150‑1 (MQ=255)
                                                                                                                           gcCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGAt                            >  2:1550120/1‑150 (MQ=255)
                                                                                                                             gggggCAAAAAGCGTGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATaa                          <  1:1638481/149‑1 (MQ=255)
                                                                                                                             cGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATaa                          <  1:4312072/150‑1 (MQ=255)
                                                                                                                             cGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATaa                          <  1:2781554/150‑1 (MQ=255)
                                                                                                                                 gCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCt                      >  2:2321491/1‑150 (MQ=255)
                                                                                                                                   aTAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGc                    <  2:3390116/150‑1 (MQ=255)
                                                                                                                                        tCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCaa               >  1:2872374/1‑150 (MQ=255)
                                                                                                                                        tCGAGCGCCGCCCGCGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCaa               >  1:3144141/1‑150 (MQ=255)
                                                                                                                                          gAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCaaaa             >  2:604734/1‑150 (MQ=255)
                                                                                                                                          gAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCaaaa             <  2:1306582/150‑1 (MQ=255)
                                                                                                                                           aGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAAc            <  1:2969688/150‑1 (MQ=255)
                                                                                                                                           aGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAAc            >  2:55751/1‑150 (MQ=255)
                                                                                                                                           aGCGCCGCCCGCGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAAc            >  2:3171251/1‑150 (MQ=255)
                                                                                                                                           aGCGCCGCCCGCGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAAc            >  2:3171161/1‑150 (MQ=255)
                                                                                                                                           aGCGCCGCCCGCGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTAAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAAc            >  1:1117061/1‑150 (MQ=255)
                                                                                                                                              gccgcGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCt         >  1:3347038/1‑150 (MQ=255)
                                                                                                                                                cgcccgCGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCTGa       >  2:509637/1‑150 (MQ=255)
                                                                                                                                                    cgggcggCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCTGACCgg   >  1:718830/1‑150 (MQ=255)
                                                                                                                                                     gggcggCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCTGACCggg  >  1:3304822/1‑150 (MQ=255)
                                                                                                                                                     gggcggCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCTGACCggg  <  2:721849/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCTGACCGGG  >  NC_002947/5748317‑5748614

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