BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,124	G→A	59.8%	A143V (GCT→GTT)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,124	0	G	A	59.8%	270.9 / inf	598	A143V (GCT→GTT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (112/127); new base A (184/170); total (299/299)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-02

Only 100 of 598 total aligned reads displayed.

AGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTC  >  CP000730/35985‑36262

   aTTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCAt                                                                                                                                        <  1:257377/141‑1 (MQ=255)

    tgcgttGCTGTATTGATACTTTGATACCTAGTGTTTCTTACTTTAATATAAGAG‑GTG‑TCTAGGTCAATGAGGTTATGCATAGATATCGTTGTACAATAGAAGTACGGTTTCAGTTT‑AAAGGAGTTAATTACTTTAACCAtt                                                                                                                                       >  2:220036/5‑141 (MQ=255)

     tttttGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTg                                                                                                                                      <  2:716214/139‑1 (MQ=255)

     tCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTg                                                                                                                                      <  1:225062/141‑1 (MQ=255)

       tttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAGCCATTGCt                                                                                                                                    <  2:90232/141‑1 (MQ=40)

        ttGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGAAAGTAAAAGTTAGGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCtg                                                                                                                                   >  1:398715/1‑140 (MQ=14)

          gccttATTGATACTTGGATATCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCtgcc                                                                                                                                 <  2:134145/137‑4 (MQ=14)

           ccgTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAGCCATTGCTACCt                                                                                                                                <  2:592981/139‑1 (MQ=255)

           cTGTATTTATACTTTGATACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCt                                                                                                                                <  2:711657/141‑1 (MQ=255)

           cTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATTTGA‑GG‑TTTTTCTTGCTTAATGAGGTTATTCAGATATTTCGTTGTATAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCt                                                                                                                                <  1:210419/141‑1 (MQ=14)

           cTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCt                                                                                                                                <  2:744875/141‑1 (MQ=255)

            tGTATTGATACTTTGATACCTTGTCTTTCTTACTTGAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAGCCATTGCTACCtt                                                                                                                               <  1:243399/141‑1 (MQ=255)

             gTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATAT‑A‑CGAGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAGCCATTGCTACCTTc                                                                                                                              >  1:338295/1‑141 (MQ=255)

               aTTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGt                                                                                                                            <  1:657855/141‑1 (MQ=255)

                ttGATACTTTGATACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGtt                                                                                                                           >  1:648697/1‑141 (MQ=255)

                 tGATACTTTGATACCTTGTCTTTCTTACTTTAATATGC‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTACCCATTTCTACCTTCGTTg                                                                                                                          >  1:367117/1‑141 (MQ=255)

                 tGATACTTTGATACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTg                                                                                                                          <  2:423570/141‑1 (MQ=255)

                     aCTTTGATACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAgg                                                                                                                      <  1:199347/141‑1 (MQ=255)

                        ttGATATCTTGTCTTTCTTACTTTAATATTA‑CT‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGc                                                                                                                   <  2:188229/141‑1 (MQ=11)

                         tgatctcctGTCGTTCTTACTCAAATAGGA‑CG‑GTGAGCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGcc                                                                                                                  <  1:15175/133‑1 (MQ=11)

                         tGATACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGcc                                                                                                                  <  1:85735/141‑1 (MQ=255)

                           aTACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAATGCTTTAATTAGTTTAACCATTGCTACCTTCGTTGAAGGTGCCTg                                                                                                                >  1:520139/1‑141 (MQ=255)

                            tatcTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGa                                                                                                               >  2:581416/4‑141 (MQ=255)

                            tACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGGTCTTGCTCAATGAGGTTATTCAGATATTTCGGTGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGa                                                                                                               >  2:562842/1‑141 (MQ=255)

                            tACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGa                                                                                                               >  1:479400/1‑141 (MQ=255)

                              tcTTGTCTTGCTTACTTTAATATGA‑CG‑GTGACCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATc                                                                                                             <  2:574678/140‑1 (MQ=255)

                                ttGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCtt                                                                                                           >  1:425639/1‑140 (MQ=255)

                                ttGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTg                                                                                                           >  1:353350/1‑141 (MQ=255)

                                ttGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTg                                                                                                           >  2:761179/1‑141 (MQ=255)

                                 tGTCTTTCTTTCTTTAATATGA‑CG‑GTGATCTTGCTTAATGGGGTTATTCTGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGt                                                                                                          <  2:244059/141‑1 (MQ=11)

                                  ttctttcgtACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGTTACGGTTTTTGAAGGTGCCTGATCTGTa                                                                                                         <  2:353407/140‑1 (MQ=11)

                                  gTCTTTGTTACTTGAATATGT‑CG‑GTGATCTTGATGAATGAGGTTATTCAGATATATCGTTGTACACTGACAGTCAGGGTTAGGTTT‑ATAAGCTTTAATTCCTTTAGCCATTGCTACCTTCGTTGAAGGTCCCGGATCTGTa                                                                                                         >  2:261176/1‑141 (MQ=11)

                                  gTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTATAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTa                                                                                                         >  2:696541/1‑141 (MQ=38)

                                  gTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTa                                                                                                         <  2:140191/141‑1 (MQ=255)

                                   tctttcttACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTaa                                                                                                        >  2:112641/1‑141 (MQ=255)

                                   tctttcttACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTaa                                                                                                        >  1:272693/1‑141 (MQ=255)

                                   tctttcttACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTaa                                                                                                        <  1:618072/141‑1 (MQ=255)

                                    ctttcttACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAAt                                                                                                       >  2:458286/1‑141 (MQ=255)

                                     tttcttACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAAtt                                                                                                      <  1:41711/141‑1 (MQ=255)

                                      ttcttACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGCTATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTaacta                                                                                                     >  1:108540/1‑138 (MQ=255)

                                      ttcttACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTa                                                                                                     <  2:120949/141‑1 (MQ=255)

                                       tcttACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTaa                                                                                                    >  1:504859/1‑140 (MQ=255)

                                       tcttACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTAc                                                                                                    >  2:366097/1‑141 (MQ=255)

                                        cttACTTTAAAATGA‑CG‑GTGATCTAGCTCAATGAGGTTATTCAGATTCATGCTCAGACAATGGATGTAAATCGGAAGTTT‑AACAGCTTTAATTACTTTGACCATTGCTACCTTCGTTGATGGTGCCTGATCTGTAATTAcc                                                                                                   >  1:560894/1‑141 (MQ=9)

                                         ttACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCt                                                                                                  >  2:359355/1‑141 (MQ=255)

                                         ttACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCt                                                                                                  >  2:81557/1‑141 (MQ=39)

                                         ttACTGTAAGAGGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAGAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCt                                                                                                  <  2:634764/141‑1 (MQ=255)

                                              ttAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGa                                                                                             >  1:128135/1‑141 (MQ=255)

                                                 atgtGA‑CG‑GGGATTTTGCTCAATGAGGTTGTTTAATTATTTCGATGTACGATGGCAGTCAGGTTTAAGTTT‑AAAAGGTTTGATTTCTTTAACCATTGCTACCTTCGCTGAAGGTGGCTGCTCTGTAATTAGCTTTTGAggg                                                                                          <  2:392111/138‑2 (MQ=11)

                                                    tGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACATTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTa                                                                                       >  1:444354/1‑141 (MQ=255)

                                                     gA‑CG‑GTGATCTTGCTCACGGTGGTTTTTCCGCTATTTCGATGTACAATTACAGTCAGGTTTAAGTTT‑AAAAGTTTTTATTATTTTAGCCATGGCTGCCGTCGTTGAAGATACCTGATCTGTAACTACCTTTTGATGCTTAc                                                                                      <  2:6891/141‑1 (MQ=2)

                                                      a‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGCGGTTTAcc                                                                                     >  2:41386/1‑141 (MQ=255)

                                                      a‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                     <  1:544995/141‑1 (MQ=255)

                                                      a‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGGTTGcc                                                                                     >  1:32035/1‑141 (MQ=255)

                                                      a‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAcc                                                                                     <  1:117079/141‑1 (MQ=11)

                                                       gcG‑GTTATCTTTCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAATTCT‑AAGAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                     <  1:178728/140‑1 (MQ=11)

                                                       gcG‑GTGATCTTGCTGAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAcc                                                                                     <  2:471672/140‑1 (MQ=11)

                                                        cG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCa                                                                                    <  1:265298/141‑1 (MQ=18)

                                                           gTGATCTTGCTCAATGAGGTTATTCAGGTATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaaa                                                                                  <  1:711484/141‑1 (MQ=12)

                                                           gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGTTCTGTAACTACCTTTTGAGGCTTACCaaa                                                                                  >  2:312583/1‑141 (MQ=11)

                                                           gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaaa                                                                                  <  2:486399/141‑1 (MQ=255)

                                                            tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                 <  1:49286/141‑1 (MQ=12)

                                                             gATCTTGCTCAATGAGGTTATTCAGATATTTCTATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAtt                                                                                <  2:650766/141‑1 (MQ=12)

                                                              aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAttt                                                                               >  1:578424/1‑140 (MQ=17)

                                                               tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                              >  2:325424/1‑141 (MQ=12)

                                                               tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGGTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGGTTACCAAATTGt                                                                              >  2:673602/1‑141 (MQ=255)

                                                                tttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGtt                                                                             <  2:736368/140‑1 (MQ=12)

                                                                cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGtt                                                                             >  2:286664/1‑141 (MQ=12)

                                                                 ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGttt                                                                            >  2:338592/1‑141 (MQ=17)

                                                                   gCTCAATGAGGTTATTCAGATATTGCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaa                                                                          <  1:194525/141‑1 (MQ=11)

                                                                       aaTGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATgag                                                                      <  2:517250/141‑1 (MQ=17)

                                                                       aaTGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTATAATTACCTTTTGAGGTTTTCCAAATTGTTTAATgag                                                                      >  1:697316/1‑141 (MQ=255)

                                                                        aTGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATgaga                                                                     >  1:59786/1‑141 (MQ=17)

                                                                        aTGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGGTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTTGTTGAAGGTGCCTGATTTGTTATTACCTTTTGAGGTTTTCGAAATGTTTTAAtagga                                                                     >  2:497477/1‑137 (MQ=11)

                                                                          gAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACg                                                                   <  1:100419/141‑1 (MQ=255)

                                                                           aGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                  <  1:662342/141‑1 (MQ=12)

                                                                           aGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                  <  2:618529/141‑1 (MQ=17)

                                                                            ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTTTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                 >  1:290662/1‑141 (MQ=12)

                                                                            ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                 >  2:271564/1‑141 (MQ=12)

                                                                                 ttCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATa                                                            >  2:531785/1‑141 (MQ=11)

                                                                                  tCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaa                                                           <  2:145989/141‑1 (MQ=17)

                                                                                     gATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAAGCTTTGATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAAc                                                         <  2:170029/141‑1 (MQ=12)

                                                                                      atatTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGGCGTTTGATAAACGc                                                       >  2:397416/1‑141 (MQ=12)

                                                                                      atatTTCGATGTACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGc                                                       <  1:306351/141‑1 (MQ=17)

                                                                                      atatTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGTGACGTTTGATAAACGc                                                       <  1:381350/141‑1 (MQ=11)

                                                                                       tatTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATCAGACGTTTGATAAACGCa                                                      >  2:447290/1‑141 (MQ=11)

                                                                                        atTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCat                                                     <  1:143626/141‑1 (MQ=11)

                                                                                           tccaTGTACAATGGCAGTCAAGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGCTTGATAAACGCATATg                                                  <  2:550032/138‑1 (MQ=11)

                                                                                            cGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                 <  1:462789/141‑1 (MQ=11)

                                                                                            cGAGGTACAATGTCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCTTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                 <  1:535716/141‑1 (MQ=11)

                                                                                                 tACAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGACGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAAt                                            >  1:678825/1‑141 (MQ=12)

                                                                                                  aCAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATg                                           >  2:659949/1‑141 (MQ=12)

                                                                                                   cAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                          >  2:240227/1‑141 (MQ=9)

                                                                                                   cAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                          >  2:302010/1‑141 (MQ=9)

                                                                                                   cAATGGCAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                          >  1:734633/1‑141 (MQ=14)

                                                                                                   aaaTGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                          <  2:153036/140‑1 (MQ=255)

                                                                                                    aaTGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                         <  1:487766/141‑1 (MQ=255)

                                                                                                                                           ttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatgt   <  1:722230/141‑1 (MQ=9)

                                                                                                                                           ttAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatgt   <  2:206782/141‑1 (MQ=14)

                                                                                                                                            tAGCCATAGCTGCCTTCGTTGAAGGTGCCTGATCTGTAAATACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAATCAAACATCTAATGTATGTc  <  1:500566/141‑1 (MQ=11)

AGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGA‑CG‑GTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTT‑AAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTC  >  CP000730/35985‑36262

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