BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,913	C→G	23.9%	S213S (TCG→TCC)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	35,913	0	C	G	23.9%	445.8 / 222.5	321	S213S (TCG→TCC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (101/141); new base G (37/39); total (138/183)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 331 total aligned reads displayed.

CTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTC  >  CP000730/35773‑36053

               aaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGCTc                                                                                                                                >  1:170029/1‑141 (MQ=255)

               aaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATc                                                                                                                                >  1:94181/1‑141 (MQ=255)

                   aaaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga                                                                                                                          <  2:325815/129‑1 (MQ=11)

                    aaatTGTTATCTAATAAAAATAGAAGTCGCTGTGGCAAAGTAACAAAATATAGCTAACCACTAATTTATCATGTCAGTGTCCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                                                                                           <  2:85341/141‑1 (MQ=14)

                    aaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                                                                                           <  1:669591/141‑1 (MQ=255)

                    aaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                                                                                                                         <  1:741271/130‑1 (MQ=11)

                      atTGTGATCTAATAAAAATAGAAGGGTTTGTTGAAAAGAAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                                                                                                                         <  2:273787/141‑1 (MQ=255)

                         gtgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                                                                                                      <  1:359196/141‑1 (MQ=255)

                         gtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                                                                                    >  2:529666/7‑141 (MQ=11)

                         gtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGATAATCACTAATATATAATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                                                                                    >  1:108390/7‑141 (MQ=11)

                           gATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAACTCAGTAGATCTAAAGAGACCTg                                                                                                                    >  2:720003/1‑141 (MQ=255)

                            ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                                                                 <  1:298128/138‑1 (MQ=11)

                            ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                                                                 <  2:54549/138‑1 (MQ=11)

                            ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCGAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                                                                 >  1:748206/4‑141 (MQ=11)

                             tCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                                                                                  <  1:19658/141‑1 (MQ=255)

                              ttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCAGtt                                                                                                               <  1:551527/140‑1 (MQ=11)

                                 ataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                                                                                                            <  1:581494/141‑1 (MQ=11)

                                                                                              cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTCAATACATTCAATACCTTTTAAAGTATGC‑TTTGCCGTATTTATAc                                                >  2:175143/1‑141 (MQ=11)

                                                                                              cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACCCCCAATACCTTTTAAAGTATTC‑TTTCCTGTATTGATAc                                                >  2:199347/1‑141 (MQ=11)

                                                                                              cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATAc                                                <  2:310403/141‑1 (MQ=11)

                                                                                              cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATAc                                                <  2:698898/141‑1 (MQ=11)

                                                                                                gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACtt                                              >  1:588571/1‑141 (MQ=17)

                                                                                                gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACtt                                              <  1:186880/141‑1 (MQ=32)

                                                                                                gggttcctttAACTTGCCAGCATAATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACtt                                              <  1:412826/133‑1 (MQ=11)

                                                                                                 tgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACttt                                             >  2:351340/1‑141 (MQ=17)

                                                                                                 tggtCGCTTAACTTGCTAGCATGATGCTAATGTCGTGGCATGGCGAAAATCCGTAGATCTTAAGAGACCTGCGGTTCTTTTTATATAGTGCGTAAATACACTCCATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACttt                                             <  2:127025/138‑1 (MQ=11)

                                                                                                  tgtCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTg                                            <  2:458465/139‑1 (MQ=12)

                                                                                                  gtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTg                                            >  1:76460/1‑141 (MQ=17)

                                                                                                  gtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTg                                            >  1:458020/1‑141 (MQ=17)

                                                                                                  gtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTg                                            <  1:234466/141‑1 (MQ=17)

                                                                                                   ttCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGa                                           <  1:335445/141‑1 (MQ=11)

                                                                                                   ttCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGa                                           <  2:687565/141‑1 (MQ=32)

                                                                                                    tcgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGAt                                          >  1:136512/5‑141 (MQ=12)

                                                                                                    tCGTTTAACTTGCTAGCATGATGCTAATTGCGTGGCATGGGGAAAAACCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGAt                                          <  1:628231/141‑1 (MQ=11)

                                                                                                    tCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGATACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGAt                                          <  1:320849/141‑1 (MQ=17)

                                                                                                     tgtttAACTTGCTAGCATGTTGCTAATTTCGTGGCATGGGGATAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATa                                         <  2:746290/138‑1 (MQ=11)

                                                                                                     cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATa                                         >  1:12204/4‑141 (MQ=17)

                                                                                                     cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATa                                         >  1:221795/4‑141 (MQ=17)

                                                                                                     cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATa                                         >  2:229078/1‑141 (MQ=17)

                                                                                                     cGCTTAACTTGCTAGCATGATACTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGCGCGTAGATACACGCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGCTACTTTGATa                                         >  1:60962/1‑141 (MQ=11)

                                                                                                     cGCTTAACTTGCTAGCATGAAGCTAATTTCGTGGCATGGCGAAAATGAGTAAGACTGGAGGAACAAGCGGTTCTTTTTATATAGAGAGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATa                                         >  2:717127/1‑141 (MQ=11)

                                                                                                      gtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATat                                        <  2:582762/139‑2 (MQ=255)

                                                                                                       tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATAcc                                       <  1:76057/140‑1 (MQ=17)

                                                                                                        ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATatct                                      >  2:146635/1‑138 (MQ=255)

                                                                                                        ttAACTGGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCt                                      <  1:273362/141‑1 (MQ=17)

                                                                                                         tAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGGATTC‑TTTGCCGTATTGATACTTTGATATCtt                                     >  1:211393/1‑141 (MQ=11)

                                                                                                         tAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACTCTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCtt                                     <  1:637615/141‑1 (MQ=17)

                                                                                                          aaCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTg                                    >  2:499624/1‑141 (MQ=17)

                                                                                                           aCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGt                                   <  2:671155/141‑1 (MQ=17)

                                                                                                            cTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGCTCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGCATTGATACTTTGATATCTTGtt                                  >  2:657855/1‑140 (MQ=11)

                                                                                                            cTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATATCTTGtc                                  <  2:605688/141‑1 (MQ=11)

                                                                                                            cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTGTTTATACTTTGATACCTTGtc                                  >  1:66273/1‑141 (MQ=12)

                                                                                                            cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGtc                                  >  2:472534/1‑141 (MQ=17)

                                                                                                            cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACAAACCATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGtc                                  <  1:148808/141‑1 (MQ=12)

                                                                                                            cTTGCTAGCATGATGCTAATGTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGtc                                  >  1:464780/1‑141 (MQ=17)

                                                                                                             ttGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGtct                                 <  2:510762/141‑1 (MQ=17)

                                                                                                             ttGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGtct                                 <  2:439767/141‑1 (MQ=17)

                                                                                                              tGCTAGCATGATTCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATATCTTGtctt                                <  2:401179/141‑1 (MQ=11)

                                                                                                              tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGtctt                                <  1:344305/141‑1 (MQ=17)

                                                                                                              tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGtctt                                <  1:433529/141‑1 (MQ=17)

                                                                                                              tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGCTTCTTTTTATATATAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGtctt                                <  2:149820/141‑1 (MQ=12)

                                                                                                                cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGtctttc                              <  2:609683/141‑1 (MQ=17)

                                                                                                                cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTACTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGGTACTTTGATACCTTGtctttc                              >  1:244059/1‑141 (MQ=12)

                                                                                                                 tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGtctttct                             <  2:77810/141‑1 (MQ=17)

                                                                                                                    cATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTAc                          >  2:379443/1‑141 (MQ=11)

                                                                                                                    cATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAGCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTA‑TTTGCAGTATTCAAACTTTGATACCTTGTCTTTCTTAc                          >  2:246468/1‑141 (MQ=11)

                                                                                                                    cATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTAc                          >  2:30319/1‑141 (MQ=17)

                                                                                                                      tgatgCTAATTTCGTGGCATGGGGAGAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACtt                        >  1:387062/1‑141 (MQ=11)

                                                                                                                      tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGCCCTACGGTTCTTTTTATATAGAGCGTAAATACATCCACTACCTTTTAAAGTATTC‑TTTCCCGTATTGATACTTGCATTTCTTGTCTTTCTTCCtt                        >  2:116820/1‑141 (MQ=11)

                                                                                                                      tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACtt                        >  1:400598/1‑141 (MQ=255)

                                                                                                                       gattgtgaTTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACttt                       <  2:6277/134‑1 (MQ=11)

                                                                                                                        atgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTa                      >  1:351935/1‑141 (MQ=25)

                                                                                                                        atgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTa                      >  2:281050/1‑141 (MQ=25)

                                                                                                                        atgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTa                      >  2:257377/1‑141 (MQ=25)

                                                                                                                        atgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTAGTCTTTCTTACTTTa                      >  2:210125/1‑141 (MQ=21)

                                                                                                                         tgCTAATTTCGTGGCATGGGGAAAATCCGTAGCTCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTaa                     >  1:16526/1‑141 (MQ=11)

                                                                                                                          gCTAATTTCGTGGCATGGGGAAAATCCGTAGATCCGAAGAGACCTACGGTTCTGTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATCTCTTGTCTTTCTTACTTTAat                    >  2:414963/1‑141 (MQ=11)

                                                                                                                           cTAATTTCGTGGCATGGCGAAAAGGCTAACATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAata                   <  2:535534/141‑1 (MQ=11)

                                                                                                                            tAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTCTTC‑TTTGCCGTATTGATACTTTGATATCTTGTATTTCTTACTTTAatat                  >  1:306942/1‑141 (MQ=11)

                                                                                                                            tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATAACTTGTCTTTCTTACTTTAatat                  >  2:650796/1‑141 (MQ=21)

                                                                                                                            tAATTTCGTGGCATGGCGAAAATCCGTAGATATGAAGCGGTTGGCGGTTATTGTAGTATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatat                  >  1:181217/1‑141 (MQ=11)

                                                                                                                             aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTGTTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                 <  2:565983/141‑1 (MQ=255)

                                                                                                                              aTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGa                <  1:255269/141‑1 (MQ=25)

                                                                                                                                 tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCCTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg             <  1:715953/141‑1 (MQ=255)

                                                                                                                                  ggtagCATGGGGAACATCCGTAGATCTGGAGAAAAATACGGTTCTGTTTATATCGCGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATATCTTGTCTCTCTTACTTTAATATGACGGt            <  1:538047/137‑1 (MQ=14)

                                                                                                                                  cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACATTTTATAG‑ATTCATGTGCTGTATTGATCCTTTTATACCTTGTCTTTCTTACTTTAATAGGACgtt            <  2:734756/141‑3 (MQ=11)

                                                                                                                                   gTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTg           <  1:733059/141‑1 (MQ=255)

                                                                                                                                   gTGGCATGGCGAAAATGCGGAGATCTGAAGAGACCTGCAGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATAC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg           >  1:429928/1‑141 (MQ=11)

                                                                                                                                   gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTGATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGaggttg           >  2:495778/1‑136 (MQ=255)

                                                                                                                                   gTGGCATGGCGAAAATCCGTAGATCCGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg           >  1:100074/1‑141 (MQ=255)

                                                                                                                                     ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAt         >  1:486993/1‑141 (MQ=255)

                                                                                                                                     ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAt         >  1:250474/1‑141 (MQ=38)

                                                                                                                                       cATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCt       >  2:219096/1‑141 (MQ=38)

                                                                                                                                       cATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGATGATCt       <  2:108390/141‑1 (MQ=255)

                                                                                                                                         tGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg     <  1:521522/141‑1 (MQ=255)

                                                                                                                                         tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg     <  1:291779/141‑1 (MQ=255)

                                                                                                                                         tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg     <  2:3328/141‑1 (MQ=38)

                                                                                                                                         tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg     <  1:648401/141‑1 (MQ=255)

                                                                                                                                          ggCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc    <  1:706818/141‑1 (MQ=39)

                                                                                                                                            cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTc  <  1:68455/141‑1 (MQ=255)

CTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTC‑TTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTC  >  CP000730/35773‑36053

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