BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,260,354	A→T	35.7%	S395S (TCA→TCT)	hemA →	glutamyl tRNA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,260,354	0	A	T	35.7%	279.7 / 292.8	258	S395S (TCA→TCT)	hemA	glutamyl tRNA reductase
Reads supporting (aligned to +/- strand): ref base A (75/91); new base T (48/44); total (123/135)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 258 total aligned reads displayed.

CGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCC  >  CP009273/1260207‑1260501

   gTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTc                                                                                                                                                 <  1:3084158/149‑1 (MQ=255)

      tCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAAc                                                                                                                                              <  2:350762/149‑1 (MQ=255)

      tCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAAc                                                                                                                                              >  2:2764680/1‑149 (MQ=255)

      tCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAAc                                                                                                                                              >  1:1909301/1‑149 (MQ=255)

      tCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAAc                                                                                                                                              >  1:2826946/1‑149 (MQ=255)

          aGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGgc                                                                                                                                          >  1:2598791/1‑149 (MQ=255)

                  agagCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTg                                                                                                                                  <  1:2682848/149‑1 (MQ=255)

                   gagCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGa                                                                                                                                 >  1:247999/1‑149 (MQ=255)

                   gagCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGa                                                                                                                                 <  2:3569623/149‑1 (MQ=255)

                   gagCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGa                                                                                                                                 >  2:3299713/1‑149 (MQ=255)

                      cAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACgg                                                                                                                              <  2:3939875/149‑1 (MQ=255)

                       aaGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACggg                                                                                                                             <  2:2245575/149‑1 (MQ=255)

                        aGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACgggg                                                                                                                            >  2:2185067/1‑149 (MQ=255)

                        aGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACgggg                                                                                                                            >  2:1808710/1‑149 (MQ=255)

                          ttCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGAt                                                                                                                          >  1:1217759/1‑149 (MQ=255)

                          ttCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGAt                                                                                                                          >  1:688839/1‑149 (MQ=255)

                               gATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATaacga                                                                                                                     >  1:1644794/1‑149 (MQ=255)

                                 tGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATaacgaac                                                                                                                   >  1:848982/1‑149 (MQ=255)

                                     ttAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCt                                                                                                               >  1:1429097/1‑149 (MQ=255)

                                       aaCCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGa                                                                                                             >  1:3595123/1‑149 (MQ=255)

                                        aCCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATaacgaa                                                                                                                    >  2:3213871/1‑141 (MQ=255)

                                        aCCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATaacgaa                                                                                                                    >  2:653228/1‑141 (MQ=255)

                                        aCCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATaacgaa                                                                                                                    >  2:1160384/1‑141 (MQ=255)

                                        aCCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGaa                                                                                                            >  2:3766727/1‑149 (MQ=255)

                                              aaaGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCt                                                                                                      >  1:3393303/1‑149 (MQ=255)

                                              aaaGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCt                                                                                                      <  1:2530540/149‑1 (MQ=255)

                                              aaaGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCt                                                                                                      <  1:459104/149‑1 (MQ=255)

                                               aaGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTg                                                                                                     <  1:2575369/149‑1 (MQ=255)

                                                aGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTgc                                                                                                    >  1:1468264/1‑149 (MQ=255)

                                                aGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTgc                                                                                                    <  2:1628340/149‑1 (MQ=255)

                                                aGCGTGAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTgc                                                                                                    <  2:2605027/149‑1 (MQ=255)

                                                    ttAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGa                                                                                                >  1:423222/1‑149 (MQ=255)

                                                    ttAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGa                                                                                                >  2:3972314/1‑149 (MQ=255)

                                                    ttAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGa                                                                                                <  1:4005718/149‑1 (MQ=255)

                                                         ggCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGcc                                                                                           >  2:2705107/1‑149 (MQ=255)

                                                              ttGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCggg                                                                                      >  2:1024924/1‑149 (MQ=255)

                                                                gAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCTTGAATATTCTGCGCGACAGCCTCGGGCt                                                                                    >  2:39125/1‑149 (MQ=255)

                                                                gAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCt                                                                                    <  1:2204607/149‑1 (MQ=255)

                                                                  gCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTgg                                                                                  <  2:1337706/149‑1 (MQ=255)

                                                                      ggcggcGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTa                                                                              <  2:596507/149‑1 (MQ=255)

                                                                      ggcggcGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCTGGCTGGAGTa                                                                              <  2:3237364/149‑1 (MQ=255)

                                                                      ggcggcGACGCGCAAGACATTATGCAGGATCTGGCATGGGAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTa                                                                              <  1:3678108/149‑1 (MQ=255)

                                                                          gcgACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAg                                                                          <  1:3555732/149‑1 (MQ=255)

                                                                          gcgACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAg                                                                          <  2:1176096/149‑1 (MQ=255)

                                                                             aCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGtcc                                                                       >  2:142647/1‑147 (MQ=255)

                                                                              cgcgCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTAca                                                                      <  1:621709/149‑1 (MQ=255)

                                                                                cgcAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTAcatc                                                                    <  2:81593/149‑1 (MQ=255)

                                                                                cgcAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTAcatc                                                                    <  2:760099/149‑1 (MQ=255)

                                                                                 gcAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTAcatca                                                                   >  1:2032972/1‑149 (MQ=255)

                                                                                    aGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCAttt                                                                >  1:1288872/1‑149 (MQ=255)

                                                                                       cATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCt                                                             >  1:3701715/1‑149 (MQ=255)

                                                                                        aTTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCtt                                                            >  1:3533605/1‑149 (MQ=255)

                                                                                        aTTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCtt                                                            <  2:3859262/149‑1 (MQ=255)

                                                                                        aTTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTCGCAGTACATCATTTTCtt                                                            >  2:508692/1‑149 (MQ=255)

                                                                                        aTTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCtt                                                            <  2:4016123/149‑1 (MQ=255)

                                                                                         ttATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCttt                                                           <  1:1570845/149‑1 (MQ=255)

                                                                                           aTGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCttttt                                                         <  1:1118338/149‑1 (MQ=255)

                                                                                           aTGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCttttt                                                         <  2:3010140/149‑1 (MQ=255)

                                                                                             gCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCttttttt                                                       >  1:1366436/1‑149 (MQ=255)

                                                                                             gCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCttttttt                                                       <  1:3157860/149‑1 (MQ=255)

                                                                                                ggATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACa                                                    >  2:3663908/1‑149 (MQ=255)

                                                                                                 gATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAg                                                   <  1:4038427/149‑1 (MQ=255)

                                                                                                  aTCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAgg                                                  >  1:3053801/1‑149 (MQ=255)

                                                                                                     tGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTg                                               >  2:2520522/1‑149 (MQ=255)

                                                                                                      ggCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGc                                              >  1:1168050/1‑149 (MQ=255)

                                                                                                      ggCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGc                                              >  1:1505498/1‑149 (MQ=255)

                                                                                                      ggCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGc                                              >  2:2815245/1‑149 (MQ=255)

                                                                                                       gCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCa                                             <  1:1621118/149‑1 (MQ=255)

                                                                                                       gCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCa                                             <  2:2357009/149‑1 (MQ=255)

                                                                                                       gCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCa                                             <  2:283928/149‑1 (MQ=255)

                                                                                                       gCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCa                                             >  1:2626744/1‑149 (MQ=255)

                                                                                                           ggAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTa                                         <  2:314465/149‑1 (MQ=255)

                                                                                                            gAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTAc                                        <  1:2037184/149‑1 (MQ=255)

                                                                                                              aaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGc                                      <  2:1598554/149‑1 (MQ=255)

                                                                                                               aCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGcc                                     >  2:3326507/1‑149 (MQ=255)

                                                                                                               aCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGcc                                     >  1:105975/1‑149 (MQ=255)

                                                                                                                 tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTa                                   >  2:3957571/1‑149 (MQ=255)

                                                                                                                  gACTAACCGCTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTAt                                  >  1:3773276/1‑149 (MQ=255)

                                                                                                                  gACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTAt                                  <  1:2800168/149‑1 (MQ=255)

                                                                                                                   aCTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATg                                 <  2:1909301/149‑1 (MQ=255)

                                                                                                                          gCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCtt                          >  2:4074888/1‑149 (MQ=255)

                                                                                                                           cTTGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTc                         <  2:102153/149‑1 (MQ=255)

                                                                                                                           cTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTc                         >  2:2082246/1‑149 (MQ=255)

                                                                                                                            ttGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCt                        >  1:3453268/1‑149 (MQ=255)

                                                                                                                            ttGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCt                        <  2:1234689/149‑1 (MQ=255)

                                                                                                                            ttGATCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCt                        >  1:1484669/1‑149 (MQ=255)

                                                                                                                                tCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCg                    <  2:423222/149‑1 (MQ=255)

                                                                                                                                tCCATGCGCCAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCg                    <  2:268038/149‑1 (MQ=255)

                                                                                                                                 ccATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGt                   <  1:2764725/149‑1 (MQ=255)

                                                                                                                                    tGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGc                <  1:3326507/149‑1 (MQ=255)

                                                                                                                                     gcgcCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGcc               <  2:1187189/149‑1 (MQ=255)

                                                                                                                                        ccAACGAAATCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCaaa            >  1:1052932/1‑149 (MQ=255)

                                                                                                                                         cAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAAc           >  2:2791746/1‑149 (MQ=255)

                                                                                                                                           aCGAAATCTCTTCAACAGGCCGCCCGTGACAGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTg         >  2:1063096/1‑149 (MQ=255)

                                                                                                                                              aaaTCTCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGaa      >  2:3134123/1‑149 (MQ=255)

                                                                                                                                              aaaTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGaa      >  2:335587/1‑149 (MQ=255)

                                                                                                                                              aaaTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGaa      >  1:1823864/1‑149 (MQ=255)

                                                                                                                                              aaaTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGaa      <  1:3215126/149‑1 (MQ=255)

                                                                                                                                              aaaTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGaa      <  2:1316828/149‑1 (MQ=255)

                                                                                                                                               aaTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAg     >  1:3791609/1‑149 (MQ=255)

CGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCC  >  CP009273/1260207‑1260501

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