Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,893,964 G→C 6.2% S365R (AGC→AGG) ‡ yegH ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,893,9640GC 6.2% 257.9 / 6.4 81S365R (AGC→AGG) ‡yegHputative transporter
Reads supporting (aligned to +/- strand):  ref base G (36/40);  new base C (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 = 2.32e-01

TGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAGGTCG  >  NC_002947/3893816‑3894111
                                                                                                                                                    |                                                                                                                                                   
tGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctgg                                                                                                                                                     <  3:294483/149‑1 (MQ=255)
tGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGCTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctgg                                                                                                                                                     <  1:316481/149‑1 (MQ=255)
   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                                                         >  4:354395/1‑147 (MQ=255)
                                                                                              gAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGa                                                         <  3:354395/147‑1 (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                                                                           <  3:241338/122‑1 (MQ=255)
                                                                                                     tcCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATcgcg                                                                           >  4:241338/1‑122 (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                              >  7:100842/1‑148 (MQ=255)
                                                                                                                        cccGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGc                              >  5:348273/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           >  7:151415/3‑149 (MQ=255)
                                                                                                                                          tctCGGGCAGCCAGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgc           >  5:377436/3‑149 (MQ=255)
                                                                                                                                            tCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcgc         <  3:627/149‑1 (MQ=255)
                                                                                                                                            tCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTgcgc         <  1:381432/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)
                                                                                                                                                    |                                                                                                                                                   
TGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAGGGCTCGTCCACGCGACCATCGCGGATCAGCGGCAGGCGCGAGTACGGCGAATTGGCCAGGGCCTGGGCAATGGCGTCGGCCGGTTGCGCCAGGTCG  >  NC_002947/3893816‑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: ‑