Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,602,238 T→A 6.9% N432I (AAC→ATC)  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,2380TA 6.9% 206.8 / 3.7 69N432I (AAC→ATC) PP_4924subtilase family serine protease
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base A (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 = 1.00e+00

CTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCATGCACCAGCA  >  NC_002947/5602090‑5602385
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                                <  2:612060/94‑1 (MQ=255)
                                                                            tcgtGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAggg                                                                                                                                >  1:612060/1‑94 (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:923458/150‑1 (MQ=255)
                                                                                           aCCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCg                                                         <  1:541284/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                                                                                                   >  2:630728/1‑100 (MQ=255)
                                                                                                   cTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCa                                                                                                   <  1:630728/100‑1 (MQ=255)
                                                                                                    ttCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAgc                                                >  1:791007/1‑150 (MQ=255)
                                                                                                           gcgcCCAGACGGCGTCTTCAACGTTGATGACGCCCTCGGCGTTGGCGTCAACAATTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCCGAGCTGCCATCACTGCCTCCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGa                                         <  1:998870/150‑1 (MQ=255)
                                                                                                           gcgcCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGa                                         <  2:183833/150‑1 (MQ=255)
                                                                                                                 aGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCg                                                                                                        >  1:19980/1‑81 (MQ=255)
                                                                                                                 aGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCg                                                                                                        <  2:19980/81‑1 (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                                                                                              <  1:1140349/77‑1 (MQ=255)
                                                                                                                               aCGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGAGCGGCGAAGCCATCCTTGAACGGCGTTTCACCaga                                                                                              >  2:1140349/1‑77 (MQ=255)
                                                                                                                                aagtgttgaCGCCAACGCCGATGGCGTCATCACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTAc                    >  2:241222/7‑150 (MQ=255)
                                                                                                                                 gTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGtt                                     >  2:557111/1‑132 (MQ=255)
                                                                                                                                 gTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGtt                                     <  1:557111/132‑1 (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)
                                                                                                                                                  cGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCATGCACCAGCa  <  1:740951/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTCGACCTGGCCGAGCTTGATCTGGCTAATGCCATGCTCCTGCAATGCCCTTAGCTCACCGGCATCGGCACGCCCGTCGTGGTTGCCATCGACCCACACCTTCAATTGCGCCCAGACGGCGTCTTCAACGTTGATGACGCCATCGGCGTTGGCGTCAACACTTGCCAGGGCGGCGAAGCCATCCTTGAACGGCGTTTCACCAGAGCTGCCATCACTGCCTGCCTTGCCGCCGTAGTATTCGGAGAACAGCTGGCTGACGTTGCTCACCTGGCCGGTACCGTTGTCATGCACCAGCA  >  NC_002947/5602090‑5602385

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: