Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,602,173 T→A 6.3% I437F (ATC→TTC)  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,1730TA 6.3% 189.4 / 3.0 63I437F (ATC→TTC) PP_4924subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base T (35/24);  new base A (2/2);  total (37/26)
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.37e-01

CCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTG  >  NC_002947/5602026‑5602321
                                                                                                                                                    |                                                                                                                                                    
ccTCTCCTACCGGCTCGACCTGTCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCGATTGCGCCCAGACGGCGTCTTCAACGTtgatg                                                                                                                                                     <  8:43387/149‑1 (MQ=255)
      cTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCATCAACGTTGAAGACGc                                                                                                                                                 <  7:226600/147‑1 (MQ=255)
       tACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAt                                                                                                                                              >  7:184187/1‑149 (MQ=255)
              tCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGAc                                                                                                                                                   >  4:132958/1‑137 (MQ=255)
                gACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGtt                                                                                                                                       <  1:64740/148‑1 (MQ=255)
                gACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGACATCGGCACGCCCGTAGTGGTTGCCATCGACCCACACATTCAATTGCGCCCAGACGGCGTCTTAAACGTTGATGACGCCAACGGCGTTgg                                                                                                                                     >  7:204276/1‑149 (MQ=255)
                   cTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGTCGTTGGCGt                                                                                                                                  <  6:28919/149‑1 (MQ=255)
                   cTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCAACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGt                                                                                                                                  <  7:148934/149‑1 (MQ=255)
                        cGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGAAGCCATCGCCCTTGGCGTCAac                                                                                                                              >  7:293303/1‑148 (MQ=255)
                         gAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAacac                                                                                                                            >  4:64746/1‑149 (MQ=255)
                            cTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTg                                                                                                                         >  8:223086/1‑149 (MQ=255)
                                        aaTGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGggcggc                                                                                                               >  2:142837/1‑147 (MQ=255)
                                        aaTGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGggcggc                                                                                                               <  1:142837/147‑1 (MQ=255)
                                         aTGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGaa                                                                                                            <  1:217952/149‑1 (MQ=255)
                                          tGCCATGCTCCTGCAATGCCCTCAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAg                                                                                                           <  1:148361/149‑1 (MQ=255)
                                            ccATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCATCAACGTTGAAGACGcc                                                                                                                                                <  5:193636/110‑1 (MQ=255)
                                            ccATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCATCAACGTTGAAGACGcc                                                                                                                                                >  6:193636/1‑110 (MQ=255)
                                              aTGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGaa                                                                                                            <  8:230092/144‑1 (MQ=255)
                                              aTGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCAt                                                                                                       >  2:31240/1‑149 (MQ=255)
                                              aTGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCAGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGaa                                                                                                            >  7:230092/1‑144 (MQ=255)
                                                  tCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCtt                                                                                                   >  4:46203/1‑149 (MQ=255)
                                                    cTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTg                                                                                                  >  7:10652/1‑148 (MQ=255)
                                                       cAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGaa                                                                                                <  8:89967/147‑1 (MQ=255)
                                                        aaTGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACgg                                                                                             >  3:155377/1‑149 (MQ=255)
                                                              cTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGttt                                                                                        >  3:258297/1‑148 (MQ=255)
                                                                 aGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTg                                                                                                                                      <  4:127095/99‑1 (MQ=255)
                                                                 aGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTg                                                                                                                                      >  3:127095/1‑99 (MQ=255)
                                                                    tCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCa                                                                                   >  1:33938/1‑147 (MQ=255)
                                                                    tCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCa                                                                                   >  4:92271/1‑147 (MQ=255)
                                                                    tCACCGGCAGCGGCACGGCCGTCGTGGGTGCGATCGACGCGCCCCTGCAATTGCGCCGAGACGGCGTCTTCAGCGTTGATGACGCGATGGGCGTTGGCGTCAGCAGTGGCCAGGGCGGCGAGGCCATCCTTGAACGGCGTTTCACCa                                                                                   >  6:241765/1‑147 (MQ=255)
                                                                         ggCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGc                                                                            <  7:223086/149‑1 (MQ=255)
                                                                          gaaTCGGCACGCCCGTCGTGGTTGCCATCGAACCAAACCTTCAATTGCGCCCAGACGTCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTTCCAGGGCGGCGAAGTCATCCTTGATCGTCGTTTCACCAGAGCTGcc                                                                           <  8:204276/147‑1 (MQ=255)
                                                                              cGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCa                                                                       >  1:324356/1‑149 (MQ=255)
                                                                               ggCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCa                                                                       <  2:127470/148‑1 (MQ=255)
                                                                                  aCGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgc                                                                   >  2:65921/1‑149 (MQ=255)
                                                                                      ccGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgcctgc                                                               >  7:241317/1‑149 (MQ=255)
                                                                                      ccGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgcctgc                                                               >  5:176853/1‑149 (MQ=255)
                                                                                           gtGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTgc                                                          >  3:205410/1‑149 (MQ=255)
                                                                                                  ccATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCgta                                                   <  3:92271/149‑1 (MQ=255)
                                                                                                  ccATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCgta                                                   <  4:317157/149‑1 (MQ=255)
                                                                                                   cATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCgtag                                                  >  4:260681/1‑149 (MQ=255)
                                                                                                    aTCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCgtag                                                  >  7:19131/1‑148 (MQ=255)
                                                                                                     tCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCt                                                                              <  8:71611/119‑1 (MQ=255)
                                                                                                     tCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCgtagta                                                <  8:293303/149‑1 (MQ=255)
                                                                                                        aCCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTc                                             <  3:46203/149‑1 (MQ=255)
                                                                                                        aCCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTc                                             <  1:65921/149‑1 (MQ=255)
                                                                                                         cccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTc                                             <  4:155377/148‑1 (MQ=255)
                                                                                                         cccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                            <  8:10652/149‑1 (MQ=255)
                                                                                                          ccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAacac                                                                                                                            >  4:118674/1‑68 (MQ=255)
                                                                                                          ccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAacac                                                                                                                            <  3:118674/68‑1 (MQ=255)
                                                                                                            acacCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                            >  7:177566/1‑146 (MQ=255)
                                                                                                            acacCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGgag                                         <  8:241317/149‑1 (MQ=255)
                                                                                                                cTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGaa                                       >  6:224806/1‑147 (MQ=255)
                                                                                                                 ttCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAg                                    >  6:316761/1‑149 (MQ=255)
                                                                                                                           gcCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGAc                           >  2:292592/1‑149 (MQ=255)
                                                                                                                             ccAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGt                         >  8:48584/1‑149 (MQ=255)
                                                                                                                               aGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTg                       >  1:35618/1‑149 (MQ=255)
                                                                                                                                    gCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCAc                  <  8:184187/149‑1 (MQ=255)
                                                                                                                                     cGTCATCAACGTTGAAGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCCCTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCAcc                 >  5:172077/1‑149 (MQ=255)
                                                                                                                                         ttCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGc             >  5:313873/1‑149 (MQ=255)
                                                                                                                                          tCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGcc            >  5:204463/1‑149 (MQ=255)
                                                                                                                                           cAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGTTGGCTGACGTTGCTCACCTGGcc            <  7:48584/148‑1 (MQ=255)
                                                                                                                                           cAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCg           >  2:65501/1‑149 (MQ=255)
                                                                                                                                               gTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTAc       >  1:306518/1‑149 (MQ=255)
                                                                                                                                                  gatgaCGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGt    >  6:128830/1‑149 (MQ=255)
                                                                                                                                                   ttgaCGCCAACGCCGATGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgcctgc                                                               >  5:159178/2‑89 (MQ=39)
                                                                                                                                                   ttgaCGCCAACGCCGATGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgcctgc                                                               <  6:159178/88‑1 (MQ=39)
                                                                                                                                                   atgaCGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCAGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGa                            >  8:282563/1‑124 (MQ=255)
                                                                                                                                                   atgaCGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCAGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGa                            <  7:282563/124‑1 (MQ=255)
                                                                                                                                                    tgaCGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGtt   <  8:19131/148‑1 (MQ=255)
                                                                                                                                                    tgaCGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTg  <  4:258297/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTCTCCTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTG  >  NC_002947/5602026‑5602321

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: