BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,479	0	T	C	12.0%	257.9 / 25.6	126	T25A (ACA→GCA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (58/51); new base C (8/7); total (67/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 130 total aligned reads displayed.

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATG  >  CP000730/36340‑36615

taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGta                                                                                                                                         <  2:226426/141‑1 (MQ=255)

taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGta                                                                                                                                         >  1:429159/1‑141 (MQ=255)

  ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGgata                                                                                                                                       >  2:630090/1‑141 (MQ=38)

     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTGACACAACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCt                                                                                                                                    >  2:519723/1‑141 (MQ=255)

        ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAg                                                                                                                                 >  1:649056/1‑141 (MQ=37)

            ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATACCTTAGAtag                                                                                                                             >  1:617848/1‑141 (MQ=25)

              tttttttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGCTAACTCA‑ATGTATATCTTAGAtagta                                                                                                                           >  2:563336/5‑141 (MQ=21)

              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGCATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGAtagta                                                                                                                           >  1:868803/1‑141 (MQ=25)

              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGAtagta                                                                                                                           >  1:84650/1‑141 (MQ=35)

              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGAtagta                                                                                                                           >  1:184727/1‑141 (MQ=35)

                gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATATTAGc                                                                                                                         >  2:921500/1‑141 (MQ=21)

                    tttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGccgtt                                                                                                                     <  1:544293/141‑4 (MQ=25)

                       cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                  >  2:489987/1‑141 (MQ=35)

                       cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                  >  2:379338/1‑141 (MQ=35)

                            aaaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAg                                                                                                             >  2:774259/1‑141 (MQ=35)

                                 tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                        >  1:795137/1‑141 (MQ=35)

                                   aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATGTCTCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAGCGGCTACAGTGATaaa                                                                                                      >  1:26528/1‑140 (MQ=255)

                                   aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAAAATGATGCAAAATGAACTTTTAGACCACGTTCCCTCAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                                                      <  1:443933/141‑1 (MQ=11)

                                          aTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                               <  1:289927/141‑1 (MQ=32)

                                            tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCCCGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                             <  2:288636/141‑1 (MQ=17)

                                             ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATCTTTCAGATATATCACGATCACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTCCAGTGATAACATCCTTGTTa                                                                                            >  2:48257/1‑141 (MQ=12)

                                             ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATGACATCCTTGTTa                                                                                            >  2:134463/1‑141 (MQ=17)

                                             ggggCATATTCTTGAACCCAACGGTAGACCGCTGAATGATGAAAGTTTACGCAACGCAACTTTAATACTTTTGTTATATGCTGATAAATTG‑ATCTAAAGCTTAAGGGGTAGCAAACTGCTACAGTTTTAACATCCTTGTTa                                                                                            >  2:785553/1‑141 (MQ=2)

                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                           <  2:84650/141‑1 (MQ=32)

                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaaa                                                                                          >  2:256505/1‑141 (MQ=32)

                                                gCATATTCTTGAACCCAACGGTAGACCGTTGATTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTAC‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                         <  1:630090/141‑1 (MQ=12)

                                                gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                         >  2:653764/1‑141 (MQ=32)

                                                  atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGCGATAACATCCTTGTTAAATTg                                                                                       <  1:286313/141‑1 (MQ=17)

                                                   tatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGt                                                                                      <  2:27748/141‑1 (MQ=32)

                                                    atTCTTGAACCCAACGGTAGGCCGTTGAATGATGAACGTGTACACCACGTTCCCTTAACATTTCAGATATATCACGATCACTCA‑ATGTATATCTTAGAGAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                                     <  1:37470/141‑1 (MQ=11)

                                                      tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCGTTAATATTTCAGATATATCACGATAACTAA‑ATTTATTTTTTAGTTAGTGGCTAACAGCTACAGTGGTAACATACTAGTTAAATTGTTta                                                                                   >  2:287431/1‑141 (MQ=11)

                                                        ttGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtata                                                                                 >  2:571846/1‑141 (MQ=32)

                                                         tGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATCTTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTatat                                                                                >  1:153761/1‑141 (MQ=17)

                                                          tgactcAATGGTAGAACCTAGAAATATGACCGTTCACACCACGTTCTCTTAACATTTCAGATATATCACGATAGATAA‑GTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                               <  1:132173/136‑1 (MQ=11)

                                                          gAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTGTATc                                                                               >  2:894903/1‑141 (MQ=17)

                                                           aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                              <  2:879713/141‑1 (MQ=32)

                                                           aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                              <  1:512041/141‑1 (MQ=32)

                                                            aCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTg                                                                             >  1:819449/1‑141 (MQ=32)

                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                            >  2:538986/1‑141 (MQ=32)

                                                              ccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaa                                                                           <  1:691818/141‑1 (MQ=32)

                                                               cAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                          <  2:367838/141‑1 (MQ=32)

                                                                 aCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATa                                                                        >  2:152787/1‑141 (MQ=32)

                                                                   ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                      <  1:530140/141‑1 (MQ=32)

                                                                     tAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTc                                                                    >  1:806273/1‑141 (MQ=32)

                                                                       gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAACTCGTTCAt                                                                  >  2:752821/1‑141 (MQ=17)

                                                                       gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                  >  2:817193/1‑141 (MQ=32)

                                                                       gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                  >  1:755353/1‑141 (MQ=32)

                                                                         tcGTTGAAGGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                <  2:701723/140‑1 (MQ=12)

                                                                         ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGATACAGTCATACCATCCTTGTTAAATTTGTTATATTTTAAATAGTTCATAc                                                                >  2:267128/1‑141 (MQ=11)

                                                                             tGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                            >  2:100273/1‑141 (MQ=32)

                                                                             tGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                            <  1:758136/141‑1 (MQ=32)

                                                                             tGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                            >  1:224502/1‑141 (MQ=32)

                                                                                atgattAACGTGTACACCAGGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAc                                                         <  2:410920/141‑1 (MQ=12)

                                                                                  gatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTc                                                       <  1:418700/141‑1 (MQ=32)

                                                                                    tgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATGTTAGATAGTAGCCAACGGCTACAGTGATAACATCCGTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAATCCt                                                     >  2:113421/1‑141 (MQ=12)

                                                                                      taCGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                   <  2:411153/140‑1 (MQ=17)

                                                                                      aaCTTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                   <  2:80749/141‑1 (MQ=17)

                                                                                        cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGATGTTGTTCAGACAGAAGACTCCttttt                                                 <  2:669330/141‑1 (MQ=255)

                                                                                        cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttttt                                                 <  2:73262/141‑1 (MQ=32)

                                                                                         gTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                                >  1:18119/1‑141 (MQ=32)

                                                                                           ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCGTACAGAAGGCTCCTTTTTGtt                                              <  1:489987/141‑1 (MQ=17)

                                                                                           ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                              <  1:379338/141‑1 (MQ=32)

                                                                                            tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCCTACAGAAGCCTCCTTTTTGTTa                                             >  2:303268/1‑141 (MQ=17)

                                                                                            tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAACCAAGTTAAATTGTTTATATCTGAACGAGTTCGTACAGAAAACTCCTTTTTGTTa                                             >  2:150523/1‑141 (MQ=11)

                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                            >  1:858728/1‑141 (MQ=32)

                                                                                                   cgttatcTTAATATTACGGATATATCACGAAAGATCAGTTGTA‑AACTTAGATTGTGGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATtat                                      >  1:452373/7‑141 (MQ=11)

                                                                                                   cGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTTTTTATATTTGACATAGTTCATAAAGACGACTACTTTTTGTTTAAATtat                                      >  1:523321/1‑141 (MQ=11)

                                                                                                       cccTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtt                                  >  1:810710/1‑140 (MQ=17)

                                                                                                        ccTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTCTATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                 >  2:617021/1‑141 (MQ=17)

                                                                                                        ccTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                 >  1:137314/1‑141 (MQ=32)

                                                                                                        ccTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                 >  1:377630/1‑141 (MQ=32)

                                                                                                          ttAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAGGACTCCTTTTTGTTAACATTATACTCTaa                               >  1:74860/1‑141 (MQ=12)

                                                                                                             atatTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                            <  2:367354/141‑1 (MQ=35)

                                                                                                             atatTTCAGATATAACACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                            <  2:120770/141‑1 (MQ=25)

                                                                                                              aatTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAGTATACTATAAATTc                           <  2:641276/140‑1 (MQ=25)

                                                                                                                  actgATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGTAGACTCCTGTTTGTTAAAATTATACTATAAATTCAACt                       <  1:187929/138‑1 (MQ=255)

                                                                                                                     gATATATCACGATAGCTCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAATTTCAACTTTg                    <  1:683885/141‑1 (MQ=255)

                                                                                                                      atatatCACGATAGCTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGCGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGc                   <  1:587874/141‑1 (MQ=255)

                                                                                                                      atatatCACGATAGCTCA‑ATGGATCTCATGGAGAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGc                   >  1:456209/1‑141 (MQ=11)

                                                                                                                        atatCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATATCATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCaa                 <  1:538986/141‑1 (MQ=25)

                                                                                                                         tatCACGATAGCTCA‑ATGTATATATTAGTTAGTAGCCAACGGCTACGGTGATAACATCCTTGTTAAATTGTTGATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCGACTTTGCAAc                <  2:184198/141‑1 (MQ=11)

                                                                                                                         tatCACGATAACTCA‑ATGTAAATTGTAGCTAGTATCCAACGGCTACAGTGATAACATCCTCGTTAAATTGTTTATACCTGAAATATTTCATACAGAAGACTCCTTTTTGTTAAAATTCTACTATACATTCAACTTTGCaca                >  2:435309/1‑139 (MQ=11)

                                                                                                                            cACGATAGCTCA‑ATGAATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGa             <  2:66054/141‑1 (MQ=12)

                                                                                                                            cACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCATTTTTGTTAAAATTATACTATAAATTCAACTTTCCacccga             >  1:116292/1‑136 (MQ=255)

                                                                                                                                 ttgctta‑aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        <  2:582/135‑1 (MQ=255)

                                                                                                                                 tagcTCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        <  2:727138/138‑1 (MQ=12)

                                                                                                                                 tagcTCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        <  1:868682/138‑1 (MQ=9)

                                                                                                                                 aagccta‑aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        <  2:324286/135‑1 (MQ=255)

                                                                                                                                  agcTCA‑ATTCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta       <  2:598123/139‑1 (MQ=255)

                                                                                                                                  agcTCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta       >  1:244209/3‑141 (MQ=9)

                                                                                                                                   aCTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat      >  2:828061/1‑141 (MQ=38)

                                                                                                                                   aCTCA‑ATGTATATCTAAAATAGTAGCCAACGGCTACAGTGATCACATCTTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat      <  1:880897/141‑1 (MQ=255)

                                                                                                                                    cTCA‑ATGCATATGTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt     >  1:665280/1‑141 (MQ=255)

                                                                                                                                    cTCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTTTAATTTTATTATAAATTCAACTTTGCAACAGAACCGtatt     >  2:882925/1‑141 (MQ=11)

                                                                                                                                     tCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta    <  2:502065/141‑1 (MQ=255)

                                                                                                                                     tCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTATTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta    <  2:502011/141‑1 (MQ=255)

                                                                                                                                     tCA‑ATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt    >  1:77696/1‑140 (MQ=16)

                                                                                                                                      cA‑ATGTATAGCGTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtattat   <  1:823997/141‑1 (MQ=255)

                                                                                                                                       a‑aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttt  >  2:469023/1‑138 (MQ=255)

                                                                                                                                       a‑aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATACAACTTTGCAACAGAACCGtatttta  >  2:583813/1‑138 (MQ=255)

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCA‑ATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATG  >  CP000730/36340‑36615

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