Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,476,108 T→C 11.1% intergenic (‑48/+132) PP_2167 ← / ← tal PhoD‑like phosphatase/transaldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,476,1080TC11.1% 245.3 / 17.2 94intergenic (‑48/+132)PP_2167/talPhoD‑like phosphatase/transaldolase
Reads supporting (aligned to +/- strand):  ref base T (11/66);  new base C (5/5);  total (23/71)
Fisher's exact test for biased strand distribution p-value = 1.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CAGCGGTTGGGTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTAT  >  NC_002947/2475959‑2476257
                                                                                                                                                     |                                                                                                                                                     
cAGCGGTTGGGTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCt                                                                                                                                                       <  1:4363509/150‑1 (MQ=255)
  gCGGTTGGGTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCTGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTa                                                                                                                                                     <  2:4405384/150‑1 (MQ=255)
  gCGGTTGGGTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTa                                                                                                                                                     <  2:1926874/150‑1 (MQ=255)
  gCGGTTGGGTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTa                                                                                                                                                     <  2:1470099/150‑1 (MQ=255)
    ggTTGGGTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACa                                                                                                                                                   <  2:1881960/150‑1 (MQ=255)
     gTTGGGTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACAt                                                                                                                                                  <  2:3907823/150‑1 (MQ=255)
         ggTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGAt                                                                                                                                              <  1:1110495/150‑1 (MQ=255)
             gTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAgt                                                                                                                                          <  2:2440407/150‑1 (MQ=255)
               gACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAgtgt                                                                                                                                        <  1:551179/150‑1 (MQ=255)
                aCGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTa                                                                                                                                       <  1:2416784/150‑1 (MQ=255)
                 cGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAg                                                                                                                                      <  1:699290/150‑1 (MQ=255)
                 cGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAg                                                                                                                                      <  1:3397278/150‑1 (MQ=255)
                  gAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCGTAACTCGGTCAGTGGAgg                                                                                                                                     >  2:871924/1‑150 (MQ=255)
                   aGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGa                                                                                                                                    <  2:2422140/150‑1 (MQ=255)
                   aGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGa                                                                                                                                    <  2:2422306/150‑1 (MQ=255)
                   aGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGa                                                                                                                                    <  2:3855381/150‑1 (MQ=255)
                     ccagacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGc                                                                                                                                  <  1:3471020/150‑1 (MQ=255)
                      cagacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGTGGAGGTCATCGACTGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCg                                                                                                                                 <  1:2276161/150‑1 (MQ=37)
                      cagacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCg                                                                                                                                 <  1:3538906/150‑1 (MQ=37)
                      cagacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCg                                                                                                                                 <  1:2276193/150‑1 (MQ=37)
                       agacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGCGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCgg                                                                                                                                <  2:2819806/150‑1 (MQ=37)
                       agacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCgg                                                                                                                                <  2:394573/150‑1 (MQ=37)
                        gacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTCCATCGATCAGTGTAGGAGCggg                                                                                                                               <  1:1742803/150‑1 (MQ=255)
                        gacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCggg                                                                                                                               <  2:3688156/150‑1 (MQ=255)
                         acggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGc                                                                                                                              <  1:3557605/150‑1 (MQ=255)
                           ggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTGGCGGGCCAGCCCCCCCCTACAAAGAGCAGGGTAGgggcgggcgt                                                                                                                            >  1:1391376/1‑148 (MQ=255)
                            gccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTg                                                                                                                           <  2:2501422/150‑1 (MQ=255)
                             ccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGc                                                                                                                          <  1:3227991/150‑1 (MQ=255)
                             ccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGc                                                                                                                          <  1:1747255/150‑1 (MQ=255)
                              cagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGcc                                                                                                                         <  1:4205187/150‑1 (MQ=255)
                              cagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGcc                                                                                                                         <  1:1936827/150‑1 (MQ=255)
                              cagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCAGTTAGCGGGCAAGCCCGCCCGTACATCGATCAGTGTAGGAGCGGGCTTGcc                                                                                                                         >  2:1203150/1‑150 (MQ=255)
                                     tGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAACCCCCCCCGTACATCGATCCGGGTACGCGCGGGCGTGCCCGCGAAc                                                                                                                  >  2:4013747/1‑150 (MQ=255)
                                       gggTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTGGCGGGCCACCCCCCCCGGACAAAGATAAGGGGGGGCGCCGGCGTGCCcgccacccc                                                                                                                >  1:2868546/1‑144 (MQ=255)
                                           tCCTGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGCCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGccc                                                                                                            <  2:4161202/150‑1 (MQ=255)
                                             cAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCga                                                                                                          <  2:2809011/150‑1 (MQ=255)
                                             cAGGCGGCGCAGGACCGGGCCGGCGAGAACGAGGGGCCACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCga                                                                                                          <  2:1473695/150‑1 (MQ=255)
                                              aGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACAACGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCgag                                                                                                         <  2:309272/150‑1 (MQ=255)
                                                gcggcgCAGTACCGGGCCGGCGAGAACGTGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGaa                                                                                                       <  2:841675/150‑1 (MQ=255)
                                                    cgcAGTACCGGGCCGGCGAGAACGAGGGGCACCGGGTTTGGGGTGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGGTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTTGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCa                                                                                                   <  1:2878322/150‑1 (MQ=39)
                                                     gcAGTACCGGGCCGGCGAGAACGAGGGTCAAGGGGTTTGGGGAGGTCATCGACAGCTGGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCaa                                                                                                  <  1:710030/150‑1 (MQ=255)
                                                     gcAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCaa                                                                                                  <  1:685913/150‑1 (MQ=255)
                                                       aGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGCTCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATa                                                                                                <  2:4409771/150‑1 (MQ=255)
                                                          aCCGGGCCGGCGAGAACGAGGGGCAACGGGTTCGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGTGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATacac                                                                                             <  1:4016644/150‑1 (MQ=255)
                                                                     gagaCCGAGGTGCAAGGGGTTAGGGGTGGTCTTCGACAGGTCGGGTGGCTAGGGCGGCGTGTTGGCGGGCACGCCGGCCCTTTCTTCTATGAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCTATACAGAACTCCGGCGt                                                                                  <  2:1246179/150‑1 (MQ=255)
                                                                     gagaACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCCACCCCGCCCGTACATCGCTCAGTGGAGGAGCGGGCGTGCCCCCCCACACGCCCGCGCAGTCAAAACAACACGAAGGCgg                                                                                  >  2:2431659/1‑149 (MQ=255)
                                                                     gagaACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCAGTTAGCGGGGAAGCCCCCCCATAAAGAGCTAAGTGTATGAGCGGGCTTGCCCCCCGACACGCCCGAGCAATCAAAAAACAACAAAGGGGt                                                                                  >  2:3784193/1‑150 (MQ=255)
                                                                         aCGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTGGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGca                                                                              <  1:647592/150‑1 (MQ=255)
                                                                           gtggggCACCGGGTTTGGGGAGTTCATCCACAGCTCGGGTGGATATGGCGGCCTGTTGGCGGGCACGCCCGCCCTAACACCGACCCGTGTAGGCGCGGGCTGGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGcaaa                                                                            <  1:1234499/148‑1 (MQ=255)
                                                                           gcggggCACCGGGTTAGGGGAGGTCATCGACCGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCCGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGcaaa                                                                            <  1:4350401/148‑1 (MQ=255)
                                                                           gAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGcaaa                                                                            <  2:4146293/150‑1 (MQ=255)
                                                                                 cAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAAc                                                                      <  2:1884800/150‑1 (MQ=255)
                                                                                 cAACCGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAACGCGTGGCAAACAAAAc                                                                      <  2:3707043/150‑1 (MQ=255)
                                                                                  aaCGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCGGTTCGCGGGCCACCCCCCCCATACAACGATAAGGGTATGAGCGGGCGTGCCCCCCCACACGCCCCACCAATAAAAACACCAACAAGGCCGGGCACACAAAAcc                                                                     >  2:783455/1‑149 (MQ=255)
                                                                                    cGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACggg                                                                   <  2:2318014/150‑1 (MQ=255)
                                                                                     ggtttAGGGGTGGTCTTCTACAGCTCGGGTGGCTATGGGGGCATGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGTGGGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGc                                                                  <  1:3432443/147‑1 (MQ=38)
                                                                                     gggTTTGGGGAGGTCATCGACAGCTCGGGGGGCTAGGGCGGCCTGTTCGCGGGAAAGCCCGCCCTTACATCTATCCGTGTAGGCGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGc                                                                  <  1:3432441/150‑1 (MQ=38)
                                                                                     gggTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGc                                                                  <  2:2442103/150‑1 (MQ=255)
                                                                                             ggcggTCATCGACAGCTCGGGTGGCTAGGGCGGCCTGTTCGGGGGAAAGCCCGCCCTTACATCTTTCCGTGTTGGAGCGGGTTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATa                                                          <  1:3387872/147‑1 (MQ=255)
                                                                                             ggAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCTATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATa                                                          <  1:2144965/150‑1 (MQ=255)
                                                                                              gAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATAc                                                         <  1:1098676/150‑1 (MQ=255)
                                                                                              gAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATAc                                                         <  1:3378350/150‑1 (MQ=255)
                                                                                              gAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATAc                                                         <  1:3378451/150‑1 (MQ=255)
                                                                                               aGGTCATCTACTGGTGGGGTGGGTTTGGGGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCCGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACt                                                        <  1:3419565/150‑1 (MQ=38)
                                                                                               aGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACt                                                        <  2:4378119/150‑1 (MQ=255)
                                                                                                ggTCATCGACTGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCATTACATCGATCAGTGTTGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTg                                                       <  1:3486829/150‑1 (MQ=255)
                                                                                                ggTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTg                                                       <  2:3467057/150‑1 (MQ=255)
                                                                                                 gTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGa                                                      <  2:684269/150‑1 (MQ=255)
                                                                                                  tCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGaa                                                     <  2:2332162/150‑1 (MQ=255)
                                                                                                  tCATCGACAGCTCGGGTGGCTAGGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCCGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGaa                                                     <  2:2336696/150‑1 (MQ=255)
                                                                                                        aCAGCTCGGGTGGCTATGCCGGCCTGTTCGCGGCCAAGCCCCCCATAACAACCATCAGTCTAGGAGCGGGCTTCCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTc                                               <  1:48668/150‑1 (MQ=255)
                                                                                                         cAGCTCGGGGGGCTATGGCGGCCGGTTCGCGGGCAAGCCCGCCCTTACAACGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTcc                                              <  1:3044427/150‑1 (MQ=255)
                                                                                                          aGCTCGGGTGGCTATGGCGGCCAGTTGGCGGGCCAGCCCCCCCTTACCTCGCTCAGGTTAGGGGCGGGCGTGCCCCCGACCAGGGCCCCGCATTAAATACACAACAAAAGCGGGGCGAACAAAAACGGCGGGAAATACAGAACAATTacc                                             >  2:3524935/1‑150 (MQ=255)
                                                                                                               gggTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTcccgcccg                                        <  1:2892762/150‑1 (MQ=255)
                                                                                                                        tGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCaaa                               <  1:2581988/150‑1 (MQ=255)
                                                                                                                           cggcCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGGGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAacac                            >  2:642880/1‑150 (MQ=255)
                                                                                                                                tGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACt                       >  2:3742773/1‑150 (MQ=255)
                                                                                                                                  ttCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGt                     <  1:2766123/150‑1 (MQ=255)
                                                                                                                                   tCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGc                                                                               <  1:147725/91‑1 (MQ=255)
                                                                                                                                   tCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGc                                                                               >  2:147725/1‑91 (MQ=255)
                                                                                                                                       gggCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAg                <  2:3977652/150‑1 (MQ=255)
                                                                                                                                        ggCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCTCCCGTCAGGCAAACACGAACTGTATCAgg               <  2:1749355/150‑1 (MQ=255)
                                                                                                                                         gCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGc              >  1:209090/1‑150 (MQ=255)
                                                                                                                                         gCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGc              >  2:180536/1‑150 (MQ=255)
                                                                                                                                           aaGCCCGCTCCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCt            >  2:172125/1‑150 (MQ=37)
                                                                                                                                           aaGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCt            >  2:2641637/1‑150 (MQ=255)
                                                                                                                                            aGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtt           >  1:1017102/1‑150 (MQ=255)
                                                                                                                                            aGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtt           >  2:895570/1‑150 (MQ=255)
                                                                                                                                            aGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtt           >  2:3805301/1‑150 (MQ=255)
                                                                                                                                             gcccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCttt          >  1:3915875/1‑150 (MQ=37)
                                                                                                                                              cccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtttt         >  2:990166/1‑150 (MQ=37)
                                                                                                                                              cccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtttt         >  1:3889975/1‑150 (MQ=37)
                                                                                                                                                   ccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCt    <  2:764894/150‑1 (MQ=255)
                                                                                                                                                    cTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTa   >  1:3176955/1‑150 (MQ=255)
                                                                                                                                                    cTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAACGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTa   <  2:1358859/150‑1 (MQ=255)
                                                                                                                                                     ttACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTAt  >  1:1456483/1‑150 (MQ=255)
                                                                                                                                                     ttACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTAt  >  2:1109562/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCGGTTGGGTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTAT  >  NC_002947/2475959‑2476257

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