BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,940	G→A	26.7%	R204R (CGC→CGT)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	35,940	0	G	A	26.7%	800.4 / inf	599	R204R (CGC→CGT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (209/224); new base A (67/92); total (279/320)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 610 total aligned reads displayed.

GTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGT  >  CP000730/35798‑36080

gtattttgataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                                                                                                              <  2:751393/133‑1 (MQ=11)

gtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                                                                                                              <  1:563105/135‑1 (MQ=11)

 tattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                                                                                                             <  1:702764/136‑1 (MQ=11)

          aaaaaTAGAAGGTTCTGTTGCAAAGTAAAAAAATGTGAATAACCACTTATTACAGATTCCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                                                                                                                                      >  1:186228/1‑141 (MQ=11)

                                                                     cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATAc                                                                           >  2:735552/1‑141 (MQ=2)

                                                                     cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAACTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATAc                                                                           >  1:215644/1‑141 (MQ=2)

                                                                     cAGTGCTGGTTTAACTTTCTAGGCTGATGCTAATTTCGTGGCATGGTGAAAATTCGTAGAGCGGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATAc                                                                           <  2:598215/141‑1 (MQ=11)

                                                                      aGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACt                                                                          >  1:385620/1‑141 (MQ=14)

                                                                       gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt                                                                         <  1:172779/141‑1 (MQ=255)

                                                                       gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt                                                                         >  2:591766/1‑141 (MQ=17)

                                                                       gtgtTCGTTTAACTTGCTAGCATGATGCTAATTGCGTGGCATGGCGAAAATAAGTAGATTTGAAGAAACCTGCTGTTCTTTTTATATAGAGCGTAAAGACACTTAATACCTTTTAAAGTATTCTTTGCTGTATTTATACtt                                                                         >  2:716532/1‑141 (MQ=11)

                                                                       acgagcacTTAACTTGCTAGCATGATGCTAATCCCGAGGCATGCCGAAAATCCGTAGATCTGAAGAGACCAGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTAAGCTGAATTGATACtt                                                                         >  2:683354/8‑141 (MQ=11)

                                                                        tgtTCGTTTAACTTGCTAGCATGATGCTAATTTGGTGGCGTGGCGGAAGTCCGTGGATTTGAGGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATAGCTTTTAAGGTTTTCTTTGGTGTATTGTTACttt                                                                        >  2:191976/1‑141 (MQ=11)

                                                                        tgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCGGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTGTTCTTTTCTATTTTGATACttt                                                                        >  2:317199/1‑141 (MQ=11)

                                                                         gtTCGTTTAACTTGCTAGCATGATGCTCATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATACAGCGTACATACATTACATACCTTTTACAGTATTATTTGCCGCATTGATACTTTg                                                                       >  1:161214/1‑141 (MQ=11)

                                                                         gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGGATTCTTTGCTGTATTGATACTTTg                                                                       <  1:575974/141‑1 (MQ=255)

                                                                          ttCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGa                                                                      <  2:618002/141‑1 (MQ=255)

                                                                          ttCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAGTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGa                                                                      >  2:11811/1‑141 (MQ=17)

                                                                           tCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAt                                                                     <  1:344342/141‑1 (MQ=17)

                                                                            cttttAACTTGCTAGCATGATGCTAATTTCGTTGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                                                    <  2:21095/138‑1 (MQ=12)

                                                                            cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGTTa                                                                    >  2:504920/1‑141 (MQ=17)

                                                                            cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                                                    >  1:50485/1‑141 (MQ=17)

                                                                            cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTATTTTTATATAGAGCGTAAATACCTTCAATACCTTTTAAAGTATTATTTGCTGTATTGAAACTTTGATa                                                                    >  2:667402/1‑141 (MQ=255)

                                                                            cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGCCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                                                    >  1:465069/1‑141 (MQ=17)

                                                                              tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATatc                                                                  <  1:634917/140‑3 (MQ=255)

                                                                              cTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                                                  <  2:546589/141‑1 (MQ=32)

                                                                              cTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                                                  <  1:424160/141‑1 (MQ=17)

                                                                              cTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACCCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                                                  >  1:434180/1‑141 (MQ=17)

                                                                               ttAACTGGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCt                                                                 <  2:101856/141‑1 (MQ=17)

                                                                                tAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCtt                                                                >  2:571044/1‑141 (MQ=17)

                                                                                  aCTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGt                                                              >  1:64482/1‑141 (MQ=255)

                                                                                  aCTTGCTAGCAGGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGCTATCTTGt                                                              <  2:694712/141‑1 (MQ=11)

                                                                                   cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtc                                                             <  1:628127/141‑1 (MQ=17)

                                                                                   cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtc                                                             <  2:307825/141‑1 (MQ=17)

                                                                                   cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtc                                                             <  1:129635/141‑1 (MQ=17)

                                                                                     tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt                                                           <  2:294747/141‑1 (MQ=17)

                                                                                     tGCTAGCAGGATGCTAATTTCGTGGCCTGGCCAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGTTGTATGGATACTCTGATACCTTGtctt                                                           <  2:357270/141‑1 (MQ=11)

                                                                                      gCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTTATAtcttttcttt                                                          >  1:370560/1‑141 (MQ=11)

                                                                                      gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtcttt                                                          <  1:226505/141‑1 (MQ=17)

                                                                                        tAGCATGATGCTAATTTCGTTGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                                                        <  1:101260/141‑1 (MQ=17)

                                                                                        tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                                                        <  2:771343/141‑1 (MQ=17)

                                                                                        caGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATAACTTTTATAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                                                        <  2:468944/140‑1 (MQ=255)

                                                                                         aGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATACCTTGtctttctt                                                       <  1:709350/141‑1 (MQ=17)

                                                                                         aGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGACCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttctt                                                       <  2:13021/141‑1 (MQ=17)

                                                                                           cATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTAc                                                     <  2:443698/141‑1 (MQ=255)

                                                                                           cATGATGCTAATTTCGTGGCATGGGGAAAACCCGTCGACATGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTAc                                                     <  1:145876/141‑1 (MQ=11)

                                                                                             tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATATCTAGGAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACtt                                                   >  1:738601/1‑141 (MQ=11)

                                                                                              gatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACttt                                                  <  1:327571/141‑1 (MQ=255)

                                                                                              gatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACttt                                                  <  1:244111/141‑1 (MQ=11)

                                                                                                tgCGAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTAGATACCTTGTCTTTCTTACTTTaa                                                >  2:494199/1‑141 (MQ=255)

                                                                                                   tAATTTCGTGGCATGGCGAAATTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatat                                             <  2:250097/141‑1 (MQ=25)

                                                                                                   tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatat                                             <  1:675118/141‑1 (MQ=25)

                                                                                                    aaTTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATg                                            >  1:451695/1‑141 (MQ=255)

                                                                                                     aTTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGa                                           >  1:317080/1‑141 (MQ=255)

                                                                                                        tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg                                        >  2:148118/1‑141 (MQ=25)

                                                                                                         cGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGt                                       <  1:531407/141‑1 (MQ=14)

                                                                                                         cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTCACTTTAATATGCCGGt                                       <  2:173158/141‑1 (MQ=255)

                                                                                                          gTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGGGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTg                                      <  2:79476/141‑1 (MQ=14)

                                                                                                          gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg                                      >  2:284112/1‑141 (MQ=38)

                                                                                                             gCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGAtt                                   >  1:259945/1‑140 (MQ=255)

                                                                                                             gCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATc                                   <  2:434961/141‑1 (MQ=255)

                                                                                                             gCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATc                                   >  1:86760/1‑141 (MQ=37)

                                                                                                             gCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATc                                   <  1:609307/141‑1 (MQ=255)

                                                                                                             gCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATc                                   <  2:384384/141‑1 (MQ=255)

                                                                                                             gCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGAGTAAAGACATTCAATAGCTTTTAAAGTATTCTCTGCTGTATTGATACTATGATACCGTGTCTTTCTTACTTTAATATGACGGTGATc                                   <  2:633601/141‑1 (MQ=11)

                                                                                                              gaTGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCCATACCTTTTATAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCt                                  <  1:271286/140‑1 (MQ=11)

                                                                                                                tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                <  1:66435/141‑1 (MQ=255)

                                                                                                                 tgggAAAATCCGTAGATCTGAAGAGACCTACTGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc                               <  1:511228/138‑1 (MQ=11)

                                                                                                                 ggggAAAATCCGTACATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATCGATACTTTGATATCTTGTCTTTCTCACTTTAATATGACGGTGATCTTGc                               <  2:655562/138‑1 (MQ=11)

                                                                                                                  ggtaaaaTCCGTAGATCTGAAGAGTCCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATCGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCt                              <  2:639112/138‑1 (MQ=11)

                                                                                                                  gCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCt                              >  1:492080/1‑141 (MQ=255)

                                                                                                                  gCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCt                              >  2:396297/1‑141 (MQ=37)

                                                                                                                  gCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACGTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCt                              >  2:406104/1‑141 (MQ=255)

                                                                                                                   cGAAAATCCGGACATCTGAAGAGACCTTCTGTTCTTTTTATAGAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTc                             <  2:219040/141‑1 (MQ=11)

                                                                                                                    gAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCa                            >  1:691296/1‑141 (MQ=255)

                                                                                                                    gAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCa                            >  1:301879/1‑141 (MQ=255)

                                                                                                                     aaaaTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCaa                           >  1:717212/1‑141 (MQ=255)

                                                                                                                       aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATg                         >  2:330121/1‑141 (MQ=38)

                                                                                                                       aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATg                         <  1:746154/141‑1 (MQ=255)

                                                                                                                         tCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAg                       <  1:789993/141‑1 (MQ=255)

                                                                                                                               gagcTGAAGAGACCAGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCGTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTAtt                 <  1:195721/138‑1 (MQ=255)

                                                                                                                                aTCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATCTGCCGGTGATCTTGCTCAATGAGGTTATTc                >  2:51674/1‑141 (MQ=255)

                                                                                                                                aTCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCCTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTc                <  2:477332/141‑1 (MQ=255)

                                                                                                                                  cTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAg              <  1:686909/141‑1 (MQ=38)

                                                                                                                                   tGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCACTGAGGTTATTCCGa             >  1:358802/1‑141 (MQ=255)

                                                                                                                                   tGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGa             <  2:266487/141‑1 (MQ=38)

                                                                                                                                   tGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATATCTTTTAACGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTCATCTTGCTCAATGAGGTTATTCAGa             <  1:581467/141‑1 (MQ=11)

                                                                                                                                    gtagagACCTGCGGTTATTTTTAGTTTGGTCGTAAATACACTCAATATCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGat            <  2:234245/139‑1 (MQ=11)

                                                                                                                                    gAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGat            >  1:776970/1‑141 (MQ=255)

                                                                                                                                     aaGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACCGTAATCTTCCTCCATGAGGTTTTTccactc           >  2:536868/1‑136 (MQ=11)

                                                                                                                                     aaGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGCCGGTGATCTTGCTCCATGAGGTTATTCACatg           >  2:786417/1‑140 (MQ=11)

                                                                                                                                      agagACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCTGatat          >  1:548129/1‑141 (MQ=255)

                                                                                                                                        agaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAttt        >  2:422015/1‑141 (MQ=255)

                                                                                                                                         gtccTTCGGTGCGTTTTCTATAGAGCGTAAATACATTCCGTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATTACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTc       <  2:344586/139‑1 (MQ=11)

                                                                                                                                         gaCCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTc       <  1:475200/141‑1 (MQ=39)

                                                                                                                                          aCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGCTACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCTGATATTTCg      >  1:284269/1‑141 (MQ=255)

                                                                                                                                          aCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCg      >  1:602022/1‑141 (MQ=255)

                                                                                                                                            cTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGAt    <  1:404615/141‑1 (MQ=255)

                                                                                                                                             tGCGGTTCTTTTTATATAGTGCCTAAATACACTCCATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATg   <  2:370737/141‑1 (MQ=11)

                                                                                                                                             tGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATg   <  1:119359/141‑1 (MQ=38)

GTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGT  >  CP000730/35798‑36080

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