BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	388,341	G→T	8.8%	S6* (TCA→TAA)	insE1 ←	IS3 transposase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	388,341	0	G	T	8.8%	362.5 / 14.6	102	S6* (TCA→TAA)	insE1	IS3 transposase A
Reads supporting (aligned to +/- strand): ref base G (40/53); new base T (7/2); total (47/55)
Fisher's exact test for biased strand distribution p-value = 7.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 105 total aligned reads displayed.

AAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCA  >  CP009273/388280‑388490

aaGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACa                                                               >  1:776638/1‑150 (MQ=11)

aaGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACa                                                               <  2:276504/150‑1 (MQ=11)

aaGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACa                                                               <  2:2849417/150‑1 (MQ=11)

aaGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACa                                                               <  2:760341/150‑1 (MQ=11)

  gCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGTTACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAAc                                                             <  1:2032068/150‑1 (MQ=17)

    ttCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCt                                                           <  1:2258583/150‑1 (MQ=33)

      cAGGGCTTCACTGCGAATTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGt                                                         <  2:1493109/150‑1 (MQ=25)

       aGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTTATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTa                                                        <  1:2283723/150‑1 (MQ=25)

       aGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTa                                                        >  1:702285/1‑150 (MQ=35)

       aGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTa                                                        >  1:3538700/1‑150 (MQ=35)

       aGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTa                                                        >  1:702163/1‑150 (MQ=35)

       aGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATAATGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTa                                                        >  1:725027/1‑150 (MQ=25)

         ggCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACg                                                      >  2:3882938/1‑150 (MQ=38)

            ttCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCa                                                   <  2:1937831/150‑1 (MQ=40)

            ttCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCa                                                   <  2:2490242/150‑1 (MQ=40)

             tCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCaa                                                  >  1:1682429/1‑150 (MQ=40)

                 tGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACt                                              <  2:597740/150‑1 (MQ=255)

                 tGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACt                                              >  1:2545766/1‑150 (MQ=255)

                  gCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTa                                             <  1:1121764/150‑1 (MQ=255)

                   cGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTTATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTAt                                            >  1:3759884/1‑150 (MQ=255)

                   cGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTTATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTAt                                            >  1:3759873/1‑150 (MQ=255)

                   cGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTAt                                            >  1:3622986/1‑150 (MQ=255)

                   cGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTAt                                            >  1:3623055/1‑150 (MQ=255)

                   cGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTAt                                            >  1:3670671/1‑150 (MQ=255)

                    gAAATTCAGGCGAATGNTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATc                                           >  2:3955487/1‑150 (MQ=255)

                    gAAATTCAGGCGAATGGTGTTTAAGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATc                                           <  2:1202806/150‑1 (MQ=255)

                    gAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATc                                           <  2:661468/150‑1 (MQ=255)

                    gAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATc                                           >  2:1313835/1‑150 (MQ=255)

                     aaaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCa                                          >  1:622287/1‑150 (MQ=255)

                     aaaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCa                                          >  1:622349/1‑150 (MQ=255)

                      aaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCAc                                         >  1:414037/1‑150 (MQ=255)

                      aaTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCAc                                         >  2:3324590/1‑150 (MQ=255)

                       aTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCATTACTATCAcc                                        <  2:3350266/150‑1 (MQ=255)

                       aTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCAcc                                        <  2:3350345/150‑1 (MQ=255)

                       aTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCAcc                                        <  2:3350324/150‑1 (MQ=255)

                       aTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCAcc                                        <  2:3350475/150‑1 (MQ=255)

                        ttCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCg                                       <  2:2102802/150‑1 (MQ=255)

                        ttCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCg                                       <  2:2102911/150‑1 (MQ=255)

                          cAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGta                                     >  2:2573308/1‑150 (MQ=255)

                           aGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaa                                    >  1:1191016/1‑150 (MQ=255)

                             gCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaact                                  <  1:2129923/150‑1 (MQ=255)

                               gAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaactaa                                >  2:13920/1‑150 (MQ=255)

                               gAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaactaa                                >  1:3732127/1‑150 (MQ=255)

                               gAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaactaa                                >  1:3732122/1‑150 (MQ=255)

                               gAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaactaa                                >  1:3732098/1‑150 (MQ=255)

                               gAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTAATGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaactaa                                >  1:3732097/1‑150 (MQ=255)

                                aaTGCTGTTTACGGGGTTTTTTACTGGTTTATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaactaac                               >  2:1186560/1‑150 (MQ=255)

                                aaTGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaactaac                               <  1:732625/150‑1 (MQ=255)

                                aaTGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGtaactaac                               <  2:1942816/150‑1 (MQ=255)

                                 aTGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAAcg                              <  1:666580/150‑1 (MQ=255)

                                 aTGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAAcg                              <  2:4118835/150‑1 (MQ=255)

                                  tGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAAcgg                             <  1:2709467/150‑1 (MQ=255)

                                  tGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAAcgg                             >  1:3713422/1‑150 (MQ=255)

                                    cTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAAcggcg                           >  2:747412/1‑150 (MQ=255)

                                       tttACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGgtg                        <  1:2333640/150‑1 (MQ=255)

                                         tACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTa                      <  2:3146848/150‑1 (MQ=255)

                                          aCGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTAt                     <  1:2071587/150‑1 (MQ=255)

                                          aCGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTAt                     <  2:133653/150‑1 (MQ=255)

                                           cGGGGTTTTTTACTGGTTTATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTAtt                    >  1:576702/1‑150 (MQ=255)

                                           cGGGGTTTTTTACTGGTTTATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTAtt                    >  1:576724/1‑150 (MQ=255)

                                           cGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTAtt                    <  1:3999077/150‑1 (MQ=255)

                                           cGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTAtt                    <  2:1148006/150‑1 (MQ=255)

                                           cGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTAtt                    <  1:3999030/150‑1 (MQ=255)

                                             gggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACTACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTgg                  <  2:365446/150‑1 (MQ=255)

                                             gggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTgg                  <  2:695629/150‑1 (MQ=255)

                                             gggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTgg                  <  2:3850167/150‑1 (MQ=255)

                                              ggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGa                 >  2:3608555/1‑150 (MQ=255)

                                              ggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGa                 >  2:304043/1‑150 (MQ=255)

                                                ttttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGAtg               <  2:499128/150‑1 (MQ=255)

                                                ttttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGAtg               <  1:300789/150‑1 (MQ=255)

                                                  ttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGAtgtg             <  1:1228029/150‑1 (MQ=255)

                                                  ttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGAtgtg             <  1:1228069/150‑1 (MQ=255)

                                                  ttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGAtgtg             <  1:1228070/150‑1 (MQ=255)

                                                  ttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGAtgtg             <  1:746015/150‑1 (MQ=255)

                                                  ttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGAtgtg             <  1:746138/150‑1 (MQ=255)

                                                    ttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGaa           >  2:2928626/1‑150 (MQ=255)

                                                    ttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGaa           >  1:3936223/1‑150 (MQ=255)

                                                    ttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGaa           <  2:3831067/150‑1 (MQ=255)

                                                    ttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGaa           <  1:3317698/150‑1 (MQ=255)

                                                     tACTGGTTTATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAAc          <  2:2038797/150‑1 (MQ=255)

                                                     tACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAAc          <  1:3814849/150‑1 (MQ=255)

                                                     tACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAAc          >  2:299499/1‑150 (MQ=255)

                                                     tACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAAc          <  1:136326/150‑1 (MQ=255)

                                                     tACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAAc          >  2:3895753/1‑150 (MQ=255)

                                                      aCTGGTTTATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACg         >  1:3732985/1‑150 (MQ=255)

                                                      aCTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACg         <  1:3339957/150‑1 (MQ=255)

                                                      aCTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACg         >  1:3956626/1‑150 (MQ=255)

                                                      aCTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACg         <  2:3874583/150‑1 (MQ=255)

                                                       cTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGt        >  2:3866373/1‑150 (MQ=255)

                                                       cTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGt        >  1:1276109/1‑150 (MQ=255)

                                                        tGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGtt       <  2:1172523/150‑1 (MQ=255)

                                                        tGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGtt       <  2:2159675/150‑1 (MQ=255)

                                                        tGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGtt       <  2:1172053/150‑1 (MQ=255)

                                                        tGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGtt       <  2:2159869/150‑1 (MQ=255)

                                                        tGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGtt       >  1:3155609/1‑150 (MQ=255)

                                                           ttGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGAt    >  2:137533/1‑150 (MQ=255)

                                                           ttGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGAt    >  2:137626/1‑150 (MQ=255)

                                                             gATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCa  <  2:227306/150‑1 (MQ=255)

                                                             gATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCa  <  2:1104531/150‑1 (MQ=255)

                                                             gATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCa  <  2:3161819/150‑1 (MQ=255)

AAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGATTACAACCTGTACGCCAATACTATCACCGTAACTAACGGCGGTGTATTGGATGTGAACGTTGATCA  >  CP009273/388280‑388490

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