BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,513,935	A→T	68.6%	Y73N (TAC→AAC)	gltX ←	glutamyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,513,935	0	A	T	68.6%	489.9 / inf	344	Y73N (TAC→AAC)	gltX	glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base A (74/34); new base T (136/100); total (210/134)
Fisher's exact test for biased strand distribution p-value = 5.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 351 total aligned reads displayed.

GTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAGGAGTAAAG  >  CP009273/2513786‑2514082

    tAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGa                                                                                                                                                 >  2:1473015/1‑150 (MQ=255)

     aaCGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCACTAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGAc                                                                                                                                                >  1:3183196/1‑150 (MQ=255)

     aaCGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGAc                                                                                                                                                <  1:128182/150‑1 (MQ=255)

       cgcgGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCt                                                                                                                                              <  1:1830109/150‑1 (MQ=255)

       cgcgGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCt                                                                                                                                              >  1:4647122/1‑150 (MQ=255)

          ggCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCa                                                                                                                                           >  1:365014/1‑150 (MQ=255)

          ggCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCa                                                                                                                                           >  2:4656434/1‑150 (MQ=255)

              tctcACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATccc                                                                                                                                       <  2:639271/150‑1 (MQ=255)

                tcACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCAc                                                                                                                                     >  1:2144734/1‑150 (MQ=255)

                 cACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACt                                                                                                                                    <  2:820884/150‑1 (MQ=255)

                    cTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCa                                                                                                                                 >  2:3193476/1‑150 (MQ=255)

                         gCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTc                                                                                                                            >  1:654/1‑150 (MQ=255)

                          ccATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCa                                                                                                                           >  2:1535311/1‑150 (MQ=255)

                             tttGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGcc                                                                                                                        <  2:2201804/150‑1 (MQ=255)

                             tttGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGcc                                                                                                                        >  2:2350913/1‑150 (MQ=255)

                              ttGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCa                                                                                                                       <  2:726143/150‑1 (MQ=255)

                              ttGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCa                                                                                                                       >  1:678932/1‑150 (MQ=255)

                               tGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAg                                                                                                                      <  1:483856/150‑1 (MQ=255)

                                 ctctTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGtt                                                                                                                    >  1:2769683/1‑150 (MQ=255)

                                       gcgcAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGcc                                                                                                              <  1:1674138/150‑1 (MQ=255)

                                        cgcAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGcca                                                                                                             <  2:4950496/150‑1 (MQ=255)

                                        cgcAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGcca                                                                                                             <  1:3549250/150‑1 (MQ=255)

                                         gcAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccat                                                                                                            >  1:5190037/1‑150 (MQ=255)

                                         gcAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccat                                                                                                            >  2:2779715/1‑150 (MQ=255)

                                          cAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatc                                                                                                           <  2:128523/150‑1 (MQ=255)

                                          cAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatc                                                                                                           >  2:2065965/1‑150 (MQ=255)

                                           aGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatcc                                                                                                          >  1:3263293/1‑150 (MQ=255)

                                           aGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatcc                                                                                                          >  2:562705/1‑150 (MQ=255)

                                            gcgcTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatcca                                                                                                         <  2:3228768/150‑1 (MQ=255)

                                            gcgcTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatcca                                                                                                         >  2:3833754/1‑150 (MQ=255)

                                            gcgcTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatcca                                                                                                         <  1:4113313/150‑1 (MQ=255)

                                            gcgcGTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGccatcca                                                                                                         <  2:1113760/150‑1 (MQ=255)

                                              gcTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATa                                                                                                       >  2:918867/1‑150 (MQ=255)

                                                ttCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAAt                                                                                                     <  1:92065/150‑1 (MQ=255)

                                                ttCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAAt                                                                                                     >  1:4689753/1‑150 (MQ=255)

                                                ttCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAAt                                                                                                     >  2:3583678/1‑150 (MQ=255)

                                                ttCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCAAAAt                                                                                                     >  2:4564486/1‑150 (MQ=255)

                                                 tCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATg                                                                                                    >  1:5145500/1‑150 (MQ=255)

                                                 tCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATg                                                                                                    >  2:988439/1‑150 (MQ=255)

                                                    aGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCt                                                                                                 <  1:1103535/150‑1 (MQ=255)

                                                    aGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCt                                                                                                 >  1:4327343/1‑150 (MQ=255)

                                                    aGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCt                                                                                                 >  2:4989651/1‑150 (MQ=255)

                                                     ggCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCtt                                                                                                <  1:275907/150‑1 (MQ=255)

                                                      gcgcTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTc                                                                                               <  2:175538/150‑1 (MQ=255)

                                                      gcgcTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTc                                                                                               <  2:3366427/150‑1 (MQ=255)

                                                      gcgcTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTc                                                                                               <  1:709535/150‑1 (MQ=255)

                                                         ctctTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGAt                                                                                            >  2:1028062/1‑150 (MQ=255)

                                                            tttAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGc                                                                                         <  1:38466/150‑1 (MQ=255)

                                                             ttAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCt                                                                                        <  1:24333/150‑1 (MQ=255)

                                                                 agCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCg                                                                                    <  1:3471595/150‑1 (MQ=255)

                                                                 agCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTGTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGATTCGATAGCTTCCg                                                                                    <  2:4687111/150‑1 (MQ=255)

                                                                   cAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGc                                                                                  >  2:4228313/1‑150 (MQ=255)

                                                                    aaTAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCg                                                                                 >  2:4488199/1‑150 (MQ=255)

                                                                    aaTAGCATTAAT‑AGCAGTGCCCTCTTCCAGCAACTGATCGATCACCGCGTTGTAGCGATAAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGGTCATGCCATCCATAATGGCGTCGATAGCTTCCGGCGt                                                                                >  1:2512663/1‑150 (MQ=255)

                                                                     aTAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGt                                                                                <  1:1952180/150‑1 (MQ=255)

                                                                         cATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGaa                                                                            >  2:3532895/1‑150 (MQ=255)

                                                                         cATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGaa                                                                            >  1:2057748/1‑150 (MQ=255)

                                                                            ttATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGc                                                                         >  2:4428/1‑150 (MQ=255)

                                                                             tataAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCt                                                                        >  2:2562415/1‑150 (MQ=255)

                                                                             tataAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGAACAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCt                                                                        >  2:1781981/1‑150 (MQ=255)

                                                                              ataAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTc                                                                       >  1:4429544/1‑150 (MQ=255)

                                                                                aaGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCaa                                                                     >  2:5177161/1‑150 (MQ=255)

                                                                                  gCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGa                                                                   >  1:23812/1‑150 (MQ=255)

                                                                                   cAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGAt                                                                  >  1:4289913/1‑150 (MQ=255)

                                                                                         ccTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGgtgt                                                            >  1:4488665/1‑150 (MQ=255)

                                                                                          ctctTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTc                                                           >  2:5278149/1‑150 (MQ=255)

                                                                                           tctTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCt                                                          <  1:2309987/150‑1 (MQ=255)

                                                                                            ctTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCtt                                                         >  2:281266/1‑150 (MQ=255)

                                                                                            ctTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCtt                                                         <  1:403062/150‑1 (MQ=255)

                                                                                                cAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAAt                                                     <  2:3989181/150‑1 (MQ=255)

                                                                                                cAGCATCTGATCGAACACCGCGTAGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTACAAt                                                     >  2:5250044/1‑150 (MQ=255)

                                                                                                 aGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCAGTGTCTTCAATa                                                    <  1:2952978/150‑1 (MQ=255)

                                                                                                 aGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATa                                                    >  2:3056904/1‑150 (MQ=255)

                                                                                                  gCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATAc                                                   <  2:2385371/150‑1 (MQ=255)

                                                                                                  gCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATAc                                                   <  1:3578463/150‑1 (MQ=255)

                                                                                                      cTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACgcag                                               >  1:842293/1‑150 (MQ=255)

                                                                                                               cACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTACCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTc                                      >  2:4841686/1‑150 (MQ=255)

                                                                                                               cACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGAAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTc                                      <  1:1699621/150‑1 (MQ=255)

                                                                                                                    cgTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCAccgc                                 >  1:4916791/1‑150 (MQ=255)

                                                                                                                     gTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCAccgcc                                >  2:2751065/1‑150 (MQ=255)

                                                                                                                     gTTGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCAccgcc                                >  2:1933277/1‑150 (MQ=255)

                                                                                                                       tGTAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGt                              >  2:3676081/1‑150 (MQ=255)

                                                                                                                         tAGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTgg                            >  1:878327/1‑150 (MQ=255)

                                                                                                                          aGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGt                           <  2:1927643/150‑1 (MQ=255)

                                                                                                                          aGCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGt                           >  1:1518280/1‑150 (MQ=255)

                                                                                                                           gCGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGtt                          <  1:3001287/150‑1 (MQ=255)

                                                                                                                            cGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTa                         >  1:2031261/1‑150 (MQ=255)

                                                                                                                            cGATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTa                         >  2:1595845/1‑150 (MQ=255)

                                                                                                                            cGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTa                         >  1:5027372/1‑150 (MQ=255)

                                                                                                                             gATCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTAc                        <  1:3700611/150‑1 (MQ=255)

                                                                                                                              aTCAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACg                       >  2:259096/1‑150 (MQ=255)

                                                                                                                              aTCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACg                       >  1:3148601/1‑150 (MQ=255)

                                                                                                                                cAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATTGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTg                     >  2:4360903/1‑150 (MQ=255)

                                                                                                                                cAAAACGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTg                     >  2:59033/1‑150 (MQ=255)

                                                                                                                                   aaCGTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCaa                  >  2:269189/1‑150 (MQ=255)

                                                                                                                                      gTTTGGTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCaaaaa               >  1:3371848/1‑150 (MQ=255)

                                                                                                                                           gTCTGGTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAg          >  1:1123229/1‑150 (MQ=255)

                                                                                                                                            tCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAgg         >  1:2689611/1‑150 (MQ=255)

                                                                                                                                               ggTAGTTCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAGGAGt      <  1:4602922/150‑1 (MQ=255)

                                                                                                                                                  agttCGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAGGAGTaaa   >  2:1450529/1‑150 (MQ=255)

GTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACGGACCTTCATCCCACTCCAGGCTCAGCCAGTTCATGCCATCCATAATGGCTTCGATAGCTTCCGGCGTGGAACGCTCAAGATCGGTGTCTTCAATACGCAGCACGAACTCACCGCCGTGGTTACGTGCAAAAAGCCAGGAGTAAAG  >  CP009273/2513786‑2514082

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