BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,997	T→C	30.7%	T185T (ACA→ACG)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	35,997	0	T	C	30.7%	1240.8 / inf	1062	T185T (ACA→ACG)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (365/363); new base C (177/146); total (548/514)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 1062 total aligned reads displayed.

CATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  CP000730/35862‑36136

caggTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGATGTAtt                                                                                                                                        >  2:984586/4‑141 (MQ=0)

     cAGGGTCCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATGCCTTTTAAAGTATTCTTTGCTGTATTGATAc                                                                                                                                   <  1:42434/141‑1 (MQ=2)

      aGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAACTATTCTTTGCTGTATTGATACt                                                                                                                                  <  2:807008/141‑1 (MQ=14)

       gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt                                                                                                                                 <  1:219355/141‑1 (MQ=255)

        tgtTCGTTTAACTTGCTAGCATGATGCAAATGTCGTGGCATGGGGAAAATCCGTAGATCTGAGGAGACCTACGGTTCTTTGAATATAGAGCGTAAATACATTCTATACCTTTTAAAGTATTCTTTGCCGTATTGATACttt                                                                                                                                >  2:1579/1‑141 (MQ=11)

        tgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGCATTGATACttt                                                                                                                                >  1:90492/1‑141 (MQ=17)

         gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTCCTGTATTGATACTTTg                                                                                                                               >  2:40449/1‑141 (MQ=12)

           tCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAt                                                                                                                             >  2:727684/1‑141 (MQ=32)

            cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                                                                                                            <  1:957263/138‑1 (MQ=17)

            cgtttAACTTGCTAGCAGGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                                                                                                            <  1:974381/138‑1 (MQ=12)

             ttttAACGTGCTAGCATGATTCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATat                                                                                                                           <  2:109242/139‑2 (MQ=11)

             gtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGTAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATat                                                                                                                           <  2:1052443/139‑2 (MQ=11)

              tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGGGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATCCCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATatc                                                                                                                          <  2:68348/140‑3 (MQ=11)

              tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                                                                                                          <  1:542492/140‑1 (MQ=17)

               ttAACTTGCTCGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATatct                                                                                                                         <  2:78845/141‑4 (MQ=11)

               ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATatct                                                                                                                         <  1:934534/141‑4 (MQ=255)

               ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCt                                                                                                                         <  2:65966/141‑1 (MQ=17)

                tAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCtt                                                                                                                        >  1:297594/1‑141 (MQ=17)

                 aaCTTGCTAGCATGATGCTAATTTCGTGGCACGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTg                                                                                                                       >  2:488609/1‑141 (MQ=17)

                  aCTTGCTAGCATGATGCAAAAACGGCGGCAGGGCGAAAATACGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGCTACTTTGATACATTGt                                                                                                                      >  1:813886/1‑141 (MQ=11)

                    ttGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtct                                                                                                                    <  1:191929/141‑1 (MQ=11)

                          gCATGATGCTAATTACGTGACATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTa                                                                                                              >  1:720022/1‑141 (MQ=11)

                             tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACtt                                                                                                           >  2:610248/1‑141 (MQ=11)

                               atgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTa                                                                                                         >  1:773849/1‑141 (MQ=25)

                                 gCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAat                                                                                                       >  1:996715/1‑141 (MQ=255)

                                 gCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAat                                                                                                       >  1:609332/1‑141 (MQ=21)

                                 gCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAat                                                                                                       >  2:593832/1‑141 (MQ=25)

                                   tAATTTCGTGGCATGGGGCAAATCACCCAATGTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATAATTTGATATCTTGTCTTTTTTACTTTAatat                                                                                                     >  2:507157/1‑141 (MQ=11)

                                   tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGGGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATAGTTTGATACATTGTCTATCTTACTTTAatat                                                                                                     >  1:210564/1‑141 (MQ=11)

                                    aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                                                                                                    >  2:452046/1‑141 (MQ=25)

                                        tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATCCCTTTCGTGGTGGTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg                                                                                                >  1:254500/1‑141 (MQ=11)

                                          gTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTTTCTTTCTTACTTTAATATGACGGTg                                                                                              >  1:172562/1‑141 (MQ=11)

                                          gTGGCATGGCGAAACTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATTACGGtt                                                                                              >  2:376450/1‑140 (MQ=255)

                                                tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                                                                        <  2:773669/141‑1 (MQ=255)

                                                 ggCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc                                                                                       <  1:409312/141‑1 (MQ=255)

                                                 ggCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc                                                                                       <  1:562917/141‑1 (MQ=255)

                                                       aaTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATg                                                                                 <  1:829626/141‑1 (MQ=255)

                                                          ccGTAGATCTGAAGAGAACTGCGGTTCTTTTTATATAGAGCGTAAATACATTCACTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGCCGGTGATCTTGCTCCATgggg                                                                              >  2:748728/1‑141 (MQ=11)

                                                           cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGt                                                                             <  1:806727/141‑1 (MQ=255)

                                                           cGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTACAGTATTCTTTGCCGTATTGCTAATTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATCAGGt                                                                             >  1:1054984/1‑141 (MQ=11)

                                                            gTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTATTACTGTAATATGGCGGTGATCTTGCTCGATGAgggt                                                                            >  1:617507/1‑139 (MQ=255)

                                                                gtCTGAAGAGGCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTc                                                                        <  2:1071606/140‑1 (MQ=255)

                                                                  ttGAAGGGACCTGCGGGTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAg                                                                      <  2:259684/140‑1 (MQ=255)

                                                                  cTGAAGATACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAg                                                                      <  2:234560/141‑1 (MQ=255)

                                                                   tGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGa                                                                     <  1:143822/141‑1 (MQ=255)

                                                                     aaGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATTAAGTTATTCAGata                                                                   >  1:281653/1‑141 (MQ=11)

                                                                       gggcccTGCCGGTCTTGGTATAGAGAGCGTAAATACACTCAATACCTTTGAAAGGCTTCGTTGCTGGATGTATACTCTGATACCTTGTCTTTCTTACGTTAATATGACTGTGATCTTGAGTAATAAGGTTATTCACATAtt                                                                 <  2:534884/137‑1 (MQ=2)

                                                                       gagaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAtt                                                                 >  1:1101521/1‑141 (MQ=255)

                                                                       gagaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAtt                                                                 <  1:878406/141‑1 (MQ=38)

                                                                          aCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGGATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCg                                                              >  1:46326/1‑141 (MQ=255)

                                                                              gctgTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTGAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGt                                                          <  2:831959/138‑1 (MQ=255)

                                                                              gCGGTTCTTTTTATATAGTGCGTAAATACACTCAATACCTTTTAATGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGt                                                          <  2:874700/141‑1 (MQ=255)

                                                                               cGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATCCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTa                                                         <  1:15668/141‑1 (MQ=255)

                                                                                 gTTTTTTTTATAGAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACa                                                       <  2:290359/141‑1 (MQ=255)

                                                                                  ttCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACaa                                                      >  2:215937/1‑141 (MQ=38)

                                                                                   tCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCGTGCTCAATGAGGTTATTCAGATATTTCGATGTACAAt                                                     <  2:182509/141‑1 (MQ=255)

                                                                                     tttttACATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCGCAATGAGGTTATTCAGATATTTCGATGTACAATgg                                                   <  1:400664/141‑2 (MQ=255)

                                                                                       tttATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACa                                                 <  2:960125/141‑1 (MQ=255)

                                                                                       tttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATggcc                                                 >  2:413303/1‑138 (MQ=255)

                                                                                       tttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAAGGAGGTTATTCAGATATTTCGATGTACAATggca                                                 <  1:486081/141‑4 (MQ=255)

                                                                                        ttATATAGAGCCTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAg                                                <  1:1079790/141‑1 (MQ=255)

                                                                                         tatataGAGCGTAACTACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGTTCAATGAAGTTATTCAGATATTTCGATGTACAATGTCAGt                                               <  1:791011/141‑1 (MQ=14)

                                                                                         tatataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGt                                               <  1:439630/141‑1 (MQ=255)

                                                                                          atataGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTc                                              >  1:1001815/1‑141 (MQ=255)

                                                                                          atataGAGCGTAAATACATTCAATACCTTTTAAAGTATTCCTTGCCGTATCGATACTTTGATATCTTGTCTTTCTTACTTGAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTc                                              <  2:116565/141‑1 (MQ=14)

                                                                                          atataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCGCGCGCAGTGAGGTTATTCAGATATTTCGATGTACCATGGCAGTc                                              >  1:853121/1‑141 (MQ=14)

                                                                                               gagCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtt                                         <  1:34450/141‑1 (MQ=255)

                                                                                               gagCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtt                                         >  2:349254/1‑141 (MQ=255)

                                                                                                   gTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaag                                     <  1:237623/141‑1 (MQ=255)

                                                                                                     aaaTACATTCAATACCTTTTAAAGTATTCTTTGCCGTATAGATCCTTTGATAACTTGCATTTGTTACTTTAACATGTTGCTGTTCTTGCTCAATGAGTATATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagtt                                   >  2:921568/1‑141 (MQ=14)

                                                                                                      aacccactCAATACCGTTTAACTGTTCCTGTGCTGGATTGATGCTTTGATACCTTGTCTTTGTTACTTTAATAGGACGGTGAGCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAAGGGCAGTCAGgtttaagttt                                  <  2:857366/134‑1 (MQ=14)

                                                                                                      aaTACACTCAAGACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                  <  2:576270/141‑1 (MQ=255)

                                                                                                       aTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGTGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttta                                 >  1:733432/1‑141 (MQ=255)

                                                                                                        tACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtttaa                                >  1:121967/1‑141 (MQ=255)

                                                                                                         acacTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTaaa                               <  1:1013421/141‑1 (MQ=255)

                                                                                                          cactCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTaaaa                              >  1:470988/4‑141 (MQ=255)

                                                                                                          cATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTaaaa                              >  1:58578/1‑141 (MQ=255)

                                                                                                            ctCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGc                            >  1:632891/2‑141 (MQ=255)

                                                                                                             tCAATACCTTTTAAATTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCt                           <  1:649845/141‑1 (MQ=255)

                                                                                                             tCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCt                           <  2:627460/141‑1 (MQ=255)

                                                                                                              cAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCtt                          >  1:508363/1‑141 (MQ=255)

                                                                                                              cAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCtt                          <  1:371930/141‑1 (MQ=255)

                                                                                                               aaTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCttt                         >  2:138006/1‑141 (MQ=255)

                                                                                                               aaTACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCttt                         >  1:321143/1‑141 (MQ=255)

                                                                                                               aaTACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTATAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCttt                         >  1:677896/1‑141 (MQ=255)

                                                                                                                aTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTa                        >  2:1029447/1‑141 (MQ=255)

                                                                                                                   ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                     >  2:934839/1‑141 (MQ=255)

                                                                                                                   ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGCACAATGGCCGTCAGGTTTAAGTTTAAAAGCTTTAAtt                     <  1:699152/141‑1 (MQ=255)

                                                                                                                   ccTTTTAAAGTATTCCTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                     <  2:903773/141‑1 (MQ=255)

                                                                                                                     ttttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTAc                   <  1:754825/141‑1 (MQ=255)

                                                                                                                     ttttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCtttatttat                   >  2:439086/1‑140 (MQ=255)

                                                                                                                       ttAATGTATTCTTGGCTGTATTGAGACTTTGATACCTTGTCTTTGTTACTTTAATCTGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACtt                 <  2:310421/141‑1 (MQ=14)

                                                                                                                        tAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACttt                <  2:647569/141‑1 (MQ=255)

                                                                                                                        tAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTACGGGTTAAGTTTAAAAGCTTTAATTACttt                >  2:698464/1‑141 (MQ=255)

                                                                                                                        tAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAGGTTTAAAAGCTTTAATTACttt                >  2:334613/1‑141 (MQ=255)

                                                                                                                         aaaGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTAttttt               >  1:1027207/1‑140 (MQ=255)

                                                                                                                            gTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAAcc            >  2:472131/1‑141 (MQ=255)

                                                                                                                                  tttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTCATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCt      <  1:359076/141‑1 (MQ=255)

                                                                                                                                   ttGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCtg     >  2:194147/1‑140 (MQ=255)

                                                                                                                                     gCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAcc   <  1:517987/141‑1 (MQ=255)

CATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  CP000730/35862‑36136

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