BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,297,485	G→C	50.2%	G44G (GGC→GGG)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,297,485	0	G	C	50.2%	‑3.4 / inf	239	G44G (GGC→GGG)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand): ref base G (56/63); new base C (60/60); total (116/123)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 243 total aligned reads displayed.

CGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCGCTTT  >  CP009273/4297338‑4297634

cgccgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGcc                                                                                                                                                     >  1:4040235/1‑150 (MQ=255)

 gccgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCg                                                                                                                                                    >  1:2622006/1‑150 (MQ=255)

 gccgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCg                                                                                                                                                    >  2:4235631/1‑150 (MQ=255)

 gccgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCAATTAGCCGATACCCGACACCAGGCCg                                                                                                                                                    <  2:2281000/150‑1 (MQ=255)

  ccgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCgg                                                                                                                                                   >  1:1105164/1‑150 (MQ=255)

  ccgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCgg                                                                                                                                                   >  1:442327/1‑150 (MQ=255)

   cgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCggg                                                                                                                                                  >  2:3721134/1‑150 (MQ=255)

          aaCGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATg                                                                                                                                           <  1:3395990/150‑1 (MQ=255)

           aCGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGa                                                                                                                                          <  1:3308408/150‑1 (MQ=255)

            cGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGAc                                                                                                                                         >  1:3798135/1‑150 (MQ=255)

             ggCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACt                                                                                                                                        <  2:688533/150‑1 (MQ=255)

               cAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTc                                                                                                                                      <  2:4144307/150‑1 (MQ=255)

               cAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTc                                                                                                                                      >  2:2343909/1‑150 (MQ=255)

                  gTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGc                                                                                                                                   <  1:98670/150‑1 (MQ=255)

                     aGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGt                                                                                                                                >  1:89255/1‑150 (MQ=255)

                        cGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCg                                                                                                                             >  1:2440886/1‑150 (MQ=255)

                            gAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCttt                                                                                                                         >  1:190142/1‑150 (MQ=255)

                               ataatGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCtttttt                                                                                                                      >  2:780203/1‑150 (MQ=255)

                                taatGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCttttttt                                                                                                                     <  1:2115840/150‑1 (MQ=255)

                                taatGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCttttttt                                                                                                                     <  1:4231872/150‑1 (MQ=255)

                                          tttATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAg                                                                                                           <  1:215092/150‑1 (MQ=255)

                                            tATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCg                                                                                                         >  1:2348030/1‑150 (MQ=255)

                                               cTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGttt                                                                                                      >  1:2087087/1‑150 (MQ=255)

                                               cTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGttt                                                                                                      <  2:1042390/150‑1 (MQ=255)

                                                tttACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTc                                                                                                     <  1:2999559/150‑1 (MQ=255)

                                                tttACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTc                                                                                                     <  2:345762/150‑1 (MQ=255)

                                                 ttACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCa                                                                                                    >  1:1368743/1‑150 (MQ=255)

                                                 ttACTGACCAGCGCCGGAAATCCCATAACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCa                                                                                                    <  2:3398274/150‑1 (MQ=255)

                                                    cTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCt                                                                                                 >  1:1761685/1‑150 (MQ=255)

                                                       aCCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTct                                                                                              <  1:3444770/150‑1 (MQ=255)

                                                         cAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTctgc                                                                                            <  1:3075949/150‑1 (MQ=255)

                                                          aGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTctgct                                                                                           <  1:1308627/150‑1 (MQ=255)

                                                          aGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTctgct                                                                                           <  1:846634/150‑1 (MQ=255)

                                                            cgcCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGt                                                                                         >  2:239752/1‑150 (MQ=255)

                                                             gcCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTc                                                                                        >  2:3165269/1‑150 (MQ=255)

                                                                ggAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGc                                                                                     >  1:3993259/1‑150 (MQ=255)

                                                                ggAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCGGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGc                                                                                     >  2:531425/1‑150 (MQ=255)

                                                                  aaaTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAg                                                                                   >  1:3194303/1‑150 (MQ=255)

                                                                  aaaTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAg                                                                                   >  1:2871862/1‑150 (MQ=255)

                                                                  aaaTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAg                                                                                   >  1:625510/1‑150 (MQ=255)

                                                                   aaTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGa                                                                                  <  1:1564656/150‑1 (MQ=255)

                                                                        cATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACaaaa                                                                             >  1:3149456/1‑150 (MQ=255)

                                                                         aTCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAg                                                                            <  2:3940652/150‑1 (MQ=255)

                                                                         aTCACAAGACCATGACAGCGAGCGTTAAAGCGCCTGATTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCATTGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAg                                                                            <  2:1753547/150‑1 (MQ=255)

                                                                              cAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGat                                                                       >  2:148028/1‑150 (MQ=255)

                                                                                 aCCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATg                                                                    <  1:1342222/150‑1 (MQ=255)

                                                                                 aCCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATg                                                                    <  1:3537171/150‑1 (MQ=255)

                                                                                 aCCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATg                                                                    <  1:3580562/150‑1 (MQ=255)

                                                                                   cATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCg                                                                  >  1:1062368/1‑150 (MQ=255)

                                                                                      gACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTg                                                               <  2:3005920/150‑1 (MQ=255)

                                                                                       aCAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGc                                                              >  1:4216314/1‑150 (MQ=255)

                                                                                        cAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGcc                                                             <  2:1841801/150‑1 (MQ=255)

                                                                                         agcgagcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGccc                                                            >  2:1650689/1‑150 (MQ=255)

                                                                                            gagcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCa                                                         <  2:1424884/150‑1 (MQ=255)

                                                                                             agcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCat                                                        <  1:169212/150‑1 (MQ=255)

                                                                                               cgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCatat                                                      <  2:387149/150‑1 (MQ=255)

                                                                                               cgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCatat                                                      <  2:642769/150‑1 (MQ=255)

                                                                                                  tAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCa                                                   <  2:892071/150‑1 (MQ=255)

                                                                                                  tAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCa                                                   <  2:3071666/150‑1 (MQ=255)

                                                                                                   aaaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCAc                                                  <  2:4203698/150‑1 (MQ=255)

                                                                                                   aaaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCAc                                                  >  1:1530498/1‑150 (MQ=255)

                                                                                                    aaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACg                                                 <  1:187536/150‑1 (MQ=255)

                                                                                                     aGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGc                                                >  1:1659400/1‑150 (MQ=255)

                                                                                                          cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCt                                           <  1:2511958/150‑1 (MQ=255)

                                                                                                          cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCt                                           <  1:2510825/150‑1 (MQ=255)

                                                                                                           tGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTa                                          >  1:2656249/1‑150 (MQ=255)

                                                                                                           tGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTa                                          >  2:2428812/1‑150 (MQ=255)

                                                                                                           tGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCGTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTa                                          >  1:3317092/1‑150 (MQ=255)

                                                                                                            gAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTac                                         >  2:3661602/1‑150 (MQ=255)

                                                                                                             aGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacg                                        <  2:1515878/150‑1 (MQ=255)

                                                                                                             aGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacg                                        <  1:794850/150‑1 (MQ=255)

                                                                                                              gTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacga                                       >  1:2424493/1‑150 (MQ=255)

                                                                                                               ttGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacgac                                      >  1:651648/1‑150 (MQ=255)

                                                                                                                tGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacgacg                                     <  1:2736784/150‑1 (MQ=255)

                                                                                                                 gCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGt                                    <  2:1922156/150‑1 (MQ=255)

                                                                                                                 gCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacgacgg                                    >  1:3343695/1‑149 (MQ=255)

                                                                                                                 gCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGt                                    <  2:2871862/150‑1 (MQ=255)

                                                                                                                  cTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGtt                                   <  2:396869/150‑1 (MQ=255)

                                                                                                                   tCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTa                                  <  1:3020991/150‑1 (MQ=255)

                                                                                                                     gTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATg                                >  1:1184630/1‑150 (MQ=255)

                                                                                                                          tCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGtt                           <  1:3864238/150‑1 (MQ=255)

                                                                                                                          tCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGtt                           >  2:2171552/1‑150 (MQ=255)

                                                                                                                          tCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGtt                           <  2:367277/150‑1 (MQ=255)

                                                                                                                           cATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGttt                          <  1:455116/150‑1 (MQ=255)

                                                                                                                           cATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACGTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGttt                          <  2:2847115/150‑1 (MQ=255)

                                                                                                                           cATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATTCTGttt                          >  1:3705181/1‑150 (MQ=255)

                                                                                                                            aTTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGtttt                         <  2:2468283/150‑1 (MQ=255)

                                                                                                                             tttCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGttttt                        >  1:1011587/1‑150 (MQ=255)

                                                                                                                              ttCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTg                       >  2:892795/1‑150 (MQ=255)

                                                                                                                                   cGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGa                  >  1:488969/1‑150 (MQ=255)

                                                                                                                                     aTACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAgg                <  2:3034361/150‑1 (MQ=255)

                                                                                                                                       aCCCGACACCATGCCGGGAGCAATGAATTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCt              <  1:1426501/150‑1 (MQ=255)

                                                                                                                                        cccGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTg             >  2:2070526/1‑150 (MQ=255)

                                                                                                                                        cccGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTg             <  2:125969/150‑1 (MQ=255)

                                                                                                                                          cGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTa           >  2:825611/1‑150 (MQ=255)

                                                                                                                                            acacCAGCCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTAcc         <  2:2970621/150‑1 (MQ=255)

                                                                                                                                               ccAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCg      <  1:3918014/150‑1 (MQ=255)

                                                                                                                                               ccAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCg      <  1:1729770/150‑1 (MQ=255)

                                                                                                                                               ccAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCATCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCAAATCCACGCCCGCTACGCCGTTATGCTGTTTTTGCATGAGGCGGTACcgcg      >  1:3720790/1‑150 (MQ=255)

                                                                                                                                                   gCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCGCttt  <  2:2139638/150‑1 (MQ=255)

CGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCGCTTT  >  CP009273/4297338‑4297634

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