Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,642,285 G→C 11.8% intergenic (‑113/+39) PP_1442 ← / ← lon‑I hypothetical protein/DNA‑binding ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,642,2850GC11.8% 188.2 / 16.2 68intergenic (‑113/+39)PP_1442/lon‑Ihypothetical protein/DNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base G (34/26);  new base C (5/3);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

AGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGG  >  NC_002947/1642144‑1642431
                                                                                                                                             |                                                                                                                                                  
aGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGt                                                                                                                                             >  1:355362/1‑149 (MQ=255)
aGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAggc                                                                                                                                             <  5:271910/149‑2 (MQ=38)
    gAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTg                                                                                                                                                  >  4:63633/1‑140 (MQ=255)
    gAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTg                                                                                                                                                  <  3:63633/140‑1 (MQ=255)
                               tggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGt                                                                                                                                             <  1:372487/118‑1 (MQ=255)
                               tggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGt                                                                                                                                             >  2:372487/1‑118 (MQ=255)
                                 gCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCt                                                                                                            <  6:48066/149‑1 (MQ=255)
                                  cTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTa                                                                                                           <  6:327216/149‑1 (MQ=255)
                                  cTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTa                                                                                                           <  2:88231/149‑1 (MQ=255)
                                   tCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGc                                                                                                             >  7:154289/1‑146 (MQ=255)
                                     cGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGTTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGCTTGCCCGCGAaga                                                                                                                          <  1:132943/130‑1 (MQ=255)
                                     cGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGTTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGCTTGCCCGCGAaga                                                                                                                          >  2:132943/1‑130 (MQ=255)
                                     cGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATATGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGaa                                                                                                        <  4:100913/149‑1 (MQ=255)
                                      gAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGaa                                                                                                        >  3:121144/1‑148 (MQ=255)
                                       aaGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                                      <  4:349805/149‑1 (MQ=255)
                                                          ggATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGc                                                                                   <  2:235312/149‑1 (MQ=255)
                                                           gATAACCGATGATCGAGGAAATTAGTGCTGCCTGTGCCGGTCCTTTCTCGAGTAAGCCAGCGCCTGCAGGGGTAGAGTACACGTGCAGGTGCGGGCTTGCCCGCGAAGAGCGCGGTGCAGGCTATCACAGCAATTTAGCCACGTCGGc                                                                                   >  6:17674/1‑148 (MQ=255)
                                                               cccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAg                                                                               <  4:121144/148‑1 (MQ=255)
                                                                  gatgatCGAGGAAATTTGTGCTGCCTGTGCTGGTCCTTTCGCGAGTAAGCCAGCGCCTGCAGGGGTAGCGTACACGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCt                                                                           >  7:113060/1‑149 (MQ=255)
                                                                   atgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCtt                                                                          <  2:31893/149‑1 (MQ=255)
                                                                   atgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCtt                                                                          >  8:28879/1‑149 (MQ=255)
                                                                       tCGAGGAAATTTGTACTGCCTGTGCCGGTCCTTTCTCGAGTAAGCCCGCTCCTACATGGTTAGCGTACCCGTGCAGGTGCGGGCTTACCCGCGAAGAGGGAGGTGCAGGATAGAACAGCACTTTATCCACTTCGGAAAAGCGCTTgac                                                                       >  8:255644/1‑146 (MQ=255)
                                                                          aGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGaaa                                                                   >  3:353903/1‑149 (MQ=255)
                                                                          aGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGaaa                                                                   >  7:9495/1‑149 (MQ=255)
                                                                                tttGTCCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCCCCTGCAGCGGTAGCGTACCCGTGCAGGTCCGGGCATGCCCCCGAAGAGGGCGCTGCAGGCTAGAACAGCACTATCGCCACGACGGCAAACCGCTTGGCGAAATGCa                                                               >  2:13461/1‑147 (MQ=255)
                                                                                    tgctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAaga                                                                                                                          >  2:289831/1‑84 (MQ=255)
                                                                                      ctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCa                                                               >  6:89752/1‑141 (MQ=255)
                                                                                      ctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCa                                                               <  5:89752/141‑1 (MQ=255)
                                                                                      ctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAgg                                                       >  7:340136/1‑149 (MQ=255)
                                                                                       tgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGc                                                      <  2:355362/149‑1 (MQ=255)
                                                                                           tgtgCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTc                                                  >  5:330655/1‑149 (MQ=255)
                                                                                             tgCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTcgc                                                >  8:111712/1‑149 (MQ=255)
                                                                                               ccGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGc                                                      <  6:341995/141‑1 (MQ=255)
                                                                                               ccGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGc                                                      >  5:341995/1‑141 (MQ=255)
                                                                                                   tCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGt                                          >  1:267859/1‑149 (MQ=255)
                                                                                                    ccTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGTCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAgg                                           >  6:401735/1‑147 (MQ=255)
                                                                                                      tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGt                                       >  2:195658/1‑149 (MQ=255)
                                                                                                          gcgAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGTCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCgg                                   <  5:401735/149‑1 (MQ=255)
                                                                                                          gcgAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCgg                                   <  2:95012/149‑1 (MQ=255)
                                                                                                           ggAGTATGCGCGGGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGGCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGAGGTTGGCGAAATGCACTGTGAGGCCTTCTCGCAGGTAGTCCGGc                                  <  5:17674/148‑1 (MQ=255)
                                                                                                            gAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCgg                                   >  4:240700/1‑147 (MQ=255)
                                                                                                              gTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGt                               <  8:340136/149‑1 (MQ=255)
                                                                                                                  gCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttctt                           >  5:117095/1‑149 (MQ=255)
                                                                                                                  gCCCGCACCTGCACGGGTACGCTACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCAGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttctt                           >  2:375270/1‑149 (MQ=38)
                                                                                                                   cccGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                          >  8:56851/1‑149 (MQ=255)
                                                                                                                    ccGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                          >  7:121590/1‑148 (MQ=255)
                                                                                                                    ccGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAAGGACCGGCACAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                          >  5:366932/1‑148 (MQ=18)
                                                                                                                     cgcaCCTGCACGGCAACCCTACCCCTGCAGGCGCGGACTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGCCAAACCGCTTGGCGAAAAGCACTGTCAGGCCTTCGCGCACGTACTCCTGCAGTTCTTCg                         >  2:401621/1‑148 (MQ=255)
                                                                                                                      gcgcCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCg                         >  8:292968/1‑147 (MQ=255)
                                                                                                                      gcaccTGCACGGGTACGCTACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCg                         >  5:219839/4‑147 (MQ=39)
                                                                                                                      gcaccTGCACGGGTACGCTACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                       <  6:219839/146‑1 (MQ=39)
                                                                                                                        ggcTGCAGGTGTATGGTACCCGTGCAGGTTCGGGCTTTCCCTAGAAGAGGGCGGTGCAGGCTAGATCAGAACTTTCGCCACGTCGTCAAAGCGCTTGGTGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                       <  3:58596/145‑1 (MQ=255)
                                                                                                                          cTGCACGGGTACGCTACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAgg                                           <  8:386848/125‑1 (MQ=38)
                                                                                                                          cTGCACGGGTACGCTACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAgg                                           >  7:386848/1‑125 (MQ=38)
                                                                                                                           tGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                          >  8:347131/1‑141 (MQ=255)
                                                                                                                           tGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                          <  7:347131/141‑1 (MQ=255)
                                                                                                                              aGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTcgccgc               >  2:71267/1‑149 (MQ=255)
                                                                                                                                gggTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCgg             <  3:240700/149‑1 (MQ=255)
                                                                                                                                gggTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCgg             >  8:50930/1‑149 (MQ=255)
                                                                                                                                 ggTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACt                                                             >  4:346441/1‑100 (MQ=255)
                                                                                                                                 ggTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACt                                                             <  3:346441/100‑1 (MQ=255)
                                                                                                                                  gtagcgtagcCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACTGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCTGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGtt           <  1:13461/149‑1 (MQ=255)
                                                                                                                                        tACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTc                                                  <  7:45961/104‑1 (MQ=255)
                                                                                                                                        tACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTc                                                  >  8:45961/1‑104 (MQ=255)
                                                                                                                                        tACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTc     >  2:113325/1‑149 (MQ=255)
                                                                                                                                        tACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTc     <  1:71267/149‑1 (MQ=255)
                                                                                                                                         aCCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGtt           >  7:341141/1‑142 (MQ=255)
                                                                                                                                         aCCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGtt           <  8:341141/142‑1 (MQ=255)
                                                                                                                                             gTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCgg  >  8:200665/1‑147 (MQ=255)
                                                                                                                                             gTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCgg  <  6:254274/147‑1 (MQ=255)
                                                                                                                                             gTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCgg  >  5:254274/1‑147 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
AGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGG  >  NC_002947/1642144‑1642431

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