Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 654,406 T→G 8.1% K114Q (AAG→CAG)  ctpL ← methyl‑accepting chemotaxis protein CtpL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947654,4060TG 8.1% 540.1 / 27.6 169K114Q (AAG→CAG) ctpLmethyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base T (101/54);  new base G (8/6);  total (109/60)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 174 total aligned reads displayed.

CAGATGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGTCGC  >  NC_002947/654257‑654555
                                                                                                                                                     |                                                                                                                                                     
 aGATGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTg                                                                                                                                                      <  2:3289468/150‑1 (MQ=255)
    tGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCcg                                                                                                                                                   <  2:3616915/150‑1 (MQ=255)
     gCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCcgg                                                                                                                                                  <  2:1143944/150‑1 (MQ=255)
        gggTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCcggcgg                                                                                                                                               <  2:1320119/150‑1 (MQ=255)
        gggTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCAAGCTGGCcggcga                                                                                                                                               <  2:3439906/150‑2 (MQ=255)
             gCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCcagc                                                                                                                                          <  2:1863504/150‑1 (MQ=255)
               gTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCcgg                                                                                                                                                  >  1:708712/1‑140 (MQ=255)
               gTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGCCGCCGGCAAGCTGGCCGGCGAcccacaa                                                                                                                                        <  1:1696218/150‑6 (MQ=255)
                tCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGt                                                                                                                                       <  2:2760231/150‑1 (MQ=255)
                  tAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTc                                                                                                                                     >  2:1794038/1‑150 (MQ=255)
                   tacctgggCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCg                                                                                                                                    >  2:2484385/6‑150 (MQ=255)
                    gCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGt                                                                                                                                   >  1:2690013/1‑150 (MQ=255)
                     cAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGtt                                                                                                                                  <  2:1941628/150‑1 (MQ=255)
                      aaGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTg                                                                                                                                 >  1:4137306/1‑150 (MQ=255)
                            agcaAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTg                                                                                                                           >  1:4988803/1‑150 (MQ=255)
                                 gTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGc                                                                                                                      >  2:1371283/1‑150 (MQ=255)
                                   aGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCt                                                                                                                    <  2:4307360/150‑1 (MQ=255)
                                    gCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCCGCAGTTCGTTGGCGGTGAAGGCCTg                                                                                                                   >  1:139162/1‑150 (MQ=255)
                                       aTTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAg                                                                                                                >  1:4070993/1‑150 (MQ=255)
                                           gCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTa                                                                                                            >  1:2598094/1‑150 (MQ=255)
                                            ccTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTAt                                                                                                           >  2:4334863/1‑150 (MQ=255)
                                            ccTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTAt                                                                                                           <  1:1230877/150‑1 (MQ=255)
                                            ccTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTAt                                                                                                           <  1:3663515/150‑1 (MQ=255)
                                            ccTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGGTCGTTGGCGGTGAAGGCCTGCAGGCTAt                                                                                                           >  2:720913/1‑150 (MQ=255)
                                             cTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATc                                                                                                          >  2:3444167/1‑150 (MQ=255)
                                             cTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATc                                                                                                          >  1:2064825/1‑150 (MQ=255)
                                             cTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATc                                                                                                          >  2:709430/1‑150 (MQ=255)
                                              tCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATcc                                                                                                         >  2:1175682/1‑150 (MQ=255)
                                                      gCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTggg                                                                                                 >  1:154219/1‑150 (MQ=255)
                                                           tGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCa                                                                                            >  2:1501792/1‑150 (MQ=255)
                                                             tCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAAc                                                                                          >  1:165516/1‑150 (MQ=255)
                                                                  ggCGGAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCGTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTtctc                                                                                     <  2:5279745/150‑1 (MQ=255)
                                                                          tGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGct                                                                             <  2:3246283/150‑1 (MQ=255)
                                                                             ccATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTa                                                                          >  2:4039219/1‑150 (MQ=255)
                                                                              cATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTAc                                                                         >  1:4758159/1‑150 (MQ=255)
                                                                                tCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCg                                                                       >  1:1672976/1‑150 (MQ=255)
                                                                                   gCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGcc                                                                    >  1:3720596/1‑150 (MQ=255)
                                                                                   gCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGcc                                                                    >  1:978286/1‑150 (MQ=255)
                                                                                    cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCt                                                                                                             >  2:894775/1‑108 (MQ=255)
                                                                                    cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCa                                                                   <  2:278303/150‑1 (MQ=255)
                                                                                    cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCa                                                                   >  1:3449127/1‑150 (MQ=255)
                                                                                    cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCa                                                                   <  1:5062841/150‑1 (MQ=255)
                                                                                     tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAg                                                                  >  2:3146389/1‑150 (MQ=255)
                                                                                     tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAg                                                                  >  2:533959/1‑150 (MQ=255)
                                                                                     tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAg                                                                  >  2:367733/1‑150 (MQ=255)
                                                                                     tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAg                                                                  >  1:1835636/1‑150 (MQ=255)
                                                                                       gAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGGCGCCGGCCAGCTTGCCGGCGACCCGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACGGGCCAGTc                                                                >  1:805719/1‑150 (MQ=255)
                                                                                           ttCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCAAGCTGGCcggcga                                                                                                                                               >  1:5170328/1‑66 (MQ=39)
                                                                                           ttCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCAAGCTGGCcggcga                                                                                                                                               <  2:5170328/67‑2 (MQ=39)
                                                                                              gCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGGTCGTTGGCGGTGAAGGCCTGCGGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTctgg                                                         >  1:2512410/1‑147 (MQ=255)
                                                                                               cccccAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTg                                                        <  1:2247297/150‑1 (MQ=255)
                                                                                                ccccAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGc                                                       >  2:1325279/1‑150 (MQ=255)
                                                                                                ccccAGCTCCCGCTCGGCTTGCAGCAGCAGGGCGTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGc                                                       >  2:5223159/1‑150 (MQ=255)
                                                                                                  ccAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCgg                                                     >  1:2333480/1‑150 (MQ=255)
                                                                                                  ccAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCgg                                                     >  2:4219424/1‑150 (MQ=255)
                                                                                                    aGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGcc                                                   >  1:3008940/1‑150 (MQ=255)
                                                                                                         ccGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGcc                                                                                                                     <  1:2213024/79‑1 (MQ=255)
                                                                                                         ccGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgc                                                                                               <  1:3180477/101‑1 (MQ=255)
                                                                                                         ccGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgc                                                                                               >  2:3180477/1‑101 (MQ=255)
                                                                                                           gCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACt                                            <  1:5319730/150‑1 (MQ=255)
                                                                                                           gCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGGTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCaccg                                            >  2:2372621/1‑147 (MQ=255)
                                                                                                             tCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGt                                          >  2:5229101/1‑150 (MQ=255)
                                                                                                                  ttGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGtt                                     >  2:1307029/1‑150 (MQ=255)
                                                                                                                  ttGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGtt                                     >  1:4526505/1‑150 (MQ=255)
                                                                                                                  ttGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGtt                                     <  2:5395679/150‑1 (MQ=255)
                                                                                                                     cagcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGc                                  >  2:407134/1‑150 (MQ=255)
                                                                                                                     cagcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGc                                  <  2:3008940/150‑1 (MQ=255)
                                                                                                                      agcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGcc                                 >  2:4989972/1‑150 (MQ=255)
                                                                                                                      agcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGcc                                 >  2:3874318/1‑150 (MQ=255)
                                                                                                                       gcagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGc                                                    >  1:3753682/1‑130 (MQ=255)
                                                                                                                        cagcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGctg                                                                                                                >  1:3173138/1‑67 (MQ=255)
                                                                                                                         agcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTcc                              <  1:179925/150‑1 (MQ=255)
                                                                                                                         agcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTACCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTcc                              >  1:3932181/1‑150 (MQ=255)
                                                                                                                          gcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGcgacaa                             <  2:723998/150‑6 (MQ=255)
                                                                                                                          gcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCa                             >  2:1053089/1‑150 (MQ=255)
                                                                                                                          gcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCa                             >  1:1830222/1‑150 (MQ=255)
                                                                                                                          gcagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCa                             >  1:4372426/1‑150 (MQ=255)
                                                                                                                           cagGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAt                            <  1:5434469/150‑1 (MQ=255)
                                                                                                                             gggCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATgg                          <  1:5282143/150‑1 (MQ=255)
                                                                                                                              ggCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAt                            <  2:549618/147‑1 (MQ=255)
                                                                                                                              ggCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggc                         >  2:1705451/1‑150 (MQ=255)
                                                                                                                              ggCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggc                         >  2:1950963/1‑150 (MQ=255)
                                                                                                                              ggCTGGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggc                         >  1:5321979/1‑150 (MQ=255)
                                                                                                                                   gctggccGCCGGCAAGCTGGCCGGCGACCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCt                                                  >  1:2146430/7‑120 (MQ=255)
                                                                                                                                   gctggccGCCGGCAAGCTGGCCGGCGACCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCt                                                  <  2:2146430/114‑1 (MQ=255)
                                                                                                                                      ggTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGtt                                     >  1:1132348/1‑130 (MQ=255)
                                                                                                                                        ccGCCGGCAAGCTGGCCGGCGACCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTa                                                                          <  1:4479282/90‑1 (MQ=39)
                                                                                                                                        ccGCCGGCAAGCTGGCCGGCGACCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTa                                                                          >  2:4479282/2‑91 (MQ=39)
                                                                                                                                         cGCCGGCAAGCTGGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTg                                        >  1:1012357/1‑124 (MQ=255)
                                                                                                                                         cGCCGGCAAGCTGGCCGGCGACCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGgc              >  2:3542965/1‑150 (MQ=255)
                                                                                                                                          gCCGGCAAGCTGGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGgcg             >  1:1260475/1‑150 (MQ=255)
                                                                                                                                             ggCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAg          <  2:1547607/150‑1 (MQ=255)
                                                                                                                                             ggCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAg          >  1:5065726/1‑150 (MQ=255)
                                                                                                                                              gCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAgg         <  2:2621976/150‑1 (MQ=255)
                                                                                                                                                   cTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGTc    >  1:77950/1‑150 (MQ=255)
                                                                                                                                                   cTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGTc    >  2:4247079/1‑150 (MQ=255)
                                                                                                                                                   cTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGTc    <  1:407134/150‑1 (MQ=255)
                                                                                                                                                   cTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGTc    <  1:1023166/150‑1 (MQ=255)
                                                                                                                                                    ttGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGTcg   >  2:5155651/1‑150 (MQ=255)
                                                                                                                                                     tGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGTcgc  >  1:1876128/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGATGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGTCGC  >  NC_002947/654257‑654555

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: