Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,456,585 Δ1 bp 6.3% coding (400/2178 nt) PP_3950 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,456,5850C. 6.3% 205.1 / 7.1 63coding (400/2178 nt)PP_3950methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base . (2/2);  total (32/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 = 5.76e-01

CAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGTTGCGCAACATCGCGCTATCGAC  >  NC_002947/4456440‑4456726
                                                                                                                                                 |                                                                                                                                             
cAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTgcgc                                                                                                                                               >  3:161112/1‑146 (MQ=255)
 aaCCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCaga                                                                                                                                           <  8:106016/149‑1 (MQ=255)
  aCCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCagag                                                                                                                                          >  4:307784/1‑149 (MQ=255)
    ccAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCt                                                                                                                                        >  7:233951/1‑149 (MQ=255)
    ccAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCt                                                                                                                                        >  3:62324/1‑149 (MQ=255)
    ccAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCt                                                                                                                                        >  3:296276/1‑149 (MQ=255)
      aGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGc                                                                                                                                      >  5:309070/1‑149 (MQ=255)
         accagccccTCATGCACGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGACCCTGCGCCAGAGCTGCCa                                                                                                                                    >  8:271323/7‑148 (MQ=255)
           tttGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCaaa                                                                                                                                  >  8:270850/1‑148 (MQ=255)
              gCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATAc                                                                                                                               <  4:296276/148‑1 (MQ=255)
                     tGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACg                                                                                                                       >  8:235190/1‑149 (MQ=255)
                      gCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGt                                                                                                                      >  7:89705/1‑149 (MQ=255)
                           aGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCg                                                                                                                  <  8:233951/148‑1 (MQ=255)
                             gggCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGc                                                                                                               >  8:279314/1‑149 (MQ=255)
                                    tAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTg                                                                                                        <  4:158518/149‑1 (MQ=255)
                                    tAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCACGCCACgtat                                                                                                        <  6:235287/149‑4 (MQ=255)
                                    tAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCG‑CAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTgg                                                                                                       >  2:221283/1‑149 (MQ=255)
                                    tAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCG‑CAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTgg                                                                                                       >  2:215904/1‑149 (MQ=255)
                                      gAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCACGCCACgtatc                                                                                                       <  6:309070/148‑5 (MQ=255)
                                        aCATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTgggg                                                                                                     <  7:279314/148‑1 (MQ=255)
                                           tCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAg                                                                                                 <  2:234195/149‑1 (MQ=255)
                                              aCGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTc                                                                                              <  7:270850/149‑1 (MQ=255)
                                              aCGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTc                                                                                              <  1:299358/149‑1 (MQ=255)
                                              aCGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCCAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTc                                                                                              <  1:303553/149‑1 (MQ=255)
                                                    cTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGc                                                                                         <  8:89705/148‑1 (MQ=255)
                                                         tCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTggg                                                                                                      >  4:246051/1‑130 (MQ=255)
                                                         tCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTggg                                                                                                      <  3:246051/130‑1 (MQ=255)
                                                           gTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGc                                                                                                                                         >  1:176890/1‑93 (MQ=255)
                                                           gTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGc                                                                                                                                         <  2:176890/93‑1 (MQ=255)
                                                           gTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGc                                                                                 >  4:53964/1‑149 (MQ=255)
                                                              gAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGc                                                                              <  2:150606/149‑1 (MQ=255)
                                                               aaGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGcc                                                                             >  1:95355/1‑149 (MQ=255)
                                                               aaGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGcc                                                                             >  6:260218/1‑149 (MQ=255)
                                                                   gcgAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGc                                                                         <  8:137458/149‑1 (MQ=255)
                                                                        gAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATc                                                                    <  1:269523/149‑1 (MQ=255)
                                                                           gACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATcgac                                                                 <  7:102389/149‑1 (MQ=255)
                                                                                gCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCACTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACa                                                             >  4:279904/1‑148 (MQ=255)
                                                                                 ccGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACaaa                                                           >  1:120895/1‑149 (MQ=255)
                                                                                 ccGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACaaa                                                           <  5:260218/149‑1 (MQ=255)
                                                                                         cAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGcc                                                                             >  6:120799/1‑123 (MQ=255)
                                                                                            ggTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCG‑CAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCTGTGCACATCGACGACAAACAGTTCGGaaaa                                               <  1:215904/149‑1 (MQ=255)
                                                                                            ggTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCG‑CAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGaaaa                                               <  1:221283/149‑1 (MQ=255)
                                                                                               cAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTAAGCCCCTGCGCCAGAGATGCAAAATACTCAATACGTGGCGTACCACGATGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGt                                             >  7:55857/1‑149 (MQ=255)
                                                                                                  cGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGtcct                                          >  1:331735/1‑149 (MQ=255)
                                                                                                   gAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGtcctc                                         <  3:53964/149‑1 (MQ=255)
                                                                                                        gCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCACTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCt                                    <  3:279904/149‑1 (MQ=255)
                                                                                                        gCTGTAGGGCCCGTGCAGAAAAGGTGCTTTAAGCCCCTGCGCCAGAGCTGCCAAATAATCGATACGTGGCGTGACACGCTGGGGCCAGGTCGAAGCGGTGACGCTGACGTGCACATCGACGACAAAGAGTTCGGAAAAGTAATCCGCCt                                    >  4:271424/1‑149 (MQ=255)
                                                                                                            tAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGTGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAATAGTTCGGAAAAGTCCTCCGCCTggg                                 <  8:55857/148‑1 (MQ=255)
                                                                                                            tAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTgg                                  <  4:161112/147‑1 (MQ=255)
                                                                                                                gCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAGTACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGc                            >  8:349385/1‑149 (MQ=255)
                                                                                                                 cccGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATcga                                                                  <  2:331486/110‑1 (MQ=255)
                                                                                                                 cccGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATcga                                                                  >  1:331486/1‑110 (MQ=255)
                                                                                                                  ccGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGaa                                                                                           <  4:17744/84‑1 (MQ=255)
                                                                                                                     tGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGATCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGtt                       <  3:271424/149‑1 (MQ=255)
                                                                                                                     tGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAATTCCTCCGCCTGGGCCAGCGTGtt                       <  7:235190/149‑1 (MQ=255)
                                                                                                                          aaaCGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGTTgcgc                   >  7:42197/1‑148 (MQ=255)
                                                                                                                           aaCGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGTTGCGCaa                 >  5:91240/1‑149 (MQ=255)
                                                                                                                                tGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGTTGCGCAACATcg            >  4:221234/1‑149 (MQ=255)
                                                                                                                                tGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGTTGCGCAACATcg            <  7:271323/149‑1 (MQ=255)
                                                                                                                                   tttCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGTTGCGCAACATCgcgc         >  8:288070/1‑149 (MQ=255)
                                                                                                                                       aGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGTTGCGCAACATCGCGCTATc     >  4:362701/1‑149 (MQ=255)
                                                                                                                                          cccTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGTTGCGCAACATCGCGCTATCGAc  >  6:225782/1‑149 (MQ=255)
                                                                                                                                           ccTGCGCCAGAGCTGCCAAATACTCGATACNTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGTTGCGCAACATCGCGCt        >  1:33455/1‑142 (MQ=255)
                                                                                                                                           ccTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGTTGCGCAACATCGCGCt        <  2:33455/142‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CAACCCAGCACTTTGCCCTCGTGCATGAGGGGCTGGTAGAACATCAACGTGACTTCATCGTGGAAGCGCGAGGAAGACGGGCCGATCTGCAAGGTCAGCGAGTCGCTGTAGGGCCCGTGCAGAAACGGTGCTTTCAGCCCCTGCGCCAGAGCTGCCAAATACTCGATACGTGGCGTGCCACGCTGGGGCCAGGTCGAAGCGGTGACGCTGCCGTGCACATCGACGACAAACAGTTCGGAAAAGTCCTCCGCCTGGGCCAGCGTGTTGCGCAACATCGCGCTATCGAC  >  NC_002947/4456440‑4456726

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