BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,269	0	C	T	57.9%	204.4 / inf	964	A95T (GCA→ACA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (210/194); new base T (237/315); total (451/513)
Fisher's exact test for biased strand distribution p-value = 5.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 981 total aligned reads displayed.

TGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGT  >  CP000730/36129‑36409

tGCTTCCTTAGTTGATGGTGCCTTATCTTTTATTACCTTTTGAGGTTTATCGGAGTTGTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGc                                                                                                                                              <  1:494199/141‑1 (MQ=11)

  cTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCAt                                                                                                                                            <  1:534818/141‑1 (MQ=32)

  atACCTTAGTTGAAGGTGCCTGATCGGTAATTACCTTTTGAGGTTTACCAAATTGTTTAGTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCAt                                                                                                                                            <  1:26839/140‑1 (MQ=12)

   tgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTACTGTATGTCCCTCTGTATc                                                                                                                                           >  1:175877/3‑141 (MQ=11)

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

      cTTCGTTGAAGGTGCCTTATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAAt                                                                                                                                        <  2:711188/141‑1 (MQ=17)

            tGAGGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGGCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCCACCAAATATCTAATGTATGTACCTCTGTAtgcctggaaca                                                                                                                                  >  2:624481/1‑131 (MQ=11)

             taaGGTGCCTGATCTGTAATTACCTTTCGCGGTGTATTAATTTCCTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGa                                                                                                                                 <  1:735417/140‑1 (MQ=11)

              aaGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGat                                                                                                                                >  1:511685/1‑141 (MQ=12)

                ggTTCCTTATCTCTAAATACCTTTTGAGGCTTACCTAATTGTTTACTGAGCCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                                                                                                                              <  1:239442/141‑1 (MQ=11)

                ggTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGatat                                                                                                                              >  2:515559/1‑141 (MQ=32)

                    ccTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt                                                                                                                          >  2:748408/1‑141 (MQ=17)

                    ccTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAAt                                                                                                                          <  1:470868/141‑1 (MQ=32)

                      tGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAg                                                                                                                        <  2:352570/141‑1 (MQ=32)

                       gATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGCTAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGc                                                                                                                       <  2:554907/141‑1 (MQ=17)

                        aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCt                                                                                                                      <  2:310238/141‑1 (MQ=32)

                        aTCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt                                                                                                                      <  2:377997/141‑1 (MQ=12)

                          cTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTcc                                                                                                                    >  2:45476/1‑141 (MQ=32)

                               actACCTTTTGAGGCTGACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATGTCATGTATGTCCCTCTGTATCAATGGCACTGTTTAAATAGCTTCAtttt                                                                                                               >  1:391407/3‑141 (MQ=11)

                                ttACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTCTCAATGGCACGATATAACTAGCTCCACTTTc                                                                                                              >  2:367592/1‑141 (MQ=255)

                                ctACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc                                                                                                              <  2:471129/140‑1 (MQ=17)

                                  tccTTTTGAAGTTTACCAAATTTTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt                                                                                                            <  1:685057/140‑1 (MQ=12)

                                    cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAt                                                                                                                  <  1:464660/133‑1 (MQ=17)

                                    cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATATCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCttt                                                                                                          >  2:13807/1‑141 (MQ=17)

                                     ttttGAGGCTTACCAAATTGTTTAATGAGACGTTTAACAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtttt                                                                                                         <  1:410155/141‑1 (MQ=255)

                                      aaagagaTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTa                                                                                                        >  2:484344/4‑141 (MQ=17)

                                          aGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGACATAAAAAGCTCCATTTTCCTTTTAtttt                                                                                                    >  2:658864/1‑141 (MQ=12)

                                           ggTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAttttt                                                                                                   >  1:780366/1‑140 (MQ=17)

                                            gtgtACAAAATGGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGa                                                                                                  <  1:435891/138‑1 (MQ=12)

                                            gcttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGa                                                                                                  >  1:594532/3‑141 (MQ=17)

                                            gTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTTCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCATTGGCACGATATAAATAGCTCCATTTTCCTTTTAttttca                                                                                                  >  1:166956/1‑139 (MQ=12)

                                               tACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATGACTGGCACGCTATAAATCGCTCCATTTTCCTTTCATTTTGAttt                                                                                               >  2:406437/1‑139 (MQ=11)

                                                 ccAAATTGTTTAATGAGACGTTAGTTATACGCATATGCTGAATGATTATCTCGTTGCTTACGCACCCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTAc                                                                                             <  2:428007/141‑1 (MQ=255)

                                                    aaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGtc                                                                                          >  2:326028/1‑141 (MQ=17)

                                                    aaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATAACTAATGTATGTCCCTCTACATCACTGGCACGATATAAAAAGCTCCAATTCACTTTTATTTTGATGGACGtc                                                                                          >  2:693368/1‑141 (MQ=11)

                                                     aTTGTTTAATGAGACGTTTGATAAAGGCATATGCTGAATGATTATCTCGTTGCTTACGCCACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATCAATAGACCCATTTTCCTTTTATTGTAATGTAcatgt                                                                                         >  2:762107/1‑137 (MQ=11)

                                                     aTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATTTAACTACCTCCATTTTCCTTTTATTTTGATGTATGtct                                                                                         >  1:607424/1‑141 (MQ=255)

                                                      ttGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTTATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTCTTTTGATTTAGGtttt                                                                                        >  1:138801/1‑138 (MQ=11)

                                                      ttGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTTTTTTGATGTACGtctc                                                                                        >  2:23867/1‑141 (MQ=17)

                                                      ttGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCGCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGtctc                                                                                        <  1:220236/141‑1 (MQ=17)

                                                        gTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcat                                                                                      <  2:732653/141‑1 (MQ=32)

                                                          ttAATTAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcatca                                                                                    <  1:91293/141‑1 (MQ=17)

                                                          ttAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcatca                                                                                    <  1:317324/141‑1 (MQ=32)

                                                          ttAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcatca                                                                                    >  1:528108/1‑141 (MQ=32)

                                                           tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTGTCCTTTTATTTTGATATACGTCTCATCaa                                                                                   >  2:443583/1‑141 (MQ=12)

                                                           tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCaa                                                                                   <  1:128021/141‑1 (MQ=32)

                                                             attagccGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCACCCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATa                                                                                 <  1:203992/135‑1 (MQ=11)

                                                              tGAGACGTTCTCTAAAAGGATATGCTGAACGACTATCTCGTTGCTTACGCACCCAAATAGCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATAc                                                                                <  2:668995/141‑1 (MQ=11)

                                                               gagaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATAcc                                                                               >  2:212893/1‑140 (MQ=12)

                                                               gagaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACg                                                                               <  2:721903/141‑1 (MQ=32)

                                                                 gaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGcc                                                                             >  2:600653/1‑141 (MQ=32)

                                                                     tttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCATTTTATTTTGATGTACGTTTAATTAATACGCCAttt                                                                         >  1:716214/1‑141 (MQ=255)

                                                                      ttaatgaaCGCGTATCCGTAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCTTTTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTg                                                                        <  2:22945/135‑1 (MQ=11)

                                                                      ttGTTAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTg                                                                        <  2:609876/141‑1 (MQ=12)

                                                                      ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGAGATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTg                                                                        >  2:85203/1‑141 (MQ=17)

                                                                        gATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATAAGCGTCTCATCAATACGCCATTTGta                                                                      >  2:283332/1‑141 (MQ=255)

                                                                        gATAAACGCATATACTGAATGATTATCTCGTTGAAATTAACACCAAATAAATTGTTTTTGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATGAATGCGCCATTTGta                                                                      >  1:379848/1‑141 (MQ=2)

                                                                          tAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGtcaa                                                                    >  2:501259/1‑138 (MQ=12)

                                                                          tAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaat                                                                    <  1:277299/141‑1 (MQ=32)

                                                                           aaaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACTAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACTCCATTTGtaata                                                                   <  1:557732/141‑1 (MQ=18)

                                                                           aaaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaata                                                                   >  1:757729/1‑141 (MQ=32)

                                                                           aaaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaata                                                                   <  2:193673/141‑1 (MQ=32)

                                                                            aaCGGATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaataa                                                                  <  1:292618/141‑1 (MQ=18)

                                                                            aaCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaataa                                                                  <  1:362924/141‑1 (MQ=33)

                                                                                    tGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCtttttt                                                          >  2:786483/1‑141 (MQ=255)

                                                                                    tGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCtttttt                                                          <  1:117733/141‑1 (MQ=255)

                                                                                       tGAATGATCATTTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATg                                                       <  1:255444/141‑1 (MQ=255)

                                                                                         aaTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCt                                                     >  1:48108/1‑141 (MQ=255)

                                                                                          aagATTATCTCGTGGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtt                                                    <  2:339400/139‑1 (MQ=255)

                                                                                             aTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCttttt                                                 >  2:573401/1‑141 (MQ=40)

                                                                                              ttattccGTTTCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTGATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttc                                                <  1:151350/135‑1 (MQ=14)

                                                                                               tATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACACCATTTGTACTCAGCTTTTTTTTGTTTTttct                                               >  1:702947/1‑141 (MQ=11)

                                                                                               tATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATACATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTGGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttct                                               <  2:107536/141‑1 (MQ=14)

                                                                                               tATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttct                                               <  2:397366/141‑1 (MQ=40)

                                                                                                ggggcgtgtcgtaagCAACGAAATATCTTATGGATGTCCCTCTGTATCAATCGCACTACATAAATTGCTCCATTTTTCTTTTATTTTGATATACGTCTCACCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttctt                                              <  2:150433/127‑1 (MQ=255)

                                                                                                aTCTCGTTGCTTACGCACCCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttctt                                              <  1:624012/141‑1 (MQ=255)

                                                                                                 tctcGTTGCTGACGCAACCAAATATCTAATGTAGGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTGTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                             <  1:103944/141‑1 (MQ=255)

                                                                                                  ctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTcc                                            >  2:159088/1‑141 (MQ=255)

                                                                                                         cTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTGATCATCAGCTTTTTTATGCTTTTTCTTCCAACTTTg                                     >  2:124833/1‑141 (MQ=11)

                                                                                                            aCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTAtt                                  >  2:177122/1‑140 (MQ=255)

                                                                                                             cGCAACCAAATATCCTGTGTATGTCCCTCTTGATCAACGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTCATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                 <  2:526183/141‑2 (MQ=14)

                                                                                                                 aCCAAATATCTAATGTATGTCCCTCTGTATGAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaaaa                             <  1:451544/141‑1 (MQ=255)

                                                                                                                 aCCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACTTCTCATAAATATGCGATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCTAATTTGATTTaaaa                             >  1:81515/1‑141 (MQ=11)

                                                                                                                  ccAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAt                            >  2:642752/1‑141 (MQ=255)

                                                                                                                   cAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTTTACGTCTCATCAATACTCAATTTGTAATAAGATTTTTTATGCTTTCTCTTCCAAGTTTGATACAACAtt                           >  1:57802/1‑141 (MQ=14)

                                                                                                                   cAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTATTCCAAATTTGATATAAAAtt                           >  1:630725/1‑141 (MQ=255)

                                                                                                                      atatCTAATGTATGTCCCTCTCCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTCATGTACGTCTCTTCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGaatgtaaattcgg                        >  2:541235/1‑129 (MQ=14)

                                                                                                                       tatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTTCGTCTCATCATTACGCCATTTGTAATAAGCTTTTTTATACTTTTTTTTCCAACTGTGATATAAAttttaca                       >  1:673300/1‑137 (MQ=11)

                                                                                                                          cTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCat                    <  1:338578/141‑1 (MQ=25)

                                                                                                                            aaTGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTTATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCatat                  >  1:659089/1‑141 (MQ=255)

                                                                                                                            aaTGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATTCTTTTTCTTCCAAATATGCTATAAAGTTGGGGCatat                  >  1:586198/1‑141 (MQ=11)

                                                                                                                            aaTGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCatat                  <  1:305171/141‑1 (MQ=25)

                                                                                                                             atgtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACTCAATTTGTAATAAGCTTTTTTATGTTTTTTCTTCCAAATTTGATGTAAAATTGGGGCATAtt                 >  1:64882/1‑141 (MQ=11)

                                                                                                                             atgtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTTTGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATAtt                 >  1:540227/1‑141 (MQ=255)

                                                                                                                             atgtaagtCCCCCTGGATCAATGGCACGATATAAATAGCTCCATTTACCTTTCCTTTTGATATACATCCCATCAATACGCCATTTGTAATACGCTTTTTCATCCTTTTTCTCTCAAACTTGATATAAGCTTAGCGCATAtt                 >  2:266322/1‑141 (MQ=2)

                                                                                                                              tgtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATAttt                >  2:250824/1‑140 (MQ=255)

                                                                                                                                     cccTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCCATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGGTTTAAAATTGGGGCATATTGTTGAAcc         >  1:191467/1‑141 (MQ=11)

                                                                                                                                      ccTCTGTATCAATGGCACGATATAAATCGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGGTATAAAATTGGGGCATTTTCTTGAAccc        >  2:298344/1‑141 (MQ=11)

                                                                                                                                      ccTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAAGTACGTCTCATCAATACGCCATTTTTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATGTTCTTGAAccc        >  1:346494/1‑141 (MQ=11)

                                                                                                                                           gCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg   >  2:585199/1‑141 (MQ=25)

TGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGT  >  CP000730/36129‑36409

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