BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,269	C→T	64.2%	A95T (GCA→ACA)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,269	0	C	T	64.2%	522.3 / inf	1018	A95T (GCA→ACA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (178/185); new base T (296/349); total (476/542)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 1026 total aligned reads displayed.

TGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGT  >  CP000730/36129‑36409

tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGc                                                                                                                                               >  2:432937/1‑141 (MQ=32)

 gCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                                                                                                                              <  2:444868/141‑1 (MQ=32)

  ctgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTACATCTGTAt                                                                                                                                             >  1:511633/4‑141 (MQ=11)

     ccTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                                                                                                                          <  2:350616/141‑1 (MQ=17)

      cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCCAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAAt                                                                                                                                         <  2:132733/141‑1 (MQ=17)

      cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAGATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                                                                                                                                         <  1:126211/141‑1 (MQ=17)

              aaGGTGCCTGATCTGTAATTACCGTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGat                                                                                                                                 <  2:614731/141‑1 (MQ=17)

                 gtggcgtaTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGAtata                                                                                                                              <  2:275268/134‑1 (MQ=17)

                 gTGCCTGATCTGTAACTTCCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGAtata                                                                                                                              <  1:271942/141‑1 (MQ=255)

                   gCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCCCTGCATCAATGGCACGAATTaaa                                                                                                                            >  1:33703/1‑141 (MQ=12)

                      tGAACTGGAATTACATTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATGTAAATat                                                                                                                         >  2:173185/1‑140 (MQ=11)

                       gATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTCTGTCCCTCTGCATCAATGGCACGATATAAATAGc                                                                                                                        >  1:110883/1‑141 (MQ=17)

                       gATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCGTCTGTATCAATGGCACGATATAAATAGc                                                                                                                        >  1:463526/1‑141 (MQ=255)

                         tCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTc                                                                                                                      >  2:203892/1‑141 (MQ=12)

                         tCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTc                                                                                                                      >  1:361806/1‑141 (MQ=12)

                          ctgatactACCTTTTGATGCTTACCAAATTGTTTAATGAGACGTTTGTAATATGAATGTTCTGAATGATTATCTCTTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTcc                                                                                                                     <  1:624470/134‑1 (MQ=11)

                          cTGTAATTCCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTcc                                                                                                                     <  2:627492/141‑1 (MQ=17)

                           tGTAACTACCTTTTCAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCa                                                                                                                    <  2:123105/141‑1 (MQ=255)

                           tGTAAATACTTTTTGAGTTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCACGTTACTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCa                                                                                                                    <  2:524309/141‑1 (MQ=11)

                            tttaTTACCTTTTTAGGTTTTCCAAATTGTTTATTGAGACGTTTTATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAt                                                                                                                   <  1:233921/138‑1 (MQ=11)

                             tAATTACCTTTTGAGGTTTACCAAATTGTTTAATGACACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGCTATAAATAGCTCCAtt                                                                                                                  <  1:369615/141‑1 (MQ=17)

                              aaTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAttt                                                                                                                 >  1:158526/1‑141 (MQ=17)

                               aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTTCATCAATGGCACGATATAAATAGCTCCAtttt                                                                                                                >  1:77742/1‑141 (MQ=17)

                               aTTACCTTTTGAGGTTTACCAAAATGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtttt                                                                                                                >  2:476854/1‑141 (MQ=17)

                                ttACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc                                                                                                               <  2:154557/141‑1 (MQ=17)

                                ttACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCCATGGCACGATATAAATAGCTCCATTTTc                                                                                                               >  1:546724/1‑141 (MQ=17)

                                ctACCTTTTGAGGCTTACTAATTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc                                                                                                               <  1:37723/140‑1 (MQ=255)

                                   ccTTTTGAGGTTTACCAAAGTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt                                                                                                            <  2:596285/141‑1 (MQ=17)

                                    cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCttt                                                                                                           <  1:353475/141‑1 (MQ=32)

                                    cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCttt                                                                                                           >  1:93805/1‑141 (MQ=32)

                                    cTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTTCGCAACCAAATATCTAATGTTTGTCCCTCTTTATCAATGGCACGCTATAAAAAGCTCCATTTTCCttt                                                                                                           >  1:259448/1‑141 (MQ=11)

                                     ttttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTGGCTTACGCAACCAAATATCTAATCTATGTCCCTCTGTATCAAGGGCACTATATAAATAGCTCCATTTTCCtttt                                                                                                          <  2:241311/141‑1 (MQ=11)

                                      tttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTTAATGATTCTATAGATTCTTACGCAACCAAATATCTACATTAGGGCCCTCCGTATCAATGGCACGATTTAAAAAGCACCATGTTCCgtttt                                                                                                         >  2:487884/1‑140 (MQ=2)

                                      tttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTa                                                                                                         <  1:401531/141‑1 (MQ=17)

                                      tttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAACGATTATCTCGTTGCTTACGCAACCAAATATCTATTTTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTa                                                                                                         <  2:366326/141‑1 (MQ=12)

                                          aggcttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGCCACGATCTAAATAGCTCCATTTTCCTTTTAtttt                                                                                                     >  2:226750/5‑141 (MQ=255)

                                            gTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCTACCAAATATCTAATGTATGTCCCTCTGCATCATTGGCACGATATAAATAGCTCGATTTTCCTTTTAtttttt                                                                                                   >  1:18013/1‑139 (MQ=12)

                                             tttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAt                                                                                                  <  2:638824/141‑1 (MQ=32)

                                             cttACCAACTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAt                                                                                                  <  1:392731/140‑1 (MQ=12)

                                              ttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAta                                                                                                 >  1:309940/1‑140 (MQ=17)

                                                     aTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAtatagttct                                                                                          >  1:594143/1‑136 (MQ=255)

                                                     aTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGtct                                                                                          <  1:616664/141‑1 (MQ=12)

                                                         ttgactgaagcGTTTGATCAACGCATCGTGTGAATGATTATGGCGTTGCTTACGCAACCAATTATCTAATGTATGGCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCtcatc                                                                                      <  2:580150/131‑1 (MQ=11)

                                                           tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATTCGTCTCATCaa                                                                                    >  1:100662/1‑141 (MQ=12)

                                                           tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCaa                                                                                    >  1:188133/1‑141 (MQ=32)

                                                              tGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑c                                                                                <  2:37196/141‑1 (MQ=17)

                                                                 gaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATG‑CGcc                                                                             >  1:638707/1‑141 (MQ=17)

                                                                  aCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCaaaaagcaa                                                                             >  1:282558/1‑134 (MQ=17)

                                                                      ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTg                                                                        >  2:166355/1‑141 (MQ=17)

                                                                         ttAATCGGATATGCTTAATGTTTATCTCCTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCCAGGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGtaa                                                                     <  2:547455/140‑1 (MQ=11)

                                                                         atatACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACAAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGtaa                                                                     <  2:53314/141‑1 (MQ=17)

                                                                         aTAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGtaa                                                                     >  2:334069/1‑141 (MQ=32)

                                                                           aaaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGtaata                                                                   <  1:248403/141‑1 (MQ=16)

                                                                             aCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAg                                                                 <  2:417585/141‑1 (MQ=255)

                                                                               gCATATGTTGAATGATTATCTCGTTGCTTACGCAACCAAATATTTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCt                                                               <  2:65039/141‑1 (MQ=255)

                                                                                 atatGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCttt                                                             <  2:571307/141‑1 (MQ=25)

                                                                                  tatGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTTTTATAAGCtttt                                                            >  1:474758/1‑141 (MQ=255)

                                                                                  tatGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCtttt                                                            >  2:460939/1‑141 (MQ=34)

                                                                                   atGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCttttt                                                           >  2:359920/1‑141 (MQ=25)

                                                                                    tGCTTAATTATTATCTGGTTGCTTACGTAACCAAATATCTAATGTATGTCCCTCGGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCtttttt                                                          <  1:587165/141‑1 (MQ=11)

                                                                                    tGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATGAGCtttttt                                                          >  1:296806/1‑141 (MQ=12)

                                                                                     gcctaaTGATTATCTCTCTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTa                                                         <  2:317812/137‑1 (MQ=255)

                                                                                      cTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCATTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCACCAATC‑CGCCATTTGTAATAAGCTTTTTTAt                                                        >  1:535239/1‑141 (MQ=11)

                                                                                       tGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATg                                                       >  2:101977/1‑141 (MQ=255)

                                                                                        gAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAGGCTTTTTTATGc                                                      <  1:522253/141‑1 (MQ=255)

                                                                                          aTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCtt                                                    <  2:382697/141‑1 (MQ=255)

                                                                                           tGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCttt                                                   <  1:441646/141‑1 (MQ=255)

                                                                                            gATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACTTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTAtttgttt                                                  >  1:14796/1‑135 (MQ=11)

                                                                                              ttATCTCGTTGCTTACTCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTttc                                                <  1:85006/141‑1 (MQ=255)

                                                                                               tATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTttct                                               <  2:555118/141‑1 (MQ=40)

                                                                                                       tGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAtt                                       <  2:114672/141‑1 (MQ=255)

                                                                                                       tGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACATCTCATCAATA‑AGCCATTTGTAATAAGCTTTATTTTGTATTttcttcttactt                                       >  1:576128/1‑135 (MQ=11)

                                                                                                       tGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATATCTCCATTTTCCTTTTATTTTTATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATTCTTTTTCTTCTAAAtt                                       >  1:663807/1‑141 (MQ=255)

                                                                                                       tGCTTACGCAACCAAAGATCTAATGTATGTCCCTCTGTATCAATGTCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAtt                                       <  1:500931/141‑1 (MQ=255)

                                                                                                         cTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAGttta                                     >  1:566319/1‑140 (MQ=255)

                                                                                                          ttACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATTCTTTTTTTTTCAAttttta                                    >  1:215665/1‑139 (MQ=14)

                                                                                                           tACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTATTCCAAAGTTGAt                                   >  1:195704/1‑141 (MQ=255)

                                                                                                            aCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                  >  2:255872/1‑141 (MQ=255)

                                                                                                            aCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                  >  2:40503/1‑141 (MQ=40)

                                                                                                             cGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                 <  1:203892/141‑2 (MQ=255)

                                                                                                              gcatcaacatgtCTAATGTATGTCCGTCTTGACCGATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATata                                <  1:650730/130‑3 (MQ=14)

                                                                                                              gCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTtatata                                >  2:498676/1‑139 (MQ=255)

                                                                                                                       tatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                       <  1:451948/141‑1 (MQ=255)

                                                                                                                       tatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                       <  1:609665/141‑1 (MQ=38)

                                                                                                                       ggtCTAATGTACGTCCCTCTGCAGCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAATTTTGATATAAAATTgggg                       <  2:538093/139‑1 (MQ=11)

                                                                                                                         tCTAATTGATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCa                     <  2:288934/141‑1 (MQ=255)

                                                                                                                          cTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGTTATAAAATTGGGGCat                    <  1:330806/141‑1 (MQ=255)

                                                                                                                           tAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTATACGTTTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCata                   >  1:383813/1‑141 (MQ=11)

                                                                                                                           tAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCata                   <  2:301793/141‑1 (MQ=255)

                                                                                                                           tAATGTATGTCCCTCTGCATCAATGGCACGATATACATAGCTCCATTTTCCTTTTATTTTGATGTACTTCTCATCAATA‑CGCCATTAGTAATAAGCTTTTTTATGCTTTTTCTTTCAAATTTGATATAAAATTGGGGCatt                   >  1:117354/1‑140 (MQ=11)

                                                                                                                             atgtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAtttt                 >  2:182924/1‑141 (MQ=255)

                                                                                                                                  tgtCCCCCTGTATCCATGGCACGATATAAATAGCTCCTTTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGa            <  1:167056/141‑1 (MQ=11)

                                                                                                                                  tgtCACTCTGTATCAATGGCACGATATAAATAGAACCATCTACCTTTTTTTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAATTTGGGGGATATTCTTGa            >  1:319659/1‑141 (MQ=11)

                                                                                                                                    tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACCACCAAGT‑TGATAAGATTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAc          >  1:484825/1‑141 (MQ=11)

                                                                                                                                     cccTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAcc         <  2:78235/141‑1 (MQ=21)

                                                                                                                                      ccTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCACATTTGATATAAAATTGTGGCATATTCTTGAAccc        >  2:391656/1‑141 (MQ=11)

                                                                                                                                      ccTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAccc        <  1:556823/141‑1 (MQ=21)

                                                                                                                                          tGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGGACCCAACg    >  1:42401/1‑141 (MQ=255)

                                                                                                                                           tcATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg   <  2:659342/140‑1 (MQ=18)

                                                                                                                                           gCATCAATGGCACGATATAAATTGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg   <  2:88435/141‑1 (MQ=18)

TGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGT  >  CP000730/36129‑36409

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`