Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,893,966 T→A 5.7% S365C (AGC→TGC) ‡ yegH ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,893,9660TA 5.7% 245.0 / 2.9 81S365C (AGC→TGC) ‡yegHputative transporter
Reads supporting (aligned to +/- strand):  ref base T (36/40);  new base A (3/2);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAGGTCG  >  NC_002947/3893819‑3894111
                                                                                                                                                   |                                                                                                                                                 
tCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGgctggctg                                                                                                                                                  >  8:217216/1‑149 (MQ=255)
 cAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGAGTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGCAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGc                                                                                                                                                 >  4:235694/1‑149 (MQ=255)
 cAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGc                                                                                                                                                 >  6:158559/1‑149 (MQ=255)
  aaCAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGcc                                                                                                                                                >  8:216629/1‑149 (MQ=255)
   aCAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGccc                                                                                                                                               >  8:96112/1‑149 (MQ=255)
      gCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCAGCCAGCCCgac                                                                                                                                            <  3:360952/149‑2 (MQ=255)
       cccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCgaga                                                                                                                                           <  4:306536/149‑1 (MQ=255)
       cccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCga                                                                                                                                             >  7:254928/1‑147 (MQ=255)
        ccGGTGAAGTCCCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGagag                                                                                                                                          <  8:200477/149‑1 (MQ=255)
         cGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGTGTCGACTGACCGCGAATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGTTCGGGCTGGATGCCCGagag                                                                                                                                          >  6:359058/1‑148 (MQ=255)
         cGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagc                                                                                                                                         <  7:216629/149‑1 (MQ=255)
            tGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagt                                                                                                                                      <  7:96112/149‑1 (MQ=255)
               aGTCACCGAACTCGTTGACCACGAAGGCAATATGGGGCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcc                                                                                                                                   <  1:286862/149‑1 (MQ=255)
                gTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcct                                                                                                                                  >  2:381432/1‑149 (MQ=255)
                gTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCAGCCAGCCCgacctgga                                                                                                                                       >  6:391697/1‑138 (MQ=255)
                gTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCAGCCAGCCCgacctgga                                                                                                                                       <  5:391697/144‑7 (MQ=255)
                   aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttca                                                                                                                               <  6:339183/149‑1 (MQ=255)
                   aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttc                                                                                                                                <  4:345114/148‑1 (MQ=255)
                   aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt                                                                                                                                 <  8:254928/147‑1 (MQ=255)
                    ccGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcag                                                                                                                              <  5:370813/149‑1 (MQ=255)
                     cGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCgag                                                                                                                                            >  1:143983/1‑134 (MQ=255)
                     cGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCgag                                                                                                                                            <  2:143983/134‑1 (MQ=255)
                          tCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagcagttc                                                                                                                        >  1:140403/1‑149 (MQ=255)
                            gTTGACCACGAAGTCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGATGATGATGAAGCTCTCCAGCAGGTTAAAAGGCGCGCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagcagttcc                                                                                                                       <  5:359058/148‑1 (MQ=255)
                                   aCGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTtgtg                                                                                                               >  1:187806/1‑149 (MQ=255)
                                     gAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGca                                                                                                             <  3:70215/149‑1 (MQ=255)
                                        ggCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGcaca                                                                                                           >  1:222175/1‑148 (MQ=255)
                                         gCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACaaa                                                                                                         <  6:215377/149‑1 (MQ=255)
                                           aaTATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGc                                                                                                       >  3:192099/1‑149 (MQ=255)
                                           aaTATGGGTAGAATGACAGCGCATTTGATCCACGGCATTCAGGATGCTGAAGCTCTCCAGCAGGTTCAAAGGCGCCCGTGCCATGCTTTCCAGGTCTGGCAGGCTGCCCGAGAGCAGTTCCTACAGCAGTTCCATCTTTAGCACaaaca                                                                                                       >  8:319976/1‑147 (MQ=255)
                                              atGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGCGCCATGCTTTCCACGACGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGAACAAAGccc                                                                                                     >  1:385688/1‑148 (MQ=255)
                                                  gTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAggg                                                                                                >  4:627/1‑149 (MQ=255)
                                                             gcgcATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCAcgc                                                                                     >  6:260608/1‑149 (MQ=255)
                                                                cATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGAc                                                                                  <  2:222175/149‑1 (MQ=255)
                                                                 aTTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGAcc                                                                                 <  3:291289/149‑1 (MQ=255)
                                                                 aTTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGAc                                                                                  >  8:122916/1‑148 (MQ=255)
                                                                       tCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCACCGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGGCCACGCGACCATcgcg                                                                           >  8:305219/1‑149 (MQ=255)
                                                                       tCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                           >  5:107815/1‑149 (MQ=255)
                                                                         cAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                           <  5:158559/147‑1 (MQ=255)
                                                                         cAGGGATTTGAGGATGCTGAAGCTCTCAAGCAGGTTCAACGGCGCCCCGGCCATGCTTTCCTGGTCGGGCTGGCTGCCGGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgc                                                                            <  7:319976/146‑1 (MQ=255)
                                                                          aGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGAt                                                                        >  8:61460/1‑149 (MQ=255)
                                                                           gggCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGAt                                                                        >  6:2984/1‑148 (MQ=255)
                                                                            ggCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCa                                                                      <  3:19923/149‑1 (MQ=255)
                                                                             gCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAg                                                                     <  2:187806/149‑1 (MQ=255)
                                                                                     gATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGgc                                                             >  3:197672/1‑149 (MQ=255)
                                                                                       tGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAGAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGgcgc                                                           <  7:217216/149‑1 (MQ=255)
                                                                                        gCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGgcgc                                                           >  4:161618/1‑148 (MQ=255)
                                                                                         cTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGgcgc                                                           >  3:69929/1‑147 (MQ=255)
                                                                                         cTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                         <  3:235694/149‑1 (MQ=255)
                                                                                          tGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAg                                                        >  8:64233/1‑149 (MQ=255)
                                                                                          tGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAg                                                        >  7:241244/1‑149 (MQ=255)
                                                                                           gAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                         <  3:354395/147‑1 (MQ=255)
                                                                                           gAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                         >  4:354395/1‑147 (MQ=255)
                                                                                            aaGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAg                                                        <  7:64233/147‑1 (MQ=255)
                                                                                            aaGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGt                                                       <  7:324560/148‑1 (MQ=255)
                                                                                             aGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTAc                                                     <  5:2984/149‑1 (MQ=255)
                                                                                                tctcCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGc                                                  >  2:47718/1‑149 (MQ=255)
                                                                                                  tcCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                           >  4:241338/1‑122 (MQ=255)
                                                                                                  tcCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                           <  3:241338/122‑1 (MQ=255)
                                                                                                      gcagGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTg                                            <  2:140403/149‑1 (MQ=255)
                                                                                                      gcagGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTg                                            <  4:192099/149‑1 (MQ=255)
                                                                                                       cagGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGGTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTgg                                           <  7:305219/149‑1 (MQ=255)
                                                                                                          gTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCa                                        <  7:122916/149‑1 (MQ=255)
                                                                                                             cAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAggg                                     >  7:255167/1‑149 (MQ=255)
                                                                                                             cAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAgg                                      <  2:385688/148‑1 (MQ=255)
                                                                                                              aaCGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGc                                    <  4:191911/149‑1 (MQ=255)
                                                                                                                 ggCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTg                                 <  8:255167/149‑1 (MQ=255)
                                                                                                                  gcgcCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTgg                                <  4:197672/149‑1 (MQ=255)
                                                                                                                     cccGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGc                              >  5:348273/1‑148 (MQ=255)
                                                                                                                     cccGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGc                              >  7:100842/1‑148 (MQ=255)
                                                                                                                      ccGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCATTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGTCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCaa                            <  8:302816/149‑1 (MQ=255)
                                                                                                                      ccGGGCCATGCTTTCCAGGTCGGGCTGGCTGCACGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCaa                            >  7:341286/1‑149 (MQ=255)
                                                                                                                          gCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGCCCAGGCCCTGGCCAAttcg                        <  6:107815/149‑4 (MQ=255)
                                                                                                                            cATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGt                      >  1:344640/1‑149 (MQ=255)
                                                                                                                                cTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCATCGGc                  <  8:241244/149‑1 (MQ=255)
                                                                                                                                 tttCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGcc                 >  2:10473/1‑149 (MQ=255)
                                                                                                                                    ccAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCgg               <  3:161618/148‑1 (MQ=255)
                                                                                                                                      aGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTg            >  3:50841/1‑149 (MQ=255)
                                                                                                                                       tctCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgc           >  5:377436/3‑149 (MQ=255)
                                                                                                                                       tctCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgc           >  7:151415/3‑149 (MQ=255)
                                                                                                                                         tCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcgc         <  1:381432/149‑1 (MQ=255)
                                                                                                                                         tCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcgc         <  3:627/149‑1 (MQ=255)
                                                                                                                                            ggctggctGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAg      >  3:39792/1‑149 (MQ=255)
                                                                                                                                            ggctggctGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAg      <  5:260608/149‑1 (MQ=255)
                                                                                                                                                ggctgCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAGGTCg  <  8:341286/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAGGTCG  >  NC_002947/3893819‑3894111

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