Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,602,235 G→C 6.8% A433G (GCC→GGC)  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,2350GC 6.8% 217.8 / 5.3 69A433G (GCC→GGC) PP_4924subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base G (31/33);  new base C (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCATGCA  >  NC_002947/5602090‑5602379
                                                                                                                                                 |                                                                                                                                                
cTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGtt                                                                                                                                              <  2:73329/150‑1 (MQ=255)
       tGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAACGCCGATGGCGTc                                                                                                                                        >  1:810871/1‑149 (MQ=255)
       tGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCAACGCCGATGGCGTc                                                                                                                                        <  2:810871/149‑1 (MQ=255)
           cGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAacac                                                                                                                                   <  2:884937/150‑1 (MQ=255)
                    tCTGGCTAATGCCATGCTCCTGCAATGGCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTTAATTGCGCCCAGACGGCGTCTTCAACGTGGATGACGCCATCGGCGTTGGCGTACACACTTGCCAggg                                                                                                                          >  2:541284/1‑150 (MQ=255)
                         cTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGTTGACGCCAACGCCGATGGCGTCatc                                                                                                                                     <  2:208146/134‑3 (MQ=255)
                         cTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGTTGACGCCAACGCCGATGGCGTCatc                                                                                                                                     >  1:208146/1‑132 (MQ=255)
                           aaTGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGaa                                                                                                                   <  1:861382/150‑1 (MQ=255)
                           aaTGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGaa                                                                                                                   <  2:211907/150‑1 (MQ=255)
                               ccATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCa                                                                                                               >  2:801940/1‑150 (MQ=255)
                                 aTGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATc                                                                                                             <  1:479808/150‑1 (MQ=255)
                                   gCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCt                                                                                                           >  1:1002697/1‑150 (MQ=255)
                                   gCTCCTGAAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCt                                                                                                           >  1:517107/1‑150 (MQ=255)
                                       cTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGaa                                                                                                       >  1:974714/1‑150 (MQ=255)
                                        tGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTg                                                                                                                                             >  2:608507/1‑111 (MQ=255)
                                        tGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTg                                                                                                                                             <  1:608507/111‑1 (MQ=255)
                                           aaTGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGc                                                                                                   <  2:894286/150‑1 (MQ=255)
                                                ccTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCAACGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTc                                                                                              >  1:66403/1‑150 (MQ=255)
                                                         aCCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCa                                                                                                                             >  2:1056750/1‑110 (MQ=255)
                                                         aCCGGCATCCGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCa                                                                                                                             <  1:1056750/110‑1 (MQ=255)
                                                          ccGGCATCGGCACGCCCGTCGTGGTTGCCACCGACTGACACCTTCAATTTCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTg                                                                                    >  1:933190/1‑150 (MQ=255)
                                                           cGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGc                                                                                   >  2:423419/1‑150 (MQ=255)
                                                              cATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCAt                                                                                <  1:801940/150‑1 (MQ=255)
                                                                 cGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCAc                                                                             <  1:423419/150‑1 (MQ=255)
                                                                  ggCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATc                                                                                                             >  1:583724/1‑117 (MQ=255)
                                                                  ggCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATc                                                                                                             <  2:583724/117‑1 (MQ=255)
                                                                   gCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCaa                                                                                                                                      >  2:426787/1‑91 (MQ=255)
                                                                   gCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCaa                                                                                                                                      <  1:426787/91‑1 (MQ=255)
                                                                    cACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgc                                                                          <  1:231255/150‑1 (MQ=255)
                                                                     aCGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgcc                                                                         >  1:117929/1‑150 (MQ=255)
                                                                      cGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgcct                                                                        <  2:517107/150‑1 (MQ=255)
                                                                        cccGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCAAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCActgcctgc                                                                      <  1:1053843/150‑1 (MQ=255)
                                                                         ccGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACtgcctgcc                                                                     >  1:183833/1‑150 (MQ=255)
                                                                            tcgtGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAggg                                                                                                                          >  1:612060/1‑94 (MQ=255)
                                                                            tcgtGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAggg                                                                                                                          <  2:612060/94‑1 (MQ=255)
                                                                                    gCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCgta                                                          >  2:845757/1‑150 (MQ=255)
                                                                                       aTCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGt                                                                                                                                         >  1:216151/1‑68 (MQ=255)
                                                                                       aTCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGt                                                                                                                                         <  2:216151/68‑1 (MQ=255)
                                                                                        tCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTAt                                                      >  2:987962/1‑150 (MQ=255)
                                                                                           aCCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                                   <  1:541284/150‑1 (MQ=255)
                                                                                           aCCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                                   <  1:923458/150‑1 (MQ=255)
                                                                                           aCCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                                   >  1:685287/1‑150 (MQ=255)
                                                                                            cccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCgg                                                  <  2:1002697/150‑1 (MQ=255)
                                                                                             ccACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCCTCCTTGACCGGCGTTTCCGAAGAACTGCCTGCACTGCCTGTCTTTCTGCCGTAGTCTTccgc                                                 >  2:11406/1‑147 (MQ=255)
                                                                                               acacCATCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCTGCGTTGGCGACAACACTTAACAAGGCGGCGAAGCCATCCTTAAACTGCGTTTCACCAGAGCTGTCATCACTGCCTGCCTTCCCGCCGTACTATTCGgaaa                                               >  2:128848/1‑148 (MQ=255)
                                                                                                   cTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCa                                                                                             <  1:630728/100‑1 (MQ=255)
                                                                                                   cTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCa                                                                                             >  2:630728/1‑100 (MQ=255)
                                                                                                    ttCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAgc                                          >  1:791007/1‑150 (MQ=255)
                                                                                                           gcgcCCAGACGGCGTCTTCAACGTTGATGACGCCCTCGGCGTTGGCGTCAACAATTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCCGAGCTGCCATCACTGCCTCCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGa                                   <  1:998870/150‑1 (MQ=255)
                                                                                                           gcgcCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGa                                   <  2:183833/150‑1 (MQ=255)
                                                                                                                 aGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCg                                                                                                  <  2:19980/81‑1 (MQ=255)
                                                                                                                 aGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCg                                                                                                  >  1:19980/1‑81 (MQ=255)
                                                                                                                     ggtgTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCAc                         <  1:72375/147‑1 (MQ=255)
                                                                                                                       cGTCTTCAACGTTGATGACGCCATCTGCGTTGGCGTCAACACTTGCCAGGGCGGCTAAGCCATGCTTGAACGGCGTTTCACCAGGGCTGCCATCACTTCTTGCCTGGCCGCCGTAGTATTCGGAGAACTGCTGGCTGACGTTGCTCACCt                       <  2:360179/150‑1 (MQ=255)
                                                                                                                           ttCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTAtt                                                     <  2:1117225/116‑1 (MQ=255)
                                                                                                                           ttCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTAtt                                                     >  1:1117225/1‑116 (MQ=255)
                                                                                                                           ttCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGcc                   <  2:94885/150‑1 (MQ=255)
                                                                                                                            tCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACCCCGTTGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCg                  <  1:11406/150‑1 (MQ=255)
                                                                                                                              aaCGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGt                >  2:214284/1‑150 (MQ=255)
                                                                                                                               aCGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGAGCGGCGAAGCCATCCTTGAACGGCGTTTCACCaga                                                                                        >  2:1140349/1‑77 (MQ=255)
                                                                                                                               aCGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGAGCGGCGAAGCCATCCTTGAACGGCGTTTCACCaga                                                                                        <  1:1140349/77‑1 (MQ=255)
                                                                                                                                aagtgttgaCGCCAACGCCGATGGCGTCATCACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTAc              >  2:241222/7‑150 (MQ=255)
                                                                                                                                 gTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGtt                               <  1:557111/132‑1 (MQ=255)
                                                                                                                                 gTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGtt                               >  2:557111/1‑132 (MQ=255)
                                                                                                                                  ttGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCg            >  1:227135/1‑150 (MQ=255)
                                                                                                                                         cGCCATCGGCGTGGGCGTCAACCCTTGCCAGCGCGGCGAAGCCATCCTTGAACGGCGTTTCACCATAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCAt     <  2:117929/150‑1 (MQ=255)
                                                                                                                                            cATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTgccgcc                                                             <  2:600013/91‑1 (MQ=255)
                                                                                                                                            cATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTgccgcc                                                             >  1:600013/1‑91 (MQ=255)
                                                                                                                                            cATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCATGCa  >  1:805347/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCATGCA  >  NC_002947/5602090‑5602379

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: