Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,360,762 G→A 8.7% G120G (GGC→GGT) ‡ 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,7620GA 8.7% 158.4 / 4.5 58G120G (GGC→GGT) ‡PP_3835two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base A (3/2);  total (29/29)
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.13e-01

TCAACATGTTGTTGAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGGGCGGCCATCCG  >  NC_002947/4360614‑4360910
                                                                                                                                                    |                                                                                                                                                    
tCAACATGTTGTTGAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTggccagg                                                                                                                                                      >  2:8071/1‑149 (MQ=255)
 cAACATGTTGTTGAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTggccaggc                                                                                                                                                     >  5:214520/1‑149 (MQ=255)
    cATGTTGTTGAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCaa                                                                                                                                                  >  2:2810/1‑149 (MQ=255)
          gttgAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCt                                                                                                                                            >  7:322194/1‑149 (MQ=255)
             gAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGc                                                                                                                                         >  5:50709/1‑149 (MQ=255)
                 tCGTGCGCCAGACGACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATTTCACGTGTCACCTTGGCCAGATCATTGGCCTGGCCaag                                                                                                                                     <  1:93478/149‑3 (MQ=255)
                   gTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACttt                                                                                                                                   >  1:189155/1‑149 (MQ=255)
                   gTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACttt                                                                                                                                   >  1:183263/1‑149 (MQ=255)
                    tGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCTCTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGATCATTGGCCTGGCCaag                                                                                                                                     <  6:334895/146‑3 (MQ=255)
                    tGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGATCATTGGCCTGGCCaag                                                                                                                                     >  5:334895/1‑144 (MQ=255)
                     gcgcCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTtca                                                                                                                                 >  2:127327/1‑149 (MQ=255)
                       gcCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTtcatc                                                                                                                               <  6:165871/149‑1 (MQ=255)
                         cAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCAc                                                                                                                             <  6:50709/149‑1 (MQ=255)
                              caccGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTgg                                                                                                                                          <  5:367495/131‑1 (MQ=255)
                              caccGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTgg                                                                                                                                          >  6:367495/1‑131 (MQ=255)
                              caccGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGAt                                                                                                                         >  6:61419/1‑148 (MQ=255)
                               accGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATca                                                                                                                       >  5:84643/1‑149 (MQ=255)
                                      aaCTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCa                                                                                                                 <  1:127327/148‑1 (MQ=255)
                                      aaCTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAg                                                                                                                >  2:318671/1‑149 (MQ=255)
                                             cGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTgccgc                                                                                                          >  1:334207/1‑148 (MQ=255)
                                                  gCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGc                                                                                                    >  7:180372/1‑149 (MQ=255)
                                                             cTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGc                                                                                          <  5:61419/148‑1 (MQ=255)
                                                             cTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGc                                                                                          <  2:334207/148‑1 (MQ=255)
                                                              tGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACtt                                                                                                                                    <  3:30362/105‑1 (MQ=255)
                                                              tGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACtt                                                                                                                                    >  4:30362/1‑105 (MQ=255)
                                                                 gcCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGcc                                                                                     >  3:102096/1‑149 (MQ=255)
                                                                 gcCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGcc                                                                                     >  2:33043/1‑149 (MQ=255)
                                                                   cTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGt                                                                                   >  8:281966/1‑149 (MQ=255)
                                                                        aCAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTg                                                                              <  1:8071/149‑1 (MQ=255)
                                                                        aCAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTg                                                                              <  1:318671/149‑1 (MQ=255)
                                                                         cAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCTCATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTgc                                                                             <  6:84643/149‑1 (MQ=255)
                                                                          aGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTgcg                                                                            >  7:186279/1‑149 (MQ=255)
                                                                          aGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTgcg                                                                            <  1:2810/149‑1 (MQ=255)
                                                                                                gctgctGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTgc                                                      <  2:189155/149‑1 (MQ=255)
                                                                                                gctgctGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTgc                                                      <  2:183263/149‑1 (MQ=255)
                                                                                                    ctgcGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTGCTTGCGCACCCGCCAGCCTTCGTCCTgcgc                                                    <  1:33043/147‑1 (MQ=255)
                                                                                                    ctgcGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTgcgt                                                    <  7:281966/147‑2 (MQ=255)
                                                                                                    ctgcGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTgcgc                                                    <  4:102096/147‑1 (MQ=255)
                                                                                                      tgaccGTCGCGAGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCcaca                                                >  1:177646/4‑149 (MQ=255)
                                                                                                       cgcCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCacac                                               <  6:214520/149‑1 (MQ=255)
                                                                                                         ccGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCt                                             >  6:190039/1‑149 (MQ=255)
                                                                                                            tctcGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGc                                          <  3:20595/149‑1 (MQ=255)
                                                                                                               cgcgACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTc                                        >  2:222989/1‑148 (MQ=255)
                                                                                                                 cgACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAg                                     >  6:16630/1‑149 (MQ=255)
                                                                                                                  gACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAgc                                    <  1:222989/149‑1 (MQ=255)
                                                                                                                      gTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCt                                <  8:322194/149‑1 (MQ=255)
                                                                                                                         aGTGATATCACGTGTCACCTTGGCCAGTCCAATGATCTCTCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTgg                              <  8:180372/148‑1 (MQ=255)
                                                                                                                           tGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGc                                          >  6:163017/1‑134 (MQ=255)
                                                                                                                           tGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGc                                          <  5:163017/134‑1 (MQ=255)
                                                                                                                             atatCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGTGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCAcag                         >  3:235161/1‑149 (MQ=255)
                                                                                                                              tatCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCa                                                                   <  2:87320/106‑1 (MQ=255)
                                                                                                                              tatCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCa                                                                   >  1:87320/1‑106 (MQ=255)
                                                                                                                                  aCGTGTCACCTTGGCCAGGCCAATGATCTGTCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGTCACAGCAGAc                    <  2:23456/149‑1 (MQ=255)
                                                                                                                                  aCGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGAc                    >  2:117332/1‑149 (MQ=255)
                                                                                                                                  aCGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGAc                    <  8:14202/149‑1 (MQ=255)
                                                                                                                                  aCGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGAc                    >  4:5632/1‑149 (MQ=255)
                                                                                                                                      gtCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTc                <  8:186279/149‑1 (MQ=255)
                                                                                                                                      atgaaagtGGCCAGATCATTGGCCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCt                 >  8:270535/8‑148 (MQ=255)
                                                                                                                                       tCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCa               <  6:26530/149‑1 (MQ=255)
                                                                                                                                        gaaagtGGCCAGATCATTGGCCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAg              >  2:352738/6‑149 (MQ=255)
                                                                                                                                                    gccaATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGGGCGGCCATCCg  >  1:96918/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCAACATGTTGTTGAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGGGCGGCCATCCG  >  NC_002947/4360614‑4360910

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