Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,287,468 C→G 6.1% G90R (GGC→CGC)  PP_3758 ← two component system response regulator/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,287,4680CG 6.1% 210.5 / 3.8 65G90R (GGC→CGC) PP_3758two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base G (2/2);  total (33/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 = 9.12e-01

GCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCCAGCAGCAACGACAAGGCAGC  >  NC_002947/4287327‑4287615
                                                                                                                                             |                                                                                                                                                   
gCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCcgcg                                                                                                                                              >  3:77267/1‑149 (MQ=255)
    gCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGa                                                                                                                                          <  4:77267/149‑1 (MQ=255)
      cAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGAc                                                                                                                                         >  4:101310/1‑148 (MQ=255)
       aaGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACg                                                                                                                                        >  7:27224/1‑148 (MQ=255)
       aaGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACg                                                                                                                                        >  4:179330/1‑148 (MQ=255)
        aGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGaa                                                                                                                                      <  3:101310/149‑1 (MQ=255)
         gCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAg                                                                                                                                     <  3:179330/149‑1 (MQ=255)
            ccTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACg                                                                                                                                  <  6:186927/149‑1 (MQ=255)
                gcgcTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCg                                                                                                                              >  8:60001/1‑149 (MQ=255)
                 cgcTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCgg                                                                                                                             >  8:1175/1‑149 (MQ=255)
                      gCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGAtgt                                                                                                                        <  4:130575/149‑1 (MQ=255)
                      gCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGAt                                                                                                                          <  1:248346/147‑1 (MQ=255)
                        ggAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCGGGCcgcgagatg                                                                                                                                         <  7:241593/130‑6 (MQ=255)
                        ggAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCGGGCcgcgagatg                                                                                                                                         >  8:241593/1‑125 (MQ=255)
                        ggAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATc                                                                                                                               <  7:99420/140‑1 (MQ=255)
                        ggAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATc                                                                                                                               >  8:99420/1‑140 (MQ=255)
                        ggAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATc                                                                                                                               >  8:95809/1‑140 (MQ=255)
                        ggAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATc                                                                                                                               <  7:95809/140‑1 (MQ=255)
                         gAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCt                                                                                                                     >  3:220019/1‑149 (MQ=255)
                         gAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCt                                                                                                                     <  8:27224/149‑1 (MQ=255)
                                  gacgaACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGTCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGt                                                                                                            >  3:204904/1‑149 (MQ=255)
                                  gacgaACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGt                                                                                                            <  7:223721/149‑1 (MQ=255)
                                            gagaCCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCGggccgc                                                                                                                                               >  5:20154/1‑104 (MQ=255)
                                            gagaCCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCGggccgc                                                                                                                                               <  6:20154/104‑1 (MQ=255)
                                                         ttACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGc                                                                                      <  5:151619/148‑1 (MQ=255)
                                                         ttACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGc                                                                                      <  7:1175/148‑1 (MQ=255)
                                                          tACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGc                                                                                      >  4:3072/1‑147 (MQ=255)
                                                          tACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGc                                                                                      <  3:3072/147‑1 (MQ=255)
                                                           aCCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAg                                                                                   >  3:250449/1‑149 (MQ=255)
                                                                cAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGtt                                                                                 <  8:229657/146‑1 (MQ=255)
                                                                cAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGtt                                                                                 >  7:229657/1‑146 (MQ=255)
                                                                      gtCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGt                                                                        >  5:123002/1‑149 (MQ=255)
                                                                         gAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTc                                                                      <  7:60001/148‑1 (MQ=255)
                                                                         gAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCa                                                                     >  2:129457/1‑149 (MQ=255)
                                                                         gAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCa                                                                     >  3:73210/1‑149 (MQ=255)
                                                                         gAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCa                                                                     >  6:255257/1‑149 (MQ=255)
                                                                                ttGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTg                                                                                                             <  1:81716/102‑1 (MQ=255)
                                                                                ttGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTg                                                                                                             >  2:81716/1‑102 (MQ=255)
                                                                                ttGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCgg                                                              >  8:218864/1‑149 (MQ=255)
                                                                                ttGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCgg                                                              >  3:231755/1‑149 (MQ=255)
                                                                                  gtgCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCa                                                            >  3:153594/1‑149 (MQ=255)
                                                                                   tgCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCgg                                                              >  4:81664/1‑146 (MQ=255)
                                                                                   tgCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCgg                                                              <  3:81664/146‑1 (MQ=255)
                                                                                   tgCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCa                                                            <  5:255257/148‑1 (MQ=255)
                                                                                     cAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATc                                                          <  5:136244/148‑1 (MQ=255)
                                                                                      aGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTg                                                        >  2:43652/1‑149 (MQ=255)
                                                                                        aaaGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGca                                                      <  4:204904/149‑1 (MQ=255)
                                                                                          aGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGcaca                                                    >  6:252568/1‑149 (MQ=255)
                                                                                          aGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGcaca                                                    <  1:134751/149‑1 (MQ=255)
                                                                                                 ggCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGaa                                             <  7:218864/149‑1 (MQ=255)
                                                                                                     ccGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGc                                         <  6:123002/149‑1 (MQ=255)
                                                                                                     ccGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGc                                         <  1:43652/149‑1 (MQ=255)
                                                                                                        ctctCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAg                                      >  1:100947/1‑149 (MQ=255)
                                                                                                             gTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTg                                    <  4:153594/146‑1 (MQ=255)
                                                                                                             gTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTg                                    <  2:100947/146‑1 (MQ=255)
                                                                                                                 cccTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTc                             <  4:73210/149‑1 (MQ=255)
                                                                                                                  ccTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCg                            >  3:262386/1‑149 (MQ=255)
                                                                                                                    ttGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCGAGTGCGAACTCGTg                          >  8:99855/1‑149 (MQ=255)
                                                                                                                            aTAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCcagc                   <  4:161472/148‑1 (MQ=255)
                                                                                                                            aTAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCcagc                   >  3:161472/1‑148 (MQ=255)
                                                                                                                               aTTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCcag                    <  6:233601/144‑1 (MQ=255)
                                                                                                                               aTTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCcag                    >  5:233601/1‑144 (MQ=255)
                                                                                                                                          ggCCCGCCGCGCTGACTATGACGATCTGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGATAATGGCCAGTGCGAACTCTTTCTCCAGCAGCAACGACAAGGCa    <  4:162688/149‑1 (MQ=255)
                                                                                                                                          ggCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCCAGCAGCAACGACAAgg      >  3:81747/1‑147 (MQ=255)
                                                                                                                                          cgCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGACTTGTCGGTACCGCGCATCAGCTCGGCCAGTTAGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCCAGCAGCAACGACAAgg      >  1:82360/2‑147 (MQ=255)
                                                                                                                                             ccgccgCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCCAGCAGCAACGACAAGGCAGc  >  4:9067/1‑148 (MQ=255)
                                                                                                                                             ccgccgCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCCAGCAGCAACGACAAGGCAGc  >  4:12761/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCCAGCAGCAACGACAAGGCAGC  >  NC_002947/4287327‑4287615

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