BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,378	C→T	89.0%	L58L (TTG→TTA)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,378	0	C	T	89.0%	1566.7 / 183.9	827	L58L (TTG→TTA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (45/46); new base T (381/347); total (434/393)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

Only 100 of 828 total aligned reads displayed.

CGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACA  >  CP000730/36238‑36516

          tatCTAATTTATGTCCCTCTGTATCAATGGTACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                                                                                                                                   <  1:5992/141‑1 (MQ=11)

          aatCTAATTTATGTCCCGTTGGATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAGTA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                                                                                                                                   <  2:278931/140‑1 (MQ=11)

          aatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                                                                                                                                   <  2:115375/140‑1 (MQ=255)

             tttatggatgtCCCTCTGCATCAATGGCACGATATAAATAGCCCCATGTTCCTTTTATTTTGATTTACGTCTCAGCTATA‑CGTCATTTGTAAGAGTTTTTTTTATGCTTTTTCGTCCAAATTTGGTTCAAAATTGGGGCat                                                                                                                                <  2:186905/138‑1 (MQ=14)

              tAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCata                                                                                                                               <  2:64087/141‑1 (MQ=255)

              tAATGTATGTCCCTCTGCATCAATGGCACGATATACATAGCTCCATTTTCCTTTTATTTTGATGTACTTCTCATCAATA‑CGCCATTAGTAATAAGCTTTTTTATGCTTTTTCTTTCAAATTTGATATAAAATTGGGGCatt                                                                                                                               >  1:117354/1‑140 (MQ=11)

              tAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCata                                                                                                                               <  1:382368/141‑1 (MQ=37)

              tAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCata                                                                                                                               <  2:353776/141‑1 (MQ=37)

                 tgtatgtCCCTCTGTAACAATCGCACGATATAAATATCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTc                                                                                                                            >  1:340794/1‑141 (MQ=11)

                   tatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCtt                                                                                                                          <  1:233163/141‑1 (MQ=37)

                    ttgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTg                                                                                                                         <  2:667541/140‑1 (MQ=255)

                     tgtgCCTCTGTATCAATGTCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGa                                                                                                                        <  1:460142/141‑1 (MQ=11)

                     tgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGa                                                                                                                        <  1:357739/141‑1 (MQ=255)

                     tgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGa                                                                                                                        >  1:613547/1‑141 (MQ=255)

                     tgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGa                                                                                                                        <  1:92546/141‑1 (MQ=255)

                     tgtCACTCTGTATCAATGGCACGATATAAATAGAACCATCTACCTTTTTTTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAATTTGGGGGATATTCTTGa                                                                                                                        >  1:319659/1‑141 (MQ=11)

                      gtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGaa                                                                                                                       <  1:206523/141‑1 (MQ=255)

                       tCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTTTTGAAc                                                                                                                      >  2:598156/1‑141 (MQ=255)

                       tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACCACCAAGT‑TGATAAGATTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAc                                                                                                                      >  1:484825/1‑141 (MQ=11)

                         ccTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAccc                                                                                                                    <  2:24331/141‑1 (MQ=18)

                          ctctGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCa                                                                                                                   >  2:427044/1‑141 (MQ=25)

                          ctctGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCa                                                                                                                   >  1:546391/1‑141 (MQ=21)

                                aTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTa                                                                                                             >  1:325317/1‑141 (MQ=21)

                                aTCAACGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTa                                                                                                             >  2:314044/1‑141 (MQ=255)

                                 tCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATTTTCTTGAACCCAACGGTAg                                                                                                            >  1:263615/1‑141 (MQ=255)

                                 tCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAg                                                                                                            <  2:204366/141‑1 (MQ=19)

                                 tCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCAAAATTTGATATAACATTGGGGGATATTCTTGaaacacaagatag                                                                                                            >  2:393798/1‑130 (MQ=11)

                                    aTGGTACGATATAAATAGCTCCATTTTCCTTTTATTTTTATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGatc                                                                                                         <  1:507493/141‑3 (MQ=11)

                                     tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTCTCTTCCAAATTTGTTATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCg                                                                                                        <  2:58920/141‑1 (MQ=255)

                                      ggCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGt                                                                                                       <  1:121937/141‑1 (MQ=255)

                                        cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTg                                                                                                     <  2:328657/141‑1 (MQ=255)

                                        cACGATATAAATAGCTCCATCTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTg                                                                                                     <  2:485776/141‑1 (MQ=255)

                                         aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTTTGTTTTTTTTTCCAAATTTGATATAAAATTGGGGTATATTCTTGAACCCAACGGTAGACCGTTGa                                                                                                    >  2:47462/1‑141 (MQ=11)

                                         aCGATATAAATAGCCCCATTTTCCTCTTATTTTGTTGTACGTATTTAAACTA‑CTCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGa                                                                                                    <  1:151884/141‑1 (MQ=11)

                                          cGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGaa                                                                                                   >  2:91864/1‑141 (MQ=255)

                                            atatAAATTGTCCCATTTTCCTTTTACTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatg                                                                                                 <  2:259350/141‑1 (MQ=11)

                                            atatAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatg                                                                                                 >  1:275155/1‑141 (MQ=255)

                                            atatAAATAGCTCCATTTTCCTTTTATTATGATATACGTCTCATCAATA‑CGCCATTTGTCATAAGTTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatg                                                                                                 <  1:661950/141‑1 (MQ=11)

                                               taAATAGCTCCATGTTCCTTTTATTTTGATGTACGTCTCATCAATG‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCATCGGTAGACCGTTGAatgatg                                                                                              <  1:468135/141‑1 (MQ=255)

                                                aaaTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTTTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTAATATAAAATTGAGGCATAATCTTGAACCCAACGTTAGCCCGTTGAAtgatga                                                                                             >  2:221585/1‑141 (MQ=11)

                                                aaaTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAAACAGTTAAAtgatga                                                                                             >  1:441836/1‑141 (MQ=255)

                                                 aaTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGaa                                                                                            >  1:1437/1‑141 (MQ=21)

                                                 aaTAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGaa                                                                                            >  1:580579/1‑141 (MQ=255)

                                                  aTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAAc                                                                                           <  1:525977/141‑1 (MQ=21)

                                                  aTAGATCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTTTTCTTTTTTTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAGCGGTAGATCGTTGAATGATGAAc                                                                                           >  2:403099/1‑141 (MQ=11)

                                                   tAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATTGTTGAATGATGAACg                                                                                          <  1:494316/141‑1 (MQ=255)

                                                     gCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAAc                                                                                           <  2:55613/138‑1 (MQ=21)

                                                      cTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTGAAAATTTGGTATAAATTTAAGGCAAATTATTGAACCCAACTGTAACAAGTTGAATGTTGaagattt                                                                                       >  1:637269/1‑136 (MQ=11)

                                                        ccATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTac                                                                                     <  2:339552/141‑1 (MQ=21)

                                                         cATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAAGGTTTaca                                                                                    >  1:40198/1‑141 (MQ=21)

                                                          tttttCCTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTaca‑c                                                                                  <  2:652041/140‑1 (MQ=255)

                                                          aTTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTaca‑c                                                                                  <  2:563621/141‑1 (MQ=21)

                                                           ttttCCTTTTATTTTGATGTACGTCTCATCCATA‑CGCCATTTGTAATAAGCTTTTTTCTGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGACCCCAACGGTAGACCGTTGAATGATGAACGTTTAca‑cc                                                                                 <  2:498823/141‑1 (MQ=255)

                                                            tttCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAca‑cca                                                                                >  1:448997/1‑141 (MQ=18)

                                                            ttcccTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAca‑cca                                                                                <  2:155372/138‑1 (MQ=255)

                                                              tCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACg                                                                              <  1:573138/141‑1 (MQ=18)

                                                               ccTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGt                                                                             >  2:37194/1‑141 (MQ=40)

                                                                cTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAAAGGTAGACCGTTGAATGATGAACGGTTACA‑CCACGtt                                                                            >  2:622020/1‑141 (MQ=255)

                                                                cTTTTATTTTGATGTACGTCACATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATAGGGGCATATTCTTGAACCCAACGGTAGACCGTTGGATGATGAACGTTTACA‑CCACGtt                                                                            >  1:419015/1‑141 (MQ=255)

                                                                cTTTTATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAAAGTTTACA‑CCACGtt                                                                            >  1:214557/1‑141 (MQ=255)

                                                                 ttttttttTGATTGACGTCTCATCAAGA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAGAATTGGTGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTc                                                                           <  2:128041/141‑1 (MQ=11)

                                                                  tttATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTct                                                                          >  1:648624/1‑140 (MQ=255)

                                                                    tATTTTGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTctct                                                                        <  1:349924/141‑4 (MQ=255)

                                                                     aTTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCtt                                                                       <  1:535492/141‑1 (MQ=21)

                                                                      ttttGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTa                                                                      >  2:515846/1‑141 (MQ=21)

                                                                         tGATATACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTCTCTTAata                                                                   <  2:105049/141‑1 (MQ=11)

                                                                          gatatACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAAGGTTTACA‑CCACGTTCTCTTAatat                                                                  >  1:283172/5‑141 (MQ=11)

                                                                          gatatACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTCTCTTAatat                                                                  >  1:390802/5‑141 (MQ=255)

                                                                          gATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAatat                                                                  >  1:218522/1‑141 (MQ=25)

                                                                          gATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGATAGACCGTTGAATGATGAACGTTTCca‑ccaccacctcctaattt                                                                  >  2:634706/1‑128 (MQ=11)

                                                                          gATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTAAACCCAACGGAAAACCGTTGAATGATGAACGTTTGCAGCCAGGTTcgtaaattt                                                                   >  2:457411/1‑133 (MQ=11)

                                                                           aTGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGGTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATAtt                                                                 >  1:177434/1‑141 (MQ=18)

                                                                                cGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGCAGCCGACTGGATTCCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGa                                                            >  1:19385/1‑141 (MQ=11)

                                                                                 ggctcAGCAATA‑CGCCATTTGTAATAAGCTGTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTGGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGat                                                           <  1:370730/139‑1 (MQ=255)

                                                                                 gTCTGATCAATA‑CGCCATTTGAAATAAGCTTTGTTATGCGTTTTATTCCAAATTTGATATAAAATTGCGGCAGCTTCTTGAACCCAACGGTAGATCGATGAGTGATGAACGTTCACA‑CCGTGTACTCATAATCTTTCAgta                                                           >  2:83555/1‑139 (MQ=2)

                                                                                 gTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATACAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTCTCTTAATTTTTCAGat                                                           >  1:166976/1‑141 (MQ=11)

                                                                                   ctcATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTCTCTTAATATTTCAGatat                                                         <  2:442069/141‑1 (MQ=255)

                                                                                   ctcATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGatat                                                         >  2:614662/1‑141 (MQ=21)

                                                                                     catcaATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTCTCTTAATATTTCAGatatat                                                       >  1:563384/1‑141 (MQ=255)

                                                                                     catcaATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTCTCTTAATATTTCAGatatat                                                       >  1:555505/1‑141 (MQ=255)

                                                                                       tcaATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGATATATCa                                                     <  2:383374/141‑1 (MQ=40)

                                                                                        caATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTCTCTTAATATTTCGGATATATCAc                                                    <  1:131111/141‑1 (MQ=255)

                                                                                         aaTA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTATTGCACCCACCGGTAGATCGTAGGAAGATGAACGTTTACA‑CCACGCTATCTTCATATTCCAGatatcataca                                                   >  2:19142/1‑135 (MQ=11)

                                                                                              cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCGTATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CAACGATCCCTTAATATTTCAAATATATCACAAtta                                               >  1:81814/1‑139 (MQ=11)

                                                                                               gCCATTTGTAATAAGCTTTTTTATGCGTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                              <  1:158398/141‑1 (MQ=18)

                                                                                               gCCATTTGTAATAAGCTTTTTTAGGCGTTTTCTTCCAAATTTGATATAAGATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTCTCTTAATATTTCAGATATATCACGATagc                                              <  1:299116/141‑3 (MQ=11)

                                                                                                ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGATATATCACGATAACt                                             >  1:587825/1‑141 (MQ=40)

                                                                                                ccATTTGTAATAAGCTTTTTTATGATTTTTCTTCCAAAGTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGAAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCGGATATATCACGATAACt                                             >  1:426839/1‑141 (MQ=11)

                                                                                                    ttaggataaGCTTGTTTCTGCTCTTGCCTCCACATTGGACATAAAATTTGCGCATATTCTTGAACCCAACGGTAGATCGTTGAATCATCCACGTGTACA‑CCACGCTCTCTTAACACTTCAGATATATCACCATAGCTCAAt                                         <  2:363399/136‑1 (MQ=2)

                                                                                                     tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAtc                                        >  2:28471/1‑140 (MQ=255)

                                                                                                     tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                        <  2:535239/141‑1 (MQ=255)

                                                                                                     tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                        <  2:580205/141‑1 (MQ=18)

                                                                                                       taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑CCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAtcca                                      >  2:297126/1‑138 (MQ=11)

                                                                                                       taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGAGATATCACGATAACTCAATgca                                      <  2:170774/141‑3 (MQ=12)

                                                                                                       taataaGCTTTTTTATGCTGTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACA‑GCACGGTGCCTTAATATTTCAGATATATCACGATAGCTCAATgca                                      >  2:138969/1‑139 (MQ=11)

                                                                                                         ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGATATATCACGATAACTCATTGtata                                    >  2:54104/1‑141 (MQ=255)

                                                                                                                tttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACAAAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGa                             >  1:625013/1‑141 (MQ=25)

                                                                                                                   ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                          <  1:100280/141‑1 (MQ=35)

                                                                                                                             tcttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACgg                >  2:547771/1‑141 (MQ=35)

                                                                                                                                       ttGATACAAAATTNGGGCATATTCTTGAACNCAGCGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTCATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt      <  1:244094/141‑1 (MQ=19)

CGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATA‑CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACA‑CCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACA  >  CP000730/36238‑36516

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