BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,315	G→A	58.3%	Y79Y (TAC→TAT)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,315	0	G	A	58.3%	417.6 / inf	1287	Y79Y (TAC→TAT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (245/278); new base A (365/374); total (624/663)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

Only 100 of 1294 total aligned reads displayed.

TTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36175‑36455

ttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATg                                                                                                                                              >  2:679299/1‑141 (MQ=32)

 tACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGt                                                                                                                                             >  2:937925/1‑141 (MQ=32)

  aCCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTa                                                                                                                                            <  2:893515/141‑1 (MQ=32)

   ccAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATAc                                                                                                                                           >  2:443782/1‑141 (MQ=17)

      aaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGAGATCAATGGCTCCATTTTCCTTTTATTTTGATGTaactc                                                                                                                                        >  2:530140/1‑137 (MQ=255)

       aTTGTTTAATGAGACATTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTTATATACGtct                                                                                                                                       >  1:702441/1‑141 (MQ=12)

         tGTTTAATGAGTCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCCACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCtca                                                                                                                                     <  2:622441/141‑1 (MQ=255)

          gTTTAATGAGACGTTTGCTAAACGCATATGCTGCATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcat                                                                                                                                    <  1:173635/141‑1 (MQ=17)

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

               aTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATa                                                                                                                               >  2:754541/1‑141 (MQ=32)

                 gagaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATTAATACg                                                                                                                             >  1:848442/1‑141 (MQ=17)

                  ggtcGTTTGATAAACGCTTATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                                                                            <  1:745620/138‑1 (MQ=255)

                  aaaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                                                                            <  1:287882/139‑1 (MQ=12)

                      gTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCAtt                                                                                                                        >  2:872445/1‑141 (MQ=17)

                       tttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCAttt                                                                                                                       >  1:263355/1‑141 (MQ=17)

                        ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATCTCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCGtttt                                                                                                                      >  1:366249/1‑140 (MQ=17)

                          gATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTGTACGTCTCATCAATTCGCCATTTGta                                                                                                                    >  1:462074/1‑141 (MQ=17)

                           aTAATCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAATTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaa                                                                                                                   <  1:101587/141‑1 (MQ=17)

                                cGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGc                                                                                                              <  1:606379/141‑1 (MQ=255)

                                cGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCC‑TTTATTTTGATATCCGTCTCATCAATACGCCATTTGtaataccct                                                                                                             >  1:426879/1‑137 (MQ=11)

                                 gCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCt                                                                                                             >  1:656243/1‑141 (MQ=25)

                                 gCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATAACCCATTTGTAATAAGCt                                                                                                             >  2:722895/1‑141 (MQ=21)

                                    tatGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACTCCATTTGTAATAAGCtttt                                                                                                          >  2:590386/1‑141 (MQ=255)

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

                                      tGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCtttttt                                                                                                        <  1:590457/141‑1 (MQ=35)

                                        cTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTAt                                                                                                      >  1:542444/1‑141 (MQ=255)

                                         tGAATGATTATCTCGTTGCTTACGCACCCAAATATCTAATGTATCTCCCACTGCATCGATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATg                                                                                                     <  1:730890/141‑1 (MQ=255)

                                         tGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTTTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCATTACGCCATTTGTAATAAGCTTTTTTATg                                                                                                     <  2:301466/141‑1 (MQ=255)

                                         ggAAGGATTATCTCGTTGCTTACGCAACCAAGTATCTAATGTATGTCCCTCTGTATCAATGTCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATg                                                                                                     <  2:5457/140‑1 (MQ=14)

                                          gAATTATTATCTCGTTGCTTACGCAACCAAAAATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  1:507597/141‑1 (MQ=255)

                                          gAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  1:718219/141‑1 (MQ=255)

                                          gAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    >  1:577855/1‑141 (MQ=39)

                                          aaaTGATTATCTCGTTGCTTACGCCACCAAACATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  2:605669/140‑1 (MQ=255)

                                              gATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtttt                                                                                                >  2:359157/1‑141 (MQ=255)

                                               aTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATTGCTCCATTTTCCTTTTATTTTGATGTACGTCTCTTCAATACGCCATTTGTCATAAACTTTTTTTTGCttttt                                                                                               >  1:197830/1‑141 (MQ=14)

                                                ttATCTCGTTGCTTACGCAACCAAATCTCACACGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAAAACGCCATTTGTAATAAGCTTTTTTATGCTTTttc                                                                                              >  2:391851/1‑141 (MQ=14)

                                                ttATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCCTCAATACGCCATTTTTAATAAGATTTTTTATGcttcttt                                                                                              >  2:781708/1‑140 (MQ=14)

                                                ttATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtttttt                                                                                              >  2:21126/1‑140 (MQ=39)

                                                ttACCTCGTTGCTTACGCAACCAAATTTCTAATGTATGTCCCTCGGCATCAATGGCACGATATAACTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttc                                                                                              <  1:20718/141‑1 (MQ=255)

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

                                                  aTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTAGTGCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttctt                                                                                            >  1:684279/1‑141 (MQ=255)

                                                     tcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCa                                                                                         <  2:196991/141‑1 (MQ=255)

                                                          gCTTACGCTACCAAATATCTAATGTATTTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAttt                                                                                    <  1:476833/141‑1 (MQ=255)

                                                           cTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAtttt                                                                                   >  2:51629/1‑140 (MQ=40)

                                                            ttACGCAGCCAAATATCTAATGCATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTCCTTCCAAATTTGa                                                                                  >  2:689089/1‑141 (MQ=14)

                                                              aCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                                                                <  2:810684/141‑1 (MQ=255)

                                                              aCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                                                                >  2:6366/1‑141 (MQ=255)

                                                              aCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCTTTTTCCTTTTATTTTGATGTTCGTCTAATTAATACGCCAGTTGTAATAGTCGTTTTTATGGTTTTTCTTGCAAATTTgctc                                                                                >  1:80749/1‑138 (MQ=11)

                                                               cGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAGATTTGATat                                                                               >  2:239623/1‑140 (MQ=255)

                                                                gCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACTTCTCATCAATACGCCATTTGTTATCAGCTTTCTTATGCTTTTTATTACAAAATTGATaaa                                                                              >  1:11951/1‑139 (MQ=11)

                                                                        atatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTggg                                                                      <  2:373709/141‑1 (MQ=255)

                                                                         tatCTAATGTATGTCGCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGAGATAAAATTgggg                                                                     <  2:476255/141‑1 (MQ=11)

                                                                         tatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                                                                     <  1:2364/141‑1 (MQ=255)

                                                                           tCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCTTTTTTTTTTGATATTTGTCTCATCAATGCGCCATTTGTAATAAGCTTTTTTATCCTTTTTCTTACAAATCTGATATGAAATTGTGGCa                                                                   >  1:378091/1‑141 (MQ=11)

                                                                            cTAATGTATGTCCCTCGTGATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCat                                                                  <  1:83509/141‑1 (MQ=11)

                                                                              aaTGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCatat                                                                >  1:556741/1‑141 (MQ=255)

                                                                              aaTGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTATGTCTCTTCAATACGCCCTTTTTTATAAGCTTTTTTATGCTTTTTCTTCCAAACTTGATATAAAGTTCGGGCatat                                                                >  1:310058/1‑141 (MQ=11)

                                                                                tgtatgtCCCTCTGCATCAATGGCAAGAAATAAATAGCTCCATTGTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTc                                                              >  2:839373/1‑141 (MQ=255)

                                                                                   atgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATATTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTg                                                           <  1:668002/141‑1 (MQ=11)

                                                                                   atgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTg                                                           <  1:618876/141‑1 (MQ=37)

                                                                                    tgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATCTAAAATTGGGGCATATTCTTGa                                                          >  2:147593/1‑141 (MQ=255)

                                                                                     gtCCCTCTGTATCAATGGCACGATATAAATAGCTCCCTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGaa                                                         <  1:784503/141‑1 (MQ=255)

                                                                                      tCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAc                                                        >  2:768628/1‑141 (MQ=255)

                                                                                      tCCCTCTGCATCAATGGCACGATATGAATAGCGCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAAc                                                        >  1:226426/1‑141 (MQ=255)

                                                                                       cccTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTCATAAGCTTTTTTATGCTTTTTCTTCCAAATTTCATATAAAATTTGGGAGTCTTCTTGAAcc                                                       >  2:459870/1‑141 (MQ=11)

                                                                                        ccTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAATATTGGGGAATATTCTTGAAccc                                                      >  2:119451/1‑141 (MQ=11)

                                                                                         ctctGTATCAATGGGACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCTATACGCTATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCTATTTTGATGTAAAATTTGGGTTTGTTCttttattcg                                                     >  2:826037/1‑134 (MQ=11)

                                                                                           ctGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTCTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAAc                                                   <  2:263355/141‑1 (MQ=255)

                                                                                             ttaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTATTTATTTTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                 <  1:410564/139‑1 (MQ=11)

                                                                                                ccattGGCACTATATAAATATCCCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAg                                              <  1:400889/137‑1 (MQ=11)

                                                                                                 cAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGa                                             <  1:153764/141‑1 (MQ=21)

                                                                                                  aaTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCACCGGTAGAc                                            >  2:79503/1‑141 (MQ=18)

                                                                                                  aaTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGAc                                            >  1:894054/1‑141 (MQ=21)

                                                                                                  aaTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGCAGat                                            >  1:928013/1‑140 (MQ=255)

                                                                                                   aTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGAcc                                           >  2:330134/1‑141 (MQ=25)

                                                                                                   aTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTTAACCCAACGGTAGatc                                           >  1:819088/1‑139 (MQ=255)

                                                                                                    tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGatcg                                          <  2:914113/141‑4 (MQ=255)

                                                                                                    tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCg                                          >  1:355605/1‑141 (MQ=255)

                                                                                                      gCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGtt                                        >  2:490359/1‑141 (MQ=255)

                                                                                                        aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTTAACCCAACGGCAGATCGTTGa                                      >  2:659181/1‑141 (MQ=11)

                                                                                                        aCGATATAAATAGCTCCATTTTCCGTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGa                                      <  2:279600/141‑1 (MQ=255)

                                                                                                           atatAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTTCATAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  1:448186/141‑1 (MQ=255)

                                                                                                           atatAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  1:185168/141‑1 (MQ=18)

                                                                                                            tataAATAGGTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatga                                  <  2:180998/141‑1 (MQ=255)

                                                                                                            tataAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatga                                  >  2:418700/1‑141 (MQ=39)

                                                                                                              taAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGActggtg                                >  2:102642/1‑141 (MQ=11)

                                                                                                                    gCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                          >  2:856249/1‑141 (MQ=255)

                                                                                                                      tCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTa                        >  1:635531/1‑141 (MQ=255)

                                                                                                                         aTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAAATATTCTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATCCTTGAACCCAACGGTAGACCGTTGAATCATGAACGTTTacac                     <  1:492935/141‑1 (MQ=11)

                                                                                                                         aTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTacac                     <  1:913769/141‑1 (MQ=255)

                                                                                                                         aTTTTCCTTTTATTTGGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTacac                     <  2:225930/141‑1 (MQ=18)

                                                                                                                          ttttCCTTTTATTTTGATGTACGTCTCATCAATACGCCGTTTGTAATAAGCTTTTTTCTGCTTTGTCTTCCAAATTTGATATAAAATTGGGGCATATGCTTGAACCCAACGGTAGACCGTTGCATGATGAACGTTTAcccc                    >  2:407242/1‑141 (MQ=11)

                                                                                                                              ccTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGt                >  2:738746/1‑141 (MQ=21)

                                                                                                                               cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGtt               >  1:12845/1‑141 (MQ=25)

                                                                                                                               cTTTTATTTTGATATACGTCTCATCATTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt               >  2:565997/1‑141 (MQ=255)

                                                                                                                                 tttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTcc             >  2:161289/1‑141 (MQ=18)

                                                                                                                                   tATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTctct           <  1:79855/141‑4 (MQ=255)

                                                                                                                                    aTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt          <  1:320834/141‑1 (MQ=11)

                                                                                                                                      tttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTaa        <  1:877550/141‑1 (MQ=11)

                                                                                                                                         gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAatat     >  1:943009/1‑141 (MQ=21)

TTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36175‑36455

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