Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 977,873 T→A 10.5% intergenic (+39/‑141) cysE → / → iscR serine acetyltransferase/DNA‑binding transcriptional regulator IscR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947977,8730TA10.5% 175.7 / 7.2 64intergenic (+39/‑141)cysE/iscRserine acetyltransferase/DNA‑binding transcriptional regulator IscR
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base A (4/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGC  >  NC_002947/977728‑978012
                                                                                                                                                  |                                                                                                                                           
tGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCt                                                                                                                                          >  6:263236/1‑149 (MQ=255)
 gACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGACATGGTGGCATGCGGGGCGTGTGCt                                                                                                                                          >  5:290524/1‑148 (MQ=255)
    cAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACg                                                                                                                                      >  6:383467/1‑149 (MQ=255)
    cAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGCACACg                                                                                                                                      >  8:352771/1‑149 (MQ=255)
    cAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCACCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGCACACg                                                                                                                                      >  3:319228/1‑149 (MQ=255)
       gATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGccc                                                                                                                                   <  2:178820/149‑1 (MQ=255)
         tGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCt                                                                                                                                 >  8:58023/1‑149 (MQ=255)
           ggCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTc                                                                                                                                <  8:293582/148‑1 (MQ=255)
           ggCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGCACACGccccgca                                                                                                                               <  7:352771/149‑4 (MQ=255)
           ggCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCACCGCACTGAAGTCATGGTGGCATGCGGGGCGTGAGCACACGccccgca                                                                                                                               <  4:319228/149‑4 (MQ=255)
            gCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGTGGGGCGTGAGCACACGccccgcat                                                                                                                              <  3:90032/149‑5 (MQ=255)
                tGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGTGCGTGTGCTCACGCCCCTCGTCTGc                                                                                                                          >  2:238744/1‑149 (MQ=255)
                tGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGc                                                                                                                          <  5:247986/149‑1 (MQ=255)
                   cTATTGTGCCAAGGAACTGCCGGCCTTGCCGGGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCt                                                                                                                         >  1:242938/1‑148 (MQ=255)
                    tATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTAc                                                                                                                       >  2:239359/1‑148 (MQ=255)
                           ccAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCgcgc                                                                                                               >  5:163959/1‑149 (MQ=255)
                            cAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgcc                                                                                                              <  7:155119/149‑1 (MQ=255)
                             aaGGAACTGCCGGCCTTGCC‑GGAGGACCACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCACCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTAACGCCCCTCGTCTGCTACAATCGcgccc                                                                                                             >  2:34878/1‑149 (MQ=255)
                                   cTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTcc                                                                                                       <  7:58023/149‑1 (MQ=255)
                                      ccGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTccc                                                                                                      <  1:156208/147‑1 (MQ=255)
                                      ccGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCt                                                                                                     >  5:101589/1‑148 (MQ=255)
                                      ccGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTg                                                                                                    <  6:302656/149‑1 (MQ=255)
                                       cGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGa                                                                                                   >  3:127021/1‑149 (MQ=255)
                                        ggCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGAt                                                                                                  <  1:239359/149‑1 (MQ=255)
                                         gCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATg                                                                                                 <  2:242938/149‑1 (MQ=255)
                                          ccTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGc                                                                                                <  6:163959/149‑1 (MQ=255)
                                            ttGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCaa                                                                                              >  7:18892/1‑149 (MQ=255)
                                                c‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATccc                                                                                          <  6:101589/149‑1 (MQ=255)
                                                c‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATccc                                                                                          >  4:83079/1‑149 (MQ=255)
                                                  ggaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgcccg                                                                                                            <  1:323769/130‑1 (MQ=255)
                                                  ggaggaCGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGcgcccg                                                                                                            >  2:323769/1‑130 (MQ=255)
                                                           cTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAg                                                                                  >  8:172105/1‑147 (MQ=255)
                                                                         cccAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGt                                                                  <  7:172105/149‑1 (MQ=255)
                                                                          ccAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGtt                                                                 >  5:389960/1‑149 (MQ=255)
                                                                            aattGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGTTGGCATGCTGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGa                                                               <  5:383467/146‑1 (MQ=255)
                                                                              gTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACt                                                             <  6:117688/149‑1 (MQ=255)
                                                                                  cccAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCTACTAAAGCAGAAGGTCTAATAGTTGACTTaaa                                                         <  1:238744/149‑1 (MQ=255)
                                                                                   ccAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAAt                                                        >  2:92339/1‑149 (MQ=255)
                                                                                      gcgcGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCCCGCCCCTCGTCTGCGACAATCGCTCCCGCCTCCCTGATGCAATCCCGCCTAAAGCAGTAGGGCTAATAGTTGACTTAAATGCt                                                     >  7:157739/1‑149 (MQ=255)
                                                                                      gcgcGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCt                                                     >  7:162734/1‑149 (MQ=255)
                                                                                        gcgACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCg                                                   >  7:294156/1‑149 (MQ=255)
                                                                                         cgACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCgg                                                  <  5:263236/149‑1 (MQ=255)
                                                                                             acTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAAt                                              >  2:237760/1‑149 (MQ=255)
                                                                                              cTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATc                                             <  1:92013/149‑1 (MQ=255)
                                                                                               tCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCg                                            >  3:323801/1‑149 (MQ=255)
                                                                                               tCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGCCGTGTGCTAACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGAAGGTCTAATAGATGACTTAAAAGCACCGGAATCg                                            >  7:268767/1‑149 (MQ=255)
                                                                                                 aGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATTCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGc                                           <  1:92339/148‑1 (MQ=255)
                                                                                                     gCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACt                                      >  4:305109/1‑149 (MQ=255)
                                                                                                     gCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACt                                      >  7:84648/1‑149 (MQ=255)
                                                                                                           aTCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGcac                                 <  3:83079/148‑1 (MQ=255)
                                                                                                               cgccgcACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATc                            <  8:42009/149‑1 (MQ=255)
                                                                                                                gccgcACTGAAGTCATGGTGGCATGGGGGGCGTGTGCTCTCGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTATAGCAGTAGGTCTAATAGTTGACTTAATTGCTCGGGAATCGCATTCTCGCACACATcc                           <  1:34878/149‑1 (MQ=255)
                                                                                                                 ccgcACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCt                          >  6:119783/1‑149 (MQ=255)
                                                                                                                           gTCATGGTTGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTTCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCACACACATCCTGAAACTCCCg                <  6:290524/149‑1 (MQ=255)
                                                                                                                           gTCATGGGTGCATGCGGGGCGTGTGCTCACGCCCCTCGGCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCg                <  8:18892/149‑1 (MQ=255)
                                                                                                                            tCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCa                                          >  3:19595/1‑122 (MQ=255)
                                                                                                                            tCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCa                                          <  4:19595/122‑1 (MQ=255)
                                                                                                                            tCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGt               <  5:119783/149‑1 (MQ=255)
                                                                                                                              aTGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTg              <  3:305109/148‑1 (MQ=255)
                                                                                                                                      aTGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCgg                                                  >  5:353857/1‑104 (MQ=255)
                                                                                                                                      aTGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCgg                                                  <  6:353857/104‑1 (MQ=255)
                                                                                                                                       tGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATa    <  1:237760/149‑1 (MQ=255)
                                                                                                                                          ggggCGTGAGCACACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGc  >  2:302339/1‑148 (MQ=255)
                                                                                                                                          ggggCGTGAGCACACGCCCCTCGTCTGCTACAATCACGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTGGCACACATCCTGAAACTCACGTGGTAGACAATAGc  >  3:283530/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
TGACCAAGATGGGCAGTGACTATTGTGCCAAGGAACTGCCGGCCTTGCC‑GGAGGACGACTTCAGCGAGGCTGCCCAAGTGGCCCAGCGCGACACTCAGCCGCATTGATCGCGCCGCACTGAAGTCATGGTGGCATGCGGGGCGTGTGCTCACGCCCCTCGTCTGCTACAATCGCGCCCGCCTCCCTGATGCAATCCCGACTAAAGCAGTAGGTCTAATAGTTGACTTAAATGCTCGGGAATCGCATACTCGCACACATCCTGAAACTCCCGTGGTACCCAATAGC  >  NC_002947/977728‑978012

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