BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,940	G→A	28.0%	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	28.0%	713.9 / inf	552	R204R (CGC→CGT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (176/216); new base A (80/74); total (256/296)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

Only 100 of 561 total aligned reads displayed.

TGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGT  >  CP000730/35799‑36080

tattttaataaAAAT‑GATGGTCCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAAGTTGCTAGCATGATGCTAATTTCGTGGCATGCAGAAAATCCGTTGATCTGGAGAGACCTGc                                                                                                                                             >  1:467471/6‑141 (MQ=11)

tattttaataaAAAT‑GATAGTTCTTTTGCAAAGT‑AAAATATATAGCTAATCACTAATATATCATGGCAGTGTTCGTTTAACTTGCTAGCATGCTGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                                                                                                                                             <  2:270311/136‑1 (MQ=11)

  ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                                                                                           >  1:843767/4‑141 (MQ=11)

   tCTAATAAAGATAGAAGTTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                                                                                                            <  2:798581/141‑1 (MQ=255)

       ataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                                                                                                                                      >  1:325594/1‑141 (MQ=11)

               aGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTCATTTATCATGACAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGCTCTATGtatg                                                                                                                                >  1:111532/1‑140 (MQ=11)

                                                     acgaaTATATCATGGCAGGGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCt                                                                                          <  1:236300/138‑1 (MQ=1)

                                                                    cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTTTATACCTTTTAAAGTATTCTTTGCCGTATTGATAc                                                                           <  2:598976/141‑1 (MQ=11)

                                                                      gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt                                                                         >  1:653364/1‑141 (MQ=255)

                                                                      gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt                                                                         <  2:889815/141‑1 (MQ=17)

                                                                      gtgtTCGCTTAACTTGCTAGCATGATGCTAATGGCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGAGACtt                                                                         <  1:620009/141‑1 (MQ=255)

                                                                       tgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGCAGAGACCTGCGGTTCTTTTTTGATAGAGCGTAAATACAGTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACttt                                                                        >  1:715855/1‑141 (MQ=255)

                                                                        gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTg                                                                       >  2:733276/1‑141 (MQ=255)

                                                                        gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTg                                                                       <  2:367376/141‑1 (MQ=255)

                                                                        gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTg                                                                       <  2:830761/141‑1 (MQ=255)

                                                                        gtTCGTTTAACTTGCTAGCACGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGTTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTg                                                                       <  1:581585/141‑1 (MQ=11)

                                                                         ttCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGa                                                                      <  1:825522/141‑1 (MQ=17)

                                                                         ttCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGa                                                                      <  1:727812/141‑1 (MQ=17)

                                                                          aacgtcaaCTTGATAGCAGCATGTTAGTTTCGAGTCATGCGGAATATCCTTAGATCTTAATAGACCTACGGTTCTTCGTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGAt                                                                     <  2:475467/135‑1 (MQ=2)

                                                                           cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATa                                                                    >  2:343640/4‑141 (MQ=255)

                                                                           cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGCTa                                                                    >  2:210135/1‑141 (MQ=12)

                                                                           cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAATGTATTCTTTGCTGTATTGATACTTTGAtc                                                                    >  2:741461/1‑140 (MQ=17)

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

                                                                                aaCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTg                                                               <  1:302645/141‑1 (MQ=17)

                                                                                 aCTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGt                                                              <  2:240296/141‑1 (MQ=11)

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

                                                                                   ttGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATGCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtct                                                            <  1:26670/141‑1 (MQ=17)

                                                                                   ttGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtct                                                            <  2:153422/141‑1 (MQ=17)

                                                                                    tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACCCTCCATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATTCCTTGtctt                                                           <  2:669501/141‑1 (MQ=255)

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

                                                                                     gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtcttt                                                          >  2:495678/1‑141 (MQ=32)

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

                                                                                       tAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctttct                                                        <  1:517616/141‑1 (MQ=11)

                                                                                       tAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctttct                                                        <  2:512056/141‑1 (MQ=11)

                                                                                        aGCATGATGCTAATGTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttctt                                                       >  2:618049/1‑141 (MQ=17)

                                                                                         gCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTa                                                      >  2:324383/1‑141 (MQ=255)

                                                                                         gCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTa                                                      >  2:424874/1‑141 (MQ=17)

                                                                                         gCATGAGGCTAATTTCGTGTCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTa                                                      <  2:72284/141‑1 (MQ=11)

                                                                                            tgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACTTTGTCTTTCTTCCtt                                                   >  2:584354/1‑141 (MQ=12)

                                                                                            tgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACtt                                                   >  1:370000/1‑141 (MQ=25)

                                                                                            tgatgCTAATTTCGGGGCATGGGGAAAATCCGTAGATCCGAAGAGACCTACGTTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATCCTTTGATATCTTGTCTTTCTTACtt                                                   <  2:450000/141‑1 (MQ=11)

                                                                                             gatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTATTACttt                                                  <  1:571367/141‑1 (MQ=255)

                                                                                               tgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTaa                                                >  2:198426/1‑141 (MQ=34)

                                                                                                tcTAATTTCGTGGCATGGCGAAAATCCGTGGATCTTAAGTGTCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTATAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAat                                               <  2:91972/140‑1 (MQ=11)

                                                                                                 cTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTGCTTTAata                                              >  1:312166/1‑141 (MQ=11)

                                                                                                 cTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAata                                              >  1:857283/1‑141 (MQ=21)

                                                                                                  tAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTCGATATCTTGTCTTTCTTACTTTAatat                                             >  2:186247/1‑141 (MQ=11)

                                                                                                   taTTTCGTGGCTTGGCGAAAATCCGTGGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                                            <  1:205320/140‑1 (MQ=255)

                                                                                                   actttCGTGGCATGGCGAAATTCCGTAGATCTGAAGGGACCTGCGGTTCTTTTTATATAGGGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                                            <  2:68008/139‑1 (MQ=255)

                                                                                                   aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTCTCTTTCTTACTTTAATATg                                            <  1:545267/141‑1 (MQ=25)

                                                                                                    aTTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCCTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGa                                           >  1:469451/1‑141 (MQ=11)

                                                                                                     tttCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGAc                                          >  2:868230/1‑141 (MQ=11)

                                                                                                      ttCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACg                                         <  1:732923/141‑1 (MQ=11)

                                                                                                       tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg                                        <  1:146447/141‑1 (MQ=38)

                                                                                                        cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTCCTTACTTTAATATGACGGt                                       <  1:483004/141‑1 (MQ=255)

                                                                                                         gTGGCATGGGGAAATTCCGTAGATCTGAAGAGACCTACGGTTCTTTTGATATAGAGCTGAAATACATTCAATACCTTTTAAAGTATTCGTTGCCGTAGTGATACTTTGATATCTTGTCTTTCTGACTTTAATATGACGGTg                                      <  1:635849/141‑1 (MQ=11)

                                                                                                         gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATCCTTTGATACGTTGTCTTTCTTACTTTATTATTACGGTg                                      >  2:566424/1‑141 (MQ=11)

                                                                                                         gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg                                      <  1:403990/141‑1 (MQ=255)

                                                                                                            gCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGAtt                                   >  2:383412/1‑140 (MQ=11)

                                                                                                            gCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATc                                   >  2:655228/1‑141 (MQ=11)

                                                                                                            gCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGGTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATc                                   <  1:248093/141‑1 (MQ=255)

                                                                                                             cATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCt                                  >  1:335462/1‑141 (MQ=37)

                                                                                                             cATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCt                                  >  2:751965/1‑141 (MQ=38)

                                                                                                               tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                <  1:77941/141‑1 (MQ=39)

                                                                                                                ggggAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGCTACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc                               >  2:359431/4‑141 (MQ=11)

                                                                                                                ggggAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc                               <  1:800899/138‑1 (MQ=255)

                                                                                                                ggggAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc                               <  1:228171/138‑1 (MQ=255)

                                                                                                                ggCGAAAATCCGTAGATCTGAAGAGCTCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc                               <  2:710639/141‑1 (MQ=255)

                                                                                                                ggCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc                               >  1:673063/1‑141 (MQ=37)

                                                                                                                  cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTc                             <  1:131394/141‑1 (MQ=255)

                                                                                                                  cGAAAATCCGTAGATCCGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTc                             <  1:152091/141‑1 (MQ=255)

                                                                                                                   gAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCa                            >  1:596991/1‑141 (MQ=38)

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

                                                                                                                      aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTCCTCAATg                         >  2:54379/1‑141 (MQ=255)

                                                                                                                       aTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGa                        <  2:309993/141‑1 (MQ=255)

                                                                                                                         ccGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAgg                      >  1:875408/1‑141 (MQ=255)

                                                                                                                          cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGt                     >  2:786893/1‑141 (MQ=255)

                                                                                                                          cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTCGCTCCATGTGGt                     >  2:64428/1‑141 (MQ=255)

                                                                                                                          cGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGt                     <  1:352581/141‑1 (MQ=255)

                                                                                                                            tAGATCTGGAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTTCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGttt                   >  1:101997/1‑140 (MQ=11)

                                                                                                                            tAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGttt                   >  1:191700/1‑140 (MQ=37)

                                                                                                                              gATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTAtt                 <  1:97712/141‑1 (MQ=255)

                                                                                                                              gATCTGAAGAGACCTACGGTTCTTTTTATATAGTGCGTACCTCTTTTGAATTACTTTTCAAGGAATCTTTGCCGTATTGATACTTTGAGATCGTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTAtt                 <  2:877330/141‑1 (MQ=14)

                                                                                                                                    aaGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGatt           >  1:690555/1‑140 (MQ=255)

                                                                                                                                    aaGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGata           >  1:434239/1‑141 (MQ=255)

                                                                                                                                    aaGAGACCTGCGGTTCTGTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATGCTTTGCTGTATTGATGCTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAGTGAGGTTATTCAGata           >  2:1335/1‑141 (MQ=11)

                                                                                                                                    aaGAGACCTGCGGTTCTATTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCCCAATGAGGTTATTCAGata           >  2:763088/1‑141 (MQ=255)

                                                                                                                                    aaGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGata           >  2:848175/1‑141 (MQ=255)

                                                                                                                                     agagACCTGCGGTTCTTTTTATCTAGCGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATATTTTTATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGatat          <  1:174283/141‑1 (MQ=11)

                                                                                                                                     agagACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGatat          >  1:381685/1‑141 (MQ=255)

                                                                                                                                     agagACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGatat          >  2:525290/1‑141 (MQ=255)

                                                                                                                                      gagaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAtt         <  1:763176/141‑1 (MQ=255)

                                                                                                                                       agaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAttt        >  2:250429/1‑141 (MQ=38)

                                                                                                                                        gaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCCGATATTTc       >  1:181565/1‑141 (MQ=255)

                                                                                                                                        gaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTc       >  1:281689/1‑141 (MQ=37)

                                                                                                                                        gaCCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGATGTTATTCAGATATTTc       <  1:451602/141‑1 (MQ=255)

                                                                                                                                        gaCCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTc       >  1:900423/1‑141 (MQ=255)

                                                                                                                                         aCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCg      >  1:553407/1‑141 (MQ=37)

                                                                                                                                          ccTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGa     <  2:809398/141‑1 (MQ=255)

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

TGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGT  >  CP000730/35799‑36080

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