Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,945,627 C→G 5.6% P318P (CCG→CCC flhB ← flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,945,6270CG 5.6% 285.4 / 4.7 89P318P (CCG→CCCflhBflagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (41/43);  new base G (3/2);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCGGGTTG  >  NC_002947/4945483‑4945774
                                                                                                                                                 |                                                                                                                                                   
cAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCgg                                                                                                                                                   <  1:369650/147‑1 (MQ=255)
  gTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGAt                                                                                                                                               >  5:363635/1‑149 (MQ=255)
     gCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCGGGCGAtcta                                                                                                                                            >  7:172273/1‑146 (MQ=255)
         tGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAg                                                                                                                                          <  7:73616/147‑1 (MQ=255)
          gggCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGa                                                                                                                                        >  1:310058/1‑149 (MQ=255)
           ggCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCGGGCGAtctatta                                                                                                                                         <  8:172273/146‑7 (MQ=255)
           ggCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCAGCGCCAGCGTCGGCGATTCGAGGATc                                                                                                                                      >  5:316744/1‑149 (MQ=255)
                  gcccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAt                                                                                                                                 <  6:303400/147‑1 (MQ=255)
                   cccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgtt                                                                                                                              <  6:321970/149‑1 (MQ=255)
                   cccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgtt                                                                                                                              <  1:78571/149‑1 (MQ=255)
                     cttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgttgt                                                                                                                            <  6:134770/149‑1 (MQ=255)
                        gccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGc                                                                                                                          <  7:259813/148‑1 (MQ=255)
                        gccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCt                                                                                                                         >  6:133283/1‑149 (MQ=255)
                              gcgGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGc                                                                                                                   >  2:50920/1‑149 (MQ=255)
                               cgGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGcc                                                                                                                  <  5:196623/149‑1 (MQ=255)
                               cgGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGcc                                                                                                                  <  2:30873/149‑1 (MQ=255)
                                ggTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGcc                                                                                                                  >  4:25078/1‑148 (MQ=255)
                                    ttGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAAttt                                                                                                             >  6:242248/1‑149 (MQ=255)
                                    ttGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAAtt                                                                                                              >  4:221345/1‑148 (MQ=255)
                                    ttGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAAtt                                                                                                              <  5:133283/148‑1 (MQ=255)
                                      ggCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCa                                                                                                           >  3:9563/1‑149 (MQ=255)
                                       gCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCa                                                                                                           <  6:2092/148‑1 (MQ=255)
                                          aaTCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGa                                                                                                                                  >  2:350912/1‑122 (MQ=255)
                                          aaTCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGa                                                                                                                                  <  1:350912/122‑1 (MQ=255)
                                          aaTCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGaa                                                                                                       >  7:32732/1‑149 (MQ=255)
                                           aTCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAAt                                                                                                      <  6:363635/149‑1 (MQ=255)
                                               ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTgg                                                                                                                                   >  2:192297/1‑116 (MQ=255)
                                               ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTgg                                                                                                                                   <  1:192297/116‑1 (MQ=255)
                                               ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCtt                                                                                                   >  4:136864/1‑148 (MQ=255)
                                                gAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCAATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTc                                                                                                  >  5:142384/1‑148 (MQ=255)
                                                  acacATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAAc                                                                                               <  5:218815/149‑1 (MQ=255)
                                                    acaTACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGc                                                                                             <  3:221345/149‑1 (MQ=255)
                                                    acaTACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGc                                                                                             >  4:403283/1‑149 (MQ=255)
                                                      aTACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCCCGCGAATTTCACGGGGCTGCAACGCCa                                                                                           >  2:244228/1‑149 (MQ=255)
                                                           ccAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGaaa                                                                                      >  5:272726/1‑149 (MQ=255)
                                                                cACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGGTCATGTTCGGTGGAGTAATAGATAGCCCGCTCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCg                                                                                   <  4:325019/148‑1 (MQ=255)
                                                                  ccTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTa                                                                                >  2:57976/1‑149 (MQ=255)
                                                                   cTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGGATAACGCCAATTTCACGAATATTCAACGCCATGAAATCGCTAc                                                                               >  1:134964/1‑149 (MQ=255)
                                                                   cTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTAc                                                                               <  1:262825/149‑1 (MQ=255)
                                                                   cTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTAc                                                                               <  1:50920/149‑1 (MQ=255)
                                                                   cTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCCCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTAc                                                                               >  3:79158/1‑149 (MQ=255)
                                                                     gcgcCACCGCCAGGTACAGCCCCGCCGGGATCCCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTAccc                                                                             <  4:79158/149‑1 (MQ=255)
                                                                      cgcCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTAccc                                                                             <  3:403283/148‑1 (MQ=255)
                                                                        ccACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTg                                                                          <  3:155854/149‑1 (MQ=255)
                                                                         cACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGa                                                                                                                                           <  4:122869/83‑1 (MQ=255)
                                                                         cACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGa                                                                                                                                           >  3:122869/1‑83 (MQ=255)
                                                                          aCCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGc                                                                        <  8:338329/149‑1 (MQ=255)
                                                                          aCCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGc                                                                        >  5:352182/1‑149 (MQ=255)
                                                                              ccAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCCATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCcag                                                                     >  7:158882/1‑148 (MQ=255)
                                                                                     cAGCGCCGCCTGGATATGCTGCTCAAGTTCGGTGGAGTAATAGATGGCGCTCGCCTGCGTCGGCGATTCGAGGATATGTATGTTGTGCTCCTCTCCAATTTGACGAATCTTCAACGCCATGAAATGGCTACCCTTTGCGAGCAGCaaa                                                              <  6:316744/148‑2 (MQ=255)
                                                                                     cAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAAc                                                              <  6:168466/148‑1 (MQ=255)
                                                                                                ggATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCg                                                  >  5:188642/1‑149 (MQ=255)
                                                                                                      ccTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGACTTtc                                            >  7:207266/1‑149 (MQ=255)
                                                                                                       cTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCt                                               <  6:352182/145‑1 (MQ=255)
                                                                                                       cTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGAAGCAACGGCGCCGCACCGCCTTtc                                            >  8:41793/1‑148 (MQ=255)
                                                                                                        tGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTtctc                                          >  1:55652/1‑149 (MQ=255)
                                                                                                          cTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACTGCGCCGAACAGCCTTTCTCCg                                        >  5:310065/1‑149 (MQ=255)
                                                                                                          cTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCg                                        <  2:55652/149‑1 (MQ=255)
                                                                                                          cTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCg                                        >  6:308237/1‑149 (MQ=255)
                                                                                                           tCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCgg                                       >  6:50601/1‑149 (MQ=255)
                                                                                                            cAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                      >  3:102947/1‑149 (MQ=255)
                                                                                                            cAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTCGTGCTCCACGCCAATTTCACGAATCTTCAACGCCAt                                                                                          >  6:52127/1‑97 (MQ=255)
                                                                                                             aaGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGt                                     >  5:393028/1‑149 (MQ=255)
                                                                                                             aaGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCAATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                      <  6:142384/148‑1 (MQ=255)
                                                                                                                   gtttgAGTAATAGATCGCCCGCGCCAGCGCCTGAGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGTGTCGTAtt                               <  6:310065/145‑1 (MQ=255)
                                                                                                                        aGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCaaa                          <  8:32732/149‑1 (MQ=255)
                                                                                                                        aGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCaaa                          <  6:393028/149‑1 (MQ=255)
                                                                                                                         gTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCaaa                          <  5:308237/148‑1 (MQ=255)
                                                                                                                           ctcgaaTCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                      <  7:368979/129‑1 (MQ=255)
                                                                                                                           ctcgaaTCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                      >  8:368979/6‑134 (MQ=255)
                                                                                                                            aTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                      >  8:251098/1‑133 (MQ=255)
                                                                                                                            aTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                      <  7:251098/133‑1 (MQ=255)
                                                                                                                            aTAGATCGCCCGCGCCAGCGCCGGCGAATCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTAAACGCCATGAAATCGCTACCCATGGCCAGCAGCAACGGCGCCGCACCGCCt                                               >  3:250768/1‑124 (MQ=255)
                                                                                                                               gATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCg                   <  1:57976/149‑1 (MQ=255)
                                                                                                                               gATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATATTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGTAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCAcc                    <  7:79331/148‑1 (MQ=255)
                                                                                                                                aacGCCCGCGCCAGCGCCGGCGATTCGAGGATGTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTa                                                                                <  4:136252/85‑1 (MQ=255)
                                                                                                                                 tCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGAGCTCCACGCCAATTTCACGAATCTTCACCGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTATCTCCGGGTCGTATTTCAAAGACACCgaa                 >  8:190396/1‑147 (MQ=255)
                                                                                                                                    cccGCGCCAGCTCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATaa              <  3:65215/149‑1 (MQ=255)
                                                                                                                                     ccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCcagca                                                                   >  1:183929/1‑95 (MQ=255)
                                                                                                                                     ccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCcagca                                                                   <  2:183929/95‑1 (MQ=255)
                                                                                                                                     ccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAAATTCACGAATCTTCAACGCCATGAACTCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAAt             >  1:97936/1‑149 (MQ=255)
                                                                                                                                      cgcgCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGcc                                                <  6:186368/113‑1 (MQ=255)
                                                                                                                                      cgcgCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTcc                                         >  7:213354/1‑120 (MQ=255)
                                                                                                                                      cgcgCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTcc                                         <  8:213354/120‑1 (MQ=255)
                                                                                                                                         gcCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGt         <  1:244228/149‑1 (MQ=255)
                                                                                                                                          ccAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTc        >  2:256989/1‑149 (MQ=255)
                                                                                                                                           tggcgcGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCg       >  1:279675/3‑149 (MQ=255)
                                                                                                                                            aGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCgg      <  4:102947/149‑1 (MQ=255)
                                                                                                                                            aGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCCTCgg      >  5:196087/1‑149 (MQ=255)
                                                                                                                                                ccGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCGGGTTg  >  7:19966/1‑149 (MQ=255)
                                                                                                                                                ccGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCGGGTTg  >  6:230820/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCGGGTTG  >  NC_002947/4945483‑4945774

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