BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,492,977	G→T	13.3%	D24Y (GAT→TAT)	ydcH →	DUF465 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,492,977	0	G	T	13.3%	679.7 / 54.1	219	D24Y (GAT→TAT)	ydcH	DUF465 family protein
Reads supporting (aligned to +/- strand): ref base G (105/84); new base T (13/17); total (118/101)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 221 total aligned reads displayed.

ATTTTTCCCTAAGTTAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAG  >  CP009273/1492830‑1493124

aTTTTTCCCTAAGTTAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCg                                                                                                                                                     <  2:882424/148‑1 (MQ=255)

 tttttCCCTAAGTTAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGa                                                                                                                                                    <  2:2095007/148‑1 (MQ=255)

  ttttCCCTAAGTTAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGAt                                                                                                                                                   >  2:793335/1‑148 (MQ=255)

  ttttCCCTAAGTTAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGAt                                                                                                                                                   >  2:793553/1‑148 (MQ=255)

        cTAAGTTAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAacac                                                                                                                                             >  2:2475485/1‑148 (MQ=255)

        cTAAGTTAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAacac                                                                                                                                             >  1:1259006/1‑148 (MQ=255)

           aGTTAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAAt                                                                                                                                          >  2:818637/1‑148 (MQ=255)

            gTTAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATa                                                                                                                                         <  1:2843916/148‑1 (MQ=255)

             ttAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATaa                                                                                                                                        <  1:1024153/148‑1 (MQ=255)

              tAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATaaa                                                                                                                                       <  1:878894/148‑1 (MQ=255)

              tAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATaaa                                                                                                                                       <  1:235756/148‑1 (MQ=255)

               aaaTGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAAc                                                                                                                                      <  1:2948881/148‑1 (MQ=255)

                  tGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTg                                                                                                                                   <  2:3015939/148‑1 (MQ=255)

                     gTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATc                                                                                                                                >  2:1586717/1‑148 (MQ=255)

                      tAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCa                                                                                                                               >  1:731588/1‑148 (MQ=255)

                           tGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAACTTGATCATGaaa                                                                                                                          <  2:2269994/148‑1 (MQ=255)

                           tGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGaaa                                                                                                                          >  1:2323807/1‑148 (MQ=255)

                            gATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAAt                                                                                                                         <  1:1107253/148‑1 (MQ=255)

                             aTGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAAtt                                                                                                                        >  1:999585/1‑148 (MQ=255)

                               ggTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGc                                                                                                                      >  1:285481/1‑148 (MQ=255)

                                gtgtgtATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGcc                                                                                                                     >  2:54876/1‑148 (MQ=255)

                                  gtgtATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAg                                                                                                                   >  1:270334/1‑148 (MQ=255)

                                   tgtATTTCAGATACACCTTGTCAGCCACTAACAGTGAGTGAGTATGTTTCCATAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGa                                                                                                                  <  1:3150641/148‑1 (MQ=255)

                                   tgtATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGa                                                                                                                  >  1:766590/1‑148 (MQ=255)

                                    gtATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGaa                                                                                                                 <  2:2685812/148‑1 (MQ=255)

                                      aTTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAcag                                                                                                               >  1:3062513/1‑148 (MQ=255)

                                       tttCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAaagg                                                                                                              <  2:3106286/148‑1 (MQ=255)

                                        ttCAGATACACCTTGTCAGCCACTAACAGGGAGTGTGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAaagga                                                                                                             >  2:612408/1‑148 (MQ=255)

                                        ttCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAaagga                                                                                                             >  1:2034697/1‑148 (MQ=255)

                                             aTACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGtt                                                                                                        >  1:2217391/1‑148 (MQ=255)

                                              tACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTc                                                                                                       <  2:803748/148‑1 (MQ=255)

                                               acacCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTcc                                                                                                      >  1:2399408/1‑148 (MQ=255)

                                                    ttGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACgg                                                                                                 <  1:1547701/148‑1 (MQ=255)

                                                     tGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGt                                                                                                >  2:1358092/1‑148 (MQ=255)

                                                     tGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGt                                                                                                >  1:1287614/1‑148 (MQ=255)

                                                       tCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACTAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCg                                                                                              >  2:2949825/1‑148 (MQ=255)

                                                       tCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCg                                                                                              <  2:1567457/148‑1 (MQ=255)

                                                        cAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGa                                                                                             <  1:587361/148‑1 (MQ=255)

                                                         aGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAg                                                                                            >  2:1241406/1‑148 (MQ=255)

                                                            cACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGt                                                                                         <  1:150492/148‑1 (MQ=255)

                                                               tAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACa                                                                                      <  2:1434068/148‑1 (MQ=255)

                                                               tAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACa                                                                                      <  2:2217391/148‑1 (MQ=255)

                                                               tAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACa                                                                                      >  2:2297653/1‑148 (MQ=255)

                                                                aaCAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACaa                                                                                     >  2:871698/1‑148 (MQ=255)

                                                                   aGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGc                                                                                  <  1:2732624/148‑1 (MQ=255)

                                                                    gggAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCg                                                                                 >  1:1856971/1‑148 (MQ=255)

                                                                     ggAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAGCTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCgg                                                                                <  1:90813/148‑1 (MQ=255)

                                                                     ggAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCgg                                                                                >  2:447172/1‑148 (MQ=255)

                                                                     ggAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCgg                                                                                >  2:1404310/1‑148 (MQ=255)

                                                                      gAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGa                                                                               <  2:1193629/148‑1 (MQ=255)

                                                                      gAGTGCGTATGTTTCCAGAATACCGAGATTTAAGATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGa                                                                               <  1:3043934/148‑1 (MQ=255)

                                                                            gTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGt                                                                         <  1:342154/148‑1 (MQ=255)

                                                                            gTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGt                                                                         <  1:1837964/148‑1 (MQ=255)

                                                                             tATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTc                                                                        >  2:578082/1‑148 (MQ=255)

                                                                             tATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTc                                                                        <  2:1305393/148‑1 (MQ=255)

                                                                                  ttCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCAt                                                                   <  1:3000143/148‑1 (MQ=255)

                                                                                     cAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGaa                                                                <  1:434528/148‑1 (MQ=255)

                                                                                      aGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGaaa                                                               >  1:2049084/1‑148 (MQ=255)

                                                                                        aaTACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGaaaaa                                                             <  1:1437376/148‑1 (MQ=255)

                                                                                         aTACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGaaaaaa                                                            <  2:3208256/148‑1 (MQ=255)

                                                                                           aCCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACa                                                          >  2:2025619/1‑148 (MQ=255)

                                                                                            ccGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACaa                                                         <  1:1950383/148‑1 (MQ=255)

                                                                                             cGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACaaa                                                        <  1:1061451/148‑1 (MQ=255)

                                                                                             cGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACaaa                                                        >  1:553705/1‑148 (MQ=255)

                                                                                                gaTTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGc                                                     <  1:2475485/148‑1 (MQ=255)

                                                                                                 aTTTAATATCCCGTCTGAAAAACTAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCt                                                    >  2:2015819/1‑148 (MQ=255)

                                                                                                  tttAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTa                                                   >  2:3295274/1‑148 (MQ=255)

                                                                                                       tatCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGtt                                              >  1:3219784/1‑148 (MQ=255)

                                                                                                       tatCCCGTCTGAAAAACGAAAATCCACGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGtt                                              >  2:2774437/1‑148 (MQ=255)

                                                                                                        atCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTa                                             >  2:1797653/1‑148 (MQ=255)

                                                                                                        atCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTa                                             >  2:202988/1‑148 (MQ=255)

                                                                                                            cGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAag                                         <  1:1914398/148‑1 (MQ=255)

                                                                                                             gTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAaga                                        <  2:740495/148‑1 (MQ=255)

                                                                                                             gTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAaga                                        <  2:2400759/148‑1 (MQ=255)

                                                                                                              tCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAagat                                       >  2:1513271/1‑148 (MQ=255)

                                                                                                                 gAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAagatgag                                    <  1:21783/148‑1 (MQ=255)

                                                                                                                 gAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAagatgag                                    <  2:908078/148‑1 (MQ=255)

                                                                                                                  aaaaaCGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCCTGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGATGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAagatgaga                                   >  1:892749/1‑148 (MQ=255)

                                                                                                                   aaaaCGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAagatgagat                                  >  1:179079/1‑148 (MQ=255)

                                                                                                                      aCGAAAATCCTCGCTTTATGTCCTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCt                               <  1:654232/148‑1 (MQ=255)

                                                                                                                       cGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTc                              >  1:2645537/1‑148 (MQ=255)

                                                                                                                            aTCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAAt                         >  1:1736421/1‑148 (MQ=255)

                                                                                                                                tCGATTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTg                     <  2:1075512/148‑1 (MQ=255)

                                                                                                                                 cGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTgc                    >  1:1709282/1‑148 (MQ=255)

                                                                                                                                   cTTTATGTCCTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTgcag                  <  1:1448999/148‑1 (MQ=255)

                                                                                                                                   cTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACATTGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTgcag                  <  1:1586717/148‑1 (MQ=255)

                                                                                                                                        tGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGga             >  2:1075806/1‑148 (MQ=255)

                                                                                                                                        tGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGga             >  2:1362873/1‑148 (MQ=255)

                                                                                                                                        tGTCCGTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGga             <  1:1448833/148‑1 (MQ=255)

                                                                                                                                         gTCCTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGgag            >  2:2586950/1‑148 (MQ=255)

                                                                                                                                          tCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGgaga           >  2:2775830/1‑148 (MQ=255)

                                                                                                                                          tCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGgaga           <  1:3243283/148‑1 (MQ=255)

                                                                                                                                           ccTTGTTCTATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGagag          >  2:680721/1‑148 (MQ=255)

                                                                                                                                            cTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGc         >  2:3266987/1‑148 (MQ=255)

                                                                                                                                             ttGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCg        <  2:1259197/148‑1 (MQ=255)

                                                                                                                                              tGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGt       <  1:3251578/148‑1 (MQ=255)

                                                                                                                                              tGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGt       <  1:2030128/148‑1 (MQ=255)

                                                                                                                                              tGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGt       <  1:2101858/148‑1 (MQ=255)

                                                                                                                                                ttCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCa     >  1:3303203/1‑148 (MQ=255)

                                                                                                                                                ttCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCa     >  1:2569011/1‑148 (MQ=255)

ATTTTTCCCTAAGTTAAATGAGTAATCTGATGGTGTGTATTTCAGATACACCTTGTCAGCCACTAACAGGGAGTGCGTATGTTTCCAGAATACCGAGATTTAATATCCCGTCTGAAAAACGAAAATCCTCGCTTTATGTCCTTGTTCGATAAACACAATAAACTTGATCATGAAATTGCCAGAAAGGAAGGTTCCGACGGTCGAGGGTACAATGCGGAAGTGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAG  >  CP009273/1492830‑1493124

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