BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,782	G→A	9.8%	D251D (GAC→GAT)	dnaK ←	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,378,782	0	G	A	9.8%	703.8 / 49.3	245	D251D (GAC→GAT)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand): ref base G (127/94); new base A (16/8); total (143/102)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 247 total aligned reads displayed.

GAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGA  >  NC_002947/5378633‑5378928

gagaTCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAgg                                                                                                                                                    >  2:6229905/1‑150 (MQ=255)

gagaTCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAgg                                                                                                                                                    >  2:1629562/1‑150 (MQ=255)

   aTCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTcc                                                                                                                                                 <  2:5715918/150‑1 (MQ=255)

      ttCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATg                                                                                                                                              <  1:5181214/150‑1 (MQ=255)

       tCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGc                                                                                                                                             >  2:6568592/1‑150 (MQ=255)

       tCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGc                                                                                                                                             >  2:6320263/1‑150 (MQ=255)

        cACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGcc                                                                                                                                            >  2:4262028/1‑150 (MQ=255)

         aCGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGccc                                                                                                                                           <  2:2521906/150‑1 (MQ=255)

           gTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGa                                                                                                                                         >  1:1762412/1‑150 (MQ=255)

           gTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGa                                                                                                                                         <  1:2471261/150‑1 (MQ=255)

             caaGGTGCTTAGGACCGGTCTGGTCAGCCGTGATGAATGGCAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGccctgtt                                                                                                                                       <  1:1601977/148‑5 (MQ=255)

                  tGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGTTTGACGTCGGTCGACTGAGCGGACGACATCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttc                                                                                                                                  <  1:767436/150‑1 (MQ=255)

                  tGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttc                                                                                                                                  <  1:6921719/150‑1 (MQ=255)

                  tGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttc                                                                                                                                  <  2:1148816/150‑1 (MQ=255)

                   gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                                                                                 >  1:1095387/1‑150 (MQ=255)

                   gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                                                                                 <  2:4615936/150‑1 (MQ=255)

                   gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                                                                                 >  2:8610004/1‑150 (MQ=255)

                    cTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttctt                                                                                                                                >  1:2846382/1‑150 (MQ=255)

                    cTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttctt                                                                                                                                <  1:3938826/150‑1 (MQ=255)

                     ttAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTg                                                                                                                               >  2:2741783/1‑150 (MQ=255)

                        ggACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAAc                                                                                                                            >  1:556197/1‑150 (MQ=255)

                         gACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCTGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACt                                                                                                                           <  1:3962401/150‑1 (MQ=255)

                         gACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTg                                                                                                                               >  2:5143285/1‑146 (MQ=255)

                           ccGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcg                                                                                                                         >  2:1437794/1‑150 (MQ=255)

                            cGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgt                                                                                                                        >  2:7189870/1‑150 (MQ=255)

                            cGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgt                                                                                                                        <  1:1057863/150‑1 (MQ=255)

                             ggTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgtc                                                                                                                       <  2:2115941/150‑1 (MQ=255)

                             ggTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgtc                                                                                                                       <  2:777677/150‑1 (MQ=255)

                              gTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgtcg                                                                                                                      >  2:4944764/1‑150 (MQ=255)

                                  cATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACaa                                                                                                                  <  1:6769282/150‑1 (MQ=255)

                                   aTCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAg                                                                                                                 >  1:6319453/1‑150 (MQ=255)

                                    tCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAgg                                                                                                                >  1:3614139/1‑150 (MQ=255)

                                      tGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTa                                                                                                              <  1:6787380/150‑1 (MQ=255)

                                       gCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAg                                                                                                             >  2:2535195/1‑150 (MQ=255)

                                         aGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTc                                                                                                           <  1:895370/150‑1 (MQ=255)

                                         aGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTc                                                                                                           >  2:8349439/1‑150 (MQ=255)

                                             aTGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATc                                                                                                       <  2:379293/150‑1 (MQ=255)

                                                tACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAgg                                                                                                    <  2:8272122/150‑1 (MQ=255)

                                                tACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTAGAACTCGTCGACAAGGTAGTCGATCAgg                                                                                                    >  2:1681487/1‑150 (MQ=255)

                                                 aCGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGgc                                                                                                   >  1:2574906/1‑150 (MQ=255)

                                                 aCGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCATGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTAGAACTCGTCGACAAGGTAGACGATCAGgc                                                                                                   >  2:3877363/1‑150 (MQ=255)

                                                  cGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGgcg                                                                                                  <  1:9059641/150‑1 (MQ=255)

                                                   ggCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGgcgc                                                                                                 <  1:2261461/150‑1 (MQ=255)

                                                     cAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCAt                                                                                               <  2:8757676/150‑1 (MQ=255)

                                                     cAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCAt                                                                                               <  1:740070/150‑1 (MQ=255)

                                                          tGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGa                                                                                          >  1:7863444/1‑150 (MQ=255)

                                                          tGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGa                                                                                          >  2:9036881/1‑150 (MQ=255)

                                                           gACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGaa                                                                                         <  2:5488054/150‑1 (MQ=255)

                                                           gACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGaa                                                                                         <  1:1622417/150‑1 (MQ=255)

                                                             tttcggtctGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGt                                                                                       >  1:4100431/3‑150 (MQ=255)

                                                             tttcggtctGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGt                                                                                       >  2:1071513/3‑150 (MQ=255)

                                                              gtcggtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTc                                                                                      >  1:8244597/1‑150 (MQ=255)

                                                                 ggtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTc                                                                                   <  2:785675/150‑1 (MQ=255)

                                                                     gACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCa                                                                               >  1:4025569/1‑150 (MQ=255)

                                                                        tGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGccg                                                                            >  2:4083978/1‑149 (MQ=255)

                                                                         gAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGccga                                                                           <  1:1500389/150‑3 (MQ=255)

                                                                             ggACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGaa                                                                       <  1:7540294/150‑1 (MQ=255)

                                                                             ggACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGaa                                                                       >  1:574934/1‑150 (MQ=255)

                                                                             ggACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGaa                                                                       <  1:246245/150‑1 (MQ=255)

                                                                                 gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCCGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgtg                                                                   >  2:363583/1‑150 (MQ=255)

                                                                                 gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgtg                                                                   >  1:3152471/1‑150 (MQ=255)

                                                                                 gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgtg                                                                   >  1:377046/1‑150 (MQ=255)

                                                                                 gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgtg                                                                   >  2:2806260/1‑150 (MQ=255)

                                                                                 gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgtg                                                                   >  2:3609455/1‑150 (MQ=255)

                                                                                 gaaaGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGAACgtg                                                                   >  2:2054099/4‑150 (MQ=255)

                                                                                 gaaaGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGAACgtg                                                                   >  1:3546664/4‑150 (MQ=255)

                                                                                 gaaaGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGAACgtg                                                                   <  1:6166721/147‑1 (MQ=255)

                                                                                   cAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTc                                                                 >  1:4584036/1‑150 (MQ=255)

                                                                                     gCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGc                                                               >  2:6495129/1‑150 (MQ=255)

                                                                                      cTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGcc                                                              >  2:5924192/1‑150 (MQ=255)

                                                                                      cTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGcc                                                              >  2:4197390/1‑150 (MQ=255)

                                                                                        cGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGt                                                            <  2:4326375/150‑1 (MQ=255)

                                                                                         gATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGtt                                                           >  1:8357220/1‑150 (MQ=255)

                                                                                                gCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCTCCACCCAGGAAGGTGTCGCCGTTGGTAGcc                                                    >  1:5986810/1‑150 (MQ=255)

                                                                                                gCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGcc                                                    >  1:6816786/1‑150 (MQ=255)

                                                                                                    tttCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTa                                                <  2:1307242/150‑1 (MQ=255)

                                                                                                          cGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGa                                          >  2:146413/1‑150 (MQ=255)

                                                                                                              ttcttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTg                                      <  1:1540882/150‑1 (MQ=255)

                                                                                                              ttcttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTg                                      >  2:8708390/1‑150 (MQ=255)

                                                                                                               tcttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTgg                                     >  1:7571038/1‑150 (MQ=255)

                                                                                                                cttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGt                                    <  1:8716791/150‑1 (MQ=255)

                                                                                                                 tttCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTg                                   >  2:5907242/1‑150 (MQ=255)

                                                                                                                      gACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCAc                              >  2:5256771/1‑150 (MQ=255)

                                                                                                                      gACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCAc                              >  2:5377681/1‑150 (MQ=255)

                                                                                                                       aCGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCAcc                             <  1:4362623/150‑1 (MQ=255)

                                                                                                                         gCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGt                           >  2:8586871/1‑150 (MQ=255)

                                                                                                                         gCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCAAGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGt                           >  2:1765015/1‑150 (MQ=255)

                                                                                                                              aGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACt                      >  2:3105507/1‑150 (MQ=255)

                                                                                                                                  ccAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGGTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCgg                  >  2:9100293/1‑150 (MQ=255)

                                                                                                                                   cAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGc                 >  1:221315/1‑150 (MQ=255)

                                                                                                                                   cAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGc                 >  1:1270929/1‑150 (MQ=255)

                                                                                                                                      agGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTATTTCGAACTGGTGCTCACCGTCGACTTCGGCGAt              >  1:1576919/1‑150 (MQ=255)

                                                                                                                                      agGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGAt              >  1:8063780/1‑150 (MQ=255)

                                                                                                                                      agGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGAt              >  1:7248892/1‑150 (MQ=255)

                                                                                                                                        gATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGAttt            >  1:1924788/1‑150 (MQ=255)

                                                                                                                                         aTCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTc           <  1:4188869/150‑1 (MQ=255)

                                                                                                                                           cGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGa         <  1:1462018/150‑1 (MQ=255)

                                                                                                                                           cGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGa         <  1:644176/150‑1 (MQ=255)

                                                                                                                                                ttcAGGTCCATGCCCGACTCTTTCTTGAACTCGTCTACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGAcc    <  1:6581083/150‑1 (MQ=255)

                                                                                                                                                 tcAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCg   >  2:453130/1‑150 (MQ=255)

GAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGA  >  NC_002947/5378633‑5378928

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