Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,094,035 A→G 5.7% intergenic (+71/+27) tpm → / ← htpX thiopurine S‑methyltransferase/protease HtpX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,094,0350AG 5.7% 210.8 / 2.6 70intergenic (+71/+27)tpm/htpXthiopurine S‑methyltransferase/protease HtpX
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base G (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

GTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCA  >  NC_002947/2093887‑2094180
                                                                                                                                                      |                                                                                                                                                  
gtgATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGaa                                                                                                                                                    >  3:305697/1‑149 (MQ=255)
        ggACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACTTCGGTCGCCCt                                                                                                                                           <  3:368105/149‑1 (MQ=255)
              cccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGttt                                                                                                                                     >  5:54871/1‑149 (MQ=255)
               ccGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGtttt                                                                                                                                    <  7:228713/149‑1 (MQ=255)
                   aGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AGAGGGCGACCGAAGTCGCCCTTTGTTTTTATc                                                                                                                                <  4:264134/149‑1 (MQ=255)
                    gTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTAt                                                                                                                                 <  6:35400/147‑1 (MQ=255)
                    gTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTAt                                                                                                                                 >  5:35400/1‑147 (MQ=255)
                     tttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAg                                                                                                                              >  5:47615/1‑149 (MQ=255)
                     tttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAACAGGTATGGCAGGCTGAACCGCTCCCACAGAAGGCTGCGTAAGGCCGGGATTAAAA‑‑AAAGGGCGCCCGAAGTCGCCCTTTGTTTTTATAAg                                                                                                                              >  1:295043/1‑149 (MQ=255)
                      ttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAACAAAGGGCGACTTCGGTCGCCCTttttttatt                                                                                                                                  >  2:232580/1‑143 (MQ=38)
                      ttGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAAACAAAGGGCGACTTCGGTCGCCCTttttttatt                                                                                                                                  <  1:232580/146‑4 (MQ=38)
                        gtgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGt                                                                                                                            <  6:47615/148‑1 (MQ=255)
                         tgCAGGCCGGGGCGTCGAGTTTTTTTGGGCGGGTTTTTACGGGGAGCCGGTAAGCAGGGATGGCCGGCGGGCCCCGACCCACAACACGGAGCGGAAGGGCCGGTATTAAA‑‑AAAAGGCGCCCGAACTGGCCCGTTTGTTTTTTAAGGTgg                                                                                                                          >  2:349000/1‑149 (MQ=255)
                         tgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGcc                                                                                                                                              <  5:269635/130‑1 (MQ=255)
                         tgCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGcc                                                                                                                                              >  6:269635/1‑130 (MQ=255)
                           cAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGa                                                                                                                         >  2:1390/1‑148 (MQ=255)
                                cGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGcc                                                                                                                   >  3:302034/1‑149 (MQ=255)
                                 ggggCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTTCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCcg                                                                                                                  <  3:88210/149‑1 (MQ=255)
                                    gCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCgcgc                                                                                                               >  4:70990/1‑149 (MQ=255)
                                      gTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCt                                                                                                              >  8:100125/1‑148 (MQ=255)
                                      gTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCt                                                                                                              >  6:313530/1‑148 (MQ=255)
                                       tCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGa                                                                                                            >  7:328891/1‑149 (MQ=255)
                                         gAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACg                                                                                                          <  6:54871/149‑1 (MQ=255)
                                          aGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGc                                                                                                         <  2:35549/149‑1 (MQ=255)
                                            tttGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAg                                                                                                       >  5:89347/1‑149 (MQ=255)
                                              tgttgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAgg                                                                                                      <  4:175889/148‑1 (MQ=255)
                                                ttgGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGc                                                                                                    <  4:204792/148‑1 (MQ=255)
                                                  ggAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTc                                                                                                 >  5:153899/1‑149 (MQ=255)
                                                         gtgtATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCgg                                                                                          >  6:159508/1‑149 (MQ=255)
                                                         gtgtATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCgg                                                                                          >  6:38223/1‑149 (MQ=255)
                                                              tCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTc                                                                                     >  6:42850/1‑149 (MQ=255)
                                                                   tCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCg                                                                                <  4:305697/149‑1 (MQ=255)
                                                                     aGCCGATAAGCAGGTATTGCATGCTGACCCGGTCCCACAGAAGGATGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggt                                                                              <  2:295043/149‑1 (MQ=255)
                                                                         gATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgggt                                                                          >  6:183988/1‑149 (MQ=255)
                                                                         gATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtggg                                                                           >  3:53893/1‑148 (MQ=255)
                                                                         gATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtggg                                                                           <  4:53893/148‑1 (MQ=255)
                                                                          aTAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCggtgggtg                                                                         <  1:1390/149‑1 (MQ=255)
                                                                           tAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGgtgggtgg                                                                        <  7:100125/149‑1 (MQ=255)
                                                                             aGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCt                                                                      <  5:313530/149‑1 (MQ=255)
                                                                              gCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGttt                                                                                                                                     <  8:32249/85‑1 (MQ=255)
                                                                              gCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGttt                                                                                                                                     >  7:32249/1‑85 (MQ=255)
                                                                                  gTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcatca                                                                 <  3:70990/149‑1 (MQ=255)
                                                                                  gTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCtcat                                                                   >  6:136091/1‑147 (MQ=255)
                                                                                    aTGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACAAGAGTCGGCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGGCTCGATTCGGTCTTCCAGCGGTGGGTGGCTCATCAg                                                                >  1:107869/1‑148 (MQ=255)
                                                                                     tGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGc                                                               <  5:159508/148‑1 (MQ=255)
                                                                                         aGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCTTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAAccc                                                          >  3:52345/1‑149 (MQ=255)
                                                                                         aGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAAccc                                                          <  1:349000/149‑1 (MQ=255)
                                                                                           gCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTATTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCgg                                                        >  6:209509/1‑149 (MQ=255)
                                                                                                 ccGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGc                                                  >  2:301146/1‑149 (MQ=255)
                                                                                                     tCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTg                                              >  2:125971/1‑149 (MQ=255)
                                                                                                         acaGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCtt                                           >  3:373575/1‑148 (MQ=255)
                                                                                                         acaGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCtt                                           >  1:191168/1‑148 (MQ=255)
                                                                                                            gAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAgg                                       <  1:301146/149‑1 (MQ=255)
                                                                                                             aaGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCg                                                                                           >  2:110596/1‑96 (MQ=255)
                                                                                                             aaGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCg                                                                                           <  1:110596/96‑1 (MQ=255)
                                                                                                               ggCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCTTGCTTGAGgccg                                    >  6:305477/1‑149 (MQ=255)
                                                                                                                     gTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCAtt                               >  5:258000/1‑148 (MQ=255)
                                                                                                                              ggAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCTTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAg                     <  4:52345/149‑1 (MQ=255)
                                                                                                                              gataaaaac‑‑aaaGGGCGACTTCGGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAg                     >  2:156041/6‑149 (MQ=255)
                                                                                                                                 aTAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGcc                  <  4:302034/149‑1 (MQ=255)
                                                                                                                                   aaaa‑‑aaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCtt                >  8:290510/1‑149 (MQ=255)
                                                                                                                                     aa‑‑aaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAg                                        <  2:330142/123‑1 (MQ=255)
                                                                                                                                     aa‑‑aaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAg                                        >  1:330142/1‑123 (MQ=255)
                                                                                                                                     aa‑‑aaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCACGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCa              >  8:79291/1‑149 (MQ=255)
                                                                                                                                         aaaGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATg            <  4:373575/149‑1 (MQ=255)
                                                                                                                                                gACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCTGCGGTGGGTGGCTCATCAGGAACCCGGCGAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCg     <  5:305477/149‑1 (MQ=255)
                                                                                                                                                gACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCg     >  2:185602/1‑149 (MQ=255)
                                                                                                                                                  ccGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGAAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGc   <  1:125971/149‑1 (MQ=255)
                                                                                                                                                   cGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAgg                                       <  7:139717/112‑1 (MQ=255)
                                                                                                                                                   cGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAgg                                       >  8:139717/1‑112 (MQ=255)
                                                                                                                                                   cGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCa  >  5:249997/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GTGATTGAGGACAGCCCGAAGTTTGTGCAGGCCGGGGCGTCGAGTTTGTTGGAGCGGGTGTATCAGGTCAGCCGATAAGCAGGTATGGCAGGCTGACCCGCTCCCACAGAAGGCTGCGTAAGGCCGGGAATAAAA‑‑AAAGGGCGACCGAAGTCGCCCTTTGTTTTTATCAGGTGGATCAGCCGCGCTGACGCAGGGCCTCGATGCGGTCTTCCAGCGGTGGGTGGCTCATCAGCAACCCGGCCAGGCCGTGCTTGAGGCCGCCATTGATGCCGAAGGCCTTCATGGTGTCCGGCA  >  NC_002947/2093887‑2094180

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: