Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,378,068 C→G 8.6% K68N (AAG→AAC)  grpE ← heat shock protein GrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,378,0680CG 8.6% 283.9 / 16.1 93K68N (AAG→AAC) grpEheat shock protein GrpE
Reads supporting (aligned to +/- strand):  ref base C (41/44);  new base G (5/3);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

TCAAGGTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCC  >  NC_002947/5377919‑5378217
                                                                                                                                                     |                                                                                                                                                     
tCAAGGTTGTAGCGCTTGGGGATTTCCTGGATCATCTTCAGGGTCAGCTCGACACCTTCGCGAACCTGCTGAACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAAc                                                                                                                                                       <  1:1971201/150‑1 (MQ=255)
tCAAGGTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAAc                                                                                                                                                       <  2:2129629/150‑1 (MQ=255)
  aaGGTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACtt                                                                                                                                                     >  1:2085567/1‑150 (MQ=255)
  aaGGTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACtt                                                                                                                                                     >  2:1433843/1‑150 (MQ=255)
   aGGTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTtc                                                                                                                                                    <  1:2150815/150‑1 (MQ=255)
     gTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTTGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTtctc                                                                                                                                                  <  1:447863/150‑1 (MQ=255)
       tGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCa                                                                                                                                                <  1:2049903/150‑1 (MQ=255)
       tGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCa                                                                                                                                                <  1:2050475/150‑1 (MQ=255)
       tGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCa                                                                                                                                                <  2:1612045/150‑1 (MQ=255)
         tAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAAGTTCTCCAgc                                                                                                                                              <  2:340541/150‑2 (MQ=255)
         tAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAgg                                                                                                                                              <  1:1777855/150‑1 (MQ=255)
         tAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAgg                                                                                                                                              <  1:1939075/150‑1 (MQ=255)
         tAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAgg                                                                                                                                              <  1:1042878/150‑1 (MQ=255)
           gcgcTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGc                                                                                                                                            >  2:206473/1‑150 (MQ=255)
           gcgcTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGc                                                                                                                                            <  1:206473/150‑1 (MQ=255)
                   gggTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTtgt                                                                                                                                    <  2:2210789/150‑1 (MQ=255)
                     gtgtCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTgg                                                                                                                                  <  1:573229/150‑1 (MQ=255)
                          cTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAAc                                                                                                                                                       <  2:8459/124‑1 (MQ=255)
                          cTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAAc                                                                                                                                                       >  1:8459/1‑124 (MQ=255)
                          cTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTAGGTCTCCGGTATCAAGTCTCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCttt                                                                                                                             >  1:534807/1‑150 (MQ=255)
                           tGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGTTTCAGGCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTtc                                                                                                                            <  1:1613263/150‑1 (MQ=255)
                           tGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTtc                                                                                                                            >  1:347478/1‑150 (MQ=255)
                           tGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTtc                                                                                                                            >  1:1980594/1‑150 (MQ=255)
                                 aTCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAAc                                                                                                                                        <  2:12/132‑1 (MQ=255)
                                 aTCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAAc                                                                                                                                        >  1:12/1‑132 (MQ=255)
                                    ttCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCca                                                                                                                        <  2:180790/145‑1 (MQ=255)
                                    ttCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCca                                                                                                                        >  1:180790/1‑145 (MQ=255)
                                    ttCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATcc                                                                                                                   >  2:1670269/1‑150 (MQ=255)
                                    ttCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATcc                                                                                                                   >  2:332291/1‑150 (MQ=255)
                                        gggTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTtctc                                                                                                                          >  2:7689/1‑139 (MQ=255)
                                        gggTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTtctc                                                                                                                          <  1:7689/139‑1 (MQ=255)
                                         ggTCAGCGCGACACCTTCGCGAAGCTGCTTGACGTGCTCATCGTCGGCGCTGGAGTGAGCCCGGGCGAGTTCCCGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGTACTTCCCCAGGGCGAGCTTCTGGGCTTTCTGCACATCCTGCTc                                                                                                              >  2:953640/1‑150 (MQ=255)
                                           tCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGACCGAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCgg                                                                                                            >  2:211081/1‑150 (MQ=255)
                                                tCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCacga                                                                                                       <  1:2117829/150‑1 (MQ=255)
                                                    cacCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACGTGTGGGCTTTCTCCACATCCTGCTCGGCacgacgac                                                                                                   <  2:2423083/150‑1 (MQ=255)
                                                     acCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAAc                                                                                                                                        >  1:991323/1‑112 (MQ=255)
                                                     acCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAAc                                                                                                                                        <  2:991323/112‑1 (MQ=255)
                                                        ttCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGacgtac                                                                                               >  1:1265631/1‑150 (MQ=255)
                                                        ttCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGacgtac                                                                                               >  1:1814952/1‑150 (MQ=255)
                                                         tCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGacgtacg                                                                                              <  2:347478/150‑1 (MQ=255)
                                                           gcgGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTGCTCCACATCCTGCTCGGCACGACGACGTACGtt                                                                                            >  2:1236041/1‑150 (MQ=255)
                                                             ggATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCt                                                                                          >  1:1248272/1‑150 (MQ=255)
                                                              gATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTg                                                                                         >  2:156428/1‑150 (MQ=255)
                                                                tCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCacgacgac                                                                                                   >  1:860377/1‑138 (MQ=255)
                                                                tCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCacgacgac                                                                                                   <  2:860377/138‑1 (MQ=255)
                                                                 cTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAg                                                                                      <  1:742265/150‑1 (MQ=255)
                                                                  tGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGacgtac                                                                                               <  1:1260496/140‑1 (MQ=255)
                                                                  tGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCAATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGacgtac                                                                                               >  2:1260496/1‑140 (MQ=255)
                                                                   gCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGt                                                                                    >  1:2399251/1‑150 (MQ=255)
                                                                   gCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGACAGTTCCAGGCTGTCGATCACCGGCAGCGAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCGCGCCACGCCGACGTACGTTCTGCAGGt                                                                                    >  2:679615/1‑150 (MQ=255)
                                                                       gACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTg                                                                                         >  2:2185156/1‑141 (MQ=255)
                                                                       gACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTg                                                                                         <  1:2185156/141‑1 (MQ=255)
                                                                          gTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGCTCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCacgacgac                                                                                                   <  1:1343100/128‑1 (MQ=255)
                                                                          gTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCacgacgac                                                                                                   >  2:1343100/1‑128 (MQ=255)
                                                                          gTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTAc                                                                             >  2:2003047/1‑150 (MQ=255)
                                                                                aTCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACtt                                                                                                                                      <  2:1051675/87‑1 (MQ=255)
                                                                                aTCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACtt                                                                                                                                      >  1:1051675/1‑87 (MQ=255)
                                                                                 tcgtcgGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGgcgcgc                                                                      >  2:166700/1‑150 (MQ=255)
                                                                                   gtcgGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAg                                                                    >  1:1937145/1‑150 (MQ=255)
                                                                                        cACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCGGCAGGTCGGCTACGGCGCGCAGCGCCt                                                               >  2:389593/1‑150 (MQ=255)
                                                                                         aCTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTg                                                              >  2:1027837/1‑150 (MQ=255)
                                                                                         aCTGGAGTGAGACAGGGCCAGTTACAGGCTGACTATTTGTGCCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTg                                                              <  2:1203810/150‑1 (MQ=255)
                                                                                          cTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTgg                                                             <  2:534807/150‑1 (MQ=255)
                                                                                                tGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTg                                                       >  2:821987/1‑150 (MQ=255)
                                                                                                tGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTg                                                       <  2:1814952/150‑1 (MQ=255)
                                                                                                      aGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGGCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCg                                                 <  2:1265631/150‑1 (MQ=255)
                                                                                                        ggCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGc                                               <  2:1457347/150‑1 (MQ=255)
                                                                                                        ggCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGc                                               <  1:1111187/150‑1 (MQ=255)
                                                                                                         gCCAGTTCCAGGCTGTCGATCACCGCCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGcc                                              >  1:713030/1‑150 (MQ=255)
                                                                                                              ttCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAgcgc                                                                 <  2:2008098/126‑1 (MQ=255)
                                                                                                              ttCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAgcgc                                                                 >  1:2008098/1‑126 (MQ=255)
                                                                                                                 cAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTc                                      <  1:804563/150‑1 (MQ=255)
                                                                                                                   ggCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCATCCCCTGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCt                                    <  1:156428/150‑1 (MQ=255)
                                                                                                                    gCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTc                                   >  2:1922650/1‑150 (MQ=255)
                                                                                                                      tGTCGATGAAAGGAAGCAAGACGCTGGAGAACTTCTCCAGGTCGAACTTGTGTGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGa                                 >  2:983856/1‑150 (MQ=255)
                                                                                                                      tGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGa                                 >  1:1022194/1‑150 (MQ=255)
                                                                                                                        tCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGc                               <  1:332291/150‑1 (MQ=255)
                                                                                                                          taTCACCCGCAGCAAGTCGCTGGAGAACTGCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTc                             <  2:969418/149‑1 (MQ=255)
                                                                                                                          gATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTc                             >  1:2301572/1‑150 (MQ=255)
                                                                                                                           aTCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTcc                            <  1:1027837/150‑1 (MQ=255)
                                                                                                                                cGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAac                       <  2:2301572/150‑1 (MQ=255)
                                                                                                                                  gcagcaAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAacac                     >  2:993113/1‑150 (MQ=255)
                                                                                                                                  gcagcaAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAacac                     >  1:1574854/1‑150 (MQ=255)
                                                                                                                                  gcagcaAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAacac                     <  1:2003047/150‑1 (MQ=255)
                                                                                                                                   cagcaAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACg                    <  1:821987/150‑1 (MQ=255)
                                                                                                                                    agcaAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACgg                   <  2:1937145/150‑1 (MQ=255)
                                                                                                                                      caAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGt                 >  1:1549808/1‑150 (MQ=255)
                                                                                                                                           cccTGGAGAAGTTCTCCAGGGCGAATGATATTTCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAg            >  1:1637733/3‑150 (MQ=255)
                                                                                                                                           cGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAg            <  2:1248272/150‑1 (MQ=255)
                                                                                                                                            ccTGGAGAAGTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCCGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGAGTATCacc           >  2:367899/2‑148 (MQ=255)
                                                                                                                                            ccTGGAGAAGTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACgg                   >  1:2179338/2‑142 (MQ=255)
                                                                                                                                            ccTGGAGAAGTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACgg                   <  2:2179338/141‑1 (MQ=255)
                                                                                                                                            ccTGGAGAAGTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTCTCAGc           >  2:1738653/2‑150 (MQ=255)
                                                                                                                                            ccTGGAGAAGTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGc           >  2:2098567/2‑150 (MQ=255)
                                                                                                                                             cTGGAGAAGTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTCGTCACCTTGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCt          <  2:1637733/150‑1 (MQ=255)
                                                                                                                                                  gaACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCAtt     <  2:2085567/150‑1 (MQ=255)
                                                                                                                                                    aCTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGc   <  2:1574854/150‑1 (MQ=255)
                                                                                                                                                     cTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTcc                                     >  1:880409/1‑115 (MQ=255)
                                                                                                                                                     cTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTcc                                     <  2:880409/115‑1 (MQ=255)
                                                                                                                                                     cTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGcc  <  1:389593/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAAGGTTGTAGCGCTTGAGGGTGTCCTGGAACATCTTCAGGGTCAGCTCGACACCTTCGCGGATCTGCTTGACGTGCTCATCGTCGGCACTGGAGTGAGCCAGGGCCAGTTCCAGGCTGTCGATCACCGGCAGCAAGTCGCTGGAGAACTTCTCCAGGGCGAACTTGTGGGCTTTCTCCACATCCTGCTCGGCACGACGACGTACGTTCTGCAGGTCGGCTACGGCGCGCAGCGCCTGGTCCTTGGCTGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCC  >  NC_002947/5377919‑5378217

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑