BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	61,741	C→T	7.4%	E265K (GAA→AAA)	USA300HOU_0055 ←	integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	61,741	0	C	T	7.4%	278.7 / 2.8	95	E265K (GAA→AAA)	USA300HOU_0055	integrase
Reads supporting (aligned to +/- strand): ref base C (38/48); new base T (7/2); total (45/50)
Fisher's exact test for biased strand distribution p-value = 7.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTATTTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATCGGCAGGTG  >  CP000730/61602‑61881

atatTATTTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCg                                                                                                                                              <  2:400124/141‑1 (MQ=255)

 tattatTTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGa                                                                                                                                             <  2:961868/141‑1 (MQ=255)

   ttattTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGatt                                                                                                                                           >  1:489073/1‑141 (MQ=255)

   ttattTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGatt                                                                                                                                           <  1:720133/141‑1 (MQ=255)

     attTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGattat                                                                                                                                         <  1:629804/141‑1 (MQ=255)

          cacaAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGAttattatta                                                                                                                                    <  2:680350/141‑1 (MQ=255)

          cacaAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGAttattatta                                                                                                                                    <  2:680414/141‑1 (MQ=255)

             aaaaaaaGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAg                                                                                                                                 >  1:104648/1‑141 (MQ=255)

              aaaaaaGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGt                                                                                                                                >  1:893194/1‑141 (MQ=255)

                 aaaGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGtttt                                                                                                                             >  2:633380/1‑141 (MQ=255)

                  aaGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGttttt                                                                                                                            >  2:792709/1‑141 (MQ=255)

                      aTGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATg                                                                                                                        <  2:982366/141‑1 (MQ=255)

                       tGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGt                                                                                                                       <  1:1016023/141‑1 (MQ=255)

                          tttttGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑tt                                                                                                                   >  1:207397/1‑141 (MQ=255)

                              tGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTTTTTTTTTTTTTTATTATTATAGTTTTTAATGTA‑TTTTCg                                                                                                               >  2:216982/1‑141 (MQ=18)

                               gTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTG‑TTTTCGa                                                                                                              <  1:858282/141‑1 (MQ=255)

                               gTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGa                                                                                                              >  1:29582/1‑141 (MQ=255)

                                ttAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTTGAc                                                                                                             >  2:763493/1‑141 (MQ=255)

                                ttAGAGGCACTGAACCATTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGAc                                                                                                             <  2:1076093/141‑1 (MQ=255)

                                  agagGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACaa                                                                                                           <  1:543119/141‑1 (MQ=255)

                                   ttttctcTGAACCCTTAAAATGGGAACTTATTAGAAACGCGACGTTCTGGGCTGCAGATTCCCTGTCTTTTACAGTTGTTGTGTTGGCTAGTTTATGTTGAATTCGAATATTATTATAGTTTTTAATGTA‑TTTTCGACAat                                                                                                          >  2:641359/7‑141 (MQ=255)

                                   gagGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAat                                                                                                          >  1:1057717/1‑141 (MQ=255)

                                   gagGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAat                                                                                                          >  2:940766/1‑141 (MQ=255)

                                   gagGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTGGTTGAATTCGGTTATTATTATAGTTTTTAATGTA‑TTTTCGACAat                                                                                                          >  1:1068156/1‑141 (MQ=255)

                                      gCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATc                                                                                                       <  2:1068456/141‑1 (MQ=255)

                                       cACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCt                                                                                                      <  2:545863/141‑1 (MQ=255)

                                          tGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTAtt                                                                                                   <  2:68137/141‑1 (MQ=255)

                                          tGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTAtt                                                                                                   <  1:745385/141‑1 (MQ=255)

                                           gAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTa                                                                                                  >  1:3572/1‑141 (MQ=255)

                                               ccTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAat                                                                                              <  1:517695/141‑1 (MQ=255)

                                                 ttAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATTGTTTTTAATGTA‑TTTTCGACAATATCTATTACAatat                                                                                            >  2:829717/1‑141 (MQ=255)

                                                    aaaTGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGAt                                                                                         <  1:173890/141‑1 (MQ=255)

                                                     aaTGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGAtt                                                                                        >  1:225249/1‑141 (MQ=255)

                                                     aaTGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGAtt                                                                                        <  2:415788/141‑1 (MQ=255)

                                                         ggAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTagag                                                                                    <  2:423837/141‑1 (MQ=255)

                                                          gAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGt                                                                                   <  2:143935/141‑1 (MQ=255)

                                                            aCTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGtta                                                                                 >  1:396346/1‑141 (MQ=255)

                                                             cTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGttat                                                                                >  1:202022/1‑141 (MQ=255)

                                                             cTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGttat                                                                                >  1:284892/1‑141 (MQ=255)

                                                             cTTAATAAAAACACCATGTTCTAGGCTATTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTTGATTATTATTATAGTTTTTAATGTACTTTTCGACAATATCTATTACAATATGATTagagcta                                                                                 <  2:1078005/141‑4 (MQ=14)

                                                              ttAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGttatt                                                                               <  2:681598/141‑1 (MQ=255)

                                                                 ataaAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAg                                                                            >  2:85395/1‑141 (MQ=255)

                                                                 ataaAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAg                                                                            <  2:424372/141‑1 (MQ=255)

                                                                 ataaAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTGCGACAATATCTATTACAATATGATTAGAGTTATTaat                                                                            <  2:752948/141‑2 (MQ=255)

                                                                  taaAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGc                                                                           <  2:742483/141‑1 (MQ=255)

                                                                    aaaaCACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTg                                                                         <  2:151718/141‑1 (MQ=255)

                                                                     aaaCACCATGTTCTAGGCTATTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTTGATTATTATTATAGTTTTTAATGTACTTTTCGACAATATCTATTACAATATGATTAGAGCTATTAAGCTg                                                                         <  1:113414/141‑1 (MQ=17)

                                                                        cacTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGatta                                                                     >  1:120029/1‑141 (MQ=255)

                                                                           tATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTg                                                                  >  2:71046/1‑141 (MQ=255)

                                                                             tGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGAt                                                                >  1:578492/1‑141 (MQ=255)

                                                                              gTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTTTTTTTTTTTATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTGCAGTTATTAAGCTGATTATTGAtt                                                               >  2:1083305/1‑140 (MQ=21)

                                                                               ttCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGt                                                              <  2:470865/141‑1 (MQ=255)

                                                                               ttCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGt                                                              <  2:508700/141‑1 (MQ=255)

                                                                                 cTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTATGTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTaa                                                            <  1:10387/141‑1 (MQ=255)

                                                                                 cTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTaa                                                            <  2:525580/141‑1 (MQ=255)

                                                                                 cTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTaa                                                            <  1:792709/141‑1 (MQ=255)

                                                                                   aGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTaaaa                                                          <  1:726139/141‑1 (MQ=255)

                                                                                        gctAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTTAATTTGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGttt                                                     >  2:223239/1‑141 (MQ=255)

                                                                                           aaTTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAg                                                  >  1:137642/1‑141 (MQ=255)

                                                                                            aTTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATACCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGa                                                 >  1:977307/1‑141 (MQ=255)

                                                                                             ttCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGAc                                                <  1:576449/141‑1 (MQ=255)

                                                                                                tctGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACttt                                             <  2:101515/141‑1 (MQ=255)

                                                                                                tctGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACttt                                             <  1:85395/141‑1 (MQ=255)

                                                                                                   gTATTTTACAGGGGATAAGTAGCCTAGTTTTTTTTTTTTTTTTTTATTATTATAGGTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAATCTGATTATTGATGTAAAAAGTTTCAGACTTTAGc                                          >  2:597962/1‑141 (MQ=18)

                                                                                                    tATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCg                                         <  1:1057550/141‑1 (MQ=255)

                                                                                                     aTTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGa                                        <  2:706492/141‑1 (MQ=255)

                                                                                                     aTTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGa                                        >  2:795883/1‑141 (MQ=255)

                                                                                                     aTTTTACAGGGGATAAGTAGCCTAGTTTTTGGTGAATTCGATTATTATTATAGTTTTTAATGTG‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGa                                        <  2:934802/141‑1 (MQ=255)

                                                                                                      ttttACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAg                                       >  1:637193/1‑141 (MQ=255)

                                                                                                      ttttACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAg                                       >  1:569589/1‑141 (MQ=255)

                                                                                                      ttttACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAg                                       >  1:703585/1‑141 (MQ=255)

                                                                                                       tttACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAgg                                      <  1:72432/141‑1 (MQ=255)

                                                                                                        ttACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGa                                     >  1:241555/1‑141 (MQ=255)

                                                                                                         tACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGaa                                    >  1:26165/1‑141 (MQ=255)

                                                                                                          aCAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAAt                                   <  2:354166/141‑1 (MQ=255)

                                                                                                           cAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATg                                  >  1:958231/1‑141 (MQ=255)

                                                                                                           cAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATg                                  >  1:843486/1‑141 (MQ=255)

                                                                                                                gATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAAc                             >  1:980441/1‑141 (MQ=255)

                                                                                                                  tAAGTAGCCTAGTTTTTTTTTTTTTTTATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACtt                           >  2:138178/1‑141 (MQ=17)

                                                                                                                  tAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACtt                           <  2:329899/141‑1 (MQ=255)

                                                                                                                    aGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTc                         <  2:489073/141‑1 (MQ=255)

                                                                                                                       aGCCTAGTTTTTTTTGTTTTTTATAATTATTAAAGTTTTTAATGTA‑TTTTTTACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTAt                      >  2:930369/1‑141 (MQ=17)

                                                                                                                         ccTAGTTTTTGTTGAATTTGATTATTATTATAGTTTTTAATGTACTTTTCGACAATATCTATTACAATATGATTAGAGCTATTAAGCTGATTATTGATGTAAAGAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATc                     >  1:734034/1‑141 (MQ=17)

                                                                                                                         ccTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCg                    <  1:164180/141‑1 (MQ=255)

                                                                                                                             gTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGc                <  1:890927/141‑1 (MQ=255)

                                                                                                                               ttttGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGt              <  1:169643/141‑1 (MQ=255)

                                                                                                                                tttGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGtt             >  1:297220/1‑141 (MQ=255)

                                                                                                                                tttGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGtt             >  1:135794/1‑141 (MQ=255)

                                                                                                                                  tgttgAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTAt           <  1:198431/141‑1 (MQ=255)

                                                                                                                                  tgttgAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTAt           <  2:225249/141‑1 (MQ=255)

                                                                                                                                   gttgATTTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATc          >  2:5252/1‑141 (MQ=255)

                                                                                                                                   gttgAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATc          <  1:835059/141‑1 (MQ=255)

                                                                                                                                   gttgAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATc          <  2:966885/141‑1 (MQ=255)

                                                                                                                                      gAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGACTGGAAACTTTCTATCGGGGCGTTATCGGc       >  1:326759/1‑141 (MQ=255)

                                                                                                                                        aTTCGATTATTATTATAGTTTTCAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATCGGCAg     >  1:238292/1‑141 (MQ=255)

                                                                                                                                           cGATTATTATTATTGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATCGGCAGgtg  <  1:829717/141‑1 (MQ=255)

                                                                                                                                           cGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATCGGCAGgtg  >  2:240385/1‑141 (MQ=255)

ATATTATTTTCACAAAAAAAGAATGCTTTTTGTTAGAGGCACTGAACCCTTAAAATGGGAACTTAATAAAAACACTATGTTCTAGGCTGCTAATTCTCTGTATTTTACAGGGGATAAGTAGCCTAGTTTTTGTTGAATTCGATTATTATTATAGTTTTTAATGTA‑TTTTCGACAATATCTATTACAATATGATTAGAGTTATTAAGCTGATTATTGATGTAAAAAGTTTCAGACTTTAGCGAGGAATGGAAACTTTCTATCGGGGCGTTATCGGCAGGTG  >  CP000730/61602‑61881

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