BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,081,697	T→A	26.8%	C19S (TGC→AGC)	galP →	D‑galactose transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,081,697	0	T	A	26.8%	527.1 / 226.7	269	C19S (TGC→AGC)	galP	D‑galactose transporter
Reads supporting (aligned to +/- strand): ref base T (124/73); new base A (57/15); total (181/88)
Fisher's exact test for biased strand distribution p-value = 1.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 269 total aligned reads displayed.

AGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCG  >  CP009273/3081550‑3081844

    aTTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTcc                                                                                                                                               >  1:702621/1‑150 (MQ=255)

    aTTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTcc                                                                                                                                               >  1:1294280/1‑150 (MQ=255)

     ttACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCT‑AAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCtt                                                                                                                                             <  2:4621998/150‑1 (MQ=255)

        cacTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGc                                                                                                                                           >  1:1211326/1‑150 (MQ=255)

               gtgATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTg                                                                                                                                    >  2:695646/1‑150 (MQ=255)

                 gATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGc                                                                                                                                  >  2:2392903/1‑150 (MQ=255)

                  aTTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCg                                                                                                                                 >  2:1157723/1‑150 (MQ=255)

                  aTTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCg                                                                                                                                 >  2:1852864/1‑150 (MQ=255)

                   tttGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCgg                                                                                                                                >  2:3928494/1‑150 (MQ=255)

                      gCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGAt                                                                                                                             <  1:4565455/150‑1 (MQ=255)

                       cTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGAtt                                                                                                                            >  1:493104/1‑150 (MQ=255)

                        ttCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTa                                                                                                                           >  2:4155467/1‑150 (MQ=255)

                          cacaTCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTAct                                                                                                                         >  1:78405/1‑150 (MQ=255)

                             aTCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCtt                                                                                                                      >  1:4137471/1‑150 (MQ=255)

                              tCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCttt                                                                                                                     <  2:344462/150‑1 (MQ=255)

                              tCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCttt                                                                                                                     <  1:3534217/150‑1 (MQ=255)

                                  tttACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGGCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGcc                                                                                                                 >  1:3270879/1‑150 (MQ=255)

                                      cgtcgtACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGTCGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGCTTACTCTTTGGCCTGGa                                                                                                             >  1:1960513/1‑150 (MQ=255)

                                         cgtACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGatat                                                                                                          >  2:1500464/1‑150 (MQ=255)

                                         cgtACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGatat                                                                                                          >  2:1232801/1‑150 (MQ=255)

                                          gtACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATc                                                                                                         >  1:1429141/1‑150 (MQ=255)

                                              tCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGgtg                                                                                                     >  1:294969/1‑150 (MQ=255)

                                               cACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGgtgt                                                                                                    <  1:1036060/150‑1 (MQ=255)

                                                aCCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTa                                                                                                   >  2:4728768/1‑150 (MQ=255)

                                                   tATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAAtt                                                                                                <  2:4658050/150‑1 (MQ=255)

                                                    aTCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTg                                                                                               >  2:2730985/1‑150 (MQ=255)

                                                    aTCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTg                                                                                               <  1:4511855/150‑1 (MQ=255)

                                                    aTCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTg                                                                                               <  1:2658673/150‑1 (MQ=255)

                                                    aTCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTg                                                                                               >  1:3989478/1‑150 (MQ=255)

                                                    aTCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTg                                                                                               >  1:3989496/1‑150 (MQ=255)

                                                     tCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGc                                                                                              >  1:686155/1‑150 (MQ=255)

                                                      cTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCt                                                                                             >  1:1509915/1‑150 (MQ=255)

                                                       ttAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTg                                                                                            >  1:2817819/1‑150 (MQ=255)

                                                        tAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTgg                                                                                           <  2:1575940/150‑1 (MQ=255)

                                                            tCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCa                                                                                       >  2:1887553/1‑150 (MQ=255)

                                                             cacaATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCAc                                                                                      >  2:3388675/1‑150 (MQ=255)

                                                             cacaATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCAc                                                                                      >  2:3586942/1‑150 (MQ=255)

                                                                  tAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCg                                                                                 >  1:2254124/1‑150 (MQ=255)

                                                                  tAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCg                                                                                 >  1:4005027/1‑150 (MQ=255)

                                                                  tAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCg                                                                                 >  2:3596825/1‑150 (MQ=255)

                                                                  tAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCg                                                                                 >  1:25312/1‑150 (MQ=255)

                                                                      aaaTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTa                                                                             >  2:2949707/1‑150 (MQ=255)

                                                                       aaTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTAt                                                                            <  2:1051366/150‑1 (MQ=255)

                                                                        aTAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTAtt                                                                           >  2:617158/1‑150 (MQ=255)

                                                                           aCCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCa                                                                        >  2:1967236/1‑150 (MQ=255)

                                                                            ccATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAg                                                                       >  2:4122800/1‑150 (MQ=255)

                                                                             cATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGa                                                                      <  1:829631/150‑1 (MQ=255)

                                                                               tatTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATg                                                                    <  2:3160501/150‑1 (MQ=255)

                                                                                atTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGa                                                                   >  1:2781504/1‑150 (MQ=255)

                                                                                 ttGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGaa                                                                  >  2:4379814/1‑150 (MQ=255)

                                                                                 ttGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGaa                                                                  >  2:1269854/1‑150 (MQ=255)

                                                                                  tGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAAt                                                                 <  1:3415729/150‑1 (MQ=255)

                                                                                   ggAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAAtt                                                                <  1:2585846/150‑1 (MQ=255)

                                                                                     aGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTcc                                                              <  2:829829/150‑1 (MQ=255)

                                                                                     aGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTcc                                                              <  2:2891144/150‑1 (MQ=255)

                                                                                     aGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTcc                                                              >  2:1683562/1‑150 (MQ=255)

                                                                                     aGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTcc                                                              >  2:4533551/1‑150 (MQ=255)

                                                                                     aGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTcc                                                              >  2:3619148/1‑150 (MQ=255)

                                                                                     aGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTcc                                                              >  2:531267/1‑150 (MQ=255)

                                                                                     aGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTcc                                                              >  1:188169/1‑150 (MQ=255)

                                                                                         catcatGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGAt                                                          >  2:3060580/1‑150 (MQ=255)

                                                                                          atcatGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGAtt                                                         >  2:2920612/1‑150 (MQ=255)

                                                                                           tcatGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGCTATCGGTGTAATTGCTGGCGCACTGCCGTGTATTGCAGATGAATTCCAGATTa                                                        >  2:1787052/1‑150 (MQ=255)

                                                                                           tcatGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTa                                                        >  2:2077098/1‑150 (MQ=255)

                                                                                           tcatGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTa                                                        <  2:2180172/150‑1 (MQ=255)

                                                                                           tcatGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTa                                                        >  2:4133372/1‑150 (MQ=255)

                                                                                             atGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACt                                                      >  2:4186695/1‑150 (MQ=255)

                                                                                              tGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACtt                                                     <  1:878852/150‑1 (MQ=255)

                                                                                                 cTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGc                                                  >  1:3737663/1‑150 (MQ=255)

                                                                                                    aCGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGcaca                                               >  1:313955/1‑150 (MQ=255)

                                                                                                     cGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCacac                                              >  1:1739467/1‑150 (MQ=255)

                                                                                                      gCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACg                                             >  1:2436463/1‑150 (MQ=255)

                                                                                                       cTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGc                                            >  2:1645502/1‑150 (MQ=255)

                                                                                                        tAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCa                                           <  2:4757305/150‑1 (MQ=255)

                                                                                                        tAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCa                                           >  1:4770184/1‑150 (MQ=255)

                                                                                                        tAAAAAACAGGGGCGGTAAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCa                                           <  2:4266576/150‑1 (MQ=255)

                                                                                                          aaaaaCAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTTCAGATGAATTACAGATTACTTAGCACACGCAAg                                         >  1:3205543/1‑150 (MQ=255)

                                                                                                          aaaaaCAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAg                                         <  2:4584013/150‑1 (MQ=255)

                                                                                                           aaaaCAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGa                                        >  1:4595831/1‑150 (MQ=255)

                                                                                                             aaCAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAAt                                      >  2:2515096/1‑150 (MQ=255)

                                                                                                               cAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATgg                                    >  2:786480/1‑150 (MQ=255)

                                                                                                                 ggggCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGt                                  >  2:3420165/1‑150 (MQ=255)

                                                                                                                      ggTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTaa                             >  2:2256839/1‑150 (MQ=255)

                                                                                                                         cAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCt                          >  2:4550276/1‑150 (MQ=255)

                                                                                                                          aaaCAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTc                         >  2:3678284/1‑150 (MQ=255)

                                                                                                                           aaCAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTcc                        <  1:1393019/150‑1 (MQ=255)

                                                                                                                           aaCAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTcc                        >  2:3997651/1‑150 (MQ=255)

                                                                                                                           aaCAAGGCAATGACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTcc                        <  1:4230371/150‑1 (MQ=255)

                                                                                                                              aaGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCatg                     >  1:3505238/1‑150 (MQ=255)

                                                                                                                               aGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTTCAGATGAATTACAGATTACTTAGCACACGCAAGAATGGGTCGTAAGCTCCaaga                    >  1:535547/1‑150 (MQ=255)

                                                                                                                               aGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCatga                    >  2:3860968/1‑150 (MQ=255)

                                                                                                                                      gACGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCgg             >  1:3068518/1‑150 (MQ=255)

                                                                                                                                       aCGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGt            <  1:239677/150‑1 (MQ=255)

                                                                                                                                        cGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTg           <  1:601819/150‑1 (MQ=255)

                                                                                                                                        cGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATGACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTg           <  1:436928/150‑1 (MQ=255)

                                                                                                                                        cGTTTTTCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTg           <  1:3113506/150‑1 (MQ=255)

                                                                                                                                            tttCGTCTGCTTCCTTGCCGCTCTTGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTACAGATTACTTCGCAAACGCAAGAAAGGGTCGTAAGCTCCATGATGTTCGGTGCGGc       >  1:689069/1‑150 (MQ=255)

                                                                                                                                             ttCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCa      <  2:3040096/150‑1 (MQ=255)

                                                                                                                                             ttCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGTTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCa      <  1:3228718/150‑1 (MQ=255)

                                                                                                                                             ttCGTCAGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCa      >  2:3188816/1‑150 (MQ=255)

AGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCG  >  CP009273/3081550‑3081844

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