Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,768,094 T→A 10.2% M255L (ATG→TTG)  PP_4218 ← lipase/esterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,768,0940TA10.2% 312.4 / 22.2 108M255L (ATG→TTG) PP_4218lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base T (50/46);  new base A (4/7);  total (55/53)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-01

Only 100 of 113 total aligned reads displayed.

ATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACC  >  NC_002947/4767945‑4768242
                                                                                                                                                     |                                                                                                                                                    
aTCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACCCCCAGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAt                                                                                                                                                      >  1:1239687/1‑150 (MQ=255)
        cctcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCac                                                                                                                                              <  1:1780013/150‑2 (MQ=255)
          tcGATGATCGAACCCATACGCATGAAGTCGTGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAggg                                                                                                                                            >  1:601621/1‑149 (MQ=255)
             atgatCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCa                                                                                                                                         <  1:87859/150‑1 (MQ=255)
                    aaCCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATGc                                                                                                                                  <  2:548064/150‑1 (MQ=255)
                    aaCCCATACGCATGAAGTCATGGGTCACTCCCGGGACCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGc                                                                                                                                  <  1:2343198/150‑1 (MQ=255)
                     aCCCATACGCATGAAGTCATGGGTCACTCCCGGGCTCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCt                                                                                                                                 <  2:2364667/150‑1 (MQ=255)
                      cccATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTg                                                                                                                                <  1:47250/150‑1 (MQ=255)
                       ccATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCACACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTgg                                                                                                                               >  1:799346/1‑150 (MQ=255)
                         aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATGCTgatt                                                                                                                             <  1:1259618/150‑4 (MQ=255)
                              cATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCt                                                                                                                        >  2:1420927/1‑150 (MQ=255)
                               aTGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTg                                                                                                                       >  2:71541/1‑150 (MQ=255)
                                tGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGc                                                                                                                      >  2:2135414/1‑150 (MQ=255)
                                 gAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTg                                                                                                                                           <  2:2215220/128‑1 (MQ=255)
                                 gAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTg                                                                                                                                           >  1:2215220/1‑128 (MQ=255)
                                    gTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAg                                                                                                                  >  2:825822/1‑150 (MQ=255)
                                       aTGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCaa                                                                                                               <  1:1886025/150‑1 (MQ=255)
                                        tGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAAc                                                                                                              <  1:1055866/150‑1 (MQ=255)
                                        tGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAAc                                                                                                              <  2:881577/150‑1 (MQ=255)
                                         gggTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAt                                                                                                                                                      <  1:731380/109‑1 (MQ=255)
                                         gggTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAt                                                                                                                                                      >  2:731380/1‑109 (MQ=255)
                                           gTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGt                                                                                                           <  2:814428/150‑1 (MQ=255)
                                            tCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAt                                                                                                                                                      <  2:2417768/106‑1 (MQ=255)
                                            tCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAt                                                                                                                                                      >  1:2417768/1‑106 (MQ=255)
                                                  ccGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                    >  1:792034/1‑150 (MQ=255)
                                                  ccGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                    >  2:1293408/1‑150 (MQ=255)
                                                  ccGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                    >  1:1626508/1‑150 (MQ=255)
                                                    gggATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCAATATGGACACCCTGATCATGCAGCGTGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAcgc                                                                                                  <  2:1775209/150‑1 (MQ=255)
                                                      gATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCgg                                                                                                <  2:1626508/150‑1 (MQ=255)
                                                         cACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATc                                                                                             <  2:725904/150‑1 (MQ=255)
                                                          aCCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTTTTGCTGGCGACGCTGCCCAGCAACGGGGAAAAACGCggctag                                                                                            >  2:2343538/1‑146 (MQ=255)
                                                          aCCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGt                                                                                                           >  1:657089/1‑135 (MQ=255)
                                                           ccGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGt                                                                                           >  2:2424709/1‑150 (MQ=255)
                                                            cGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTg                                                                                          >  2:2161776/1‑150 (MQ=255)
                                                             gCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGa                                                                                         <  2:2219732/150‑1 (MQ=255)
                                                             gCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGa                                                                                         <  2:2440597/150‑1 (MQ=255)
                                                               aaGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACg                                                                                       >  2:1635014/1‑150 (MQ=255)
                                                                   tGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTgatg                                                                                                                                   >  1:1260605/1‑99 (MQ=255)
                                                                   tGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATg                                                                                                                                   <  2:1260605/102‑1 (MQ=255)
                                                                    gCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCt                                                                                                                                 >  2:638918/1‑103 (MQ=255)
                                                                    gCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCt                                                                                                                                 <  1:638918/103‑1 (MQ=255)
                                                                    gCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAcgcg                                                                                                 >  2:1683245/1‑135 (MQ=255)
                                                                    gCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAcgcg                                                                                                 <  1:1683245/135‑1 (MQ=255)
                                                                      aCCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGcc                                                                                <  1:1293408/150‑1 (MQ=255)
                                                                      aCCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATGc                                                                                                                                  >  1:1000431/1‑100 (MQ=38)
                                                                              cccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGc                                                                                 <  1:1688561/141‑1 (MQ=255)
                                                                              cccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGAGAAAAACGCGGATCGTGACGGTCTGc                                                                                 >  2:1688561/1‑141 (MQ=255)
                                                                               ccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTg                                                                       <  2:480385/150‑1 (MQ=255)
                                                                                 ggCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGc                                                                     <  2:1113452/150‑1 (MQ=255)
                                                                                   cTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCACCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGt                                                                   >  2:515491/1‑150 (MQ=255)
                                                                                      gCTCCAAGCGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTAtt                                                                <  1:18628/150‑1 (MQ=255)
                                                                                        tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGcc                                                                                <  1:748377/132‑1 (MQ=255)
                                                                                        tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGcc                                                                                >  2:748377/1‑132 (MQ=255)
                                                                                        tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCCCCGTGGCGTATTGc                                                              >  1:973782/1‑150 (MQ=255)
                                                                                        tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGc                                                                     >  1:952291/1‑143 (MQ=255)
                                                                                        tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGc                                                              >  1:752887/1‑150 (MQ=255)
                                                                                        tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGGGACGGTGTGCCGGCACCGTGGCGTATTGc                                                              >  2:29881/1‑150 (MQ=255)
                                                                                         ccAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCt                                                             <  1:1635014/150‑1 (MQ=255)
                                                                                          cAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTTGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTg                                                            <  2:601621/150‑1 (MQ=255)
                                                                                           aaGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                           <  1:2111390/150‑1 (MQ=255)
                                                                                           aaGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                           <  2:2124532/150‑1 (MQ=255)
                                                                                                ccGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCACACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGaa                                                      <  2:799346/150‑1 (MQ=255)
                                                                                                 cGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAAc                                                     <  1:1385218/150‑1 (MQ=255)
                                                                                                   tGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCa                                                   >  2:661561/1‑150 (MQ=255)
                                                                                                    gCATAAGCGACACCCTGGTGATGCAGGGGGTCGCACTCGGCAATCAGGATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAGAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCTTATTGCTGGTAGAACCAc                                                  <  2:826269/150‑1 (MQ=255)
                                                                                                      aTAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGatcggaag                                                >  1:1717338/1‑143 (MQ=255)
                                                                                                       tAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTcc                                               >  1:987359/1‑150 (MQ=255)
                                                                                                          gCGACACCCTGGTCATGCAGCGGGTCGCATTCGGCGCCCGGCTTCAGTGTATGTGCACTGTTGCTGGCGCAGTATTCTAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAgg                                            <  1:532504/150‑1 (MQ=255)
                                                                                                          gCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAgg                                            >  2:1125/1‑150 (MQ=255)
                                                                                                            gACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGt                                          >  1:1031466/1‑150 (MQ=255)
                                                                                                             acacCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGccc                                                                                                                    <  1:1797858/75‑1 (MQ=255)
                                                                                                             acacCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTc                                         <  1:1543961/150‑1 (MQ=255)
                                                                                                             acacCCTGGTCATGACGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGccc                                                                                                                    >  2:1797858/1‑75 (MQ=255)
                                                                                                               acCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTg                                                               <  2:902311/126‑1 (MQ=255)
                                                                                                               acCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTg                                                               >  1:902311/1‑126 (MQ=255)
                                                                                                                 ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                    >  1:2251268/1‑55 (MQ=255)
                                                                                                                 ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                    <  2:2251268/55‑1 (MQ=255)
                                                                                                                 ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                     <  2:752887/150‑1 (MQ=255)
                                                                                                                 ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                     >  2:385193/1‑150 (MQ=255)
                                                                                                                  cTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                           >  2:831866/1‑127 (MQ=255)
                                                                                                                  cTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                           <  1:831866/127‑1 (MQ=255)
                                                                                                                     gTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGGAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTc                                 <  1:2343538/150‑1 (MQ=255)
                                                                                                                     gTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGaa                                                      >  1:2109134/1‑129 (MQ=255)
                                                                                                                        aTGCAGCGGGACCCACCCGGCATTCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAg                              >  2:2112276/1‑150 (MQ=255)
                                                                                                                           cAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                           <  1:2161776/150‑1 (MQ=255)
                                                                                                                                    gCACTCGGCAATCAGCATCAGTGCAGGTGCGCTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCt                  >  2:77970/1‑150 (MQ=255)
                                                                                                                                        tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCa              >  2:890322/1‑150 (MQ=255)
                                                                                                                                        tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCa              >  2:1124438/1‑150 (MQ=255)
                                                                                                                                        tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGGGAAAAACGCGGAGCGTGACGGTCAGGAGGCACCTTGGCTTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAGCCGCTGCCa              >  1:1750394/1‑150 (MQ=255)
                                                                                                                                           gCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAg           >  2:520186/1‑150 (MQ=255)
                                                                                                                                           gCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAg           >  2:520878/1‑150 (MQ=255)
                                                                                                                                               tCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCt       >  2:1758436/1‑150 (MQ=255)
                                                                                                                                               ccAGCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGc                             <  1:2396069/127‑1 (MQ=255)
                                                                                                                                               ccAGCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGc                             >  2:2396069/2‑128 (MQ=255)
                                                                                                                                                 aGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATc                                                                                             <  1:1974471/62‑1 (MQ=255)
                                                                                                                                                 aGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATc                                                                                             >  2:1974471/1‑62 (MQ=255)
                                                                                                                                                 aGCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCaa                               <  1:1304754/124‑1 (MQ=255)
                                                                                                                                                 aGCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGAGGCGTATTGCTGGGAGAACCACTCCAGGGTCTGCGCTTCaa                               >  2:1304754/1‑124 (MQ=255)
                                                                                                                                                   cATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCAc   <  1:2135414/150‑1 (MQ=255)
                                                                                                                                                    aTCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCAcc  >  2:2233327/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACC  >  NC_002947/4767945‑4768242

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