Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,360,293:1 +G 10.1% coding (829/1161 nt) PP_3835 ← two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,360,2931.G10.1% 249.6 / 20.8 69coding (829/1161 nt)PP_3835two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base . (29/33);  new base G (4/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

AAGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGCGTTAAC  >  NC_002947/4360145‑4360437
                                                                                                                                                     |                                                                                                                                                
aaGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATc                                                                                                                                                   <  5:84014/149‑1 (MQ=255)
       gcgGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCa                                                                                                                                           >  8:314986/1‑149 (MQ=255)
       gcgGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCa                                                                                                                                           >  2:52263/1‑149 (MQ=255)
          gCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACGGGCAtgg                                                                                                                                        <  3:194719/149‑3 (MQ=255)
            aGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGaa                                                                                                                                      <  1:52263/149‑1 (MQ=255)
            aGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACGGGCAtggcc                                                                                                                                      <  3:301824/149‑5 (MQ=255)
                    tgtgGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAacaaca                                                                                                                                 <  7:186544/146‑1 (MQ=255)
                    tgtgGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGt                                                                                                                              >  3:352325/1‑149 (MQ=255)
                        ggTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCTCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGa                                                                                                                           <  5:155516/148‑1 (MQ=255)
                         gTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGc                                                                                                                         <  4:46150/149‑1 (MQ=255)
                           aTGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCaa                                                                                                                       >  3:145817/1‑149 (MQ=255)
                            tGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCaa                                                                                                                       <  5:162326/148‑1 (MQ=255)
                            tGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAac                                                                                                                      <  4:352325/149‑1 (MQ=255)
                               cTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGa                                                                                                                           >  8:288927/1‑141 (MQ=255)
                               cTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGa                                                                                                                           <  7:288927/141‑1 (MQ=255)
                               cTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAacac                                                                                                                    <  3:171504/148‑1 (MQ=255)
                               cTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAacac                                                                                                                    >  4:171504/1‑148 (MQ=255)
                                ttCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCa                                                                                                                  >  7:98642/1‑149 (MQ=255)
                                        tGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGt                                                                                                          >  4:61080/1‑149 (MQ=255)
                                               cACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCg                                                                                                    <  7:314986/148‑1 (MQ=255)
                                               cACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCg                                                                                                    <  4:289874/148‑1 (MQ=255)
                                               cACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCg                                                                                                    <  4:123277/148‑1 (MQ=255)
                                                aCCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCATCACCATCTGCAGTTGCCCgg                                                                                                   >  5:285630/1‑148 (MQ=255)
                                                aCCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCgg                                                                                                   >  2:330808/1‑148 (MQ=255)
                                                aCCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCgg                                                                                                   >  5:290618/1‑148 (MQ=255)
                                                 ccTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCt                                                                                                               >  7:96140/1‑135 (MQ=255)
                                                 ccTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCt                                                                                                               <  8:96140/135‑1 (MQ=255)
                                                 ccTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCa                                                                                                 <  1:330808/149‑1 (MQ=255)
                                                    cGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATc                                                                                               >  1:354727/1‑148 (MQ=255)
                                                     gAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATcc                                                                                              <  1:295034/148‑1 (MQ=255)
                                                        gcgcACGAAGCTGCCCTTCAGCCCTTCATCTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCt                                                                                          <  8:134229/149‑1 (MQ=255)
                                                           cACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGa                                                                                       >  6:360337/1‑149 (MQ=255)
                                                            aCGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGAt                                                                                      <  8:98642/149‑1 (MQ=255)
                                                            aCGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGAt                                                                                      <  4:138957/149‑1 (MQ=255)
                                                                   tGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGACTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGccc                                                                                                     <  4:234413/127‑1 (MQ=255)
                                                                   tGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGACTGATGCGGCCATGCCCGACGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGccc                                                                                                     >  3:234413/1‑127 (MQ=255)
                                                                     cccttcagcccttcaACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAg                                                                             >  4:122813/1‑149 (MQ=255)
                                                                     cccttcagcccttcaACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAg                                                                             >  7:43/1‑149 (MQ=255)
                                                                      ccttcagcccttcaACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAgg                                                                            <  5:360337/149‑1 (MQ=255)
                                                                             cccttcaACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAg                                                                                                                          >  4:242359/1‑96 (MQ=255)
                                                                             cccttcaACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAg                                                                                                                          <  3:242359/96‑1 (MQ=255)
                                                                                  caACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTgcag                                                                 >  6:49604/1‑148 (MQ=255)
                                                                                      ttCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTc                                                                                         <  6:95668/120‑1 (MQ=255)
                                                                                      ttCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTc                                                                                         >  5:95668/1‑120 (MQ=255)
                                                                                      ttCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAaca                                                            >  1:65570/1‑149 (MQ=255)
                                                                                          ccGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCtt                                                        <  4:145817/149‑1 (MQ=255)
                                                                                           cGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAATAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCa                                                                                                 <  2:350796/107‑1 (MQ=255)
                                                                                           cGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAATAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCa                                                                                                 >  1:350796/1‑107 (MQ=255)
                                                                                            gAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTc                                                       >  2:192896/1‑148 (MQ=255)
                                                                                              acaacTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAg                                                    <  2:354727/149‑1 (MQ=255)
                                                                                               caacTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGt                                                   <  3:61080/149‑1 (MQ=255)
                                                                                                  cTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTgca                                                                  >  8:282413/1‑131 (MQ=255)
                                                                                                  cTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTgca                                                                  <  7:282413/131‑1 (MQ=255)
                                                                                                     gACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAAt                                               <  3:122813/147‑1 (MQ=255)
                                                                                                     gACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCg                                             <  5:49604/149‑1 (MQ=255)
                                                                                                          caTCGACAGTCAGGCTGATGCTGCCTTGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCGACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCtt                                        <  5:325757/149‑1 (MQ=255)
                                                                                                           aTCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAg                                                                                                                          <  2:300702/66‑1 (MQ=255)
                                                                                                           aTCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAg                                                                                                                          >  1:300702/1‑66 (MQ=255)
                                                                                                                 aGTCAGGCTGATGCGCCCATGCCCGTCGATGGCATCACGGGCAATGAACCACAGCTTGAGCAACCCCATCTGCAGTTGACCGGCAACAACGGCGATCAACGGCAGGGCCGGTTGCAGCAACACCTACAAGACAATCGAGCTTGGCAAtc                                 >  6:267378/1‑148 (MQ=255)
                                                                                                                   tCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATCCa                               >  2:303462/1‑149 (MQ=255)
                                                                                                                              gttcaaTGCCCGTGATGCCATCGACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAg                    >  2:269554/6‑149 (MQ=255)
                                                                                                                                 caaTGCCCGTGATGCCATCGACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTgg                                      <  2:164503/126‑1 (MQ=255)
                                                                                                                                 caaTGCCCGTGATGCCATCGACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTgg                                      >  1:164503/3‑128 (MQ=255)
                                                                                                                                 caaTGCCCGTGATGACATCGACGGGCATTGAACAAAAGGTTGAGCAACACCATCTTCAGTTGCCAGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTgg                                      >  1:169883/3‑128 (MQ=255)
                                                                                                                                  aaTGCCCGTGATGCCATCGACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCa                >  3:6794/2‑149 (MQ=255)
                                                                                                                                       ccGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCAcgcg            >  5:135104/1‑148 (MQ=255)
                                                                                                                                            gATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGCg       <  7:111100/148‑1 (MQ=255)
                                                                                                                                                gCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGCGTTAAc  <  1:192896/149‑1 (MQ=255)
                                                                                                                                                 cATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAAc                                                                                  <  1:71460/68‑1 (MQ=255)
                                                                                                                                                 cATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAAc                                                                                  >  2:71460/1‑68 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AAGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGCGTTAAC  >  NC_002947/4360145‑4360437

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: