BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,269	C→T	58.3%	A95T (GCA→ACA)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,269	0	C	T	58.3%	115.7 / inf	585	A95T (GCA→ACA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (97/145); new base T (153/184); total (252/333)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

Only 100 of 589 total aligned reads displayed.

TGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGT  >  CP000730/36129‑36409

  ctgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTAt                                                                                                                                            >  1:251027/4‑141 (MQ=255)

   tatcTTTGTTGAAGGTGCTTTATCTGTTATTACCTTTTGTGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGGATATGCTGAATGATTATCTCGTTGCTTACGGAACCAACTATCTAATGTATGTCCCTCTGTATc                                                                                                                                           <  1:113688/138‑1 (MQ=11)

   tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATc                                                                                                                                           <  2:279414/141‑1 (MQ=32)

   catcgtCGGTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATc                                                                                                                                           <  2:194458/136‑1 (MQ=255)

     ccTTCGTTGCAGGTGACTTATCTGCAATTACCTTTTGAAGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATATCGTTGCTTACGCAACCAAATAACTCATGTATGTCCCGCTGCATTaa                                                                                                                                         >  2:124268/1‑141 (MQ=11)

      ctacGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAAt                                                                                                                                        <  1:21375/138‑1 (MQ=17)

       ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                                                                                                                                       >  1:315331/1‑141 (MQ=17)

       aacGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATg                                                                                                                                       <  2:254114/139‑1 (MQ=17)

         cGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGc                                                                                                                                     <  2:16523/141‑1 (MQ=32)

          gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATTTATGTCCCTCTGCATCAATGGCa                                                                                                                                    >  2:36669/1‑141 (MQ=17)

          gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCa                                                                                                                                    <  2:8531/141‑1 (MQ=32)

           tggAAGGTGCCTGCTCTGTTATTACCTTTTTAGGTTTACTAAATTTTTTAATCAGTCGTTTGATAAACGCATATGCTGAATTATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                                                                                                                   <  1:140894/139‑1 (MQ=11)

             gAAGTTGCCTTATCTTTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGa                                                                                                                                 <  2:71767/141‑1 (MQ=12)

             gAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGa                                                                                                                                 >  1:6710/1‑141 (MQ=32)

               aGGTGCCTGATCTGTAATTGCCTTTTGAGGTTTACCAGATTGTTTAATGTGACGTTTGATAAATGCATATGCTGAATGATTATGTGGTTGGTTACGCAGCTAAATTTGTGATGTATGGTCCTCTGCATCAATGGAACgggt                                                                                                                               >  2:140241/1‑138 (MQ=11)

                    ccTGATCTGTAATTACTCTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAAt                                                                                                                          <  1:864/141‑1 (MQ=17)

                       taTCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc                                                                                                                       <  1:255400/140‑1 (MQ=255)

                             tAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTATGCAAAAAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtt                                                                                                                 >  1:74440/1‑141 (MQ=255)

                                  aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCTTATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCt                                                                                                            <  1:298355/141‑1 (MQ=17)

                                  aCCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt                                                                                                            <  1:153677/141‑1 (MQ=17)

                                   ccTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt                                                                                                           <  2:351846/141‑1 (MQ=17)

                                    cTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATCTAAATAGCTCCATTTTCCttt                                                                                                          <  1:205939/141‑1 (MQ=12)

                                        tGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGACACGATATAAATAGCACCATTTTCCTTTTAtt                                                                                                      >  1:1942/1‑141 (MQ=12)

                                        tGAGGCTTACCAAATTGTTTAATGAGGCGTTTGATAAACGCATATTCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtt                                                                                                      <  1:183741/141‑1 (MQ=255)

                                          tggTTTACCAAATTTGTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtttt                                                                                                    <  1:286298/140‑1 (MQ=17)

                                           ggTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTTGTTGCTTTCGCAACCAATTATCTAATTTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTg                                                                                                   >  1:3227/1‑141 (MQ=255)

                                             tttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAt                                                                                                 <  2:376015/141‑1 (MQ=17)

                                               tACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTGAGTCCCTCTGTATCAATGGCACGATATAAATAACTCCAGTTTCCTTTTATTTTGAtat                                                                                               >  2:135109/1‑139 (MQ=11)

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

                                                aCCAAATTGTTCAATGGGACGTTTGATAAACGCATATGCTGAATGATTATATCGTTGCTTACCCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAtata                                                                                              <  2:42566/141‑4 (MQ=11)

                                                 ccAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTAc                                                                                             <  1:373640/141‑1 (MQ=32)

                                                  cAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACg                                                                                            <  1:96408/141‑1 (MQ=32)

                                                     aTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAAACAAATATCAAATGTATGTGCCTCTCCATCAATAGCAAGATATAACTAGCAACATATCCCTTTACTTTTAAAGTatctcc                                                                                         >  2:174269/1‑136 (MQ=11)

                                                      ttGTTTAATGTGACGTTTGATAAACGCATATGCGGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGtctc                                                                                        <  2:200802/141‑1 (MQ=255)

                                                       tGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtca                                                                                       <  1:351307/141‑1 (MQ=32)

                                                        gTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcat                                                                                      <  1:327969/141‑1 (MQ=32)

                                                        gTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcat                                                                                      <  1:200971/141‑1 (MQ=32)

                                                         tttAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCtcatc                                                                                     <  1:222472/141‑1 (MQ=32)

                                                           tAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCaa                                                                                   <  2:29302/141‑1 (MQ=32)

                                                            tagtagaCGTTTGACAACCGCATATGCTGAACGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAACGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAAt                                                                                  <  2:189981/137‑1 (MQ=11)

                                                             aTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCtcattcata                                                                                 >  2:262786/1‑136 (MQ=255)

                                                               gagaCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACg                                                                               <  1:162997/141‑1 (MQ=32)

                                                                   cGTTTGATAAACGATTACTCTGAATGATTATCTCGTTGCTTACGCACCTAAATATTGAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTGTCATCAATACGCCAt                                                                           <  1:291638/141‑1 (MQ=11)

                                                                    gTTTGATAAACGCATATGCTGAATGATTACCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAACGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCAtt                                                                          <  2:203014/141‑1 (MQ=255)

                                                                       tGATAACCGCAGTTGCGGAATGATTCTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGAGATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGt                                                                       <  1:344292/141‑1 (MQ=11)

                                                                        gATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGta                                                                      >  1:360692/1‑141 (MQ=17)

                                                                        gATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGta                                                                      >  1:284345/1‑141 (MQ=32)

                                                                          tAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTAGGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGtaat                                                                    <  1:178846/141‑1 (MQ=12)

                                                                            atcGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGtaataa                                                                  <  2:297681/139‑1 (MQ=21)

                                                                              cGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATTTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCACTTGTAATAAGc                                                                >  2:211108/1‑141 (MQ=255)

                                                                              cGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGc                                                                <  1:259413/141‑1 (MQ=255)

                                                                              cGCATATGCTGAATGATTATCTCGTTGCTTACGCAAACAAATATCTAAAGTATGTCCATCTGCATCAATGGCACGATATAAATAAGTCCATTTTCATTTTATTTTGATGTACGTCTCATGACTAACACATTTGAAATAAGc                                                                >  2:242418/1‑141 (MQ=11)

                                                                               gCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCt                                                               <  2:316548/141‑1 (MQ=255)

                                                                                 atatGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATGTAAATAGCTCCATTTTCCTTTAATTTTGTTATGCGTCTCATCAATCCGCCATTTGTTATAAGCttt                                                             >  2:118063/1‑141 (MQ=11)

                                                                                    ggCTGAATGCTTATCTCGTTGCTCACCCAACCAAATATCTAATGTATTTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCtttttt                                                          <  2:95613/140‑1 (MQ=11)

                                                                                     gCTGAATGATTATCTCGCTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTa                                                         <  2:304893/141‑1 (MQ=255)

                                                                                      cTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCTCGTTTTAAATAGCTCAATTTTCATTTTATTTTGATATAAGTCTCATCAATACGCCAATAGTAATAAGCTTTTTTAt                                                        >  1:270513/1‑141 (MQ=11)

                                                                                      cTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTAt                                                        <  2:220507/141‑1 (MQ=38)

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

                                                                                       cgAATGATTATCTCGTTGCTTACGCAACCAAATATCTTATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATg                                                       <  1:305987/140‑1 (MQ=38)

                                                                                          aTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtt                                                    <  1:384695/141‑1 (MQ=255)

                                                                                          aTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtt                                                    <  1:149812/141‑1 (MQ=40)

                                                                                            gATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCtttt                                                  >  2:314492/1‑141 (MQ=40)

                                                                                             aTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATTCttttt                                                 >  2:33017/1‑141 (MQ=255)

                                                                                             aTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCttttt                                                 >  2:394462/1‑141 (MQ=40)

                                                                                              ttATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTCTATTTTGATATACGTCTCATCAATACGCCATTTGTAATGAGCTTTTTTATGCTTTttc                                                <  1:237573/141‑1 (MQ=255)

                                                                                                aTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttctt                                              >  1:243835/1‑141 (MQ=255)

                                                                                                 tctcGTGGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTttcttc                                             <  2:102889/141‑1 (MQ=255)

                                                                                                  ttcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTcc                                            <  1:48450/140‑1 (MQ=255)

                                                                                                  ctcGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTcc                                            <  1:101410/141‑1 (MQ=255)

                                                                                                   tcGTTGCTTGCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCa                                           <  2:61067/141‑1 (MQ=255)

                                                                                                    cGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTGATAAGCTTTTTTATGCTTTTTCTTCCaa                                          <  2:78047/141‑1 (MQ=255)

                                                                                                         cTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTGTTTTTTGCTTTTTCTTCCAAATTTg                                     >  1:218890/1‑141 (MQ=255)

                                                                                                            aCGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATa                                  <  1:168512/141‑1 (MQ=40)

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

                                                                                                             cGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                 >  2:146806/1‑140 (MQ=38)

                                                                                                             cGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATCTGTAATCACCTTTTTTCTGCTTTTTCTTTCAACTTTGATat                                 >  1:12493/1‑140 (MQ=11)

                                                                                                             cGCAACCAAATATCTAATGTATGTCCCACTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATat                                 <  2:386618/141‑2 (MQ=255)

                                                                                                              gCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTCTTTTATGCTTTTTCTTCCAAATTTGATata                                >  1:388003/1‑139 (MQ=255)

                                                                                                              gCAACCAAATATCTAATGTATGGCCCTCTGCATCTATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATata                                <  1:266596/141‑3 (MQ=255)

                                                                                                               cAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTTCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATataa                               <  2:270513/141‑4 (MQ=255)

                                                                                                                aaCCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaaa                              <  1:138941/141‑1 (MQ=255)

                                                                                                                     aaTATCTAATGTATGTCCTACGTTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgg                         >  2:382390/1‑141 (MQ=11)

                                                                                                                     aaTATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTGGTAATAAGCTTTTTTATGCGTTTTCTTCCAAATTTGGTATTAATTTgg                         >  1:60003/1‑141 (MQ=11)

                                                                                                                       tatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                       <  1:122114/141‑1 (MQ=255)

                                                                                                                       tatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggg                       <  2:248554/141‑1 (MQ=255)

                                                                                                                         tCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTAAACGATATGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgggggt                     >  1:188353/1‑139 (MQ=11)

                                                                                                                           tAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTCTTTTGATATACGTCTCATCAATACGACATTTTTAATAAGATTTTTTATGCTTTTTATTACAAATTTTATTTAAAATTGAGGAata                   >  1:105474/1‑141 (MQ=11)

                                                                                                                           tAATGTAAGTCCCTCTGTAACAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCata                   <  1:40599/141‑1 (MQ=11)

                                                                                                                            aaTGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGTatat                  >  1:334899/1‑141 (MQ=255)

                                                                                                                              tgtatgtCCCTCAGTATCAATGGCACGATATAAAAAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTCATAAGCTTTTTTCTGCTTTTTCTTCCACATCTGATCTCAAATTGGCGCACATTc                >  2:91112/1‑141 (MQ=11)

                                                                                                                               ttatgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCt               <  2:356395/140‑1 (MQ=255)

                                                                                                                               gtctgtCCCTCTGCATCTATGGCACGGTATAAATAGCTCCATTTTCCTTTTGTTTTGATGTACGTCTCATCAATACGCCTTTTGTATTAAGCTTTTTTCTCCTTTTTCTTCCGTATTTGATGTAAAATTGGGGCGTATGCt               >  2:324037/4‑141 (MQ=11)

                                                                                                                                 atgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTg             <  1:36669/141‑1 (MQ=25)

                                                                                                                                  tgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGa            <  2:48060/141‑1 (MQ=25)

                                                                                                                                   gtCCCTTTGGATCAATGGCCCGATATAAATAGCCCCCTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATATGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGaa           <  2:237472/141‑1 (MQ=11)

                                                                                                                                    tCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTGTTTCTTCCAACTTTGCTACAAAATTGGGGCTTATTCTTTAAc          >  2:337162/1‑141 (MQ=14)

                                                                                                                                      ccTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTCTTTATGCCTTTTCTTCCAAATTTGATATAAAATTGGGACATATTCTTGAAccc        >  2:41352/1‑141 (MQ=11)

                                                                                                                                       ctctGTATCAATGGCACGATATAAATCGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAACTTTGATGTAAAATTTGGCTCTATTCTTGAAcccc       >  2:56577/1‑140 (MQ=11)

                                                                                                                                           gCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGCGGCATATTCTTGAACCCAACgg   >  2:155238/1‑141 (MQ=21)

TGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGT  >  CP000730/36129‑36409

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