Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,143,484 G→C 5.2% L88L (CTG→CTC)  PP_1899 → ExbD/TolR family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,143,4840GC 5.2% 257.2 / 3.1 77L88L (CTG→CTC) PP_1899ExbD/TolR family protein
Reads supporting (aligned to +/- strand):  ref base G (35/38);  new base C (2/2);  total (37/40)
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.65e-01

ACGCAACTTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCG  >  NC_002947/2143338‑2143629
                                                                                                                                                  |                                                                                                                                                  
aCGCAACTTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAg                                                                                                                                                  >  7:26224/1‑149 (MQ=255)
         cgcgTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGcc                                                                                                                                           <  1:191507/147‑1 (MQ=255)
          gcgTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCAt                                                                                                                                         <  5:109143/148‑1 (MQ=255)
          gcgTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATa                                                                                                                                        <  1:60436/149‑1 (MQ=255)
               gAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgc                                                                                                                                    <  5:26125/148‑1 (MQ=255)
                aGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                                  <  1:156854/149‑1 (MQ=255)
                  cTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                                  >  2:216625/1‑147 (MQ=255)
                  cTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGa                                                                                                                                <  3:197357/149‑1 (MQ=255)
                   tGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCTCGCTCAGGGTGGcc                                                                                                                                           >  8:147726/1‑137 (MQ=255)
                   tGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCTCGCTCAGGGTGGcc                                                                                                                                           <  7:147726/137‑1 (MQ=255)
                    gCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGt                                                                                                                              >  3:104836/1‑149 (MQ=255)
                     cccGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTc                                                                                                                             >  8:356946/1‑149 (MQ=255)
                      ccGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCg                                                                                                                            >  6:93633/1‑149 (MQ=255)
                        gAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCggg                                                                                                                          <  2:90227/149‑1 (MQ=255)
                          aGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGggcg                                                                                                                        <  4:242554/149‑1 (MQ=255)
                            ccaccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGggcggc                                                                                                                      <  7:6450/149‑1 (MQ=255)
                               ccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Gac                                                                                                                  <  6:381585/149‑1 (MQ=255)
                               ccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Ga                                                                                                                   >  6:328486/1‑148 (MQ=255)
                               ccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Ga                                                                                                                   >  6:334207/1‑148 (MQ=255)
                                caGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCggg                                                                                                                          >  5:161740/1‑141 (MQ=255)
                                caGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCggg                                                                                                                          <  6:161740/141‑1 (MQ=255)
                                  gcgcCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Gacaa                                                                                                                >  5:333221/1‑148 (MQ=255)
                                    gcCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGGGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Gacaac                                                                                                               >  1:336324/1‑147 (MQ=255)
                                       gAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGtt                                                                                                          <  8:26224/149‑1 (MQ=255)
                                         ctgtgcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GAcaacaa                                                                                                             >  4:206454/5‑144 (MQ=255)
                                            ggcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGc                                                                                                     >  2:315860/1‑149 (MQ=255)
                                            ggcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGc                                                                                                     >  4:130395/1‑149 (MQ=255)
                                            ggcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGc                                                                                                     <  4:104836/149‑1 (MQ=255)
                                             gcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCt                                                                                                    >  6:216614/1‑149 (MQ=255)
                                              cgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTg                                                                                                   <  1:315860/149‑1 (MQ=255)
                                              cgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTg                                                                                                   <  3:80277/149‑1 (MQ=255)
                                              cgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTg                                                                                                   <  5:93633/149‑1 (MQ=255)
                                               gccggccgACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCGTGAGCGAGGCCATTGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTgg                                                                                                  <  2:178944/149‑1 (MQ=255)
                                                    ccgACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATc                                                                                             <  1:216625/149‑1 (MQ=255)
                                                     cgACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCa                                                                                            <  6:398347/149‑1 (MQ=255)
                                                         ccGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCAGGCCATCAgcgc                                                                                        <  7:356946/149‑1 (MQ=255)
                                                            ggCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGa                                                                                     <  2:336324/149‑1 (MQ=255)
                                                                aGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGc                                                                                 <  3:130395/149‑1 (MQ=255)
                                                                aGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGc                                                                                 <  5:216614/149‑1 (MQ=255)
                                                                   tggtggAAATCACCTTCAGCGCCGCAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAg                                                                              <  7:76099/149‑1 (MQ=255)
                                                                       ggAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGAccc                                                                          >  1:357301/1‑149 (MQ=255)
                                                                        gAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                                  <  1:195909/93‑1 (MQ=255)
                                                                        gAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                                  >  2:195909/1‑93 (MQ=255)
                                                                          aaTCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGc                                                                       <  6:333221/149‑1 (MQ=255)
                                                                               ccATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAgcg                                                                                         >  8:219072/1‑126 (MQ=255)
                                                                               ccATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAgcg                                                                                         <  7:219072/126‑1 (MQ=255)
                                                                               ccATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAg                                                                  >  2:213117/1‑149 (MQ=255)
                                                                                cATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGc                                                                                                                                            >  5:228272/1‑75 (MQ=255)
                                                                                cATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGc                                                                                                                                            <  6:228272/75‑1 (MQ=255)
                                                                                cATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAgg                                                                 >  2:334457/1‑149 (MQ=255)
                                                                                 aTCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCAGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAgg                                                                 <  1:101671/149‑1 (MQ=255)
                                                                                 aTCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGc                                                                >  4:186232/1‑149 (MQ=255)
                                                                                        ccGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGt                                                          >  2:315051/1‑148 (MQ=255)
                                                                                        ccGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTc                                                         <  1:213117/149‑1 (MQ=255)
                                                                                                gtgtATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAAt                                                 >  6:162473/1‑149 (MQ=255)
                                                                                                         gTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCAGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGAtgctgc                                        >  2:129030/1‑149 (MQ=255)
                                                                                                           gAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCg                                      >  6:194844/1‑149 (MQ=255)
                                                                                                           gAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCg                                      >  6:187935/1‑149 (MQ=255)
                                                                                                               aacCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGc                                    >  7:162558/1‑147 (MQ=255)
                                                                                                                     ctgctgCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCgg                            >  5:15706/1‑149 (MQ=255)
                                                                                                                        ctgcCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGttt                         >  7:316022/1‑149 (MQ=255)
                                                                                                                           cccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                       <  6:15706/148‑1 (MQ=255)
                                                                                                                           cccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                       <  5:194844/148‑1 (MQ=255)
                                                                                                                            ccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGc                                                                                                 >  4:97459/1‑73 (MQ=255)
                                                                                                                            ccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGc                                                                                                 <  3:97459/73‑1 (MQ=255)
                                                                                                                            ccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGc                     >  7:34785/1‑149 (MQ=255)
                                                                                                                             cAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGa                                              <  7:351045/123‑1 (MQ=255)
                                                                                                                             cAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGa                                              >  8:351045/1‑123 (MQ=255)
                                                                                                                             cAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGc                     >  5:93459/1‑148 (MQ=255)
                                                                                                                              aaaaGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCa                   <  6:93459/149‑1 (MQ=255)
                                                                                                                                    tctatGGCCTCGCTCAGGGTGGCCATAGAGCGCGAGTCGGGCGGC‑GGCAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                       >  5:114946/5‑139 (MQ=255)
                                                                                                                                    tctatGGCCTCGCTCAGGGTGGCCATAGAGCGCGAGTCGGGCGGC‑GGCAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                       <  6:114946/135‑1 (MQ=255)
                                                                                                                                     aCCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTgcgc            <  5:187935/149‑1 (MQ=255)
                                                                                                                                     aCCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTgcgc            <  1:334457/149‑1 (MQ=255)
                                                                                                                                          gCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGa        <  1:315051/148‑1 (MQ=255)
                                                                                                                                               ccTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCg  >  3:17842/1‑149 (MQ=255)
                                                                                                                                               ccTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCg  >  1:383841/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACGCAACTTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCG  >  NC_002947/2143338‑2143629

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑