Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,426,755 A→T 6.5% intergenic (+52/+16) PP_4765 → / ← rhlE‑I hypothetical protein/ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,426,7550AT 6.5% 188.6 / 2.9 61intergenic (+52/+16)PP_4765/rhlE‑Ihypothetical protein/ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base T (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 = 9.91e-01

GGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTC  >  NC_002947/5426610‑5426903
                                                                                                                                                 |                                                                                                                                                    
ggAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCg                                                                                                                                                   <  8:296850/149‑1 (MQ=255)
 gAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCgg                                                                                                                                                  >  3:140862/1‑149 (MQ=255)
  aGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGc                                                                                                                                                 >  3:17567/1‑149 (MQ=255)
      gcCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCTTTGCGGccctt                                                                                                                                             >  5:267674/1‑147 (MQ=255)
      gcCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCg                                                                                                                                                   <  8:140730/143‑1 (MQ=255)
      gcCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCg                                                                                                                                                   >  7:140730/1‑143 (MQ=255)
        cGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCTTTGCGGccctttc                                                                                                                                           <  6:267674/149‑5 (MQ=255)
         gACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCaa                                                                                                                                           >  5:299481/1‑148 (MQ=255)
          aCCGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCaaaa                                                                                                                                         <  6:30575/149‑1 (MQ=255)
           ccGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGCAAAGGCCGCGTTGCGGcccgttagc                                                                                                                                         >  8:140347/1‑142 (MQ=255)
           ccGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCTTTGCGGccctttcg                                                                                                                                          <  7:140347/147‑6 (MQ=255)
               ggggTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCa                                                                                                                                    <  3:150549/149‑1 (MQ=255)
                 ggTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACg                                                                                                                                  >  2:239810/1‑149 (MQ=255)
                  gTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACgg                                                                                                                                 <  4:140862/149‑1 (MQ=255)
                   tATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGc                                                                                                                                <  6:299481/149‑1 (MQ=255)
                          tCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCt                                                                                                                         >  6:8364/1‑149 (MQ=255)
                            aCGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCtt                                                                                                                        <  5:8364/148‑1 (MQ=255)
                            aCGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTg                                                                                                                       >  8:92662/1‑149 (MQ=255)
                                  tGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGGGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGc                                                                                                                 <  1:35175/149‑1 (MQ=255)
                                   gCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGcc                                                                                                                >  5:115533/1‑149 (MQ=255)
                                    cTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGccc                                                                                                               <  1:239810/149‑1 (MQ=255)
                                          cAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCaaaa                                                                                                                                         <  4:7198/117‑1 (MQ=255)
                                          cAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCAtt                                                                                                         >  8:259906/1‑149 (MQ=255)
                                          cAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCAtt                                                                                                         <  4:17567/149‑1 (MQ=255)
                                           aGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTg                                                                                                        <  7:259906/149‑1 (MQ=255)
                                                   gcgcgcCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTc                                                                                                  <  7:92662/147‑1 (MQ=255)
                                                   gcgcgcCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCg                                                                                                 <  3:207274/148‑1 (MQ=255)
                                                            aCGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCg                                                                                       >  1:302976/1‑149 (MQ=255)
                                                                 cAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAg                                                                                  >  1:81725/1‑149 (MQ=255)
                                                                     ggcggATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCtt                                                                              >  1:86435/1‑149 (MQ=255)
                                                                     ggcggATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCtt                                                                              >  1:81470/1‑149 (MQ=255)
                                                                        ggATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCATGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCAc                                                                           <  3:146146/149‑1 (MQ=255)
                                                                              ggTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCgg                                                                     >  2:86744/1‑149 (MQ=255)
                                                                                  cAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTcc                                                                 <  6:115533/149‑1 (MQ=255)
                                                                                   aGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCg                                                                >  7:90469/1‑149 (MQ=255)
                                                                                          cTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGc                                                         <  2:86435/149‑1 (MQ=255)
                                                                                          cTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGc                                                         <  2:81470/149‑1 (MQ=255)
                                                                                           tGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCTCGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGcc                                                        >  6:291316/1‑149 (MQ=255)
                                                                                           tGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGcc                                                        >  3:239455/1‑149 (MQ=255)
                                                                                           tGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGcc                                                        >  1:288157/1‑149 (MQ=255)
                                                                                                 aTTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCa                                                  <  1:86744/149‑1 (MQ=255)
                                                                                                             aTGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAGGCTCGGCACCTTGCGTACAGCCTTCACCGCCgg                                                                     >  8:273086/1‑118 (MQ=255)
                                                                                                             aTGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAGGCTCGGCACCTTGCGTACAGCCTTCACCGCCgg                                                                     <  7:273086/118‑1 (MQ=255)
                                                                                                             aTGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGtccatc                                      >  6:228745/1‑149 (MQ=255)
                                                                                                               gCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGtccatcc                                     >  6:112310/1‑148 (MQ=255)
                                                                                                                     gCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGccc                              <  1:323702/149‑1 (MQ=255)
                                                                                                                         ccTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCa                            <  6:231762/147‑1 (MQ=255)
                                                                                                                         ccTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCa                            >  5:231762/1‑147 (MQ=255)
                                                                                                                         ccTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGc                          >  3:201803/1‑149 (MQ=255)
                                                                                                                         ccTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGc                          >  6:141942/1‑149 (MQ=255)
                                                                                                                         ccTAGCCTCGCGAAAGGGCCGCAAAGCGGCACCCAAAATCTCACCGCTTACGCTTCCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCAGCGCCGGTACCCCACTGGCCGCTTCAGTGCTGTCCATCCAACCCCCCAGc                          >  2:16784/1‑149 (MQ=255)
                                                                                                                          cTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGcc                         >  5:294070/1‑149 (MQ=255)
                                                                                                                             gCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGccc                        <  4:169819/147‑1 (MQ=255)
                                                                                                                               gTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCgcgc                    <  3:246890/149‑1 (MQ=255)
                                                                                                                                  gcgAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGcc                               >  8:313163/1‑135 (MQ=255)
                                                                                                                                  gcgAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGcc                               <  7:313163/135‑1 (MQ=255)
                                                                                                                                   cgAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTgg                <  8:90469/149‑1 (MQ=255)
                                                                                                                                    gAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCTTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGgc               <  4:201803/149‑1 (MQ=255)
                                                                                                                                     aaaGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGgcg              <  2:237647/149‑1 (MQ=255)
                                                                                                                                            cGCAAAGCGGCCCCCAAAATCTCAGGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGGACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCCGCCCGCGCTTGGCGGTGtttt       <  2:283334/149‑1 (MQ=255)
                                                                                                                                            cGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGCCACACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGCTTCCGCACTCGCCGCTTCAGTGCAGTCCAAGCAACGCCCCCGCCAGCGCTTGGCGCTGttt        >  1:283334/1‑148 (MQ=255)
                                                                                                                                               aaaGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTtttctt    <  2:81725/149‑1 (MQ=255)
                                                                                                                                                 aGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCCTTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTtttctttc  <  5:73259/149‑1 (MQ=255)
                                                                                                                                                 aGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTtttctttc  <  5:228745/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCAGCGGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCGGCCCCCAAAATCTCACGGCTTACGCTTGCGCGGCCCACCATTGAAGCTCGGCACCTTGCGTACAGCCTTCACCGCCGGTTCCGCACTGGCCGCTTCAGTGCTGTCCATCCAACGCCCCAGCCCGCGCTTGGCGCTGTTTTCTTTC  >  NC_002947/5426610‑5426903

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