Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,367,848 G→A 8.4% F226F (TTC→TTT ftsH ← integral membrane ATP‑dependent zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,367,8480GA 8.4% 591.5 / 31.6 202F226F (TTC→TTTftsHintegral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (95/90);  new base A (6/11);  total (101/101)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 208 total aligned reads displayed.

GCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCG  >  NC_002947/5367700‑5367997
                                                                                                                                                    |                                                                                                                                                     
gCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGa                                                                                                                                                      >  1:784520/1‑150 (MQ=255)
 cTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGaa                                                                                                                                                     <  2:3949210/150‑1 (MQ=255)
     tcgtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATCATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGtaa                                                                                                                                                 <  2:3088046/150‑3 (MQ=255)
        tGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGaa                                                                                                                                              >  2:806113/1‑150 (MQ=255)
            cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCgg                                                                                                                                          >  2:2169525/1‑150 (MQ=255)
            cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCgg                                                                                                                                          >  1:2008691/1‑150 (MQ=255)
            cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCgg                                                                                                                                          >  2:2186088/1‑150 (MQ=255)
            cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGAGGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACAGGTCGCGCACACGGCTGGCACCGACACCCACGAACATGTCCACGAAGTCCGAAGCgg                                                                                                                                          >  2:5437862/1‑150 (MQ=255)
            cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTGCGTCGAGGACGATGATGCACGGGGCGTGCGGCGTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCgg                                                                                                                                          >  2:3450888/1‑150 (MQ=255)
             gccgccAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGa                                                                                                                                         >  1:3451261/1‑150 (MQ=255)
                gccCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAAt                                                                                                                                      >  1:4984955/1‑150 (MQ=255)
                 cccATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATg                                                                                                                                     >  2:2691911/1‑150 (MQ=255)
                   cATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGt                                                                                                                                   <  1:2132463/150‑1 (MQ=255)
                    aTGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTg                                                                                                                                  <  1:2253841/150‑1 (MQ=255)
                     tGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTga                                                                                                                                 <  1:4462083/150‑1 (MQ=255)
                     gaccGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTga                                                                                                                                 <  1:4572227/148‑1 (MQ=255)
                     gaccGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTga                                                                                                                                 <  1:3400275/148‑1 (MQ=255)
                      gCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                                                                                                <  1:4508804/150‑1 (MQ=255)
                      gCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                                                                                                <  1:283497/150‑1 (MQ=255)
                      gCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                                                                                                <  1:2676662/150‑1 (MQ=255)
                      gCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                                                                                                <  1:2147941/150‑1 (MQ=255)
                         ggCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaagaa                                                                                                                             <  1:2477805/150‑1 (MQ=255)
                          gcgcCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAAc                                                                                                                            >  1:5228603/1‑150 (MQ=255)
                            gcCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACgg                                                                                                                          >  2:2310542/1‑150 (MQ=255)
                              cACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTa                                                                                                                        <  2:4984955/150‑1 (MQ=255)
                                     gCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTttcgc                                                                                                                 <  2:1726420/150‑1 (MQ=255)
                                       gACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTttcgctt                                                                                                               >  1:4213585/1‑150 (MQ=255)
                                        aCCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTttcgcttc                                                                                                              <  1:3146877/150‑1 (MQ=255)
                                        aCCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTttcgcttc                                                                                                              <  2:4229784/150‑1 (MQ=255)
                                            aCGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGccc                                                                                                          <  1:5380843/150‑1 (MQ=255)
                                            aCGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACATTGGCTTCGcca                                                                                                          <  1:3950419/150‑2 (MQ=255)
                                            aCGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACATTGGCTTCGcca                                                                                                          <  1:2560415/150‑2 (MQ=255)
                                            aCGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACATTGGCTTCGcca                                                                                                          <  2:3769776/150‑2 (MQ=255)
                                              ggCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCcgc                                                                                                        <  1:5416588/150‑1 (MQ=255)
                                                 gTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACATTGGCTTCGCCagcaac                                                                                                     >  1:2044072/1‑147 (MQ=255)
                                                  tCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATg                                                                                                    <  2:1259437/150‑1 (MQ=255)
                                                  tCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATg                                                                                                    <  1:3988993/150‑1 (MQ=255)
                                                   cGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATgg                                                                                                   <  1:4456846/150‑1 (MQ=255)
                                                   cGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATgg                                                                                                   >  1:3633302/1‑150 (MQ=255)
                                                        tcgtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTg                                                                                              <  1:3859653/150‑1 (MQ=255)
                                                          gtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGc                                                                                            <  2:68653/150‑1 (MQ=255)
                                                          gtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGc                                                                                            >  1:68647/1‑150 (MQ=255)
                                                               tGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCcagca                                                                                       <  1:4624576/150‑1 (MQ=255)
                                                               tGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCcagca                                                                                       <  1:3249980/150‑1 (MQ=255)
                                                               tGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCAACGAACATTTCCACGAAGTCCGAACCGGAAATGGTg                                                                                                                                  >  2:3561861/1‑107 (MQ=255)
                                                                   gatgatGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGACCCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGt                                                                                   >  2:1635959/1‑150 (MQ=255)
                                                                        tgCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                                                                                                <  2:4058489/100‑1 (MQ=255)
                                                                        tgCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                                                                                                >  1:4058489/1‑100 (MQ=255)
                                                                        tgCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtac                                                                              >  2:4876099/1‑150 (MQ=255)
                                                                        tgCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtac                                                                              >  2:3914652/1‑150 (MQ=255)
                                                                         gCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtacc                                                                             >  1:1588760/1‑150 (MQ=255)
                                                                         gCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtacc                                                                             <  2:3663537/150‑1 (MQ=255)
                                                                         gCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtacc                                                                             <  2:3869947/150‑1 (MQ=255)
                                                                          cACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccg                                                                            >  1:544632/1‑150 (MQ=255)
                                                                          cACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccg                                                                            >  2:4011163/1‑150 (MQ=255)
                                                                             ggggCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccggta                                                                         >  2:765015/1‑150 (MQ=255)
                                                                               ggCGTGCCTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccggtacc                                                                       >  2:3416577/1‑150 (MQ=255)
                                                                                   tGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGt                                                                   <  1:637096/150‑1 (MQ=255)
                                                                                   tGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGt                                                                   <  2:2484175/150‑1 (MQ=255)
                                                                                   tGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGt                                                                   <  1:3519665/150‑1 (MQ=255)
                                                                                     cttcttGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCggg                                                                    >  2:3346832/1‑147 (MQ=255)
                                                                                     cttcttGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTgg                                                                 <  2:1677157/150‑1 (MQ=255)
                                                                                     cttcttGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTgg                                                                 <  1:2202634/150‑1 (MQ=255)
                                                                                      ttcttGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTggg                                                                >  1:4606195/1‑150 (MQ=255)
                                                                                      ttcttGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTggg                                                                <  2:4963186/150‑1 (MQ=255)
                                                                                       tcttGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGc                                                               >  2:1583950/1‑150 (MQ=255)
                                                                                          tGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTACACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCCGGGTCTTACCGGTACCGGGTGGGCCGa                                                            <  1:1077757/150‑1 (MQ=255)
                                                                                           ggCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGAc                                                           <  2:1534990/150‑1 (MQ=255)
                                                                                              cTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCAt                                                        >  2:3578280/1‑150 (MQ=255)
                                                                                               tGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATc                                                       <  2:3828248/150‑1 (MQ=255)
                                                                                               tGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATc                                                       <  2:474135/150‑1 (MQ=255)
                                                                                               tGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATc                                                       >  2:3394344/1‑150 (MQ=255)
                                                                                               tGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATc                                                       >  1:4876055/1‑150 (MQ=255)
                                                                                                 cTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGAc                                                           >  1:111929/1‑144 (MQ=255)
                                                                                                  tCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGc                                                    >  1:3006713/1‑150 (MQ=255)
                                                                                                    gAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGcac                                                  >  1:3483138/1‑150 (MQ=255)
                                                                                                         tGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCTCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACcgc                                             >  1:5034584/1‑150 (MQ=255)
                                                                                                         tGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACcgc                                             >  1:3677803/1‑150 (MQ=255)
                                                                                                               gcACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTGGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGa                                       <  1:366128/150‑1 (MQ=255)
                                                                                                                   acGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATCCg                                   >  1:954476/1‑150 (MQ=255)
                                                                                                                       cTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCa                               >  2:4029096/1‑150 (MQ=255)
                                                                                                                        tGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCAc                              >  1:1259984/1‑150 (MQ=255)
                                                                                                                            aCCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAg                          >  1:702156/1‑150 (MQ=255)
                                                                                                                            aCCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAg                          >  2:1873800/1‑150 (MQ=255)
                                                                                                                            aCCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAg                          >  1:2721803/1‑150 (MQ=255)
                                                                                                                            aCCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAg                          >  2:4523358/1‑150 (MQ=255)
                                                                                                                                  acCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTg                    <  1:429799/150‑1 (MQ=255)
                                                                                                                                  acCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTg                    <  1:163914/150‑1 (MQ=255)
                                                                                                                                    ccACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGa                  >  1:2182059/1‑150 (MQ=255)
                                                                                                                                        gAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACtt              <  1:1760281/150‑1 (MQ=255)
                                                                                                                                            aTTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGccc          <  2:4333694/150‑1 (MQ=255)
                                                                                                                                             tttCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCg         >  2:3202716/1‑150 (MQ=255)
                                                                                                                                               tCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTTCCCGGa       >  2:101137/1‑150 (MQ=255)
                                                                                                                                                 cACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATc     <  2:5228603/150‑1 (MQ=255)
                                                                                                                                                 cACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATc     >  2:1220829/1‑150 (MQ=255)
                                                                                                                                                  aCGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcg    >  2:4942414/1‑150 (MQ=255)
                                                                                                                                                    gAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcgcg  >  2:3601565/1‑150 (MQ=255)
                                                                                                                                                    gAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcgcg  <  2:2034275/150‑1 (MQ=255)
                                                                                                                                                    gAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcgcg  <  2:59818/150‑1 (MQ=255)
                                                                                                                                                    gAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTGACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcgcg  <  2:1439953/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCG  >  NC_002947/5367700‑5367997

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: