BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,315	G→A	54.6%	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	54.6%	71.7 / inf	963	Y79Y (TAC→TAT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (188/236); new base A (244/272); total (446/517)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

Only 100 of 971 total aligned reads displayed.

TTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36175‑36455

ttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCACTTTCCTTTTATTTTGATg                                                                                                                                              >  1:742884/1‑141 (MQ=17)

     aaaTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGt                                                                                                                                         >  2:398264/1‑141 (MQ=32)

     aaaTTGTTTAATGACACGTTTGATAAAGGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGt                                                                                                                                         <  1:913286/141‑1 (MQ=255)

       aTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACGATATAGCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGtct                                                                                                                                       >  2:850920/1‑141 (MQ=11)

       aTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGCATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATTTACTtcc                                                                                                                                       >  1:186793/1‑140 (MQ=255)

          gTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCtcat                                                                                                                                    <  2:749731/141‑1 (MQ=12)

             tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCaa                                                                                                                                 >  2:973779/1‑141 (MQ=32)

                  agaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                                                                            <  1:667602/141‑1 (MQ=17)

                    aCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTATCATCAATACACCa                                                                                                                          >  2:1027799/1‑141 (MQ=17)

                       tttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCTTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATAAGcacttt                                                                                                                       >  1:163330/1‑136 (MQ=11)

                       tttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAGATACGTCTCATCAATACGCCAttt                                                                                                                       >  2:1013989/1‑141 (MQ=12)

                        ttGATAAACGCATATGCTGATTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTg                                                                                                                      <  2:569510/141‑1 (MQ=12)

                         tGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGt                                                                                                                     >  1:1038910/1‑141 (MQ=17)

                           aTAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGtaa                                                                                                                   <  1:597799/141‑1 (MQ=17)

                            tAAATGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTTTGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTTATCAATACTCCCtttttttt                                                                                                                  >  1:881725/1‑136 (MQ=11)

                            tAAACGCAACTGCTTAATGAGTATCTCGTTGCTTACGCAACCATATATCAATTGTATGTCCCTCTGCATCAATGGCACGATATAAATCGCTCCATTTGCCTTTTATTTTGATGGACGTCTCATCAATACGCCATTTGtaat                                                                                                                  <  2:129961/141‑1 (MQ=11)

                                  cATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATTTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCtt                                                                                                            <  2:693845/141‑1 (MQ=255)

                                  cATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCACTATGGTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAAAAAGCtt                                                                                                            >  2:907062/1‑141 (MQ=11)

                                      tGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCCTCACTACGCCATAGGGAATAAGCttttat                                                                                                        >  1:296580/1‑139 (MQ=11)

                                      tGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCtttttt                                                                                                        >  1:214691/1‑141 (MQ=35)

                                      tGCTGAATGATTACCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCtttttt                                                                                                        <  2:246665/141‑1 (MQ=25)

                                        cTGAATGATTATCTCGTTGCTTACGCAACCAATTATCTAATGTATTGCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTAt                                                                                                      <  1:455328/141‑1 (MQ=255)

                                        cTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTAt                                                                                                      <  1:476953/141‑1 (MQ=255)

                                         tGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCTATGGCACGATCTAAATAGCTCCATTTTCCTTTTATTTTGTTGTACGTCTGCTCAATGAGCCATTTGTATGAATCTTTTTTAtt                                                                                                     >  1:310443/1‑140 (MQ=11)

                                          gAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  1:395174/141‑1 (MQ=39)

                                          gAATGATTATCTCGTTCCTTACGCAACTAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTGTAGTTGGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGc                                                                                                    <  1:615145/141‑1 (MQ=14)

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

                                                 tATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACTTCTCATCAATACGCCATTTGTAGTAAGCTTTTTTATGCTTTttct                                                                                             >  1:94643/1‑141 (MQ=255)

                                                 tATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttct                                                                                             <  1:49096/141‑1 (MQ=255)

                                                 tATCTCGTCGCTTACGCATCCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttct                                                                                             <  1:1058962/141‑1 (MQ=255)

                                                  aTCTGGTTGCTTACGCAACGAAGTGTCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATAGACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTGttctt                                                                                            <  2:162078/141‑1 (MQ=14)

                                                  aTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTTATCAATACGCCATTTGTAATAAGCTTTTTTTTGCTTTttctt                                                                                            >  1:628028/1‑141 (MQ=255)

                                                  aTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATATTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttctt                                                                                            >  2:951353/1‑141 (MQ=255)

                                                   tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                                                                           <  1:795800/141‑1 (MQ=255)

                                                   tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                                                                           <  1:83228/141‑1 (MQ=255)

                                                   tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                                                                           <  1:272140/141‑1 (MQ=39)

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

                                                         tGCTTACGCAACCAAATATCTAATGTATATCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAtt                                                                                     <  1:898517/141‑1 (MQ=255)

                                                         tGCTTACGCAACCAAATATCTAAAGTATGTCCCTCTGCATCAATGGAACGATTTAAATAGACACATTGTACTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTTTTCCAAttt                                                                                     >  2:434942/1‑141 (MQ=14)

                                                              aCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACCTCAATTTTAATAAAATTTTTTTTGCTTTATATTCCAAAtttatga                                                                                >  1:539144/1‑137 (MQ=14)

                                                              aCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                                                                <  1:783659/141‑1 (MQ=40)

                                                               cGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                                                               >  1:424691/1‑140 (MQ=255)

                                                                gCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCACCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATata                                                                              <  1:558265/141‑3 (MQ=255)

                                                                ccAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATata                                                                              <  2:704112/140‑3 (MQ=255)

                                                                  aaCCAAATATCTAAGGTATGTCCCGCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaaa                                                                            <  2:455850/141‑1 (MQ=11)

                                                                    ccAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCGTTTGTAATAATATTTTTTATGCTTTTTCTTCCAAATTTGATATTAAAt                                                                          >  1:370139/1‑141 (MQ=11)

                                                                     cAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCACTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTTTTCCAAGGTTGCCAAAAAAtt                                                                         >  1:712715/1‑141 (MQ=14)

                                                                     cAAATATCTAATCTATGTCCCTCTGCATCAATTGCACGATATAAATAGCACCATTTTCCTTTTATTTTGATGTACGTTTCATCAATAGGCCATTTGTAATAAGCTTTTTTATTGTTTTTTTTCTACATGTTTTACAACAtt                                                                         >  2:35865/1‑141 (MQ=14)

                                                                         tatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATTGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTgggg                                                                     <  2:416705/141‑1 (MQ=255)

                                                                         gagcTAATGTATGTCCCTCTGGATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTgggg                                                                     <  2:327800/138‑1 (MQ=255)

                                                                           tCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTTTGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCa                                                                   >  2:281844/1‑141 (MQ=255)

                                                                           tCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATTTACGTCTCATCTATACGCCATTTGTATTAAGCTTTTTTATGCCTTTTAATCCAAAATGGATATAAAATGGTGGCa                                                                   >  1:81147/1‑141 (MQ=11)

                                                                           tCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATTTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGGTTTTTCTTTCAGATTGGGTATAAAATTGGGGCa                                                                   >  1:352958/1‑141 (MQ=11)

                                                                           tCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCa                                                                   >  2:20723/1‑141 (MQ=38)

                                                                             tAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCACCTTCCTTTTATTTTGATGTACGCCTCATCAATACCCCATTTGTAATAAGATTTTTAATGCTTTTACTTCCCACAATGATATAAAAtagagcaata                                                                 >  2:107257/1‑132 (MQ=11)

                                                                               atgtatgtTCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAATTTTGATACAAATTTGGGGCATAtt                                                               <  2:514305/141‑1 (MQ=255)

                                                                               atgtatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATAtt                                                               <  1:522753/141‑1 (MQ=255)

                                                                                 gtatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCt                                                             <  1:137149/141‑1 (MQ=255)

                                                                                   atgtCCCTCTGCATCAATGGGACGATATAAATAGCTCCATTTTCATTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAACACGCTTTTTTATGCTTTTTCTGCCAAATTTGATACAAAATTGGGGCATATTCttt                                                           >  2:580186/1‑140 (MQ=14)

                                                                                    cctacctctctaacAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGa                                                          <  2:789667/128‑1 (MQ=11)

                                                                                          tctGTATCAATGGCACGATAGAAATAGCTCCATTTTCCTTTTATTTTGATATACTTCTCATCAATACGCCCTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATGCGTGAACCCaa                                                    >  1:687789/1‑141 (MQ=11)

                                                                                          tctGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATTCGCCCTTTGTACTAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCaa                                                    >  1:85666/1‑141 (MQ=255)

                                                                                             gtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTTCAAATTTCATATAAAATTGGGGCTTATTATTGAAGCCAACgg                                                 >  1:115938/3‑141 (MQ=11)

                                                                                             gCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                 <  1:511530/141‑1 (MQ=21)

                                                                                             ccATCAATGGCACGAAATAAAAAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                 <  1:1048970/140‑1 (MQ=255)

                                                                                              cATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATCTTCTTGAACCCAACGGt                                                >  1:441751/1‑141 (MQ=255)

                                                                                                tCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAg                                              <  1:964317/141‑1 (MQ=21)

                                                                                                tCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAg                                              <  2:862093/141‑1 (MQ=21)

                                                                                                tCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTATTCCATACGCCATGTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAATCCAACGGTAg                                              <  1:770981/141‑1 (MQ=11)

                                                                                                acAATGGCACTATATAAACAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAAGTGGGGCATATTCTTGAACCCAACGGTAg                                              <  2:774013/140‑1 (MQ=255)

                                                                                                   aTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCATCGGTAGAcc                                           >  2:1042739/1‑141 (MQ=255)

                                                                                                    tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGatcg                                          <  1:752423/141‑4 (MQ=255)

                                                                                                     ggCAAGATATTAATAGCTCCATTTTCCTTTCATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTGATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAACCAACGGTAGACCGt                                         <  1:227432/141‑1 (MQ=11)

                                                                                                      gCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAGATTTGATATAAAATTGGGGCATATTCTTGAATCTAACGGTAGACCGtt                                        >  1:585018/1‑141 (MQ=255)

                                                                                                      gCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGtt                                        >  2:661745/1‑141 (MQ=40)

                                                                                                      gCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGtt                                        <  2:381220/141‑1 (MQ=255)

                                                                                                       cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAGTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTAGAACCCTACGGTAGACCGTTg                                       <  2:726516/141‑1 (MQ=255)

                                                                                                       cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCGATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTg                                       <  2:616546/141‑1 (MQ=21)

                                                                                                         cGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAGATTCTTGAACCCAACGGTAGATCGtgtaa                                     >  2:673703/1‑137 (MQ=11)

                                                                                                          gATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTGGAat                                    >  1:1021522/1‑141 (MQ=25)

                                                                                                           atatAAATTGCTCCTATTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatg                                   <  1:260287/141‑1 (MQ=255)

                                                                                                                 aTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGTTAGACCGTTGAATGATGAAc                             >  1:72705/1‑141 (MQ=25)

                                                                                                                   aGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGt                           >  1:92862/1‑141 (MQ=21)

                                                                                                                   aGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGt                           <  2:588058/141‑1 (MQ=21)

                                                                                                                   aGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGt                           >  1:338362/1‑141 (MQ=255)

                                                                                                                    gCTCCATTTTCCTTTTATTTTGATATGCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                          <  1:643101/141‑1 (MQ=255)

                                                                                                                    gCTCCACTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGACCGtt                          <  2:831601/141‑1 (MQ=255)

                                                                                                                      tCCATTTTCCTTTTATTTAGAAGGAGGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATACCAAATTGGGGCATATTCTTGAACCCAACGGTAGGCCGTTGGATGATGAACGTTTa                        >  2:489085/1‑141 (MQ=14)

                                                                                                                         aTTTTCCTTTTTTTTGTATATACTTCTCATCTATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTTAACCCAACGGTAGATCGTTGAATGATGAACGTTTacac                     <  2:63219/141‑1 (MQ=11)

                                                                                                                         aTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATCTAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTacac                     >  2:523637/1‑141 (MQ=255)

                                                                                                                          ttttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacc                    <  2:179148/141‑1 (MQ=21)

                                                                                                                          ttttCCTTTTATTTTGATGTACGTCTAATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacc                    <  1:1002229/141‑1 (MQ=21)

                                                                                                                           tttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacca                   <  2:144493/141‑1 (MQ=18)

                                                                                                                             tCCTTTTATTTTGACGTACTTCTCATCAATACGCCATTTGTAATAAGCTTTTTTACGCCTTTTCTTCCAAATTTCATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGTCCGTTGAATGACGAACGTTTACACCACg                 <  2:547353/141‑1 (MQ=11)

                                                                                                                              ccTTTGATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGACCCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGt                <  2:193882/141‑1 (MQ=11)

                                                                                                                                ttttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTc              >  2:886615/1‑141 (MQ=21)

                                                                                                                                   tATTTTGATATACGTCTCATCTATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTctct           <  1:693965/141‑4 (MQ=11)

                                                                                                                                    ttttcGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt          <  2:925336/140‑1 (MQ=11)

                                                                                                                                      tttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTaa        <  1:981727/141‑1 (MQ=18)

                                                                                                                                      tttGATGTACCTCTCATCAATACGCCACTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGCGCATATTCTTGAACCCAACCGTAGACCGTTGAATTATGAACGTTTACACCACGTTCCCTTaa        <  2:808891/141‑1 (MQ=14)

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` `< 27 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`