BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,378	C→T	85.3%	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	85.3%	1291.6 / 247.3	785	L58L (TTG→TTA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (52/64); new base T (317/344); total (377/408)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 788 total aligned reads displayed.

CGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATC  >  CP000730/36238‑36518

   aagcAAATATCTAATGTATGTCCCTCTGCATCAATGGCACCATCTAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaaa                                                                                                                                           <  2:171269/138‑1 (MQ=11)

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

    aCCAAATATCTAATCTTTGTCCCTCTGCATCAATGGCACGATATAACTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaaaa                                                                                                                                          <  2:811881/141‑1 (MQ=255)

      cAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTAGGCTTTTTCTTCCATATTTTATATAAAAtt                                                                                                                                        >  2:206814/1‑141 (MQ=11)

      cAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAACTTTTTTTTTTTTTTTCTTACAAATTTGATACAAAAtt                                                                                                                                        >  2:560022/1‑141 (MQ=14)

         atatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATTAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGTTATAAACTTggg                                                                                                                                     >  1:510864/1‑141 (MQ=11)

         atatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATACTTTTTCTTCCAAATTTGATATAAAATTggg                                                                                                                                     >  1:326719/1‑141 (MQ=255)

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

             atAATGGATGTCCCTCTGGAGCAATGGCACGATATAAATAGCTCCATTTTGCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCat                                                                                                                                 <  2:836660/140‑1 (MQ=11)

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

                 tgtatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATAttt                                                                                                                             >  2:66892/1‑140 (MQ=255)

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

                    ttgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCACTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTCTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTg                                                                                                                          <  2:380744/140‑1 (MQ=255)

                    aggttccTCTGGATCAACGGAACTATATAAATAGCTCTATTTTCCTTTTATTTTGATATACGTCTCATCAATATGCCATTTGTAATAAGCTTTGTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTg                                                                                                                          <  1:571038/136‑1 (MQ=11)

                     tgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAGTTGGGGCATATTCTTGa                                                                                                                         <  2:1058476/141‑1 (MQ=255)

                     tgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGa                                                                                                                         <  1:380463/141‑1 (MQ=255)

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

                      ggcccTCTGTATCATTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTGTTTCTTGCAAATTTGATATAAAATTGGGGCATATTCTTGaa                                                                                                                        <  2:518667/139‑1 (MQ=11)

                       tCCCTCTGTATCAATGGCAGGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTATTCCAAATTTGATATAAAAGTGGGGCATATTCTTAAAc                                                                                                                       >  1:688242/1‑141 (MQ=11)

                       tCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTCATATAAAATGGGGGCATATTCTTGAAc                                                                                                                       >  1:191367/1‑141 (MQ=11)

                       tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAGAATTGTGGCATATTCTTGAAc                                                                                                                       >  1:851553/1‑141 (MQ=255)

                            ttGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCTATACGTCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAAc                                                                                                                  <  1:164981/140‑1 (MQ=255)

                             tgtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATCCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATAGAAAGTTGGGGCATTGGCTTGAACCCAACg                                                                                                                 >  1:883420/4‑141 (MQ=11)

                              gtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                                                                                <  2:656645/139‑1 (MQ=255)

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

                               cATCAATGGCACGATATAAATAGCTCAATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGACATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGGTATAAAATTGGGGCATCTTCTTGAACCCCACGGt                                                                                                               >  2:19821/1‑141 (MQ=11)

                                aTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTa                                                                                                              >  1:490590/1‑141 (MQ=25)

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

                                 tCAATGGCAAGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTAGGGCATATTCTGGAACCCAAAGGtcg                                                                                                             >  1:777507/1‑139 (MQ=255)

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

                                   ggattcACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTTGGGCATATTATTGAACCCAACGGTAGat                                                                                                           >  2:535670/6‑140 (MQ=11)

                                   aaTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAAAc                                                                                                           >  1:339208/1‑141 (MQ=21)

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

                                      ggCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGt                                                                                                        >  1:452689/1‑141 (MQ=25)

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

                                       gCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGtt                                                                                                       >  2:189738/1‑141 (MQ=255)

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

                                        cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGCTATTAAATTGGGGCCTATTATTAGACGCACAGGTAGAcagggg                                                                                                      >  1:349552/1‑136 (MQ=11)

                                        cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATCCTTGAACCCAACGGTAGACCGTTg                                                                                                      >  1:422360/1‑141 (MQ=21)

                                         aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACACAACGGTAGACCGTTGa                                                                                                     >  1:3021/1‑141 (MQ=21)

                                         aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGa                                                                                                     <  1:226775/141‑1 (MQ=255)

                                         aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAAATCGTTGa                                                                                                     >  2:576011/1‑141 (MQ=255)

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

                                             tataAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatga                                                                                                 >  1:694348/1‑141 (MQ=21)

                                               taAATAGCTCCATTTTCCTTTTATTTTGATATGCGCTTTATTAGTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatgatg                                                                                               <  2:479960/141‑1 (MQ=11)

                                                 aaTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCCTCATCTTGAACCCAACGATAGACCGTTGAATGATGaa                                                                                             >  1:743428/1‑141 (MQ=11)

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

                                                  aTAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAAc                                                                                            <  2:922852/141‑1 (MQ=255)

                                                     gCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATAGAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                                                                                         >  2:582139/1‑141 (MQ=255)

                                                       tCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGAATATTTTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTa                                                                                       >  1:801680/1‑141 (MQ=255)

                                                        ccATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTac                                                                                      >  2:403497/1‑141 (MQ=21)

                                                         cATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTaca                                                                                     >  1:752787/1‑141 (MQ=21)

                                                         cATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTaca                                                                                     <  1:102807/141‑1 (MQ=255)

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

                                                          aTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTacac                                                                                    <  2:852277/141‑1 (MQ=21)

                                                           ttttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCCTATTCCTGAACCCAACGGTAGACCGTTGAATGATTAACGTTTAcacc                                                                                   <  1:206514/141‑1 (MQ=255)

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

                                                            tttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCATCGGTAGACCGTTGAATGATGAACGTTTAcacca                                                                                  <  1:944499/141‑1 (MQ=255)

                                                            tttCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcacca                                                                                  <  2:119460/141‑1 (MQ=255)

                                                            tttCCTTTTAGTTTGATGTACGTCCCATCAATACGCCATTTGTAATAAGCTCTTTTATGCTTGTTCTTCCAAAGTTGTTATAAAATTGTGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacca                                                                                  <  2:30344/141‑1 (MQ=11)

                                                             ttCCTTTTATTTTGATCTATGTCTTATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCATTGAATGATTAACGTTTAcaccac                                                                                 >  2:1025207/1‑141 (MQ=11)

                                                             ttCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcaccac                                                                                 >  1:956235/1‑141 (MQ=255)

                                                              tCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACg                                                                                <  1:596207/141‑1 (MQ=18)

                                                              tCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACg                                                                                <  1:557486/141‑1 (MQ=18)

                                                              tCCTTTTATTTTGATCTACGTCTCATCAATCACCCGTTTGTCTTACGAGTTTTTATGCTTTTTCTTGATGCATTTTGATAAACCTAGGCCGAACTCCGGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACg                                                                                <  1:1047136/141‑1 (MQ=2)

                                                              tCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAGTGATGAACGTTTACACCACg                                                                                <  1:647777/141‑1 (MQ=255)

                                                               ccTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                                                                                         <  2:693318/131‑1 (MQ=255)

                                                                cTTTTATTTTGATGTACGTCTCATCAATACGACATTTGTATTAAGTTTTTTTCTGCTTTTTATTCCAAATTTTTTATAAAATTGCGTCATATTTTTGAACCCAACGGTGGACCGTTGAATGATGAACGTTTACACCACGtt                                                                              >  2:509585/1‑141 (MQ=11)

                                                                cTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt                                                                              >  1:473547/1‑141 (MQ=255)

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

                                                                 ttttATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTc                                                                             >  2:246767/1‑141 (MQ=255)

                                                                     aTTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                         <  1:105657/141‑1 (MQ=18)

                                                                      ttttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTa                                                                        <  2:742695/141‑1 (MQ=11)

                                                                      ttctGATTTACGTCTCATCAATACGCTATTTTTAATAAGCTTTTTTATGCGTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCTTTGAATGATGAACGTTTACACCACGTTCTCTTa                                                                        <  1:830307/138‑1 (MQ=11)

                                                                        ttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAat                                                                      <  2:62173/141‑1 (MQ=21)

                                                                          ggcaaggatctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCCTAatat                                                                    >  2:275609/9‑141 (MQ=11)

                                                                          gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAGCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTCatat                                                                    >  2:418027/5‑141 (MQ=11)

                                                                          gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGTTCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                    >  1:619626/5‑141 (MQ=11)

                                                                              ttcGTCTCATCAATACGCCATTTATAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCa                                                                <  1:712820/139‑1 (MQ=11)

                                                                              tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATGTTTCa                                                                >  1:523714/1‑141 (MQ=255)

                                                                              tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCa                                                                >  1:976893/1‑141 (MQ=40)

                                                                               aCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAg                                                               <  2:80253/141‑1 (MQ=21)

                                                                               aCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTATGCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAg                                                               <  2:1000901/141‑1 (MQ=255)

                                                                                cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGa                                                              >  2:972092/1‑141 (MQ=21)

                                                                                cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATCTTTCAGa                                                              <  1:661745/141‑1 (MQ=255)

                                                                                 gTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGat                                                             <  1:904550/141‑1 (MQ=255)

                                                                                     catcaATACGCCATTTGTAATAATCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCATGTTCTCTTAATATTTCCGatatat                                                         <  1:674521/141‑1 (MQ=11)

                                                                                      atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                        <  2:442966/141‑1 (MQ=18)

                                                                                        caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATccc                                                      >  1:394322/1‑139 (MQ=255)

                                                                                             cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATag                                                 <  2:94643/141‑2 (MQ=255)

                                                                                             cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                 <  1:1016137/141‑1 (MQ=21)

                                                                                               ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGACGATCGTTGAAAGATGACCGTTTACACCACCATCTCTTAAACTTTCAAAATTATCACGAtggct                                               >  1:239280/1‑137 (MQ=11)

                                                                                                cATTTGTCATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCAGGTCCCCTTACTGTTTCTGATATATCACGATCACTc                                              >  2:596130/1‑141 (MQ=11)

                                                                                                  tctGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                            <  1:281844/139‑1 (MQ=11)

                                                                                                    tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATCTTTCCGATATATCACGATAGCTCAATg                                          >  2:345732/1‑141 (MQ=11)

                                                                                                    tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                          <  2:366991/141‑1 (MQ=21)

                                                                                                     gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                         <  1:1001540/141‑2 (MQ=11)

                                                                                                     gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATTACGATAACTCAATgc                                         <  1:789018/141‑2 (MQ=12)

                                                                                                     gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATA‑TAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAGCTCAATgca                                        >  1:858045/1‑139 (MQ=11)

                                                                                                      taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                        <  1:120540/141‑3 (MQ=11)

CGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATC  >  CP000730/36238‑36518

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