BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,130,671	T→C	D186G (GAT→GGT)	yghR ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,130,671	0	T	C	83.8%	709.1 / 134.3	315	D186G (GAT→GGT)	yghR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (26/25); new base C (144/120); total (170/145)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 318 total aligned reads displayed.

TGCCGTTATAGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCAGAAACATGCTGATT  >  NC_000913/3130579‑3130760

   cGTTATAGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCaa                                                                                        <  2:3585344/93‑1 (MQ=255)

   cGTTATAGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCaa                                                                                        >  1:1535082/1‑93 (MQ=255)

    gTTATAGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTAAATACCAAt                                                                                       >  1:3136827/1‑93 (MQ=255)

     ttATAGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATg                                                                                      >  2:5317690/1‑93 (MQ=255)

     ttATAGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATg                                                                                      >  2:4477937/1‑93 (MQ=255)

     ttATAGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATg                                                                                      <  1:2618977/93‑1 (MQ=255)

       ataGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGcc                                                                                    <  1:5019509/93‑1 (MQ=255)

        taGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGccc                                                                                   <  1:1033742/93‑1 (MQ=255)

        taGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGccc                                                                                   <  1:2809942/93‑1 (MQ=255)

         aGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACTAATGCCCa                                                                                  >  1:4512094/1‑93 (MQ=255)

         aGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCa                                                                                  <  2:7297544/93‑1 (MQ=255)

         aGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCa                                                                                  >  1:866007/1‑93 (MQ=255)

          gCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAg                                                                                 >  1:6124290/1‑93 (MQ=255)

          gCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAg                                                                                 >  1:2972275/1‑93 (MQ=255)

           ccAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAgg                                                                                <  2:3666034/93‑1 (MQ=255)

           ccAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAgg                                                                                >  2:3005373/1‑93 (MQ=255)

            cAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGc                                                                               >  2:7521751/1‑93 (MQ=255)

            cAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGc                                                                               >  2:2740539/1‑93 (MQ=255)

            cAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGc                                                                               >  1:2267742/1‑93 (MQ=255)

             aaTCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCg                                                                              >  2:2831713/1‑93 (MQ=255)

             aaTCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCg                                                                              >  2:2434191/1‑93 (MQ=255)

              aTCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGa                                                                             >  1:6794749/1‑93 (MQ=255)

              aTCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGa                                                                             >  1:3549344/1‑93 (MQ=255)

              aTCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGa                                                                             >  1:1261124/1‑93 (MQ=255)

               tCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGaa                                                                            >  2:3124332/1‑93 (MQ=255)

                 ttCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGaata                                                                          >  2:224724/1‑93 (MQ=255)

                  tCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGaataa                                                                         >  2:405672/1‑93 (MQ=255)

                   cGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGaataat                                                                        <  2:2457422/93‑1 (MQ=255)

                   cGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGaataat                                                                        <  1:6334033/93‑1 (MQ=255)

                    gACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAAtaatt                                                                       <  2:861774/93‑1 (MQ=255)

                      cATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAAtaattaa                                                                     <  1:1168812/93‑1 (MQ=255)

                      cATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAAtaattaa                                                                     >  1:4087653/1‑93 (MQ=255)

                       aTGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAAtaattaat                                                                    <  1:4686539/93‑1 (MQ=255)

                       aTGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAAtaattaat                                                                    <  1:5862385/93‑1 (MQ=255)

                       aTGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAAtaattaat                                                                    <  2:5213166/93‑1 (MQ=255)

                        tGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAATaattaatt                                                                   <  2:6658391/93‑1 (MQ=255)

                        tGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATaattaatt                                                                   >  2:3559666/1‑93 (MQ=255)

                         gACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTc                                                                  <  1:5309422/93‑1 (MQ=255)

                         gACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTc                                                                  <  2:6921586/93‑1 (MQ=255)

                         gACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTc                                                                  <  1:3304635/93‑1 (MQ=255)

                          acacCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTcc                                                                 >  2:5610236/1‑93 (MQ=255)

                          acacCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTcc                                                                 >  1:6612250/1‑93 (MQ=255)

                           cacCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCg                                                                >  1:1579144/1‑93 (MQ=255)

                           cacCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCg                                                                <  2:4825562/93‑1 (MQ=255)

                            acCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCgg                                                               >  1:6015550/1‑93 (MQ=255)

                            acCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCgg                                                               >  1:6745072/1‑93 (MQ=255)

                            acCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCgg                                                               <  1:1688876/93‑1 (MQ=255)

                            acCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCgg                                                               >  2:1268374/1‑93 (MQ=255)

                             ccGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCggg                                                              >  2:1037833/1‑93 (MQ=255)

                               gATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAATAATTAATTCCGGGCg                                                            >  1:3500529/1‑93 (MQ=255)

                               gATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCg                                                            >  1:4948919/1‑93 (MQ=255)

                                 ttttGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGat                                                          <  2:3336047/93‑1 (MQ=255)

                                  tttGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGata                                                         <  1:3961540/93‑1 (MQ=255)

                                  tttGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCATGCGAATAATTAATTCCGGGCGata                                                         <  2:7362972/93‑1 (MQ=255)

                                  tttGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGata                                                         <  2:5969229/93‑1 (MQ=255)

                                  tttGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGata                                                         <  2:3755411/93‑1 (MQ=255)

                                   ttGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGatat                                                        <  1:2380127/93‑1 (MQ=255)

                                       ccTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTgcg                                                    >  2:2072933/1‑93 (MQ=255)

                                         tGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGcc                                                  <  1:2501669/93‑1 (MQ=255)

                                         tGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGcc                                                  >  1:3799351/1‑93 (MQ=255)

                                          gCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCa                                                 >  1:2007113/1‑93 (MQ=255)

                                           cAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCAt                                                <  1:5124479/93‑1 (MQ=255)

                                            aGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCAtt                                               >  2:1651586/1‑93 (MQ=255)

                                             gCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCAttt                                              <  1:2129581/93‑1 (MQ=255)

                                              cTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATttg                                             <  1:6419526/93‑1 (MQ=255)

                                               tCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATttgt                                            <  2:4744714/93‑1 (MQ=255)

                                                cGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATttgtt                                           <  2:5148347/93‑1 (MQ=255)

                                                 ggCATCGTCATGGGCAGGCTTGCGGGCGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGGTATTGCGCCATttgttg                                          >  1:5698352/1‑93 (MQ=255)

                                                 ggCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATttgttg                                          <  1:5143595/93‑1 (MQ=255)

                                                 ggCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATttgttg                                          >  2:450508/1‑93 (MQ=255)

                                                  gCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTgg                                         <  2:6372510/93‑1 (MQ=255)

                                                   cATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGt                                        <  1:3653785/93‑1 (MQ=255)

                                                     tAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGtat                                      <  2:352639/93‑1 (MQ=255)

                                                     tAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGtat                                      >  1:807446/1‑93 (MQ=255)

                                                     tAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGtat                                      >  2:4944097/1‑93 (MQ=255)

                                                      aGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGtata                                     <  1:4282440/93‑1 (MQ=255)

                                                       gTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATaa                                    >  2:4163141/1‑93 (MQ=255)

                                                        tCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATaaa                                   <  1:5463716/93‑1 (MQ=255)

                                                        tCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATaaa                                   <  2:5536194/93‑1 (MQ=255)

                                                         cATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAAc                                  <  1:7227385/93‑1 (MQ=255)

                                                               cAGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGtt                            <  2:4004906/93‑1 (MQ=255)

                                                                aGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTc                           <  1:4814084/93‑1 (MQ=255)

                                                                aGGCTTGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTc                           <  2:6268339/93‑1 (MQ=255)

                                                                    ttGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAATCGCCGTTCgcgc                       <  2:1925827/93‑1 (MQ=255)

                                                                    ttGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCgcgc                       >  1:3014846/1‑93 (MQ=255)

                                                                     tGCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCt                      <  2:5657967/93‑1 (MQ=255)

                                                                      gCGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTg                     <  1:7139443/93‑1 (MQ=255)

                                                                       cGGGAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTgc                    <  2:6020824/93‑1 (MQ=255)

                                                                        gggAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTgcg                   >  2:6229400/1‑93 (MQ=255)

                                                                         ggAGATGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTgcgc                  >  1:6536576/1‑93 (MQ=255)

                                                                           agaTGGCAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCa                <  2:7123639/93‑1 (MQ=255)

                                                                                 cAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCAGAAACa          <  2:6526488/93‑1 (MQ=255)

                                                                                 cAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCAGAAACa          >  2:1115019/1‑93 (MQ=255)

                                                                                 cAGTCTCAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCAGAAACa          <  1:7445197/93‑1 (MQ=255)

                                                                                   gTCTCAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCAGAAACATg        >  2:1589320/1‑93 (MQ=255)

                                                                                     ctcAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCAGAAACATGCt      >  1:4442496/1‑93 (MQ=255)

                                                                                      tcAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCAGAAACATGCTg     >  1:5619794/1‑93 (MQ=255)

                                                                                      tcAATACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCAGAAACATGCTg     >  2:138183/1‑93 (MQ=255)

                                                                                        aaTACCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCAGAAACATGCTGAt   >  2:2366109/1‑93 (MQ=255)

                                                                                         atacCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCAGAAACATGCTGAtt  <  2:3818712/93‑1 (MQ=255)

TGCCGTTATAGCCAATCTTCGACATGACACCGATTTTGTCCTGCAGCTCGGCATAGTCATGGTCAGGCTTGCGGGAGATGGCAGTCTCAATATCAATGCCCAGGCGAATAATTAATTCCGGGCGATATTGCGCCATTTGTTGGTATAAACGCCGTTCGCGCTGCGCCAGAAACATGCTGATT  >  NC_000913/3130579‑3130760

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