Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,623,383 G→A 6.1% intergenic (‑27/+23) 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,3830GA 6.1% 241.4 / 3.6 82intergenic (‑27/+23)PP_1422/PP_1423hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/34);  new base A (3/2);  total (46/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

TTCCGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCACCGCGA  >  NC_002947/1623237‑1623530
                                                                                                                                                  |                                                                                                                                                    
ttCCGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGcc                                                                                                                                                    >  5:126913/1‑149 (MQ=255)
   cGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCcaca                                                                                                                                                 <  4:55759/149‑1 (MQ=255)
   cGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCcac                                                                                                                                                  <  3:23129/148‑1 (MQ=255)
    gCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAg                                                                                                                                                <  1:49803/149‑1 (MQ=255)
      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                                                                                               <  3:387481/149‑1 (MQ=255)
                                                    aGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGtt                                                                                               >  1:84998/1‑149 (MQ=255)
                                                    aGGTCAGCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGtt                                                                                               <  6:230229/149‑1 (MQ=255)
                                                      gTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAgcg                                                                                                                          >  3:273601/1‑121 (MQ=255)
                                                      gTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAgcg                                                                                                                          <  4:273601/121‑1 (MQ=255)
                                                           ccTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTaa                                                                                          <  3:232862/147‑1 (MQ=255)
                                                           ccTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTaa                                                                                          >  4:232862/1‑147 (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                                                                                                                        >  8:21750/1‑103 (MQ=255)
                                                                          cGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCt                                                                                                                        <  7:21750/103‑1 (MQ=255)
                                                                               aaTTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGgc                                                                    >  5:99400/1‑149 (MQ=255)
                                                                                     gCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCgg                                                                <  2:51969/147‑1 (MQ=255)
                                                                                     gCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCgg                                                                >  1:51969/1‑147 (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                                                       <  8:399501/120‑1 (MQ=255)
                                                                                                                         tCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAg                                                       >  7:399501/1‑120 (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    >  8:146629/1‑149 (MQ=255)
                                                                                                                                               cccGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCACcgc    >  2:407919/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)
                                                                                                                                                  |                                                                                                                                                    
TTCCGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGCCAGCTTGTTCATGGGCATGTCTACCGCAGTGGTGGTGCGGCCCCGCCACAGGGGTGGGGCAGGCTGGTCAGCGCTGGATTTCCCGGTCGCGGTCCAGTTCGTAAGGCAGGGTCAGCACCTGCAGGCGCGGGCCTTCGAGGCTGCCCAGGTGCAGGTGGCCGTCTTCCACCGCCTCGGCGCTGAGCACCGCGA  >  NC_002947/1623237‑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: