Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,135,888 G→C 6.6% L228V (CTC→GTC)  PP_4521 ← aerotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,135,8880GC 6.6% 198.7 / 3.9 61L228V (CTC→GTC) PP_4521aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base C (2/2);  total (30/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 = 4.36e-01

GCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGA  >  NC_002947/5135750‑5136024
                                                                                                                                          |                                                                                                                                        
gcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCGCtt                                                                                                                                 >  8:118440/1‑148 (MQ=255)
 ccGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTc                                                                                                                                <  4:116415/148‑1 (MQ=255)
 ccGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCa                                                                                                                               <  5:243503/149‑1 (MQ=255)
       gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGc                                                                                                                                         >  1:125436/1‑133 (MQ=255)
       gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGc                                                                                                                                         <  2:125436/133‑1 (MQ=255)
       gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCt                                                                                                                          <  5:134077/148‑1 (MQ=255)
        aGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTgg                                                                                                                        <  8:123762/149‑1 (MQ=255)
         gCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCTACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGATGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTtgg                                                                                                                       <  6:114644/149‑1 (MQ=255)
           aTCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTgg                                                                                                                        >  1:192862/1‑146 (MQ=255)
           aTCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTgg                                                                                                                        <  2:192862/146‑1 (MQ=255)
             cTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCt                                                                                                                                        <  4:243923/128‑1 (MQ=255)
             cTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCt                                                                                                                                        >  3:243923/1‑128 (MQ=255)
                   cTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTc                                                                                                                  >  3:134195/1‑144 (MQ=255)
                   cTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTc                                                                                                                  <  4:134195/144‑1 (MQ=255)
                   cTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTg                                                                                                              >  5:241004/1‑148 (MQ=255)
                   cTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGt                                                                                                             <  5:146544/149‑1 (MQ=255)
                     gggCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTc                                                                                                                                <  5:18770/128‑1 (MQ=255)
                     gggCATCCGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGACGAGCTGGCCCTTc                                                                                                                                >  6:18770/1‑128 (MQ=255)
                                aactCGGAAGATTTCGCCGCCCGCGACTTTACCCCGACCCCCCCATCCACCACGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTCTAGGTCAGCGccac                                                                                                <  6:51029/146‑1 (MQ=255)
                                aaGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGccac                                                                                                <  7:242418/149‑1 (MQ=255)
                                          aTTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTg                                                                                      <  3:249986/149‑1 (MQ=255)
                                               gccgccCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTa                                                                                   >  1:178931/1‑147 (MQ=255)
                                                ccgcccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAgcgc                                                                                                   >  3:179050/1‑130 (MQ=255)
                                                ccgcccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAgcgc                                                                                                   <  4:179050/130‑1 (MQ=255)
                                                ccgcccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCACGTCGAGCTGGCCCTTCAGGCCTGGGTTCTCGCTGTAGGTCAGCGCCACCAACGGAATGGTAAAg                                                                                >  6:134893/1‑149 (MQ=255)
                                                 cgcccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGGCAGCGCCACCAACGGACTGGTAAAgg                                                                               >  1:133762/1‑149 (MQ=255)
                                                     cgcgACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCtt                                                                           <  2:178931/149‑1 (MQ=255)
                                                        gACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTg                                                                          >  8:178008/1‑147 (MQ=255)
                                                        gACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTg                                                                          <  7:178008/147‑1 (MQ=255)
                                                           tttACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTggg                                                                                                                       <  3:103215/99‑1 (MQ=255)
                                                           tttACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTggg                                                                                                                       >  4:103215/1‑99 (MQ=255)
                                                           tttACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATg                                                                     >  3:71057/1‑149 (MQ=255)
                                                           tttACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGccatc                                                                                                                                 >  1:113849/1‑86 (MQ=255)
                                                           tttACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGccatc                                                                                                                                 <  2:113849/89‑4 (MQ=255)
                                                                cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGt                                                                                    <  7:181709/129‑1 (MQ=255)
                                                                cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGt                                                                                    >  8:181709/1‑129 (MQ=255)
                                                                 ccgcccccgGCATGCACCAGGCGGGTCAGGGACGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGGCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGc                                                               <  1:180864/149‑1 (MQ=255)
                                                                      cccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACg                                                          >  7:187183/1‑149 (MQ=255)
                                                                         ggCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGc                                                       <  3:248261/149‑1 (MQ=255)
                                                                           cATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGCCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGcga                                                     <  5:185940/149‑1 (MQ=255)
                                                                                accaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAAt                                                <  8:187183/149‑1 (MQ=255)
                                                                                accaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAAt                                                >  3:131041/1‑149 (MQ=255)
                                                                                 ccaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATg                                               <  4:131041/149‑1 (MQ=255)
                                                                                    ggCGGGGCAGCGCCGTTTGCAGCCGGGCTTCTTCGCTGATGATGGCCAGGTCGCGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCCGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGcc                                            <  4:71057/149‑1 (MQ=255)
                                                                                       gggTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTgg                                         >  1:107522/1‑149 (MQ=255)
                                                                                          tCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTc                                      >  8:163624/1‑149 (MQ=255)
                                                                                                              gCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGGCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTa                  <  5:50915/149‑1 (MQ=255)
                                                                                                                  cttcGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTAcc                <  3:168612/147‑1 (MQ=255)
                                                                                                                  cttcGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTAcc                >  3:121033/1‑147 (MQ=255)
                                                                                                                  cttcGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTCCCGTTGATCGACCCAGTACCaa              >  6:137426/1‑149 (MQ=255)
                                                                                                                   ttcGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAg             >  6:271254/1‑149 (MQ=255)
                                                                                                                    tcGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGc            >  6:100984/1‑149 (MQ=255)
                                                                                                                      gCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTaaa                                                                                 >  4:193775/1‑78 (MQ=255)
                                                                                                                      gCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTaaa                                                                                 <  3:193775/78‑1 (MQ=255)
                                                                                                                        tgatgaTGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcag         >  4:29460/1‑148 (MQ=255)
                                                                                                                        tgatgaTGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcag         >  4:26617/1‑148 (MQ=255)
                                                                                                                        tgatgaTGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcag         >  5:26553/1‑148 (MQ=255)
                                                                                                                           tgatgGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcagc     <  6:241004/149‑1 (MQ=255)
                                                                                                                            aagggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGc            <  2:203693/138‑1 (MQ=255)
                                                                                                                            aagggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGc            >  1:203693/4‑141 (MQ=255)
                                                                                                                              tgGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGa  >  1:6172/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                        
GCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGA  >  NC_002947/5135750‑5136024

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: ‑