Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,602,181 T→A 5.2% D434V (GAT→GTT)  PP_4924 ← subtilase family serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,602,1810TA 5.2% 236.8 / 2.3 77D434V (GAT→GTT) PP_4924subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base T (41/32);  new base A (2/2);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

ACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCATGCAC  >  NC_002947/5602034‑5602329
                                                                                                                                                   |                                                                                                                                                    
aCCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATc                                                                                                                                                     >  1:128018/1‑149 (MQ=255)
 ccGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCg                                                                                                                                                    >  8:185139/1‑149 (MQ=255)
      tCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGtt                                                                                                                                               >  5:119149/1‑149 (MQ=255)
            tGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCCCACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGt                                                                                                                                          >  1:153365/1‑148 (MQ=255)
              gCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAACGCCGATGGCGTCat                                                                                                                                       <  2:313690/149‑2 (MQ=255)
                   gCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACtt                                                                                                                                  >  8:30202/1‑149 (MQ=255)
                   gCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACtt                                                                                                                                  <  7:185139/149‑1 (MQ=255)
                   gCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCGTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACtt                                                                                                                                  >  2:259445/1‑149 (MQ=255)
                    cTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTg                                                                                                                                 >  3:120537/1‑149 (MQ=255)
                      tGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGc                                                                                                                                >  3:48562/1‑148 (MQ=255)
                       gATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGTATGCCCGTCGTTGTTGCCTTCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGc                                                                                                                                <  1:104692/147‑1 (MQ=255)
                         tCTGGCTAATGCCATGCTCCTGCAATGCGCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCa                                                                                                                                        >  1:183088/1‑137 (MQ=255)
                         tCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCa                                                                                                                                        <  2:183088/137‑1 (MQ=255)
                            ggCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGggcg                                                                                                                         >  7:185115/1‑149 (MQ=255)
                             gCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGggcgg                                                                                                                        >  3:118120/1‑149 (MQ=255)
                                   gCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGc                                                                                                                  >  3:255729/1‑149 (MQ=255)
                                   gCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGc                                                                                                                  >  8:97835/1‑149 (MQ=255)
                                     cATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGc                                                                                                                  >  8:137823/1‑147 (MQ=255)
                                      aTGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCAt                                                                                                               >  5:123205/1‑149 (MQ=255)
                                        gCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATcc                                                                                                             >  4:260691/1‑149 (MQ=255)
                                          tCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCtt                                                                                                           >  1:129508/1‑149 (MQ=255)
                                            cTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCtt                                                                                                           <  4:255729/147‑1 (MQ=255)
                                               cAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACg                                                                                                      <  4:177782/149‑1 (MQ=255)
                                                    cccTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGt                                                                                                  <  6:195147/148‑1 (MQ=255)
                                                      cTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGttt                                                                                                >  5:32433/1‑148 (MQ=255)
                                                       ttAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAACGCCGATGGCGTc                                                                                                                                         <  8:120114/106‑1 (MQ=255)
                                                       ttAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAACGCCGATGGCGTc                                                                                                                                         >  7:120114/1‑106 (MQ=255)
                                                              aCCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCaga                                                                                         <  4:120537/147‑1 (MQ=255)
                                                                cGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTg                                                                                     >  8:241532/1‑149 (MQ=255)
                                                                cGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTg                                                                                     <  3:260691/149‑1 (MQ=255)
                                                                cGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTg                                                                                     >  3:300395/1‑149 (MQ=255)
                                                                cGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTg                                                                                     >  6:210614/1‑149 (MQ=255)
                                                                    aTCGGCACGCCCGTCGTGGTTGCCATCGCCCCACCCCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCAt                                                                                 >  8:72532/1‑149 (MQ=255)
                                                                       ggCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGaa                                                                                                        >  1:177607/1‑123 (MQ=255)
                                                                       ggCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGaa                                                                                                        <  2:177607/123‑1 (MQ=255)
                                                                       ggCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCAc                                                                              <  8:185115/149‑1 (MQ=255)
                                                                          aCGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgc                                                                           >  4:204391/1‑149 (MQ=255)
                                                                          aCGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgc                                                                           >  1:273046/1‑149 (MQ=255)
                                                                               cgtcgtGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACtgcctgcc                                                                      >  6:54122/1‑149 (MQ=255)
                                                                                    tGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTgcc                                                                 <  6:119149/149‑1 (MQ=255)
                                                                                      gTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTgccg                                                                <  4:48562/148‑1 (MQ=255)
                                                                                        tGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGccgccg                                                             >  5:296762/1‑149 (MQ=255)
                                                                                            aTCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCgtag                                                          >  4:193260/1‑148 (MQ=255)
                                                                                                aCCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTc                                                     <  4:300395/149‑1 (MQ=255)
                                                                                                 cccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                                    <  4:167154/149‑1 (MQ=255)
                                                                                                  ccACACCTTCAATTGCGCCCAGACGGCGTGTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                                    >  2:261849/1‑148 (MQ=255)
                                                                                                  ccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCgg                                                   <  2:273046/149‑1 (MQ=255)
                                                                                                    acacCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGgag                                                 <  3:204391/149‑1 (MQ=255)
                                                                                                      acCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGaa                                               <  7:241532/149‑1 (MQ=255)
                                                                                                       ccTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAAc                                              <  7:30202/149‑1 (MQ=255)
                                                                                                        cTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTc                                                     <  5:312781/141‑1 (MQ=255)
                                                                                                        cTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTc                                                     >  6:312781/1‑141 (MQ=255)
                                                                                                        cTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGaa                                               >  4:210514/1‑147 (MQ=255)
                                                                                                           cAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAgct                                          >  6:312705/1‑149 (MQ=255)
                                                                                                           cAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAgc                                           <  5:44541/148‑1 (MQ=255)
                                                                                                             aTTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGggc                                                                                                                          >  8:248219/1‑67 (MQ=255)
                                                                                                             aTTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGggc                                                                                                                          <  7:248219/67‑1 (MQ=255)
                                                                                                             aTTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAgctgg                                        >  3:14664/1‑149 (MQ=255)
                                                                                                             aTTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAgctgg                                        <  7:97835/149‑1 (MQ=255)
                                                                                                              ttGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAgctggc                                       >  7:120712/1‑149 (MQ=255)
                                                                                                                gcgcCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAgctggct                                      <  7:137823/148‑1 (MQ=255)
                                                                                                                gcgcCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAgctggct                                      <  6:32433/148‑1 (MQ=255)
                                                                                                                  gcCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAg                                            >  8:34628/1‑140 (MQ=255)
                                                                                                                  gcCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAg                                            <  7:34628/140‑1 (MQ=255)
                                                                                                                  gcCCAGACGGCGTATTCAACGTTGATGACGCAATCGGCGTTGGCGTCTCCACTTGCCAGGGAGGCGAAGCCATCCTTGAACGGCGGTACACAAGAGCTGCCAGGACTGCCTGCCGAGCCGGCGTAGTATTCGCAGAACAGCTGGCAGAc                                   >  6:97141/1‑149 (MQ=255)
                                                                                                                              tCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTc                                                                                               >  1:277165/1‑77 (MQ=255)
                                                                                                                              tCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTc                                                                                               <  2:277165/77‑1 (MQ=255)
                                                                                                                                 tCAACGTTGATGACGCCATCGGGGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGcc                    >  7:201373/1‑149 (MQ=255)
                                                                                                                                   aaCGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCg                   <  4:118120/148‑1 (MQ=255)
                                                                                                                                     cGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTa                <  6:123205/149‑1 (MQ=255)
                                                                                                                                         gatgaCGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGt            >  6:13036/1‑149 (MQ=255)
                                                                                                                                            gaCGCCAACGCCGATGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTg          >  4:56215/1‑148 (MQ=255)
                                                                                                                                              cGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCa       <  5:312705/149‑1 (MQ=255)
                                                                                                                                               gCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCAt      <  2:153365/149‑1 (MQ=255)
                                                                                                                                                  aTCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTg                               >  1:79160/1‑121 (MQ=255)
                                                                                                                                                  aTCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTg                               <  2:79160/121‑1 (MQ=255)
                                                                                                                                                  aTCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCATGCa   <  5:210614/149‑1 (MQ=255)
                                                                                                                                                   tCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCATGCAc  <  1:261849/149‑1 (MQ=255)
                                                                                                                                                   tCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCATGCAc  <  3:210514/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCATGCAC  >  NC_002947/5602034‑5602329

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: