BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,444	C→T	46.5%	R36R (AGG→AGA)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,444	0	C	T	46.5%	101.5 / inf	316	R36R (AGG→AGA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (89/79); new base T (64/82); total (154/162)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 321 total aligned reads displayed.

TTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTAT  >  CP000730/36304‑36584

tttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTcc                                                                                                                                              >  1:677631/1‑141 (MQ=21)

 ttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc                                                                                                                                             <  1:169602/141‑1 (MQ=21)

  tATTTTGATATACGTCTTATCAATACGCCATCTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATGTGATCTAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTctct                                                                                                                                            <  2:664153/141‑1 (MQ=11)

   tttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                                                           <  2:772140/140‑1 (MQ=21)

   aTTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGTGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACTTTTACACCACGTTCCCtt                                                                                                                                           <  2:29644/141‑1 (MQ=255)

   aTTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                                                           <  2:454947/141‑1 (MQ=21)

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

   aTTTTGATATACGTCTCAGCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                           <  1:82213/141‑1 (MQ=11)

   aTTGTGATGTACGTCGCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                                                           <  2:631671/141‑1 (MQ=255)

     tttGATGGACGTCGCTCTAAGACGCGTTTTGTAATAATGTGTTTTATGTGTTTTCTGCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTaa                                                                                                                                         <  2:592907/141‑1 (MQ=11)

     tttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTaa                                                                                                                                         <  2:364919/141‑1 (MQ=11)

     tttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGACGAACGTTTAGACCACGTTCTCTTaa                                                                                                                                         <  1:421359/141‑1 (MQ=11)

       tGATGTACCTCTCATCACTACGCCTTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAata                                                                                                                                       <  1:743249/141‑1 (MQ=255)

       tGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAata                                                                                                                                       <  1:727343/141‑1 (MQ=255)

       tGATATACGTCTAATCAATACGTCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAata                                                                                                                                       <  2:139969/141‑1 (MQ=11)

        gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAatat                                                                                                                                      <  1:563075/141‑1 (MQ=40)

        gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAatat                                                                                                                                      <  1:603263/141‑1 (MQ=40)

        gATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCCTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACAACGTTCTCTTAatat                                                                                                                                      >  1:78220/1‑141 (MQ=11)

         atatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGATCTCTTAATAtt                                                                                                                                     >  1:691497/4‑141 (MQ=11)

         aTGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAtt                                                                                                                                     <  2:708341/141‑1 (MQ=39)

            tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGATCCCTTAATATTTCa                                                                                                                                  >  1:51341/1‑141 (MQ=21)

             ttttttcatcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATTTTTTTTAACCCTACGGTAGACCGTTGAATTATGAACGTTTACACCACGTTCCCTTACTATTTAAg                                                                                                                                 >  2:535818/6‑141 (MQ=11)

             aCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAg                                                                                                                                 >  2:382129/1‑141 (MQ=255)

              cgtctacccaATACGCTATTTGTAATAAGCTTTGTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGTATATTCTTGAACCCAACGGTAGATCGCTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCCGa                                                                                                                                <  1:636782/133‑1 (MQ=11)

              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGa                                                                                                                                <  2:713009/141‑1 (MQ=255)

              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGa                                                                                                                                <  1:312332/141‑1 (MQ=21)

              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGa                                                                                                                                <  2:511577/141‑1 (MQ=18)

               gTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGCACCCACCGGTAGACCGTTGAATGATGAACGCTTACACTACGCTCACTTAATATTTCAGat                                                                                                                               <  2:830313/141‑1 (MQ=14)

               gTCTCATCAATACGCCATTGGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGGAATCAAAGTGGGGCATATTCTTGAACCCAACGGTAGAGCGTTGAATGATGAACGTTTCCACCCCGTTCCCTTGATATTTCAGat                                                                                                                               >  1:87059/1‑141 (MQ=11)

                tctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGata                                                                                                                              >  1:769212/1‑141 (MQ=255)

                 ctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatat                                                                                                                             <  1:274556/141‑1 (MQ=255)

                  tcatcaATACGCCATTTGTAATAAGCTTTTTTTTTTTTTTTTTTACCAATTTAATACAAAATTGAGGCATATTCTTGAACCCAACGGTAGACCGTTGTATTCTGAACGTTTATACCACGTTCCCTTAATATTTCAGatata                                                                                                                            >  2:57402/1‑141 (MQ=14)

                  tcatcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAACTTGAGGCATCTTCTTGAAGCCCACCGTAGAACCTTGAATGATGAACCTTTACACCACGACCCCTTAATCTTTCAGttata                                                                                                                            >  2:90053/1‑141 (MQ=11)

                  tcatcaATACGCCATTTGTAATAAGCTTATTTATGCTTTTTCTTCCAAATTTGACACAAAATTGGGGCATATTCTTGAAGCCAACGGTAGACCGTTGAATGATGAACGTTTACNCCANGNTCCCTTAATATTTCAAatata                                                                                                                            >  1:541288/1‑141 (MQ=14)

                   catcaATACGCCATTCGTAATAAGCTATTTTATGTTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCATCGGTAGACCGTTGTATGATGAACGTTTACACCACGTTTCCTTAATATTTCAGatatat                                                                                                                           <  1:508669/141‑1 (MQ=11)

                    atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCCAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATAGTTGAATTATGAACGTTTACACCACGGTCTCTTAATATTTTAGatatata                                                                                                                          >  2:105160/1‑140 (MQ=11)

                    atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                                                          <  2:135360/141‑1 (MQ=21)

                     tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                                                         >  2:736945/1‑141 (MQ=255)

                     tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                                                         >  1:254203/1‑141 (MQ=255)

                     cctaTACGCCAGTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTCCCACCCAGTTCTATTACTATTTCAGATATATCa                                                                                                                         <  2:701125/138‑1 (MQ=11)

                      caATCCCCCATTTGTAATAAGCTTTTGTATGCGTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCAc                                                                                                                        >  2:489743/1‑141 (MQ=11)

                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCAc                                                                                                                        >  1:267673/1‑141 (MQ=21)

                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCAc                                                                                                                        >  1:717764/1‑141 (MQ=21)

                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCAc                                                                                                                        >  1:756834/1‑141 (MQ=21)

                        aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAGTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGa                                                                                                                      <  1:358641/141‑1 (MQ=255)

                        aTACGCCATTTGTAACAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGa                                                                                                                      <  1:226734/141‑1 (MQ=21)

                          aCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTATTATTTCAGATATATCACGATa                                                                                                                    >  1:95206/1‑141 (MQ=255)

                          aCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATa                                                                                                                    <  1:130524/141‑1 (MQ=255)

                          aCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATa                                                                                                                    >  1:795335/1‑141 (MQ=40)

                          aCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATa                                                                                                                    >  2:685115/1‑141 (MQ=40)

                           ctccAATTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                   <  2:161173/139‑1 (MQ=255)

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

                            gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATagc                                                                                                                  <  1:43232/141‑3 (MQ=255)

                            cccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACCCCACGTTCCCTTAATATTTCAGGTATATCACGATAAc                                                                                                                  <  1:449238/140‑1 (MQ=255)

                             ccATTTGTAATAAGCTTTTTTATTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTTAGATATATCACgtttgct                                                                                                                 >  2:10875/1‑135 (MQ=11)

                             ccATTTGTAATAAGCTTTTTTATGCTTTTTTTTTTTTTTTTGATATAAAATTGGGGCATATTCTTTAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGTTCCCTTAATATTTCAGATATATCCCGATAACt                                                                                                                 >  2:782907/1‑141 (MQ=11)

                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACAACACGGTCTCTTCATATTTCAGATATATCACGATAGCt                                                                                                                 >  1:723407/1‑141 (MQ=11)

                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACt                                                                                                                 >  1:258046/1‑141 (MQ=21)

                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTCCAGATATATCACCATAAca                                                                                                                 >  1:103718/1‑140 (MQ=255)

                              cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                <  2:22411/141‑1 (MQ=18)

                              cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                >  1:387106/1‑141 (MQ=21)

                               aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATTTTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTACTATTTCAGATATATCACGCTAACTCa                                                                                                               >  2:540692/1‑141 (MQ=255)

                                tttGTAATCAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATGTATCACGATAACTCaa                                                                                                              <  1:416217/141‑1 (MQ=12)

                                tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                              >  2:314489/1‑141 (MQ=255)

                                tttGTAATAAGCTCATTTATGCTTCTTCTTCCAAATTTGATATAAAATTGGTGCATATTCTTGAACCCAACGGTAGACCGATGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                              >  2:75087/1‑141 (MQ=11)

                                 ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                             >  1:203371/1‑141 (MQ=11)

                                 ttGTAATAAGCTTGTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                             <  1:550557/141‑1 (MQ=17)

                                  tgcaataaGCTTTTTTATGCCTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                            <  1:156535/138‑1 (MQ=11)

                                  tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCACCGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                            >  2:258692/1‑141 (MQ=11)

                                  tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGCTGAACGTTTACACCACGGTCTCTTAATATTTCGGATATATCACGATAGGTCaagg                                                                                                            >  1:335451/1‑139 (MQ=11)

                                  tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                            <  2:331475/141‑1 (MQ=11)

                                  tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                            <  1:125845/141‑1 (MQ=18)

                                  tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTGATATTTCAGATATATCACGATAACTCAATg                                                                                                            >  2:522802/1‑141 (MQ=255)

                                  tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAGTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGCTGAACGTTTTCATCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                            <  1:704584/141‑1 (MQ=11)

                                   gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                           >  2:1219/1‑140 (MQ=255)

                                   gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTATTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                           <  1:463114/141‑2 (MQ=11)

                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                          >  2:484158/1‑139 (MQ=255)

                                    taataaGCTTTTTGATGCTTTTTCTTCCAAATTTCATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGTACTTTTACACCACGATCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                          <  1:325140/141‑1 (MQ=14)

                                     aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGCTAACTCAATGaat                                                                                                         >  2:150208/1‑141 (MQ=12)

                                      ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGaata                                                                                                        >  2:5688/1‑141 (MQ=17)

                                               tttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAta                                                                                               >  1:188222/1‑141 (MQ=35)

                                                  ttttttttttttACAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGAACGTTGAATGATGAACGTATCCACTACGTTCCCTTAATTTTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                            >  2:61471/5‑141 (MQ=11)

                                                             tcCAAATTTGATACAAACTTGGAAATGCTATTTTAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTa                                                                                 >  1:830863/1‑141 (MQ=11)

                                                                 aaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGt                                                                             <  2:556495/141‑1 (MQ=35)

                                                                    ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                          <  1:295031/141‑1 (MQ=35)

                                                                                 ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTTCAGTGATAACATCCTTGTTa                                                             <  2:595522/141‑1 (MQ=17)

                                                                                        atTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTTCTTCCCGTTGTCTTAATATTTCTGATCTCTTACGATAACAAAAACTACATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                      <  1:401720/141‑1 (MQ=11)

                                                                                          tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                    <  1:136086/141‑1 (MQ=32)

                                                                                          tCTGGAACCCAACGTGAGGCCCTTGAATGATGAACGTTTACTCCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAGTATTTATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                    <  1:650565/141‑1 (MQ=11)

                                                                                              gAACCCAACGGGAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACATTGATAACATCCTTGTTAAATGGTTTATATc                                                >  2:284525/1‑141 (MQ=12)

                                                                                                aCCCAAACGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTg                                              <  1:179927/141‑1 (MQ=17)

                                                                                                 cccAACGGTAGACCGTTGAATTCTGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTAGATCTGa                                             <  2:464436/141‑1 (MQ=12)

                                                                                                        gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                      >  2:645567/1‑141 (MQ=32)

                                                                                                        gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                      >  1:470750/1‑141 (MQ=32)

                                                                                                                gtGAATGATGAACGGTTAAACCCAGTTCTCTTAAAATTACCGGTATATCACGATCGTTCAATGTCGATCTGGGATATTAGACGTGGGCCACAGGTAGAATAGCCTTGGTATATTGGTTATATCTGAAATAGTTCATACaga                              <  2:198023/140‑1 (MQ=2)

                                                                                                                          aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                    <  1:614025/141‑1 (MQ=32)

                                                                                                                           aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACGTAGTGGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtttt                   >  1:478394/1‑141 (MQ=255)

                                                                                                                           aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtttt                   >  2:691288/1‑141 (MQ=32)

                                                                                                                                 tcaccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa             <  2:349532/140‑1 (MQ=17)

                                                                                                                                  caccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaa            >  1:211831/1‑141 (MQ=32)

TTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTAT  >  CP000730/36304‑36584

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