Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,486,703 C→G 37.8% A159A (GCC→GCG yhdZ → amino acid ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,486,7030CG37.8% 390.8 / inf 272A159A (GCC→GCGyhdZamino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (84/85);  new base G (40/62);  total (124/148)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 275 total aligned reads displayed.

AACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAA  >  NC_002947/1486554‑1486850
                                                                                                                                                     |                                                                                                                                                   
aaCTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGcc                                                                                                                                                     <  1:1655746/150‑1 (MQ=255)
  cTGCACCCTGGCGCCGATGTGGGTGAGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCAt                                                                                                                                                   >  1:4733426/1‑150 (MQ=255)
   tGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCAtt                                                                                                                                                  >  1:1742070/1‑150 (MQ=255)
     cACCCTGGCACCGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTATCAGGCGGCCAGCAGCAGCGTGTGGCGATTGc                                                                                                                                                <  2:3200571/150‑1 (MQ=255)
        ccTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCcg                                                                                                                                             >  2:1540187/1‑150 (MQ=255)
        ccTGGCACCGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcg                                                                                                                                             >  2:1617093/1‑150 (MQ=255)
        ccTGGCACCGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcg                                                                                                                                             >  2:615671/1‑150 (MQ=255)
             caccGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcgtgcgc                                                                                                                                        >  2:4158773/3‑149 (MQ=255)
              accGATTTGGGTACGCAAGATGCCAAAGAAGAAGCCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcgtgcgct                                                                                                                                       >  2:3299194/2‑148 (MQ=255)
                cGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTtgt                                                                                                                                     >  1:1069578/1‑150 (MQ=255)
                cGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcgtgcgcttt                                                                                                                                     <  1:4149224/150‑5 (MQ=255)
                 gATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTtgtg                                                                                                                                    >  2:2548601/1‑150 (MQ=255)
                  ctgtgGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTATGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGc                                                                                                                                   <  2:4278518/149‑1 (MQ=255)
                    gtgGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCAt                                                                                                                                 <  1:2021056/150‑1 (MQ=255)
                        gTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAg                                                                                                                             >  2:1357745/1‑150 (MQ=255)
                           cgcAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGccc                                                                                                                          <  2:2381735/150‑1 (MQ=255)
                                gATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGAt                                                                                                                     >  1:1569988/1‑150 (MQ=255)
                                gATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGAt                                                                                                                     >  1:4363884/1‑150 (MQ=255)
                                   gCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCAt                                                                                                                  <  1:2374593/150‑1 (MQ=255)
                                      aaagaagaaGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCt                                                                                                               >  1:2150254/8‑150 (MQ=255)
                                            gaaGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTcga                                                                                                         >  1:4294175/2‑150 (MQ=255)
                                            gaaGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTcga                                                                                                         >  1:2068573/2‑150 (MQ=255)
                                            gaaGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCACGCTGTTcga                                                                                                         <  1:1421035/149‑1 (MQ=255)
                                                gCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAg                                                                                                     >  1:1286812/1‑150 (MQ=255)
                                                   gAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCAt                                                                                                                                 <  1:45977/119‑1 (MQ=255)
                                                   gAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCg                                                                                                  >  1:5061104/1‑150 (MQ=255)
                                                    aagAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGa                                                                                                 <  2:3920086/148‑1 (MQ=12)
                                                     tgAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGAc                                                                                                <  1:3355875/149‑1 (MQ=255)
                                                         atcgCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCg                                                                                            >  2:4679217/4‑150 (MQ=12)
                                                         aTTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCg                                                                                            >  1:395364/1‑150 (MQ=255)
                                                          tcgCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATCTCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCgg                                                                                           <  2:1840044/148‑1 (MQ=12)
                                                          tcgCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCgg                                                                                           >  2:717697/3‑150 (MQ=12)
                                                             ccATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGcgc                                                                                        >  1:4766939/1‑148 (MQ=12)
                                                              cATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACt                                                                                       <  2:3224485/150‑1 (MQ=255)
                                                               aTGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGt                                                                                                             >  1:2214403/1‑127 (MQ=255)
                                                               aTGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTg                                                                                      <  1:5292179/150‑1 (MQ=255)
                                                               aTGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTg                                                                                      >  1:4998781/1‑150 (MQ=255)
                                                                tGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTgg                                                                                     >  1:1872455/1‑150 (MQ=255)
                                                                     tACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCg                                                                                >  1:1935590/1‑150 (MQ=255)
                                                                       ccTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGa                                                                              >  1:3683626/1‑150 (MQ=255)
                                                                        cTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGaa                                                                             >  2:3727401/1‑150 (MQ=255)
                                                                          ggAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAAt                                                                           >  2:3403146/1‑150 (MQ=255)
                                                                          ggAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAAt                                                                           >  1:2851776/1‑150 (MQ=255)
                                                                           gAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATg                                                                          <  2:3613275/150‑1 (MQ=255)
                                                                              cGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTc                                                                       <  2:1651751/150‑1 (MQ=255)
                                                                                tgtgCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCaa                                                                     <  2:4989395/150‑1 (MQ=255)
                                                                                tgtgCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCaa                                                                     <  1:2726966/150‑1 (MQ=255)
                                                                                cgtgaaaaTTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGTAAATGGTCaa                                                                     <  2:488983/143‑1 (MQ=9)
                                                                                cgtgaaaaTTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCaa                                                                     <  1:244296/143‑1 (MQ=9)
                                                                                 gtgaaaaTTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAg                                                                    >  1:996997/7‑150 (MQ=9)
                                                                                    cGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGaa                                                                 <  2:5351942/150‑1 (MQ=255)
                                                                                    cGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGaa                                                                 >  2:4826658/1‑150 (MQ=255)
                                                                                     ggATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAgtttg                                                                >  1:4488125/1‑146 (MQ=255)
                                                                                     ggATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAg                                                                <  1:483757/150‑1 (MQ=255)
                                                                                     aaaTTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAg                                                                <  2:2787372/148‑1 (MQ=12)
                                                                                         tCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCt                                                            <  1:780459/150‑1 (MQ=17)
                                                                                          ccGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTc                                                           <  1:1137743/150‑1 (MQ=255)
                                                                                          ccGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCtg                                                           <  1:3299194/150‑2 (MQ=12)
                                                                                          ccGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCtg                                                           <  1:3572698/150‑2 (MQ=12)
                                                                                           cGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCtgg                                                          <  1:862253/150‑3 (MQ=9)
                                                                                            ggAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGa                                                         >  1:5419962/1‑150 (MQ=255)
                                                                                            ggAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGa                                                         >  2:2106395/1‑150 (MQ=255)
                                                                                            ggAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGa                                                         <  1:1675353/150‑1 (MQ=9)
                                                                                             gaacAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGAt                                                        <  1:1760732/147‑1 (MQ=9)
                                                                                              aacAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCATCAGCAGCGTGTGGCTATTTCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAATTGGTCAAGGAAGTGCTGGATa                                                       <  1:3146286/148‑1 (MQ=12)
                                                                                                  ggccaacaAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGAt                                                   >  1:3525585/6‑150 (MQ=9)
                                                                                                   gCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCAtg                                                  >  2:3080033/1‑150 (MQ=255)
                                                                                                       acaAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGAtggtga                                              <  1:292953/150‑2 (MQ=17)
                                                                                                         aaGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGc                                            >  2:5032299/1‑150 (MQ=255)
                                                                                                               ccgggCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGcc                                      <  1:900089/147‑1 (MQ=17)
                                                                                                                 tGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGa                                    >  2:5133311/1‑150 (MQ=255)
                                                                                                                 tGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGa                                    >  2:198388/1‑150 (MQ=255)
                                                                                                                 tGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGa                                    >  2:2744634/1‑150 (MQ=255)
                                                                                                                 tGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGa                                    >  2:231463/1‑150 (MQ=255)
                                                                                                                 gggCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGa                                    >  2:349921/2‑150 (MQ=17)
                                                                                                                  ggCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAg                                   <  2:5305897/150‑1 (MQ=255)
                                                                                                                  ggCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAg                                   <  1:5246679/150‑1 (MQ=255)
                                                                                                                  ggCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGat                                   <  1:4973802/150‑2 (MQ=17)
                                                                                                                   gCCAGCTGTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGatg                                  <  2:4812631/150‑3 (MQ=17)
                                                                                                                      agctttCTGGCGGGCAGCAGCAGCGTGTAGCCATCGCGCGTTCACTCTGCGTTAAACCATCGGTCCTGCTATGCGACGAGCCAACCTCCGCACTCGACCCGGAAATGGTCAAGGAGGTGCTGGAAACGTTGGCAGACCTGGCGCAGGacg                               <  2:1137841/145‑3 (MQ=255)
                                                                                                                      aGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGaa                                                                 >  2:4750375/1‑116 (MQ=255)
                                                                                                                      aGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATg                               <  1:2785874/150‑1 (MQ=255)
                                                                                                                          ttcaggcggcCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCTGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCAt                           <  2:50576/146‑1 (MQ=12)
                                                                                                                          gTCCGGCGGGCAGCAGCAGCGTGTCGCCATAGCGCGTTCCTTGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGATGGTGTCGCTCGCCGAGGATGGTAt                           >  2:2828821/1‑150 (MQ=255)
                                                                                                                             tGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGAc                        >  2:3870417/1‑150 (MQ=255)
                                                                                                                             tGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATAGCATGAc                        >  2:1821061/1‑150 (MQ=255)
                                                                                                                              ggcggcCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGAcc                       >  1:362780/1‑150 (MQ=255)
                                                                                                                               gcggcCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGACGAGGATGGCATGACCa                      >  2:2583821/1‑150 (MQ=255)
                                                                                                                                     agcagcagcGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTAAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGctct                <  1:2314757/150‑1 (MQ=255)
                                                                                                                                       cagcagcGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTgt              <  2:5229061/150‑2 (MQ=17)
                                                                                                                                       cagcagcGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTgt              <  1:2621891/150‑2 (MQ=17)
                                                                                                                                       cagcagcGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCg                                                                                            >  1:4685318/1‑72 (MQ=255)
                                                                                                                                       cagcagcGTGTCGCCATAGCGCGTTCCTTGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGATGGTGTCGCTCGCCGAGGATGGTATGACGATGATCTgt              >  2:3882208/1‑149 (MQ=255)
                                                                                                                                           agcGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGa          <  1:4038360/150‑1 (MQ=17)
                                                                                                                                           agcGTGTCGCCATAGCGCGTTCCTTGTGCATGCAGCCCGACATCATCCTGTTCGACGAGCCGACCTCAGCCCTCGACCCTGAACTGGTGCGCGAAGTTCTCGATACGATGGTGTCGCTCGCCGAGGATGGTATGACGATGATCTgtgtaa          <  2:2968739/150‑3 (MQ=255)
                                                                                                                                             cGTGTGGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGAcg        <  1:493291/150‑2 (MQ=12)
                                                                                                                                             cGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGAcc        <  1:2899261/150‑1 (MQ=255)
                                                                                                                                                gtgGCGATTGCCCGTGCGCTTTGCATGAACCCGAAGGTCATGCTGTTCGACGAGCCGACCTCGGCGCTCGATCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCAc     <  1:408399/150‑1 (MQ=12)
                                                                                                                                                   gCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGaa  >  1:4445459/1‑150 (MQ=255)
                                                                                                                                                     cATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTc                                                           <  2:4530668/91‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGTTCGACGAGCCGACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAA  >  NC_002947/1486554‑1486850

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: