BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,315	G→A	56.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	56.6%	107.4 / inf	639	Y79Y (TAC→TAT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (126/145); new base A (175/181); total (309/330)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-02

Only 100 of 644 total aligned reads displayed.

TTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTC  >  CP000730/36175‑36455

 tACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGt                                                                                                                                             <  2:518693/141‑1 (MQ=32)

   ccAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTGGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATAc                                                                                                                                           >  1:22287/1‑141 (MQ=12)

   ccAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGTTCCATTTTCCTTTTATTTTGATATAc                                                                                                                                           >  2:634015/1‑141 (MQ=12)

   ccAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATAc                                                                                                                                           <  2:105359/141‑1 (MQ=17)

   ccAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTCTTATTTTGATATAc                                                                                                                                           <  2:456273/141‑1 (MQ=12)

   ccAAATGGTTTAATGAGACGTTTGATAAACGCATATGCTGAACGATTATCTTGGTGCTGACGCTACCAAATATCTACTGTATGCCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTTCTTTTATTTTGATATAc                                                                                                                                           <  2:346271/141‑1 (MQ=11)

    cAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAGATACg                                                                                                                                          >  1:40651/1‑141 (MQ=12)

     aaaTTTTTTACGGAGACGTTTTATAAACGTATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATCTCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGt                                                                                                                                         <  2:214997/141‑1 (MQ=11)

        ttGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGGACCTCTGCATCGATGGCACGATATAAATAGGTCCATTTTCCTTTTATTTTGATGtgaggttc                                                                                                                                      >  2:553854/1‑134 (MQ=11)

         tGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTCCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtca                                                                                                                                     <  1:523885/141‑1 (MQ=17)

           tttAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTATGATATACGTCtcatc                                                                                                                                   <  1:332577/141‑1 (MQ=12)

           tttAATGAGACGGTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCCAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcatc                                                                                                                                   <  1:508653/141‑1 (MQ=17)

              aaTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAAt                                                                                                                                <  1:537363/141‑1 (MQ=32)

               aTGAGACGTTTTATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATa                                                                                                                               <  1:115743/141‑1 (MQ=255)

                tGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATAc                                                                                                                              <  1:30506/141‑1 (MQ=32)

                tGAGACGTTTGATAAACGCATAGTCTGAATGCTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATAc                                                                                                                              <  2:440868/141‑1 (MQ=12)

                  ggactgttGATTAACGAATATGCTGAATTATTAGCCCGTTGCTTACGCAAACAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGc                                                                                                                            <  2:153765/135‑1 (MQ=11)

                    aCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCa                                                                                                                          >  1:269275/1‑141 (MQ=17)

                     cGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCAt                                                                                                                         <  2:541047/141‑1 (MQ=12)

                     cGTTTGAAAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGCTATACGTCTCTTCAATACGCCAt                                                                                                                         >  1:416648/1‑141 (MQ=11)

                       tttGATAAACTCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCAttt                                                                                                                       <  1:320609/141‑1 (MQ=12)

                       tttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTGATTAATACGCCAttt                                                                                                                       >  1:68161/1‑141 (MQ=255)

                        ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTAATTCTTTTACACGTCTCATCAATACGCCATTTg                                                                                                                      <  1:661603/141‑1 (MQ=11)

                        ttGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATAACTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGGCTCATCATTAcacatttta                                                                                                                      >  1:231203/1‑133 (MQ=11)

                         tGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGt                                                                                                                     <  1:377793/141‑1 (MQ=17)

                          gATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATTCGTCTTATCAATACGCCATTTGta                                                                                                                    <  1:633200/141‑1 (MQ=255)

                          gATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATCCGCCCTTTGta                                                                                                                    >  2:118441/1‑141 (MQ=17)

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

                            tAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGtaat                                                                                                                  <  1:584037/141‑1 (MQ=12)

                                cGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCCCATCAATACGCCATTTGTAATAAGc                                                                                                              >  1:567901/1‑141 (MQ=25)

                                  cATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCtt                                                                                                            <  1:402369/141‑1 (MQ=255)

                                      ttcTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATTGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCtttttt                                                                                                        <  2:660182/139‑1 (MQ=25)

                                       gCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAGGCttttttt                                                                                                       >  2:371503/1‑140 (MQ=255)

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

                                         tGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTTCTTTTATTTTGATATATGTTTCATCAATACGCGCGTTGTAATAAGATTTTTTAtt                                                                                                     >  1:13714/1‑140 (MQ=11)

                                         tGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATg                                                                                                     <  2:420768/141‑1 (MQ=38)

                                           aaTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGGCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGACGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCt                                                                                                   <  2:356502/141‑1 (MQ=255)

                                              gATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtttt                                                                                                >  2:255970/1‑141 (MQ=39)

                                               aTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCttttt                                                                                               >  2:528629/1‑141 (MQ=255)

                                                ttATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTTATATACGTCTCATCACTACGCCCTTTGTAATAAGCTTTTTTATGttttttt                                                                                              >  1:258505/1‑140 (MQ=14)

                                                 tttcgcgctGCTTACGCAACCAAATAGCTAATGGATGGTGCTCTGTATCAATGGCACGATATAAATAGCTGCATTTTCCTTGTATTTTGATATACGTCTCATCAATACGCCATGTGTAATAAGCTTTTTTATGCTTTttct                                                                                             <  2:325658/133‑1 (MQ=14)

                                                   tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                                                                           >  2:397335/1‑141 (MQ=40)

                                                   tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                                                                           <  2:237656/141‑1 (MQ=40)

                                                   tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                                                                           <  1:596825/141‑1 (MQ=40)

                                                   tctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                                                                           <  1:199154/141‑1 (MQ=40)

                                                    ctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATAGGTGTCAATAATATGCCATATGTAATAAGCTTTTTTATGCTTTTTCTTcc                                                                                          >  1:622426/1‑141 (MQ=14)

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

                                                         ggCTTACGCACCCAAATATCTAATGTATGTCCCTCTGTATCATTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTGGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAtt                                                                                     <  2:311361/140‑1 (MQ=14)

                                                          gCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTGTACGTCTCATCAATACGCCATTTGTAATAAGTTGTTTTATGCTTTTTCTTCCAAAttt                                                                                    <  2:643424/141‑1 (MQ=255)

                                                           cTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAACTTTg                                                                                   >  2:191259/1‑141 (MQ=255)

                                                              aCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTGTATGCTTTTTCTTCCAAATTTGATa                                                                                <  2:473346/141‑1 (MQ=255)

                                                               cGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                                                               <  1:100285/141‑2 (MQ=255)

                                                               cGCAACCAAATATCTAATGTATGTCCCCCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                                                               <  2:14153/141‑2 (MQ=255)

                                                               agctacaaaaTATCAAATGTATGTCCCTCTTTATCAATGTCACGATATAAATAGCTTCATTTTCCTTTTATTTTGTTATACGTCTCATCAATACTCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                                                               <  2:404976/134‑2 (MQ=14)

                                                                 cAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATataa                                                                             <  1:457663/141‑4 (MQ=255)

                                                                   aCCTAATTTCTAATGTATTTCCCTCTGCATCAATGGCACTATATATATTGTTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaaaa                                                                           <  1:288651/141‑1 (MQ=11)

                                                                   aCCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaaaa                                                                           <  1:179078/141‑1 (MQ=255)

                                                                     cAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTTTTTTGTTATACTTCTCATCAATACGCCAGTTGTATTATGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAttt                                                                         >  1:36917/1‑141 (MQ=11)

                                                                     cAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATAATTCATCATTAAGAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAtt                                                                         >  1:549401/1‑141 (MQ=11)

                                                                     cAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAtt                                                                         >  1:559627/1‑141 (MQ=37)

                                                                        atatCTAATGTATGTTCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTggg                                                                      >  1:495732/1‑141 (MQ=255)

                                                                        atatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTCTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCACATTTGCTATAAAATTggg                                                                      >  2:72600/1‑141 (MQ=11)

                                                                        agatCTAATGGATGTCCTTGTGTATCAATGGCACGATATAAATAGCTCCATTGTCCTTTTCTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTTCAAATTTGATATAAAATTggg                                                                      <  1:344665/139‑1 (MQ=11)

                                                                         tatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCATTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                                                                     >  1:432402/1‑141 (MQ=255)

                                                                          atCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGc                                                                    <  1:169410/141‑1 (MQ=37)

                                                                           tCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTACAAATTTGATATAAAATTGGGGCa                                                                   >  1:189113/1‑141 (MQ=255)

                                                                            ttAACGTATTTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCTATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCat                                                                  <  1:275589/140‑1 (MQ=11)

                                                                             tAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCata                                                                 <  2:378064/141‑1 (MQ=255)

                                                                               atgtatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCTATTAGTAATAAGCTTTTTTTTGCTATTAATTCAACATTGGATACAAAATTGGGGCATAtt                                                               >  1:88456/1‑141 (MQ=14)

                                                                                 gtatgtCCCTGTGTACGATTGGCACGCTATACATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATACGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCt                                                             <  2:169317/141‑1 (MQ=11)

                                                                                  ggtgtCCGTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCtt                                                            <  2:625525/139‑1 (MQ=11)

                                                                                            tgtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAACTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACg                                                  >  1:384291/4‑141 (MQ=255)

                                                                                            tGCATCAATGTCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTTTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACg                                                  <  2:12606/141‑1 (MQ=255)

                                                                                             gtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTGTTCCAAATTTGATATAAAAGTGGGGCATATATTTGAACCCAACgg                                                 >  1:237974/3‑141 (MQ=11)

                                                                                             gtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTTAACCCAACgg                                                 >  1:103931/3‑141 (MQ=255)

                                                                                             gtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTACTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                 <  1:14471/139‑1 (MQ=255)

                                                                                             gCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGAAGATGCTTGAACCCAACgg                                                 >  1:604582/1‑141 (MQ=255)

                                                                                             gCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACgg                                                 <  1:577629/141‑1 (MQ=39)

                                                                                                   aTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGatc                                           >  2:604924/1‑139 (MQ=255)

                                                                                                    tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCg                                          <  1:295886/141‑1 (MQ=40)

                                                                                                    tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAAATTCTTGAACCCAACGGTAGATCg                                          >  2:230715/1‑141 (MQ=255)

                                                                                                     ggCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGaa                                                         <  2:560486/125‑1 (MQ=25)

                                                                                                         cGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAAGGGTAGATGGTTGaa                                     >  1:35738/1‑141 (MQ=11)

                                                                                                         cGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCg                                          <  2:467876/136‑1 (MQ=255)

                                                                                                         cGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGAATATTCTTGAACCCACCGGTAGATCg                                          >  1:467876/1‑136 (MQ=11)

                                                                                                          gATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAat                                    <  1:531579/141‑1 (MQ=255)

                                                                                                            tataAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatga                                  <  1:642242/141‑1 (MQ=21)

                                                                                                              taAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatgatg                                <  2:298501/141‑1 (MQ=255)

                                                                                                                   aGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGt                           <  2:269275/141‑1 (MQ=255)

                                                                                                                        ttttttCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGTTTTTTCTTCCAAAGTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTaca                      <  2:508381/139‑1 (MQ=11)

                                                                                                                        cATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTaca                      >  1:438965/1‑141 (MQ=255)

                                                                                                                          ttttCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcacc                    <  1:208404/141‑1 (MQ=255)

                                                                                                                          ttttCCTGTTATTTTGATATACGTCTCATCAATACGCCCTTTGTAATAAGCTTTTTTATGCTTTTGCTTCCACATTTGATATACCCTCGGAGCTGATTCTTGACCCCAACGCTAGATCGTTGAATGATGAACGTTTAcacc                    <  1:95117/141‑1 (MQ=11)

                                                                                                                              tcTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGt                <  1:245119/140‑1 (MQ=255)

                                                                                                                              ccTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGGATATTTTTGAAGCCATAGGTTGACCGTTGAATTATGAAGGTTTGCACCAAGt                >  2:436305/1‑141 (MQ=11)

                                                                                                                                ctttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACCCCATGTTc              >  1:31324/2‑141 (MQ=255)

                                                                                                                                 ttttttttGATATACGTCTCATCAATATGCCATTTTTAATAATCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGTTCATATTCTTGAACCCCACAATATATCGGTTAATGATGAACGTTTACACCACTTTct             <  1:132020/141‑2 (MQ=11)

                                                                                                                                    aTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTTGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt          <  2:576944/141‑1 (MQ=11)

                                                                                                                                      tttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCTATTTTTTATATAAAATTGGGGCATATTCTTTAACCCAACGGTAGACGGTTGAATGATGAACGTTTCCACCACGtcgattttt        >  2:283483/1‑133 (MQ=11)

                                                                                                                                      tttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTaa        <  1:490478/141‑1 (MQ=255)

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