Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,067,254 G→C 5.8% intergenic (+14/‑160) phoR → / → PP_5322 two‑component system sensor histidine kinase/metal ion transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,067,2540GC 5.8% 269.4 / 6.4 87intergenic (+14/‑160)phoR/PP_5322two‑component system sensor histidine kinase/metal ion transporter
Reads supporting (aligned to +/- strand):  ref base G (44/38);  new base C (3/2);  total (47/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 = 6.67e-01

GGGGCTCGCAATCGTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTA  >  NC_002947/6067111‑6067402
                                                                                                                                                |                                                                                                                                                    
ggggCTCGCAATCGTCAAGCATGTGCTGATGCGCCACCAGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGc                                                                                                                                                <  2:1360599/150‑1 (MQ=255)
    cTCGCAATCGTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCtt                                                                                                                                            >  1:2074015/1‑150 (MQ=255)
        cAATCGTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTcgcg                                                                                                                                        <  1:176616/150‑1 (MQ=255)
         aaTCGTTAAGCATGTTCTTCTGCGCCACCGGGGCAAGCTTGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACGTGTTCCTTTGCGCGAGCACAAGTGGCAGACGAAAAGGCCTGAGATTGATCCTGCGGGGCAGCGCTTcactg                                                                                                                                       >  2:648777/1‑146 (MQ=255)
         aaTCGTCAAGCATGTGCTGATGCTCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATCGGGCTTCCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAACAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCgg                                                                                                                                       <  1:1955543/150‑1 (MQ=255)
         aaTCGTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCgg                                                                                                                                       >  2:350662/1‑150 (MQ=255)
         aaTCGTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCgg                                                                                                                                       >  2:1816774/1‑150 (MQ=255)
          aTCGTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCggg                                                                                                                                      <  2:2055240/150‑1 (MQ=255)
              tCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGGAATCAGCAGTGTACCGGGGCATGGGAGCACCTTTACCTGTCACTTTGCGCCAGCGCAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATg                                                                                                                                  >  1:203548/1‑150 (MQ=255)
                   cATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGc                                                                                                                             <  2:2268784/150‑1 (MQ=255)
                   cATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGc                                                                                                                             <  1:350662/150‑1 (MQ=255)
                       tgCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTGACATGGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGcc                                                                                                                         >  1:781137/1‑150 (MQ=255)
                           gATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTg                                                                                                                     <  2:2101590/150‑1 (MQ=255)
                            aTGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTgg                                                                                                                    >  1:2196552/1‑150 (MQ=255)
                              gcgcCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGc                                                                                                                  <  2:2196552/150‑1 (MQ=255)
                               cgcCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGcc                                                                                                                 <  1:2335176/150‑1 (MQ=255)
                                gcCACCGGGGCCCGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGccc                                                                                                                <  1:1100064/150‑1 (MQ=255)
                                gcCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTcgcg                                                                                                                                        >  1:1451765/1‑126 (MQ=255)
                                gcCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTcgcg                                                                                                                                        <  2:1451765/126‑1 (MQ=255)
                                   aCCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTgg                                                                                                                    >  1:626159/1‑143 (MQ=255)
                                   aCCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTgg                                                                                                                    <  2:626159/143‑1 (MQ=255)
                                       gggCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGa                                                                                                         <  2:1754300/150‑1 (MQ=255)
                                           aaGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTgg                                                                                                                    >  2:660012/1‑135 (MQ=255)
                                           aaGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTgg                                                                                                                    <  1:660012/135‑1 (MQ=255)
                                           aaGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTcc                                                                                                     >  2:2135283/1‑150 (MQ=255)
                                               tGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCACATCGCGGGCATGCCGGCATCCTCTTGGCCCGCAATGAGTCCGCTg                                                                                                 >  2:265252/1‑150 (MQ=255)
                                                       agcagTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGc                                                                                                                  >  2:1671377/1‑125 (MQ=255)
                                                       agcagTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGc                                                                                                                  <  1:1671377/125‑1 (MQ=255)
                                                       agcagTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGc                                                                                         <  1:265252/150‑1 (MQ=255)
                                                        gcagTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCgg                                                                                                                                       <  1:940045/103‑1 (MQ=255)
                                                        gcagTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCgg                                                                                                                                       >  2:940045/1‑103 (MQ=255)
                                                        gcagTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCt                                                                                        >  2:1155914/1‑150 (MQ=255)
                                                        gcagTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCt                                                                                        >  1:1698941/1‑150 (MQ=255)
                                                           gtgtACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTCAGCGCCGCGGGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCtttt                                                                                     >  1:1621307/1‑150 (MQ=255)
                                                           gtgtACCGGGGCATGGCAGCACCTTGACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCtttt                                                                                     <  2:640691/150‑1 (MQ=255)
                                                            tgtACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTa                                                                                    >  1:1494080/1‑150 (MQ=255)
                                                             gtACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTAc                                                                                   <  2:781137/150‑1 (MQ=255)
                                                                 cGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCAt                                                                                                                                   >  1:396383/1‑98 (MQ=255)
                                                                 cGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCAt                                                                                                                                   <  2:396383/98‑1 (MQ=255)
                                                                 cGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGcc                                                                               >  2:880232/1‑150 (MQ=255)
                                                                      cATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGc                                                                          >  2:2008368/1‑150 (MQ=255)
                                                                          gcagcaCCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCt                                                                      >  1:102462/1‑150 (MQ=255)
                                                                               aCCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACAtt                                                                 >  2:1623975/1‑150 (MQ=255)
                                                                                 cTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGcc                                                                         >  2:404295/1‑140 (MQ=255)
                                                                                 cTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGcc                                                                         <  1:404295/140‑1 (MQ=255)
                                                                                       cTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATcccc                                                         >  2:1027784/1‑150 (MQ=255)
                                                                                        tGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGc                                                                                                                  <  1:2340860/92‑1 (MQ=255)
                                                                                        tGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGc                                                                                                                  >  2:2340860/1‑92 (MQ=255)
                                                                                        tGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATGGGATCCCCt                                                        >  1:371251/1‑150 (MQ=255)
                                                                                         gTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTgg                                                               >  2:653084/1‑142 (MQ=255)
                                                                                         gTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTgg                                                               <  1:653084/142‑1 (MQ=255)
                                                                                            aCTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCCAAATTAAATcccaggtgt                                                    >  2:136472/1‑144 (MQ=255)
                                                                                              tttGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCa                                                  >  2:2060057/1‑150 (MQ=255)
                                                                                                    ccAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCaaaa                                            <  1:1995046/150‑1 (MQ=255)
                                                                                                    ccAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCaaaa                                            >  1:390751/1‑150 (MQ=255)
                                                                                                     cAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTa                                                                                    <  2:960937/109‑1 (MQ=255)
                                                                                                     cAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTa                                                                                    >  1:960937/1‑109 (MQ=255)
                                                                                                      aGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACAAACGGTTGCCTGTGCCAGCAAAAGc                                          >  1:125944/1‑150 (MQ=255)
                                                                                                      aGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCATTCGCGGGCATGCCGGCTGCCACTGGGCCAGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGc                                          <  1:1425845/150‑1 (MQ=255)
                                                                                                       gCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCt                                         <  2:1494080/150‑1 (MQ=255)
                                                                                                         ccagttGGCCGACGCAAAGGCTTGATATTGCGGCTGCGGCGCAACCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCCACGAGTCCGCTGCCCTTTGCTTCTCCGGCCCCAGCCGCTACATTGGCTCCCCTGTGCCAGCAAAAGCttg                                       <  2:203548/146‑3 (MQ=255)
                                                                                                         acaATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTgg                                       <  1:1589997/150‑1 (MQ=255)
                                                                                                         acaATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTgg                                       <  1:136472/150‑1 (MQ=255)
                                                                                                           aaTTGGCCAACGAAAAGGTTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCa                                     >  1:1903549/1‑150 (MQ=255)
                                                                                                               ggCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCCCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTCTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCTTTTGCTGGCAcagg                                 <  2:2074015/150‑4 (MQ=255)
                                                                                                                 ccAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCt                                                        >  2:1255353/1‑125 (MQ=255)
                                                                                                                 ccAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCt                                                        <  1:1255353/125‑1 (MQ=255)
                                                                                                                   aaCGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACttttttt                             >  1:173963/1‑150 (MQ=255)
                                                                                                                   aaCGAAAAGGCTGGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCTCTGAGCCCGCAATGAGGCCGTTGCCCTTTGCTTTTACGGCCCCAGTCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACttttttt                             <  1:1608930/150‑1 (MQ=255)
                                                                                                                        aaaGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGc                                                                                                                             >  2:739827/1‑49 (MQ=255)
                                                                                                                        aaaGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGc                                                                                                                             <  1:739827/49‑1 (MQ=255)
                                                                                                                         aaGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTtctct                         >  1:481280/1‑148 (MQ=255)
                                                                                                                         aaGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTtctct                         <  2:481280/148‑1 (MQ=255)
                                                                                                                           ggCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTAt                     <  2:173963/150‑1 (MQ=255)
                                                                                                                            gCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCt                                                                                        <  1:2081390/82‑1 (MQ=255)
                                                                                                                            gCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCt                                                                                        >  2:2081390/1‑82 (MQ=255)
                                                                                                                              ttGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCtt                                                                                       <  2:1735716/81‑1 (MQ=255)
                                                                                                                              ttGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCtt                                                                                       >  1:1735716/1‑81 (MQ=255)
                                                                                                                              ttGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTAtctc                  >  2:1365396/1‑150 (MQ=255)
                                                                                                                              ttGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTAtctc                  >  1:1652643/1‑150 (MQ=255)
                                                                                                                               tGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCa                 >  2:2210212/1‑150 (MQ=255)
                                                                                                                                 agaTTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCaaa               >  2:930934/1‑150 (MQ=255)
                                                                                                                                      agggCTGCGCCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACAtt                                                                 >  2:5212/3‑95 (MQ=255)
                                                                                                                                      agggCTGCGCCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACAtt                                                                 <  1:5212/93‑1 (MQ=255)
                                                                                                                                       gggCTGCGCCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGa         >  1:1896675/2‑150 (MQ=255)
                                                                                                                                        ggCTGCGCCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCaa                                                                                                            <  1:1753711/50‑1 (MQ=255)
                                                                                                                                        ggCTGCGCCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCaa                                                                                                            >  2:1753711/1‑50 (MQ=255)
                                                                                                                                             gcggcgCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCg    >  2:1639560/1‑150 (MQ=255)
                                                                                                                                               ggcgAAACCCCTCGCGGGCATCCCGGCTGCCACAGGGCACGCAATGAGAGCGCTGCCCTTTGCTGTTATGGCCCCAGACGCCACATTGGATCCCCTGTGCCAGCACCAGCGGGCACTTTTTTTCCCTCTATCTCAAACACGGGCCTCGTa  >  2:888025/1‑150 (MQ=255)
                                                                                                                                                gcgcAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTgg                                       >  1:2273875/1‑112 (MQ=255)
                                                                                                                                                gcgcAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTgg                                       <  2:2273875/112‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGGGCTCGCAATCGTCAAGCATGTGCTGATGCGCCACCGGGGCAAGCTGGAAATCAGCAGTGTACCGGGGCATGGCAGCACCTTTACCTGTCACTTTGCGCCAGCACAATTGGCCAACGAAAAGGCTTGAGATTGAGGCTGCGGCGCAGCCCTTCGCGGGCATGCCGGCTGCCACTGGGCCCGCAATGAGTCCGCTGCCCTTTGCTTTTACGGCCCCAGCCGCTACATTGGATCCCCTGTGCCAGCAAAAGCTGGCACTTTTTTTCTCTCTATCTCAAAGACGGACCCCGTA  >  NC_002947/6067111‑6067402

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: