Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,019,127 T→G 5.4% L120R (CTG→CGG)  PP_1797 → HlyD family secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,019,1270TG 5.4% 605.3 / 15.8 177L120R (CTG→CGG) PP_1797HlyD family secretion protein
Reads supporting (aligned to +/- strand):  ref base T (86/81);  new base G (7/3);  total (93/84)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

GGCATCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAACCG  >  NC_002947/2018982‑2019274
                                                                                                                                                 |                                                                                                                                                   
ggCATCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGt                                                                                                                                                 <  2:3044018/150‑1 (MQ=255)
   aTCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGcc                                                                                                                                              >  2:2440896/1‑150 (MQ=255)
    tCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCAGGCGTGcct                                                                                                                                             <  1:4496555/150‑2 (MQ=255)
         aaaGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGt                                                                                                                                        >  1:3716562/1‑150 (MQ=255)
         aaaGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGt                                                                                                                                        >  2:2576972/1‑150 (MQ=255)
           agagCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGGGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGa                                                                                                                                      <  2:969095/150‑1 (MQ=255)
           agagCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGa                                                                                                                                      >  2:2889294/1‑150 (MQ=255)
           agagCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGa                                                                                                                                      >  2:654018/1‑150 (MQ=255)
                 gCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTTTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCtt                                                                                                                                <  1:2916644/150‑1 (MQ=255)
                 gCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCtt                                                                                                                                >  1:4002270/1‑150 (MQ=255)
                   aTGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTgg                                                                                                                              >  1:4856043/1‑150 (MQ=255)
                     gCCCGCGAAGGGGACCTGGTGAAGAAGGGCGATCCGCAGATTTTGCTGGATGACTCGCAGCTCAGCGCCGAGAATGAGACCACCCGCAACCAACTGATTGTTGCCCGTTAAAAAGAGGCACGACTGCGAGCCGAGCGTGACGGCTTGGAc                                                                                                                            >  2:1807109/1‑150 (MQ=255)
                       ccGCGATGGGGACCTGGTGAAGAATGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTc                                                                                                                          >  1:4215225/1‑150 (MQ=255)
                                ggACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTc                                                                                                                 >  2:1631388/1‑150 (MQ=255)
                                 gACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGAGTGACGGCTTGGACTCGATCCCGTCg                                                                                                                <  2:486426/150‑1 (MQ=255)
                                  aCCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCgg                                                                                                               <  1:4573446/150‑1 (MQ=255)
                                   ccTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGt                                                                                                              <  1:3806161/150‑1 (MQ=255)
                                      ggTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGAt                                                                                                           <  1:5182695/150‑1 (MQ=255)
                                       gTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATa                                                                                                          >  2:4577922/1‑150 (MQ=255)
                                        tGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATaa                                                                                                         <  2:3645153/150‑1 (MQ=255)
                                         gaagaaGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAAt                                                                                                        >  2:4987443/1‑150 (MQ=255)
                                          aagaagGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATg                                                                                                       >  2:4281479/1‑150 (MQ=255)
                                             aagGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGaa                                                                                                    <  2:1935782/150‑1 (MQ=255)
                                             aagGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGaa                                                                                                    <  2:1570200/150‑1 (MQ=255)
                                             aagGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGaa                                                                                                    <  1:3445983/150‑1 (MQ=255)
                                              agGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAg                                                                                                   <  1:5130613/150‑1 (MQ=255)
                                               gggCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAgg                                                                                                  <  2:4002270/150‑1 (MQ=255)
                                                ggCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGt                                                                                                 >  1:3596210/1‑150 (MQ=255)
                                                  cGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTAc                                                                                               <  1:4678245/150‑1 (MQ=255)
                                                        gCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTc                                                                                         <  1:1270802/150‑1 (MQ=255)
                                                        gCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTc                                                                                         <  1:4822555/150‑1 (MQ=255)
                                                         cTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCt                                                                                        <  1:3815144/150‑1 (MQ=255)
                                                            aTTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGaa                                                                                                    <  1:5386774/135‑1 (MQ=255)
                                                              tgtgCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACcg                                                                                   >  2:3657295/1‑150 (MQ=255)
                                                                tgCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTg                                                                                                                                                <  2:3171313/87‑1 (MQ=255)
                                                                tgCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTg                                                                                                                                                >  1:3171313/1‑87 (MQ=255)
                                                                  cTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACcgcgcg                                                                               >  2:526268/1‑150 (MQ=255)
                                                                   tGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCgcgcgc                                                                              >  2:73971/1‑150 (MQ=255)
                                                                    ggATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCa                                                                             >  1:1399912/1‑150 (MQ=255)
                                                                     gATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCaa                                                                            >  2:4888362/1‑150 (MQ=255)
                                                                     gATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCaa                                                                            >  1:2661652/1‑150 (MQ=255)
                                                                        gAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCaagaa                                                                         <  1:4770829/150‑1 (MQ=255)
                                                                             gCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCt                                                                    <  2:3882629/150‑1 (MQ=255)
                                                                             gCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCt                                                                    <  2:1236241/150‑1 (MQ=255)
                                                                             gCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCt                                                                    <  2:4854181/150‑1 (MQ=255)
                                                                              cAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTg                                                                   >  2:1534919/1‑150 (MQ=255)
                                                                                  tCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAg                                                               >  2:2429236/1‑150 (MQ=255)
                                                                                  tCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAg                                                               <  1:2105455/150‑1 (MQ=255)
                                                                                     gcgcCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTg                                                            >  1:2929609/1‑150 (MQ=255)
                                                                                       gcCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGag                                                          >  1:967706/1‑150 (MQ=255)
                                                                                        ccGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGagc                                                         <  2:1177974/150‑1 (MQ=255)
                                                                                          gAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATCATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGagcag                                                       >  2:3436151/1‑150 (MQ=255)
                                                                                           aGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGagcagc                                                      <  1:3024654/150‑1 (MQ=255)
                                                                                            gTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAgcagca                                                     <  2:613607/150‑1 (MQ=255)
                                                                                             tATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGcagcag                                                    <  1:5267360/150‑1 (MQ=255)
                                                                                               tGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGgt                                                  >  1:4918933/1‑150 (MQ=255)
                                                                                               tGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGgt                                                  <  2:2901338/150‑1 (MQ=255)
                                                                                                gAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATACTGGAAGGTACCGATTCTGACCGCGCGCaagaag                                                                        >  1:3552193/1‑127 (MQ=255)
                                                                                                gAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCaagaag                                                                        <  2:3552193/127‑1 (MQ=255)
                                                                                                gAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGgtg                                                 >  1:107583/1‑150 (MQ=255)
                                                                                                gAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGgtg                                                 >  1:2071881/1‑150 (MQ=255)
                                                                                                  gTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGtt                                               >  1:1454097/1‑150 (MQ=255)
                                                                                                   tCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTc                                              <  2:1724780/150‑1 (MQ=255)
                                                                                                     cACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACcgcgc                                                                                >  1:856145/1‑114 (MQ=255)
                                                                                                     cACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCaa                                            <  1:2298376/150‑1 (MQ=255)
                                                                                                          gCAACCAACTGATTGGTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGccc                                       <  2:4158361/150‑1 (MQ=255)
                                                                                                          gCAACCAACTGATGGGTGCCCGGTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCATGGACTCGATCCCGTCGGTGATAATTGAAGGTACCGATTCTGACCGCGCGCTAGAAGCGCTGGCAGGTGAGGAGCAGGTGTTCAAGGcca                                       <  2:3813343/150‑2 (MQ=255)
                                                                                                              ccaacTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCa                                   <  2:5321995/150‑1 (MQ=255)
                                                                                                              ccaacTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTAAAGGCCCGCCa                                   >  2:2150581/1‑150 (MQ=255)
                                                                                                               caacTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCAc                                  >  1:4415700/1‑150 (MQ=255)
                                                                                                                aacTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACg                                 >  2:3179389/1‑150 (MQ=255)
                                                                                                                aacTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACg                                 >  1:2929553/1‑150 (MQ=255)
                                                                                                                   tGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAg                                           <  1:3518605/137‑1 (MQ=255)
                                                                                                                   tGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAg                                           >  2:3518605/1‑137 (MQ=255)
                                                                                                                   tGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGAtt                              <  1:3645704/150‑1 (MQ=255)
                                                                                                                    gATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTc                             >  2:1596343/1‑150 (MQ=255)
                                                                                                                    gATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCATCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTc                             <  2:1887705/150‑1 (MQ=255)
                                                                                                                        gttgCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTg                         <  2:4330742/150‑1 (MQ=255)
                                                                                                                           gCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGag                                                          <  1:5427598/114‑1 (MQ=255)
                                                                                                                           gCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGag                                                          >  2:5427598/1‑114 (MQ=255)
                                                                                                                           gCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTg                      <  2:5344817/150‑1 (MQ=255)
                                                                                                                                 tataAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGaa                <  2:2661652/150‑1 (MQ=255)
                                                                                                                                  ataAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGaaa               <  2:2295728/150‑1 (MQ=255)
                                                                                                                                   acgctcggCACGCAGGCGTGCCTAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCCCGATTCGTTGCTGGGCGAAAt              >  1:2293507/7‑150 (MQ=255)
                                                                                                                                    aaaGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAAtc             >  2:1035681/1‑150 (MQ=255)
                                                                                                                                    aaaGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAAtc             >  2:1115888/1‑150 (MQ=255)
                                                                                                                                      ctcggCACGCAGGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCCAGGCCCGCCACGATTCGTTGCTGGGCGAAAcctc           >  2:2289041/4‑150 (MQ=255)
                                                                                                                                       tcggCACGCAGGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTg                      <  2:1704453/136‑1 (MQ=255)
                                                                                                                                       tcggCACGCAGGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTg                      >  1:1704453/3‑138 (MQ=255)
                                                                                                                                       tcggCACGCAGGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCg          >  2:1914358/3‑150 (MQ=255)
                                                                                                                                       gagGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCg          <  1:2573045/150‑1 (MQ=255)
                                                                                                                                        agGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGt                                                             <  1:3251000/98‑1 (MQ=255)
                                                                                                                                        agGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGt                                                             >  2:3251000/1‑98 (MQ=255)
                                                                                                                                          gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCaag                                                                           >  2:5395292/1‑82 (MQ=255)
                                                                                                                                          gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTc       >  1:1651680/1‑150 (MQ=255)
                                                                                                                                          gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTc       >  2:4648234/1‑150 (MQ=255)
                                                                                                                                          gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTc       >  2:3881700/1‑150 (MQ=255)
                                                                                                                                           cACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCa      >  1:4994862/1‑150 (MQ=255)
                                                                                                                                            aCGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGt        >  1:627561/1‑147 (MQ=255)
                                                                                                                                               ccTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAACCg  >  2:611050/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGCATCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAACCG  >  NC_002947/2018982‑2019274

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: