BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,137,957	C→A	5.3%	A14A (GCC→GCA)	tolC →	agglutination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,137,957	0	C	A	5.3%	1124.7 / 34.9	332	A14A (GCC→GCA)	tolC	agglutination protein
Reads supporting (aligned to +/- strand): ref base C (152/162); new base A (16/2); total (168/164)
Fisher's exact test for biased strand distribution p-value = 9.93e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 346 total aligned reads displayed.

GAATGTACTTAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCGT  >  NC_002947/5137808‑5138106

     tACTTAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGc                                                                                                                                                  <  1:1914082/150‑1 (MQ=255)

      aCTTAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCg                                                                                                                                                 >  2:7233160/1‑150 (MQ=255)

       cTTAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGa                                                                                                                                                <  1:4642812/150‑1 (MQ=255)

        ttAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAg                                                                                                                                               >  2:5824962/1‑150 (MQ=255)

        ttAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTGGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAg                                                                                                                                               <  2:7421238/150‑1 (MQ=255)

          aGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTg                                                                                                                                             >  1:4063956/1‑150 (MQ=255)

             tgtgTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCa                                                                                                                                          <  2:864147/150‑1 (MQ=255)

                 ttgttAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGt                                                                                                                                      >  1:4213311/1‑150 (MQ=255)

                 ttgttAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGt                                                                                                                                      <  1:2970888/150‑1 (MQ=255)

                   gttAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTAc                                                                                                                                    <  1:7866924/150‑1 (MQ=255)

                    ttAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACa                                                                                                                                   >  1:3333102/1‑150 (MQ=255)

                    ttAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACa                                                                                                                                   >  1:6711417/1‑150 (MQ=255)

                      aaaaaCCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAgg                                                                                                                                 >  1:6302736/1‑150 (MQ=255)

                       aaaaCCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGc                                                                                                                                <  1:1451080/150‑1 (MQ=255)

                        aaaCCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCg                                                                                                                               <  1:2664990/150‑1 (MQ=255)

                          aCCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGAt                                                                                                                             <  2:5647193/150‑1 (MQ=255)

                            cTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGt                                                                                                                           <  2:2680604/150‑1 (MQ=255)

                             tGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTc                                                                                                                          >  2:609237/1‑150 (MQ=255)

                              gCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCg                                                                                                                         <  2:1224906/150‑1 (MQ=255)

                               cAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGa                                                                                                                        <  2:2307769/150‑1 (MQ=255)

                                aaCACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGAt                                                                                                                       <  2:3672235/150‑1 (MQ=255)

                                aaCACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGAt                                                                                                                       <  1:3110001/150‑1 (MQ=255)

                                aaCACTATACGCGATAAGTTTATCGTCA‑TTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGAtt                                                                                                                      >  2:113951/1‑150 (MQ=255)

                                    cTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTAcc                                                                                                                   >  2:8319380/1‑150 (MQ=255)

                                    cTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTAcc                                                                                                                   <  1:5298768/150‑1 (MQ=255)

                                      ataCGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGa                                                                                                                 <  2:1744574/150‑1 (MQ=255)

                                       taCGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAg                                                                                                                >  1:2458433/1‑150 (MQ=255)

                                        aCGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGAGTACCGAgg                                                                                                               >  2:8173471/1‑150 (MQ=255)

                                          gcgATCAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCt                                                                                                             >  1:2701137/1‑150 (MQ=255)

                                             aTAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTc                                                                                                          <  2:7207199/150‑1 (MQ=255)

                                                 gTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCaaa                                                                                                      <  2:6781838/150‑1 (MQ=255)

                                                 gTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCaaa                                                                                                      <  1:8711546/150‑1 (MQ=255)

                                                     aTCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAgcgc                                                                                                  <  1:5239848/150‑1 (MQ=255)

                                                        gTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGt                                                                                               <  2:5296570/150‑1 (MQ=255)

                                                        gTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGt                                                                                               <  1:7490221/150‑1 (MQ=255)

                                                        gTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGt                                                                                               >  1:497365/1‑150 (MQ=255)

                                                         tCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTg                                                                                              <  2:2592907/150‑1 (MQ=255)

                                                           aTTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGa                                                                                            <  2:8107742/150‑1 (MQ=255)

                                                            tttttttGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGAc                                                                                           >  2:4868487/1‑150 (MQ=255)

                                                             ttttttGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACt                                                                                          <  2:2958994/150‑1 (MQ=255)

                                                               ttttGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTac                                                                                        <  2:4007048/150‑1 (MQ=255)

                                                                 ttGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTacca                                                                                      <  2:4485642/150‑1 (MQ=255)

                                                                  tGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTaccac                                                                                     <  1:7266105/150‑1 (MQ=255)

                                                                   ggCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTaccacc                                                                                    <  1:1768256/150‑1 (MQ=255)

                                                                       aaCAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCcaca                                                                                >  2:2016593/1‑150 (MQ=255)

                                                                        aCAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACaa                                                                               <  1:5417209/150‑1 (MQ=255)

                                                                            cagGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCGGTCCAAAGCGCCGTGGACTACCACCCACAAGTCa                                                                           <  2:5190092/150‑1 (MQ=255)

                                                                              ggATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGc                                                                         <  2:2195944/150‑1 (MQ=255)

                                                                                aTCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTc                                                                       >  2:3648487/1‑150 (MQ=255)

                                                                                aTCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTc                                                                       >  1:7563404/1‑150 (MQ=255)

                                                                                 tCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTcc                                                                      >  1:6553928/1‑150 (MQ=255)

                                                                                   aaTTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGCCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCaa                                                                    <  2:6261448/150‑1 (MQ=255)

                                                                                    aTTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAAc                                                                   >  2:4763059/1‑150 (MQ=255)

                                                                                     ttCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAAcc                                                                  <  1:3565175/150‑1 (MQ=255)

                                                                                     ttCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAAcc                                                                  >  1:680588/1‑150 (MQ=255)

                                                                                         aGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTaa                                                              <  1:8186548/150‑1 (MQ=255)

                                                                                           gATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACa                                                            >  1:5932335/1‑150 (MQ=255)

                                                                                             tCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGc                                                          >  2:40579/1‑150 (MQ=255)

                                                                                                gggAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAg                                                       <  1:3292713/150‑1 (MQ=255)

                                                                                                gggAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAg                                                       <  1:6368166/150‑1 (MQ=255)

                                                                                                        tACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGc                                               >  2:7141017/1‑150 (MQ=255)

                                                                                                           cATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGa                                            >  1:5785760/1‑150 (MQ=255)

                                                                                                            aTGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGAt                                           >  1:2731153/1‑150 (MQ=255)

                                                                                                            aTGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGAt                                           <  1:6911005/150‑1 (MQ=255)

                                                                                                            aTGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGACGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGAt                                           <  2:4305334/150‑1 (MQ=255)

                                                                                                               cgcgTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGaa                                        <  1:3812113/150‑1 (MQ=255)

                                                                                                                gcgTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAg                                       >  1:2303438/1‑150 (MQ=255)

                                                                                                                  gTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGAc                                     >  1:7126112/1‑150 (MQ=255)

                                                                                                                  gTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGAc                                     <  1:5555108/150‑1 (MQ=255)

                                                                                                                   ttttGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACg                                    >  2:3109006/1‑150 (MQ=255)

                                                                                                                          ccccATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGtt                             >  1:6619248/1‑150 (MQ=255)

                                                                                                                           cccATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTc                            >  2:3868194/1‑150 (MQ=255)

                                                                                                                           cccATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTc                            >  2:7865920/1‑150 (MQ=255)

                                                                                                                           cccATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTc                            >  2:454679/1‑150 (MQ=255)

                                                                                                                           cccATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTc                            >  1:3473171/1‑150 (MQ=255)

                                                                                                                             cATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTcgc                          >  1:8383807/1‑150 (MQ=255)

                                                                                                                              aTCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTcgcg                         >  2:6609723/1‑150 (MQ=255)

                                                                                                                               tCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGATcgcgc                        >  2:7732006/1‑150 (MQ=255)

                                                                                                                                 aCCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGAGAAGTTcgcgcgg                      >  2:5454895/1‑149 (MQ=255)

                                                                                                                                    aGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCgtggt                   <  2:7601958/150‑1 (MQ=255)

                                                                                                                                    aGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCgtggt                   <  1:2159085/150‑1 (MQ=255)

                                                                                                                                    aGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCgtggt                   <  2:163770/150‑1 (MQ=255)

                                                                                                                                         tttACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGcta              <  2:6714955/150‑1 (MQ=255)

                                                                                                                                         tttACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGcta              >  1:1086930/1‑150 (MQ=255)

                                                                                                                                              tGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTAcc         >  2:5143221/1‑150 (MQ=255)

                                                                                                                                              tGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTAcc         >  1:808805/1‑150 (MQ=255)

                                                                                                                                              tGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTAcc         >  1:4810638/1‑150 (MQ=255)

                                                                                                                                              cgTTGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTAcc         >  2:7998751/2‑150 (MQ=255)

                                                                                                                                                ttGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCt       >  2:2867055/1‑150 (MQ=255)

                                                                                                                                                ttGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAAGAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCt       >  1:1247652/1‑150 (MQ=255)

                                                                                                                                                ttGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTCCTACCCt       >  2:7964217/1‑150 (MQ=255)

                                                                                                                                                ttGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCt       >  1:4333975/1‑150 (MQ=255)

                                                                                                                                                   gCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTcc    <  2:6487993/150‑1 (MQ=255)

                                                                                                                                                   gCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTcc    >  1:2365195/1‑150 (MQ=255)

                                                                                                                                                    cccTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCg   >  1:1978257/1‑150 (MQ=255)

                                                                                                                                                    cccTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCg   >  2:70183/1‑150 (MQ=255)

                                                                                                                                                    cccTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCg   <  2:239553/150‑1 (MQ=255)

                                                                                                                                                    cccTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCg   >  2:1256083/1‑150 (MQ=255)

                                                                                                                                                     ccTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCGt  <  1:1018578/150‑1 (MQ=255)

                                                                                                                                                     ccTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCGt  <  1:7198140/150‑1 (MQ=255)

GAATGTACTTAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAAGCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCGT  >  NC_002947/5137808‑5138106

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: `‑`