Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,050,221 G→A 5.9% intergenic (+42/+64) PP_0905 → / ← PP_0906 hypothetical protein/RND family multidrug transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,050,2210GA 5.9% 1085.3 / 37.1 340intergenic (+42/+64)PP_0905/PP_0906hypothetical protein/RND family multidrug transporter
Reads supporting (aligned to +/- strand):  ref base G (138/182);  new base A (14/6);  total (152/188)
Fisher's exact test for biased strand distribution p-value = 2.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 357 total aligned reads displayed.

GCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGTGC  >  NC_002947/1050072‑1050370
                                                                                                                                                     |                                                                                                                                                     
gcgTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCgg                                                                                                                                                       <  1:1039194/150‑1 (MQ=255)
   ttGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGc                                                                                                                                                    <  1:2376395/150‑1 (MQ=255)
     gCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAAAGGCCCCGCAATCGCGAGGCCt                                                                                                                                                  <  2:2859932/150‑1 (MQ=255)
       cgcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCttt                                                                                                                                                <  2:3398018/150‑1 (MQ=255)
       cgcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAAAGGCCCCGCAATCGCGAGGCCtct                                                                                                                                                <  2:5228493/150‑3 (MQ=255)
        gcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTg                                                                                                                                               >  1:2519777/1‑150 (MQ=255)
           ggATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtt                                                                                                                                            <  2:425306/150‑1 (MQ=255)
           ggATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtt                                                                                                                                            <  2:4218723/150‑1 (MQ=255)
            gATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGttt                                                                                                                                           >  2:696008/1‑150 (MQ=255)
              tGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGttttt                                                                                                                                         <  2:5346211/150‑1 (MQ=255)
               gCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtttttt                                                                                                                                        <  1:6714915/150‑1 (MQ=255)
               gCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtttttt                                                                                                                                        <  1:2082739/150‑1 (MQ=255)
                 aaaaCCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGc                                                                                                                                      <  2:6210508/150‑1 (MQ=255)
                 aaaaCCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGc                                                                                                                                      >  1:6000413/1‑150 (MQ=255)
                     ccGCAGCGGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCt                                                                                                                                  <  2:6879255/150‑1 (MQ=255)
                     ccGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCt                                                                                                                                  >  2:4020047/1‑150 (MQ=255)
                     ccGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCt                                                                                                                                  <  2:201045/150‑1 (MQ=255)
                          gcagcagGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCtt                                                                                                                             <  2:5965349/150‑1 (MQ=255)
                          gcagcagGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCtt                                                                                                                             >  1:4201869/1‑150 (MQ=255)
                            agcagGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGc                                                                                                                           <  2:2129247/150‑1 (MQ=255)
                               agGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCgg                                                                                                                        >  1:1516743/1‑150 (MQ=255)
                                 gCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGcc                                                                                                                      >  1:878976/1‑150 (MQ=255)
                                        gcgagcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTcgc                                                                                                               <  2:5705980/150‑1 (MQ=255)
                                         cgagcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTcgcg                                                                                                              <  2:4441811/150‑1 (MQ=255)
                                           agcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCggg                                                                                                            <  2:2709258/150‑1 (MQ=255)
                                             ccAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCa                                                                                                          >  2:6101928/1‑150 (MQ=255)
                                             ccAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCa                                                                                                          >  1:555734/1‑150 (MQ=255)
                                             ccAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCa                                                                                                          <  1:1665444/150‑1 (MQ=255)
                                              cAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAg                                                                                                         >  1:6608827/1‑150 (MQ=255)
                                              cAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAg                                                                                                         >  1:2463069/1‑150 (MQ=255)
                                                 aCCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgccc                                                                                                      >  2:6049939/1‑150 (MQ=255)
                                                    aGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccgcc                                                                                                   <  1:364715/150‑1 (MQ=255)
                                                     ggACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccgccc                                                                                                  >  2:31336/1‑150 (MQ=255)
                                                       aCCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCa                                                                                                <  2:4358590/150‑1 (MQ=255)
                                                        ccTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCaa                                                                                               <  2:1944140/150‑1 (MQ=255)
                                                        ccTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCaa                                                                                               <  1:7312092/150‑1 (MQ=255)
                                                         cTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAg                                                                                              <  2:6950722/150‑1 (MQ=255)
                                                          ttttGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAgg                                                                                             >  1:7043877/1‑150 (MQ=255)
                                                              gCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCgg                                                                                         <  1:6892953/150‑1 (MQ=255)
                                                               cTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGt                                                                                        <  1:5476030/150‑1 (MQ=255)
                                                                   tCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCaa                                                                                    <  1:2817432/150‑1 (MQ=255)
                                                                    ccTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAg                                                                                   <  2:6280277/150‑1 (MQ=255)
                                                                      tGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCg                                                                                 >  1:1091439/1‑150 (MQ=255)
                                                                       gTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCgg                                                                                         <  1:4145198/141‑1 (MQ=255)
                                                                       gTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCgg                                                                                         >  2:4145198/1‑141 (MQ=255)
                                                                       gTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGa                                                                                <  1:2703610/150‑1 (MQ=255)
                                                                          tGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGAGTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTg                                                                             <  1:5277588/150‑1 (MQ=255)
                                                                              cTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAgg                                                                         >  1:3527501/1‑150 (MQ=255)
                                                                                gCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccg                                                                                                     <  1:4794098/120‑1 (MQ=255)
                                                                                gCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccg                                                                                                     >  2:4794098/1‑120 (MQ=255)
                                                                                   tCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCt                                                                    >  2:1273814/1‑150 (MQ=255)
                                                                                    cAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCtt                                                                   >  1:4196426/1‑150 (MQ=255)
                                                                                    cAAGGGTTGAGGTCTCAGTGATGAACACTGAGTTAGCTGTAAAAGCCAGAGGCCTCGCGATTGCGGTGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCtt                                                                   <  2:4599389/150‑1 (MQ=255)
                                                                                      aGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGAGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTgg                                                                 >  2:2386360/1‑150 (MQ=255)
                                                                                           tGCGGTCCCAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTc                                                            <  1:6123445/150‑1 (MQ=255)
                                                                                           tGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTc                                                            <  2:4676204/150‑1 (MQ=255)
                                                                                             cGGTCACAGGGATGAACCCTGCGTTAGCTGGAATAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGc                                                          <  2:2998245/150‑1 (MQ=255)
                                                                                             cGGTCACAGGGATGAACCCTGCGTGAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGc                                                          <  2:7089718/150‑1 (MQ=255)
                                                                                               gTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCt                                                        <  2:2068878/150‑1 (MQ=255)
                                                                                               gTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCt                                                        >  1:3242472/1‑150 (MQ=255)
                                                                                                tCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCtt                                                       >  2:5857495/1‑150 (MQ=255)
                                                                                                  acaGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTtgt                                                     <  1:216377/150‑1 (MQ=255)
                                                                                                   caGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTtgtg                                                    >  1:3148621/1‑150 (MQ=255)
                                                                                                        aTGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCACCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCca                                               <  1:6438992/150‑1 (MQ=255)
                                                                                                          gAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCcagc                                             >  2:3565057/1‑150 (MQ=255)
                                                                                                          gAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCcagc                                             >  2:21947/1‑150 (MQ=255)
                                                                                                             cccTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAgg                                          <  2:4855472/150‑1 (MQ=255)
                                                                                                             cccTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAgg                                          <  2:5189035/150‑1 (MQ=255)
                                                                                                              ccTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAggg                                         >  2:5499061/1‑150 (MQ=255)
                                                                                                                tGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgtg                                       >  1:4299480/1‑150 (MQ=255)
                                                                                                                tGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCGTTGGTTTTGTGCGTCTGCCGTGCCGGCCTCTTCGCGGGCAGCCCCCCCCAAGGGCGGTTCAAGCGAGGGCAGGTGCCGAGGCGCTCGCCGGGTGCGCCCGCCggggg                                       >  2:5595853/1‑148 (MQ=255)
                                                                                                                  cGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTa                                     <  1:5093763/150‑1 (MQ=255)
                                                                                                                    ttAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGa                                   >  1:5220260/1‑150 (MQ=255)
                                                                                                                         tGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCAc                              <  2:1070381/150‑1 (MQ=255)
                                                                                                                         tGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCAc                              <  1:2131763/150‑1 (MQ=255)
                                                                                                                         tGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCAc                              <  1:4706238/150‑1 (MQ=255)
                                                                                                                           gAAAAGCCAGAGGCCTCGTGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgt                                        >  2:808730/1‑138 (MQ=255)
                                                                                                                           gAAAAGCCAGAGGCCTCGTGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgt                                        <  1:808730/138‑1 (MQ=255)
                                                                                                                           gAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACgg                            <  1:138140/150‑1 (MQ=255)
                                                                                                                              aaGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCaa                         <  1:7357986/150‑1 (MQ=255)
                                                                                                                              aaGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCaa                         >  2:6497874/1‑150 (MQ=255)
                                                                                                                               aGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCCAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAac                        >  2:5279837/1‑150 (MQ=255)
                                                                                                                                gCCAGAGGCCTCGCGCTGGCGGCGCCTTAGGCATTTTGCGTCGGCCGTGCCGGCCTCGCCGCGGGCAGCCCGCCCCAACGTCGGTTCAAGCGATCGCAGGCGCCTCGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAaca                       <  2:1515734/150‑1 (MQ=255)
                                                                                                                                gCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAaca                       <  2:7316163/150‑1 (MQ=255)
                                                                                                                                gCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAaca                       <  1:3224950/150‑1 (MQ=255)
                                                                                                                                 ccAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAacac                      <  1:5903409/150‑1 (MQ=255)
                                                                                                                                  caaagGCCCCGCAATCGCGAGGCCTCTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAAcaca                     >  1:4416594/4‑150 (MQ=255)
                                                                                                                                   agagGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCcag                                              >  2:1870780/1‑124 (MQ=255)
                                                                                                                                   agagGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCcag                                              <  1:1870780/124‑1 (MQ=255)
                                                                                                                                   agagGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaa                    >  2:3701877/1‑150 (MQ=255)
                                                                                                                                    aagGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaaa                   >  1:2249225/2‑150 (MQ=255)
                                                                                                                                       gCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaaacaa                >  1:7172261/1‑150 (MQ=255)
                                                                                                                                       gCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaaacaa                <  2:6603869/150‑1 (MQ=255)
                                                                                                                                        ccTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAg               <  1:3902004/150‑1 (MQ=255)
                                                                                                                                          tCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGt             >  1:4467003/1‑150 (MQ=255)
                                                                                                                                          tCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGt             >  2:3232060/1‑150 (MQ=255)
                                                                                                                                             cgATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTaaa          <  1:2638402/150‑1 (MQ=255)
                                                                                                                                                   cGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaa                    >  1:3070786/1‑134 (MQ=255)
                                                                                                                                                   cGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaa                    <  2:3070613/134‑1 (MQ=255)
                                                                                                                                                    ggggCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGTg   <  1:347670/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTAAACAACGTGC  >  NC_002947/1050072‑1050370

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