Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,623,390 G→T 6.5% intergenic (‑34/+16) PP_1422 ← / ← PP_1423 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,623,3900GT 6.5% 245.4 / 6.2 77intergenic (‑34/+16)PP_1422/PP_1423hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (42/30);  new base T (3/2);  total (45/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-01

TGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCACCGCGA  >  NC_002947/1623243‑1623530
                                                                                                                                                   |                                                                                                                                             
tGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAggg                                                                                                                                              <  1:247087/149‑1 (MQ=255)
       cgGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCACCCCTGTGGCGGGGc                                                                                                                                       <  6:277751/149‑1 (MQ=255)
       cgGATGCGCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGc                                                                                                                                       <  2:377715/149‑1 (MQ=255)
              tCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTgg                                                                                                                                >  8:331509/1‑149 (MQ=255)
              tCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTgg                                                                                                                                >  4:279245/1‑149 (MQ=255)
              tCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTgg                                                                                                                                >  4:74416/1‑149 (MQ=255)
              tCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCACCCCTGTGGCGGGGc                                                                                                                                       <  3:420895/142‑1 (MQ=255)
              tCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCACCCCTGTGGCGGGGc                                                                                                                                       >  4:420895/1‑142 (MQ=255)
                     aaCCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAgcgc                                                                                                                         >  6:229638/1‑149 (MQ=255)
                      aCCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCt                                                                                                                        <  1:182500/149‑1 (MQ=255)
                      aCCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCt                                                                                                                        >  2:182500/1‑149 (MQ=255)
                        cTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGCCATGTCTACCGCAGTTGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTgg                                                                                                                     <  5:298804/149‑1 (MQ=255)
                          cGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGAtt                                                                                                                  >  4:283820/1‑149 (MQ=255)
                             gCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGc                                                                                                                                   >  5:200239/1‑131 (MQ=255)
                             gCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGc                                                                                                                                   <  6:200239/131‑1 (MQ=255)
                             gCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTg                                                                                                                      >  5:222741/1‑143 (MQ=255)
                             gCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTg                                                                                                                      <  6:222741/143‑1 (MQ=255)
                             gCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTc                                                                                                                <  4:224733/149‑1 (MQ=255)
                             gCAGGGTTGGCAGGACCAGGTCATACTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTc                                                                                                                >  6:213176/1‑149 (MQ=255)
                              cAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTcc                                                                                                               <  6:126913/149‑1 (MQ=255)
                                gggTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCGACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCg                                                                                                             >  1:220722/1‑149 (MQ=255)
                                 ggTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCgg                                                                                                            <  5:108597/149‑1 (MQ=255)
                                    tGGCAGGACCAGGTCATACTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTcg                                                                                                         >  2:65621/1‑149 (MQ=255)
                                      gCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTg                                                                                                                                 <  2:14221/124‑1 (MQ=255)
                                      gCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTg                                                                                                                                 >  1:14221/1‑124 (MQ=255)
                                      gCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTcgcg                                                                                                       <  8:43141/149‑1 (MQ=255)
                                         ggACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTc                                                                                                    >  2:351384/1‑149 (MQ=255)
                                           aCCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCa                                                                                                  <  5:159741/149‑1 (MQ=255)
                                           aCCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCACCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTc                                                                                                    >  4:387481/1‑147 (MQ=255)
                                             cAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGt                                                                                                <  2:195903/149‑1 (MQ=255)
                                              aGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGtt                                                                                               >  1:84998/1‑149 (MQ=255)
                                              aGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGtt                                                                                               <  3:387481/149‑1 (MQ=255)
                                              aGGTCAGCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGtt                                                                                               <  6:230229/149‑1 (MQ=255)
                                                gTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAgcg                                                                                                                          >  3:273601/1‑121 (MQ=255)
                                                gTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAgcg                                                                                                                          <  4:273601/121‑1 (MQ=255)
                                                     ccTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTaa                                                                                          >  4:232862/1‑147 (MQ=255)
                                                     ccTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTaa                                                                                          <  3:232862/147‑1 (MQ=255)
                                                     ccTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAgg                                                                                        >  2:277390/1‑149 (MQ=255)
                                                       ttttCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAgg                                                                                        >  4:327941/1‑147 (MQ=255)
                                                          tCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAggg                                                                                   <  3:283820/148‑1 (MQ=255)
                                                           cGATGGCGTCGTGCAATTGTGCTTGAAACAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTCGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAgg                                                                                    <  1:65621/148‑1 (MQ=255)
                                                           cGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGCGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGt                                                                                  >  2:401495/1‑149 (MQ=255)
                                                            gATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGt                                                                                  <  3:279245/148‑1 (MQ=255)
                                                              tGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAg                                                                               <  3:327941/149‑1 (MQ=255)
                                                                    cGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCt                                                                                                                        <  7:21750/103‑1 (MQ=255)
                                                                    cGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCt                                                                                                                        >  8:21750/1‑103 (MQ=255)
                                                                         aaTTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGgc                                                                    >  5:99400/1‑149 (MQ=255)
                                                                               gCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCgg                                                                >  1:51969/1‑147 (MQ=255)
                                                                               gCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCgg                                                                <  2:51969/147‑1 (MQ=255)
                                                                                cTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGcc                                                             >  1:375912/1‑149 (MQ=255)
                                                                                   gAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCt                                                            <  7:331509/147‑1 (MQ=255)
                                                                                      ccAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAg                                                       <  2:25943/149‑1 (MQ=255)
                                                                                       cAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGAGGTCCAGTTCGTAAGGAAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAgg                                                      >  4:224469/1‑149 (MQ=255)
                                                                                         aCTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCt                                                    >  7:165242/1‑149 (MQ=255)
                                                                                         aCTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCt                                                    >  4:182458/1‑149 (MQ=255)
                                                                                                      gTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCcaggtgcag                                        <  2:84998/148‑1 (MQ=255)
                                                                                                       ttCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCGACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCcaggtgcaggt                                      <  2:220722/149‑1 (MQ=255)
                                                                                                            gggCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCg                                 >  3:95437/1‑149 (MQ=255)
                                                                                                              gCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTc                               >  4:414051/1‑149 (MQ=255)
                                                                                                               cATGTCTACCGCAGTGGTGGGGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGGCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCt                              <  3:23932/149‑1 (MQ=255)
                                                                                                                   tCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAg                                                       >  7:399501/1‑120 (MQ=255)
                                                                                                                   tCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAg                                                       <  8:399501/120‑1 (MQ=255)
                                                                                                                     tACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCAcc                        >  7:333780/1‑149 (MQ=255)
                                                                                                                       ccGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGc                      >  7:314624/1‑149 (MQ=255)
                                                                                                                        cGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGc                      >  2:159306/1‑148 (MQ=255)
                                                                                                                                  gtgCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCtta           <  5:213176/149‑3 (MQ=255)
                                                                                                                                   tgCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAg          <  2:375912/149‑1 (MQ=255)
                                                                                                                                    gCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAg          >  1:350381/1‑148 (MQ=255)
                                                                                                                                    cctgCCCCACCCCTGTGGCGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAg          >  5:404312/4‑148 (MQ=255)
                                                                                                                                     ctgCCCCACCCCTGTGGCGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCa        >  4:49894/3‑149 (MQ=255)
                                                                                                                                       gCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCCGCACCTGCAGGCGCGGGCCTTCTATGCTGCCCAGGTGCAGGTGGCCGTCTTACACCGCCTCTGCCCTGAGCAcc      >  3:225384/1‑149 (MQ=255)
                                                                                                                                       gCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCAcc      >  6:35995/1‑149 (MQ=255)
                                                                                                                                        ccccGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCACcg     >  2:172989/1‑149 (MQ=255)
                                                                                                                                         cccGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCACcgc    >  2:407919/1‑149 (MQ=255)
                                                                                                                                         cccGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCACcgc    >  8:146629/1‑149 (MQ=255)
                                                                                                                                         cccGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGAATTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCACcgc    >  7:146067/1‑149 (MQ=255)
                                                                                                                                          ccGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCACcgc    >  3:37921/1‑148 (MQ=255)
                                                                                                                                           cGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCACCGCGa  >  3:393795/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCACCGCGA  >  NC_002947/1623243‑1623530

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: