Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,602,173 T→A 5.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 5.3% 226.2 / 2.1 75I437F (ATC→TTC) PP_4924subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base T (33/38);  new base A (2/2);  total (35/40)
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.85e-01

CTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTT  >  NC_002947/5602032‑5602320
                                                                                                                                             |                                                                                                                                                    
cTACCGGCTCGCCCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGcc                                                                                                                                                >  2:113606/1‑148 (MQ=255)
         cGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTg                                                                                                                                      <  7:128156/149‑1 (MQ=255)
          gACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGtt                                                                                                                                       <  6:171715/147‑1 (MQ=255)
                gCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTc                                                                                                                                 <  8:143196/147‑1 (MQ=255)
                  cGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAac                                                                                                                              >  8:110410/1‑148 (MQ=255)
                  cGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAACGCCGATGGCGtcatca                                                                                                                             <  1:113606/149‑1 (MQ=255)
                         gATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGt                                                                                                                                  <  4:5065/137‑1 (MQ=255)
                         gATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGt                                                                                                                                  >  3:5065/1‑137 (MQ=255)
                           tCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGACCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAgg                                                                                                                    >  5:264428/1‑149 (MQ=255)
                            cTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAggg                                                                                                                   >  3:19021/1‑149 (MQ=255)
                            cTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCGCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAgg                                                                                                                    <  8:113010/148‑1 (MQ=255)
                               gCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGc                                                                                                                                    <  1:26341/129‑1 (MQ=255)
                               gCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGc                                                                                                                                    >  2:26341/1‑129 (MQ=255)
                               gCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGggcgg                                                                                                                >  4:280726/1‑149 (MQ=255)
                               gCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGggc                                                                                                                  <  6:264428/147‑1 (MQ=255)
                                cTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGggcggc                                                                                                               <  7:163690/149‑1 (MQ=255)
                                   aTGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGaa                                                                                                            <  7:119790/149‑1 (MQ=255)
                                    tGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAg                                                                                                          <  4:231317/149‑1 (MQ=255)
                                     gCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGc                                                                                                         >  5:286071/1‑149 (MQ=255)
                                        aTGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCAGGGCGGCGAAGCCATc                                                                                                     <  6:113699/149‑1 (MQ=255)
                                            tCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCCCACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGAGTCCACACTTTCCAGGGCGGAGAAGCCGTCCt                                                                                                   >  8:285672/1‑148 (MQ=255)
                                             ccTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCtt                                                                                                  >  6:62588/1‑148 (MQ=255)
                                              cTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTAGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTCCCACGGCGGCGAAGCCATCCTTg                                                                                                 >  2:84057/1‑148 (MQ=255)
                                                    tGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCg                                                                                          <  1:119528/149‑1 (MQ=255)
                                                      cccTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAAAAGCGCCCAGACGGCGTCGTCAACGTTGATGACGCCATCGGCGTGGGGGTCAACACTTGCCAGGGCGGCGAAGCCATCCGTGAACGTCGtt                                                                                        >  4:216249/1‑149 (MQ=255)
                                                           aGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCAc                                                                                    <  1:243981/148‑1 (MQ=255)
                                                           aGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCAc                                                                                    <  4:19021/148‑1 (MQ=255)
                                                            gCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCa                                                                                  <  2:230554/149‑1 (MQ=255)
                                                             cTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACTCTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCag                                                                                 <  3:280726/149‑1 (MQ=255)
                                                              tCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCagc                                                                                >  4:247859/1‑148 (MQ=255)
                                                              tCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCa                                                                                  >  3:107314/1‑147 (MQ=255)
                                                              tCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCa                                                                                  >  7:232558/1‑147 (MQ=255)
                                                                aCCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGc                                                                              >  7:41684/1‑149 (MQ=255)
                                                                    gCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGcc                                                                          <  7:273028/149‑1 (MQ=255)
                                                                         ggCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATTACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCAGCAGAGCTGCCATCAc                                                                     <  8:193484/149‑1 (MQ=255)
                                                                          gCACGCCCGTCGTGGTTGCCATCGACCCTCACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCAct                                                                    >  1:172938/1‑149 (MQ=255)
                                                                                      tGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTgcc                                                        <  3:2080/149‑1 (MQ=255)
                                                                                        gTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTgccgc                                                      >  7:1904/1‑149 (MQ=255)
                                                                                        gTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTgccg                                                       <  3:247859/148‑1 (MQ=255)
                                                                                            ccATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTGAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGGGAATGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGccgccg                                                    <  1:84057/148‑1 (MQ=255)
                                                                                            ccATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCgta                                                  <  5:253324/149‑1 (MQ=255)
                                                                                            ccATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCgta                                                  <  5:62588/149‑1 (MQ=255)
                                                                                              aTCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCgtag                                                 >  1:28172/1‑148 (MQ=255)
                                                                                                   cccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                           >  7:79017/1‑149 (MQ=255)
                                                                                                   cccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                           <  8:232558/149‑1 (MQ=255)
                                                                                                    ccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCgg                                          <  2:28172/149‑1 (MQ=255)
                                                                                                    ccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                           >  4:132915/1‑148 (MQ=255)
                                                                                                     cacaCCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGga                                         <  7:110410/149‑1 (MQ=255)
                                                                                                      acacCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                           >  2:223868/1‑146 (MQ=255)
                                                                                                         ccTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGaa                                      >  3:152346/1‑148 (MQ=255)
                                                                                                          cTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGaa                                      >  7:81621/1‑147 (MQ=255)
                                                                                                          cTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGaa                                      >  2:61402/1‑147 (MQ=255)
                                                                                                           ttCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAg                                   >  2:5219/1‑149 (MQ=255)
                                                                                                           ttCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAg                                   >  2:116571/1‑149 (MQ=255)
                                                                                                           ttCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAACgc                                                                                                                                           >  1:143087/1‑45 (MQ=255)
                                                                                                           ttCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAACgc                                                                                                                                           <  2:143087/46‑2 (MQ=255)
                                                                                                                ttGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTc                                            <  6:217154/135‑1 (MQ=255)
                                                                                                                ttGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTc                                            >  5:217154/1‑135 (MQ=255)
                                                                                                                    gcCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAgctggctg                            >  7:198553/1‑147 (MQ=255)
                                                                                                                     cccAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACg                         <  4:107314/149‑1 (MQ=255)
                                                                                                                       cAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGtt                       <  1:5219/149‑1 (MQ=255)
                                                                                                                        cgatggCGTCATCAACGTTGAAGACGCCGTCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgcct                                                                <  3:48013/103‑1 (MQ=38)
                                                                                                                        cgatggCGTCATCAACGTTGAAGACGCCGTCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgcct                                                                >  4:48013/5‑107 (MQ=38)
                                                                                                                        aGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACg                         >  7:81210/1‑146 (MQ=255)
                                                                                                                        aGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACg                         <  8:81210/146‑1 (MQ=255)
                                                                                                                        aGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTg                      <  6:286071/149‑1 (MQ=255)
                                                                                                                         gatggCGTCATCAACGTTGAAGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTc                                                                                      <  2:89846/81‑1 (MQ=38)
                                                                                                                         gatggCGTCATCAACGTTGAAGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTc                                                                                      >  1:89846/4‑84 (MQ=38)
                                                                                                                         gACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGc                     <  8:41684/149‑1 (MQ=255)
                                                                                                                            ggCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCa                  <  1:116571/149‑1 (MQ=255)
                                                                                                                               gTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCt               <  3:66365/149‑1 (MQ=255)
                                                                                                                                   tCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGcc           >  5:310987/1‑149 (MQ=255)
                                                                                                                                     aaCGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCg          <  1:223868/148‑1 (MQ=255)
                                                                                                                                     aaCGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCgg         >  5:90230/1‑148 (MQ=255)
                                                                                                                                         ttGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTAcc     <  2:172938/149‑1 (MQ=255)
                                                                                                                                            atgaCGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGtt  >  8:264515/1‑149 (MQ=255)
                                                                                                                                             tgaCGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGtt  <  8:1904/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CTACCGGCTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTT  >  NC_002947/5602032‑5602320

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: