BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,502,324	C→A	9.6%	T429T (ACC→ACA)	ydcP →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,502,324	0	C	A	9.6%	928.1 / 55.5	251	T429T (ACC→ACA)	ydcP	putative peptidase
Reads supporting (aligned to +/- strand): ref base C (117/110); new base A (14/10); total (131/120)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 255 total aligned reads displayed.

ATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTG  >  CP009273/1502176‑1502471

aTTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCAcc                                                                                                                                                     <  1:1740771/149‑1 (MQ=255)

 tttGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCa                                                                                                                                                    >  1:3109392/1‑149 (MQ=255)

 tttGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCa                                                                                                                                                    <  2:5149471/149‑1 (MQ=255)

  ttGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAg                                                                                                                                                   <  1:2188261/149‑1 (MQ=255)

     cacaAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACACTCACAAGTGa                                                                                                                                                >  2:3860207/1‑149 (MQ=255)

      acaAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGaa                                                                                                                                               <  1:986575/149‑1 (MQ=255)

        aaaaTTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAaga                                                                                                                                             <  1:3617329/149‑1 (MQ=255)

          aaTTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAgg                                                                                                                                           <  1:4838675/149‑1 (MQ=255)

             tCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGgtg                                                                                                                                        >  1:2318104/1‑149 (MQ=255)

              cGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGgtgt                                                                                                                                       <  1:21253/149‑1 (MQ=255)

                 ccACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAg                                                                                                                                    <  1:2535939/149‑1 (MQ=255)

                   acaTCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCa                                                                                                                                  >  1:5214076/1‑149 (MQ=255)

                   acaTCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCa                                                                                                                                  >  1:2592812/1‑149 (MQ=255)

                   acaTCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCa                                                                                                                                  <  1:1618567/149‑1 (MQ=255)

                   acaTCACCCAATAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAATACCTCCAGCGTACGTCGGGTGGCGGTAGACATTGAACTGGGAGGCTGGGAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCa                                                                                                                                  <  2:4081915/149‑1 (MQ=255)

                    caTCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCAt                                                                                                                                 <  2:4634522/149‑1 (MQ=255)

                     aTCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCTGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATc                                                                                                                                <  2:1393726/149‑1 (MQ=255)

                          ccACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCa                                                                                                                           >  1:2683388/1‑149 (MQ=255)

                          ccACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCa                                                                                                                           >  1:1286328/1‑149 (MQ=255)

                           cACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCAt                                                                                                                          >  1:1052967/1‑149 (MQ=255)

                             cTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATAc                                                                                                                        <  1:4918140/149‑1 (MQ=255)

                              tAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACg                                                                                                                       >  1:509510/1‑149 (MQ=255)

                               aaaCCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGc                                                                                                                      >  2:2322768/1‑149 (MQ=255)

                                 aCCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTg                                                                                                                    >  1:4512309/1‑149 (MQ=255)

                                   cGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCTGGAACAACTGTTTCTGACCATCACCAGTGAAGAGGGTGTCAGCATCACGCCTACGCTGGa                                                                                                                  >  2:2847217/1‑149 (MQ=255)

                                   cGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGa                                                                                                                  <  2:3963515/149‑1 (MQ=255)

                                    gTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGAc                                                                                                                 <  2:4164795/149‑1 (MQ=255)

                                     tAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACg                                                                                                                >  2:5173119/1‑149 (MQ=255)

                                      aaTCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACgg                                                                                                               <  2:2379941/149‑1 (MQ=255)

                                       aTCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACggg                                                                                                              >  1:4646497/1‑149 (MQ=255)

                                        tCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGc                                                                                                             <  1:3993189/149‑1 (MQ=255)

                                            gATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGtt                                                                                                         <  1:1590855/149‑1 (MQ=255)

                                               cATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTcga                                                                                                      <  1:671952/149‑1 (MQ=255)

                                               cATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTcga                                                                                                      <  1:671714/149‑1 (MQ=255)

                                                aTAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTcgac                                                                                                     >  2:925688/1‑149 (MQ=255)

                                                aTAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACAAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTcgac                                                                                                     <  2:5243556/149‑1 (MQ=255)

                                                 tAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTcgacg                                                                                                    <  2:4302514/149‑1 (MQ=255)

                                                  aaCTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTcgacga                                                                                                   >  1:4613967/1‑149 (MQ=255)

                                                   aCTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGaa                                                                                                  >  1:3013030/1‑149 (MQ=255)

                                                       gcagcaGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCa                                                                                              <  2:1130986/149‑1 (MQ=255)

                                                         agcagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCTGTTCGACGAAGCCAAt                                                                                            <  2:2853785/149‑1 (MQ=255)

                                                         agcagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAAt                                                                                            <  2:1240564/149‑1 (MQ=255)

                                                         agcagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAAt                                                                                            >  1:1383410/1‑149 (MQ=255)

                                                           cagGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATaa                                                                                          <  2:3153282/149‑1 (MQ=255)

                                                             ggCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACg                                                                                        >  1:1012790/1‑149 (MQ=255)

                                                              gCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGc                                                                                       <  2:2525416/149‑1 (MQ=255)

                                                                 cTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGa                                                                                    >  2:3296568/1‑149 (MQ=255)

                                                                 cTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGa                                                                                    >  1:2513717/1‑149 (MQ=255)

                                                                    aCAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGaaaa                                                                                 >  1:3892626/1‑149 (MQ=255)

                                                                     cAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGaaaaa                                                                                <  2:2915191/149‑1 (MQ=255)

                                                                     cAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACAAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGaaaaa                                                                                <  1:4559389/149‑1 (MQ=255)

                                                                       aaaaCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGc                                                                              <  2:2409270/149‑1 (MQ=255)

                                                                        aaaCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCa                                                                             >  1:4790093/1‑149 (MQ=255)

                                                                        aaaCCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCa                                                                             >  1:156333/1‑149 (MQ=255)

                                                                           ccTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATg                                                                          >  1:2532422/1‑149 (MQ=255)

                                                                             tCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGaa                                                                        <  1:3171919/149‑1 (MQ=255)

                                                                               cAGCGAACGTCGGGTGGCGGTAGACATTGACCTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACa                                                                      <  2:2631010/149‑1 (MQ=255)

                                                                               cAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACa                                                                      <  2:4422653/149‑1 (MQ=255)

                                                                                  cGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATc                                                                   <  2:74063/149‑1 (MQ=255)

                                                                                   gAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCt                                                                  <  2:3781323/149‑1 (MQ=255)

                                                                                     aCGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGACGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGa                                                                >  2:1049531/1‑149 (MQ=255)

                                                                                     aCGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGa                                                                >  2:3527849/1‑149 (MQ=255)

                                                                                     aCGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGa                                                                >  1:288448/1‑149 (MQ=255)

                                                                                        tCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAgg                                                             >  2:4815595/1‑149 (MQ=255)

                                                                                         cGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGa                                                            >  2:27907/1‑149 (MQ=255)

                                                                                         cGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGa                                                            <  2:4895665/149‑1 (MQ=255)

                                                                                           ggTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATg                                                          <  2:1861894/149‑1 (MQ=255)

                                                                                            gTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATgg                                                         <  2:1980857/149‑1 (MQ=255)

                                                                                                cGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTg                                                     >  1:2101450/1‑149 (MQ=255)

                                                                                                  gTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGc                                                   >  2:1517936/1‑149 (MQ=255)

                                                                                                    aGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCaa                                                 >  1:1002765/1‑149 (MQ=255)

                                                                                                       cATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAAc                                              <  1:186863/149‑1 (MQ=255)

                                                                                                       cATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACAAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAAc                                              <  1:2937236/149‑1 (MQ=255)

                                                                                                          tGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTgg                                           >  2:1315718/1‑149 (MQ=255)

                                                                                                           gAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTggg                                          <  1:80114/149‑1 (MQ=255)

                                                                                                           gAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTggg                                          <  1:2502740/149‑1 (MQ=255)

                                                                                                            aaCTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTgggg                                         >  2:1849851/1‑149 (MQ=255)

                                                                                                            aaCTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCAAGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTgggg                                         >  1:3932301/1‑149 (MQ=255)

                                                                                                                 ggcggcTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAAc                                    <  1:2737804/149‑1 (MQ=255)

                                                                                                                 ggcggcTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAAc                                    >  2:2113822/1‑149 (MQ=255)

                                                                                                                 ggcggcTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAAc                                    >  1:4769175/1‑149 (MQ=255)

                                                                                                                 ggcggcTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAAc                                    <  1:2737034/149‑1 (MQ=255)

                                                                                                                 ggcggcTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAAc                                    <  1:4142151/149‑1 (MQ=255)

                                                                                                                  gcggcTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAAcc                                   <  1:1135882/149‑1 (MQ=255)

                                                                                                                    ggcTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCct                                 >  1:4180975/1‑149 (MQ=255)

                                                                                                                     gcTGGCAGGAACAACTGATTCTGACCCTCACAAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCctc                                <  1:345077/149‑1 (MQ=255)

                                                                                                                       tGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAATCCCTCTa                              >  2:3412078/1‑149 (MQ=255)

                                                                                                                        ggCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTAt                             >  2:2152363/1‑149 (MQ=255)

                                                                                                                        ggCAGGAACAACTGATTCTGACCCTCACAAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTAt                             >  2:4022395/1‑149 (MQ=255)

                                                                                                                            ggAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACg                         >  1:4439312/1‑149 (MQ=255)

                                                                                                                             gAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGc                        >  2:1514733/1‑149 (MQ=255)

                                                                                                                                   cTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGa                  <  1:2624306/149‑1 (MQ=255)

                                                                                                                                   cTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGa                  <  2:666477/149‑1 (MQ=255)

                                                                                                                                    tGATTCTGACCCTCACAAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGAt                 >  1:4151901/1‑149 (MQ=255)

                                                                                                                                         cTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCa            >  2:4543345/1‑149 (MQ=255)

                                                                                                                                              ccTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAaatta       <  2:4215456/149‑1 (MQ=255)

                                                                                                                                                tCACAAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAaattaat     >  2:181265/1‑149 (MQ=255)

                                                                                                                                                 cACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAaattaatt    <  2:367148/149‑1 (MQ=255)

                                                                                                                                                 cACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAaattaatt    >  2:2026229/1‑149 (MQ=255)

                                                                                                                                                  aCCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAAttt   >  2:653827/1‑149 (MQ=255)

ATTTGCACAAAATTCGTCCACATCACCCACTAAACCGTAATCTTGATCATAACTGGCAGCAGGCACTGACAAAAACCTCCAGCGAACGTCGGGTGGCGGTAGACATTGAACTGGGCGGCTGGCAGGAACAACTGATTCTGACCCTCACCAGTGAAGAGGGTGTCAGCATCACGCATACGCTGGACGGGCAGTTCGACGAAGCCAATAACGCCGAAAAAGCAATGAACAATCTGAAGGATGGTCTGGCAAAACTGGGGCAAACCCTCTATTACGCCCGCGATGTGCAAATTAATTTG  >  CP009273/1502176‑1502471

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