Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 979,599 Δ1 bp 6.1% coding (1049/1215 nt) iscS‑I → cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947979,5990A. 6.1% 209.9 / 7.2 66coding (1049/1215 nt)iscS‑Icysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base . (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CTGAACCTGAGCTTCAACTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGACATCTCCAAG  >  NC_002947/979451‑979744
                                                                                                                                                     |                                                                                                                                                  
ctgaacctgaGCTTCAACTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGca                                                                                                                                                    <  3:235090/149‑1 (MQ=255)
   accctgaTCTTTCACTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGcaca                                                                                                                                                 <  1:129620/146‑1 (MQ=255)
   aacctgaGCTTCAACTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAg                                                                                                                                                <  8:348753/149‑1 (MQ=255)
         aGCTTCAACTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGTGCGCCAGCTCgt                                                                                                                                           <  4:368315/148‑2 (MQ=255)
          gCTTCAACTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTCCGTACTGCGTGCTCTGGGCCGCAACGACGAGGTGTGCGCCAGCTCgtcg                                                                                                                                         <  1:87794/149‑4 (MQ=255)
             tCAGCTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGc                                                                                                                                      >  3:278834/1‑149 (MQ=255)
             tCAACTACGTCGAAGGCGAATCGCTGCTTATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGc                                                                                                                                      >  5:236944/1‑149 (MQ=255)
             tCAACTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGc                                                                                                                                      >  6:379680/1‑149 (MQ=255)
                aCTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTtc                                                                                                                                   <  5:47310/149‑1 (MQ=255)
                aCTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTtc                                                                                                                                   <  3:334928/149‑1 (MQ=255)
                   aCGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTtctc                                                                                                                                 <  4:188421/148‑1 (MQ=255)
                       cGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCg                                                                                                                            >  4:62844/1‑149 (MQ=255)
                         aaGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGc                                                                                                                          <  3:3569/149‑1 (MQ=255)
                                       ctgATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAAGCTCGATCCGCTTCTCCTTCGGCCGTTTCaccacc                                                                                                              <  1:61595/148‑1 (MQ=255)
                                            gTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCgaag                                                                                                          >  6:226357/1‑146 (MQ=255)
                                               gCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGaagaagaag                                                                                                    >  7:70212/1‑149 (MQ=255)
                                               gCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGaagaagaag                                                                                                    >  3:291009/1‑149 (MQ=255)
                                               gCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGaagaagaag                                                                                                    >  8:184524/1‑149 (MQ=255)
                                                cTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGt                                                                                                   <  7:184524/149‑1 (MQ=255)
                                                   aaaGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGa                                                                                                <  5:223197/149‑1 (MQ=255)
                                                         aTCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTAcgc                                                                                          >  3:338924/1‑149 (MQ=255)
                                                         aTCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTAcgc                                                                                          >  2:353938/1‑149 (MQ=255)
                                                          tCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCTCACAGCTCGATCCGCTTCTCCTTCGGCCGttt                                                                                                                     <  7:55665/121‑1 (MQ=255)
                                                          tCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGttt                                                                                                                     >  8:55665/1‑121 (MQ=255)
                                                            gCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTAcgccgc                                                                                       <  6:236944/149‑1 (MQ=255)
                                                                    gTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGaaa                                                                                <  5:339460/148‑1 (MQ=255)
                                                                    gTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGaaa                                                                                <  5:283987/148‑1 (MQ=255)
                                                                    gTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGa                                                                                  >  6:283987/1‑146 (MQ=255)
                                                                    gTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAg                                                                               >  2:195097/1‑149 (MQ=255)
                                                                         gTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGc                                                                          <  7:39148/149‑1 (MQ=255)
                                                                         gTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGc                                                                          <  3:62844/149‑1 (MQ=255)
                                                                          ttCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGc                                                                          >  1:15488/1‑148 (MQ=255)
                                                                             ggCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGACTTTCTGCgcg                                                                       <  2:186508/148‑3 (MQ=255)
                                                                             ggCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGATTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGatg                                                                      <  8:223207/149‑3 (MQ=255)
                                                                               ccTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGcc                                                                    <  1:306259/149‑1 (MQ=255)
                                                                                  gTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGt                                                                  >  7:4693/1‑148 (MQ=255)
                                                                                   tACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGt                                                                  <  1:353938/147‑1 (MQ=255)
                                                                                      cTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAAt                                                              >  1:11941/1‑148 (MQ=255)
                                                                                         cgcggggCTCGAGCCGTCGTAAGTACTGCTTGCTGTGGGCCGCTACGACGAGCTGGCGCTCAGCTCGATCCGCTTCTCCTTCGGCCGTTTCTCCACCGAAGAAGAAGTCGACTACGCCGCGCAGATAGTCTGCGAGGCCGTCAATaat                                                           <  1:74370/142‑2 (MQ=255)
                                                                                         cgcgTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAAc                                                          >  2:40608/1‑149 (MQ=255)
                                                                                             tCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCg                                                      >  8:116732/1‑149 (MQ=255)
                                                                                               gCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCg                                                      >  3:345334/1‑147 (MQ=255)
                                                                                               gCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCg                                                      <  4:345334/147‑1 (MQ=255)
                                                                                               gCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGt                                                     <  4:278834/148‑1 (MQ=255)
                                                                                                cTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGa                                                   <  8:70212/149‑1 (MQ=255)
                                                                                                cTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGa                                                   <  8:4693/149‑1 (MQ=255)
                                                                                                  cGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGc                                                 >  7:200522/1‑149 (MQ=255)
                                                                                                   gAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTGTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCt                                                >  6:3920/1‑149 (MQ=255)
                                                                                                      ccGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTc                                             <  5:226357/149‑1 (MQ=255)
                                                                                                      ccGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTc                                             >  1:37285/1‑149 (MQ=255)
                                                                                                        gTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTcgc                                           <  1:40608/149‑1 (MQ=255)
                                                                                                               gTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCtgtg                                    >  5:108586/1‑149 (MQ=255)
                                                                                                                tACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTgggg                                   <  6:108586/149‑1 (MQ=255)
                                                                                                                 aCTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCtgtg                                    <  5:217531/147‑1 (MQ=255)
                                                                                                                 aCTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCtgtg                                    >  6:217531/1‑147 (MQ=255)
                                                                                                                     cGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGata                               <  5:379680/148‑1 (MQ=255)
                                                                                                                          tctGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGatat                              <  6:96902/144‑1 (MQ=255)
                                                                                                                          tctGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGatat                              >  5:96902/1‑144 (MQ=255)
                                                                                                                          tctGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTAc                          >  8:256793/1‑148 (MQ=255)
                                                                                                                            tGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACa                         <  7:116732/147‑1 (MQ=255)
                                                                                                                              ggCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGa                     >  3:357666/1‑149 (MQ=255)
                                                                                                                                   cAACCACCAGCTGGCCCACAGCTCGATCCGCTTCTCCTTCGGCCCTTTCACCACCCAAGAAGAACACGACTACGCCGCGCAGAAAGTCACCGAGGCCGGCAATAAACTGCGTGAGCTCTCGCAGCTGTGGCATATGTACAAACACGGCg                >  8:277315/1‑149 (MQ=255)
                                                                                                                                    aaCGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGt               >  3:68677/1‑149 (MQ=255)
                                                                                                                                      gatcgaGCTGTGCGCCAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCTAGAAGAAAGTCTGCGAGTAAGTCAATAAAATGCGTGAGATGTCGACGCTGTGGTATATGTACAAAGACGGAGtt              >  5:185734/4‑148 (MQ=255)
                                                                                                                                      gatcgaGCTGTGCGCCAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGtt              >  7:155761/4‑148 (MQ=255)
                                                                                                                                      cgacgaGCTGGCGCACAGGTGGTTCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGTAGAAGTCGACTACGCCGGGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTGGCCGCTGTGGGATATGTACAAAGACGGCGTTg             <  8:200522/149‑1 (MQ=255)
                                                                                                                                                  cgcACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGACATCTCCAAg  <  2:11941/149‑1 (MQ=255)
                                                                                                                                                   gcACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGACATCTcc     >  3:364791/1‑145 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTGAACCTGAGCTTCAACTACGTCGAAGGCGAATCGCTGCTGATGTCGCTGAAAGACATCGCCGTATCGTCCGGTTCGGCCTGTACCTCCGCGTCGCTCGAGCCGTCCTACGTACTGCGTGCTCTGGGCCGCAACGACGAGCTGGCGCACAGCTCGATCCGCTTCTCCTTCGGCCGTTTCACCACCGAAGAAGAAGTCGACTACGCCGCGCAGAAAGTCTGCGAGGCCGTCAATAAACTGCGTGAGCTGTCGCCGCTGTGGGATATGTACAAAGACGGCGTTGACATCTCCAAG  >  NC_002947/979451‑979744

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: