BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,049,149	G→A	G75S (GGC→AGC)	yihI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,049,149	0	G	A	99.2%	1372.4 / ‑5.0	377	G75S (GGC→AGC)	yihI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (2/1); new base A (180/194); total (182/195)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 383 total aligned reads displayed.

CGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACG  >  NC_000913/4049059‑4049238

 gcgcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCg                                                                                        <  1:2001394/93‑1 (MQ=255)

  cgcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGt                                                                                       <  1:4513101/93‑1 (MQ=255)

  cgcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGt                                                                                       <  2:5112868/93‑1 (MQ=255)

   gcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTg                                                                                      >  1:3152063/1‑93 (MQ=255)

    cAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGa                                                                                     <  1:5316050/93‑1 (MQ=255)

     aGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGAc                                                                                    <  2:278957/93‑1 (MQ=255)

        gggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                                 <  2:3780219/93‑1 (MQ=255)

        gggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                                 >  1:6799560/1‑93 (MQ=255)

           cggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaa                                                                              >  2:5721318/1‑93 (MQ=255)

            ggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaaa                                                                             >  1:898200/1‑93 (MQ=255)

            ggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaaa                                                                             <  1:7161031/93‑1 (MQ=255)

             gAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAg                                                                            <  2:5813414/93‑1 (MQ=255)

             gAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTAAGCGTGACTGAAAAAg                                                                            >  2:1902109/1‑93 (MQ=255)

              aaaCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAgc                                                                           >  1:6738139/1‑92 (MQ=255)

                acacCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCa                                                                         <  2:2317853/93‑1 (MQ=255)

                acacCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCa                                                                         >  1:3563930/1‑93 (MQ=255)

                  accacGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCAcc                                                                       >  1:2960342/1‑93 (MQ=255)

                   ccacGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCa                                                                      <  1:2075854/93‑1 (MQ=255)

                     acGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCaaa                                                                    <  1:5435209/93‑1 (MQ=255)

                     acGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCaaa                                                                    >  2:3891120/1‑93 (MQ=255)

                     acGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCaaa                                                                    <  1:2704737/93‑1 (MQ=255)

                       gTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACa                                                                  <  1:1573035/93‑1 (MQ=255)

                       gTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACa                                                                  <  2:5940341/93‑1 (MQ=255)

                       gTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACa                                                                  >  1:5991454/1‑93 (MQ=255)

                        tCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAg                                                                 >  2:2550092/1‑93 (MQ=255)

                          aGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGca                                                               <  2:1163800/93‑1 (MQ=255)

                          aGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGca                                                               <  1:6902290/93‑1 (MQ=255)

                            gcagcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAACCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGcaca                                                             <  2:7175909/93‑1 (MQ=255)

                              agcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACaaa                                                           <  1:99704/93‑1 (MQ=255)

                              agcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACaaa                                                           >  2:7108965/1‑93 (MQ=255)

                              agcaAAGGCCAGAACGCACAAAAAGATCCACGCATTGGCAGGAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACaaa                                                           >  2:3938400/1‑93 (MQ=255)

                               gcaAAGTCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAc                                                          <  1:5771670/93‑1 (MQ=255)

                               gcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAc                                                          >  1:7907133/1‑93 (MQ=255)

                               gcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAc                                                          <  2:6332813/93‑1 (MQ=255)

                                caAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAcc                                                         >  2:4433682/1‑93 (MQ=255)

                                caAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAcc                                                         <  1:6606450/93‑1 (MQ=255)

                                  aaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGa                                                       >  2:1368856/1‑93 (MQ=255)

                                  aaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGa                                                       <  2:6179950/93‑1 (MQ=255)

                                  aaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGa                                                       >  1:7483199/1‑93 (MQ=255)

                                   aGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                                                      <  2:3620812/93‑1 (MQ=255)

                                     gCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAACCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAaga                                                    <  2:3330791/93‑1 (MQ=255)

                                     gCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAaga                                                    <  1:3652240/93‑1 (MQ=255)

                                     gCCAGAACGCACCAAAAGATCCACGTATGGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAaga                                                    <  1:5128108/93‑1 (MQ=255)

                                      ccAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAagag                                                   <  2:6527846/93‑1 (MQ=255)

                                       cAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGt                                                  <  1:1754516/93‑1 (MQ=255)

                                       cAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGt                                                  <  1:3244389/93‑1 (MQ=255)

                                       cAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGt                                                  <  1:4855921/93‑1 (MQ=255)

                                         gAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTga                                                <  2:3526056/93‑1 (MQ=255)

                                         gAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTga                                                >  2:7394735/1‑93 (MQ=255)

                                          aaCGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTgag                                               >  2:3001725/1‑93 (MQ=255)

                                             gCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGaaa                                            >  2:5386145/1‑93 (MQ=255)

                                             gCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGaaa                                            >  1:6867669/1‑93 (MQ=255)

                                              cACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAc                                           >  2:4706943/1‑93 (MQ=255)

                                                ccAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCt                                         >  1:4442999/1‑93 (MQ=255)

                                                 cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                                        <  2:7497331/93‑1 (MQ=255)

                                                 cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                                        >  2:2710243/1‑93 (MQ=255)

                                                  aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGa                                                       >  1:1214477/1‑77 (MQ=255)

                                                  aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                                       >  1:5528940/1‑93 (MQ=255)

                                                  aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                                       >  2:4957310/1‑93 (MQ=255)

                                                    aaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGc                                     <  2:6083680/93‑1 (MQ=255)

                                                      gATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCtt                                   >  1:220275/1‑93 (MQ=255)

                                                       aTCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCttt                                  >  2:1435785/1‑93 (MQ=255)

                                                         ccACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACTGCACAAACCGAAGAGTGAGAAACCTATGCTTTCa                                <  1:6414504/93‑1 (MQ=255)

                                                         ccACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCa                                >  2:6029916/1‑93 (MQ=255)

                                                         ccACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCa                                <  1:3828507/93‑1 (MQ=255)

                                                          cACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCAc                               <  2:5259883/93‑1 (MQ=255)

                                                          cACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCAc                               <  1:5619365/93‑1 (MQ=255)

                                                          cACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCAc                               >  1:534033/1‑93 (MQ=255)

                                                          cACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCAc                               <  2:45251/93‑1 (MQ=255)

                                                             gTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGc                            >  2:3192558/1‑93 (MQ=255)

                                                              tATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCa                           <  1:4747402/93‑1 (MQ=255)

                                                              tATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCa                           >  1:1102670/1‑93 (MQ=255)

                                                               aTTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAg                          <  2:5632790/93‑1 (MQ=255)

                                                               aTTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAg                          >  2:4224861/1‑93 (MQ=255)

                                                                ttGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAgg                         <  2:829848/93‑1 (MQ=255)

                                                                ttGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAgg                         >  1:7068244/1‑93 (MQ=255)

                                                                  ggCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCg                       >  1:5028915/1‑93 (MQ=255)

                                                                    cAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGa                     >  1:1129801/1‑93 (MQ=255)

                                                                    cAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGa                     <  1:4243053/93‑1 (MQ=255)

                                                                     aGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAg                    >  2:2484154/1‑93 (MQ=255)

                                                                      gTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGt                   >  1:623127/1‑93 (MQ=255)

                                                                       tAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGtt                  <  2:4132499/93‑1 (MQ=255)

                                                                       tAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGtt                  >  2:6692535/1‑93 (MQ=255)

                                                                         aacccccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTgg                <  1:6937023/90‑1 (MQ=255)

                                                                          aaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGa               >  1:2370203/1‑93 (MQ=255)

                                                                          aaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGa               >  2:1196860/1‑93 (MQ=255)

                                                                            ccccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAt             >  1:25245/1‑93 (MQ=255)

                                                                            ccccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTCAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAt             <  2:6871637/93‑1 (MQ=255)

                                                                             cccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAtt            >  2:5334526/1‑93 (MQ=255)

                                                                             cccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAtt            <  2:7625258/93‑1 (MQ=255)

                                                                             cccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAAtt            >  2:5553844/1‑93 (MQ=255)

                                                                              ccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTa           >  1:3338895/1‑93 (MQ=255)

                                                                               ctgttCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTAc          >  1:1921601/4‑93 (MQ=255)

                                                                               ctgttCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTAc          >  2:173150/4‑93 (MQ=255)

                                                                                tgttCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACt         <  2:2912394/91‑1 (MQ=255)

                                                                                    ccATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGaa     >  1:2253411/1‑93 (MQ=255)

                                                                                    ccATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGaa     >  2:1547977/1‑93 (MQ=255)

                                                                                    ccATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGAGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGaa     >  1:2806257/1‑93 (MQ=255)

                                                                                     cATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGaaa    <  1:2656318/93‑1 (MQ=255)

                                                                                       ttGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACg  >  2:2199875/1‑93 (MQ=255)

CGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACG  >  NC_000913/4049059‑4049238

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