BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,431,925	A→G	I79V (ATC→GTC)	b4212 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,431,925	0	A	G	99.7%	1135.0 / ‑5.7	365	I79V (ATC→GTC)	b4212	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (180/184); total (180/185)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 381 total aligned reads displayed.

AGCGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAC  >  NC_000913/4431835‑4432014

 gCGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTcg                                                                                        <  2:6868355/93‑1 (MQ=255)

  cGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTcgc                                                                                       <  1:5677676/93‑1 (MQ=255)

   gACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTcgcg                                                                                      >  2:3591266/1‑93 (MQ=255)

   gACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTcgcg                                                                                      >  2:6288259/1‑93 (MQ=255)

    aCCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGt                                                                                     >  2:2031865/1‑93 (MQ=255)

       tGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATc                                                                                  <  2:4521385/93‑1 (MQ=255)

        gcggcgCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATcc                                                                                 <  1:2375881/93‑1 (MQ=255)

         cggcgCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCg                                                                                >  2:7395989/1‑93 (MQ=255)

          ggcgCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCgg                                                                               <  2:5100752/93‑1 (MQ=255)

           gcgcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggt                                                                              <  2:139014/93‑1 (MQ=255)

             gcAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCggtgg                                                                            <  2:1752015/93‑1 (MQ=255)

                   tCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTGGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgccgc                                                                      <  1:5942197/93‑1 (MQ=255)

                   tCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTgccgc                                                                      >  1:7070855/1‑93 (MQ=255)

                     ggTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGccgccg                                                                    >  1:1088970/1‑93 (MQ=255)

                       tGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCa                                                                  <  1:6670058/93‑1 (MQ=255)

                         gggTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAtg                                                                >  2:7571394/1‑93 (MQ=255)

                          ggTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCAtgt                                                               >  2:4746632/1‑93 (MQ=255)

                            tgagtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTgg                                                             >  1:6967500/1‑93 (MQ=255)

                            tgagtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTgg                                                             >  2:2416099/1‑93 (MQ=255)

                             gagtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGa                                                            >  2:3618259/1‑93 (MQ=255)

                              agtgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGaa                                                           <  1:5486135/93‑1 (MQ=255)

                               gtgaAAAGATGCTTGCGCAGTCGTTACAGGGGGTAGAACAGGATGGTTGTCTTAACCGTGTCGCGTATCCGGGGGTGCCGCCGCATGTGCaag                                                          >  1:5261568/1‑92 (MQ=255)

                                tgaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAata                                                         <  1:5134487/93‑1 (MQ=255)

                                 gaAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAatat                                                        <  1:7263533/93‑1 (MQ=255)

                                   aaaGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAg                                                      >  2:53887/1‑93 (MQ=255)

                                   aaaGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAg                                                      >  2:55912/1‑93 (MQ=255)

                                    aaGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGc                                                     >  1:2197822/1‑93 (MQ=255)

                                     aGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGcc                                                    <  2:1487328/93‑1 (MQ=255)

                                      gATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCt                                                   >  1:2910063/1‑93 (MQ=255)

                                      gATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCt                                                   >  1:6388455/1‑93 (MQ=255)

                                       aTGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTc                                                  >  2:1137609/1‑93 (MQ=255)

                                        tGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCa                                                 >  2:6057709/1‑93 (MQ=255)

                                        tGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCa                                                 >  2:6171353/1‑93 (MQ=255)

                                        tGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCa                                                 >  2:7281119/1‑93 (MQ=255)

                                         tcTTGCGCAGTCGTTACAGGCGTTAGAACAGGAGTGTTTTCTTAACCGTGTCGCGTAGCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAc                                                <  1:7622635/92‑1 (MQ=255)

                                         gCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAc                                                <  1:3940430/93‑1 (MQ=255)

                                         gCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAc                                                <  2:5448806/93‑1 (MQ=255)

                                           ttGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAcgc                                              <  2:5959165/93‑1 (MQ=255)

                                            tGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAcgcc                                             <  1:1714534/93‑1 (MQ=255)

                                            tGCGCAGTCGTTACAGGCGTTAGAACAGGAGGGCTTTCTTAACCGAGTCGCGTATCCGGTGTTGCCGCCGCATGTGAAATAAAGCCTCGcgcc                                             <  2:7438623/93‑1 (MQ=255)

                                              cgcAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAcgccgc                                           <  2:1865959/93‑1 (MQ=255)

                                              cgcAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAcgccgc                                           <  2:5179215/93‑1 (MQ=255)

                                              cgcAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAcgccgc                                           >  2:4168269/1‑93 (MQ=255)

                                              cgcAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAcgccgc                                           >  1:5258355/1‑93 (MQ=255)

                                               gcAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCAcgccgcg                                          >  1:825634/1‑92 (MQ=255)

                                               gcAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCt                                          <  1:6606081/93‑1 (MQ=255)

                                                cAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTg                                         >  1:7173231/1‑93 (MQ=255)

                                                  gTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTggg                                       <  1:4302999/93‑1 (MQ=255)

                                                  gTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTggg                                       >  1:1486559/1‑93 (MQ=255)

                                                  gTCGTTACAGGCGTTAGAACAGGATGGTTTTCTGAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTggg                                       <  1:3420615/93‑1 (MQ=255)

                                                  gTCGTTACAGGCGGTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCAGGTGGTGCCGCCGCATGTGGAATATAGCCTACAGCCGCTggg                                       >  2:7748955/1‑93 (MQ=255)

                                                   tCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGc                                      <  1:4682687/93‑1 (MQ=255)

                                                     gTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGa                                    >  1:3237757/1‑93 (MQ=255)

                                                     gTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGa                                    >  1:3373540/1‑93 (MQ=255)

                                                      ttACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAg                                   >  2:1612798/1‑93 (MQ=255)

                                                      ttACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAg                                   >  1:1868701/1‑93 (MQ=255)

                                                       tACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGc                                  >  1:3824010/1‑93 (MQ=255)

                                                        aCAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCa                                 >  1:6212712/1‑93 (MQ=255)

                                                        aCAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCa                                 >  1:4944366/1‑93 (MQ=255)

                                                        aCAGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCa                                 <  2:5862044/93‑1 (MQ=255)

                                                          aGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAgg                               >  2:4877324/1‑93 (MQ=255)

                                                          aGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAgg                               >  2:6426172/1‑93 (MQ=255)

                                                          aGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAgg                               <  1:2646919/93‑1 (MQ=255)

                                                          aGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAgg                               <  1:5738363/93‑1 (MQ=255)

                                                          aGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAgg                               >  1:6804146/1‑93 (MQ=255)

                                                          aGGCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAgg                               >  1:7579691/1‑93 (MQ=255)

                                                          aGGCGTTAGAACAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAgg                               <  2:6476871/93‑1 (MQ=255)

                                                            gCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTg                             <  2:2485359/93‑1 (MQ=255)

                                                            gCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTg                             >  2:6696238/1‑93 (MQ=255)

                                                            gCGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTg                             >  2:3151856/1‑93 (MQ=255)

                                                             cGTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTgg                            >  2:403866/1‑92 (MQ=255)

                                                              gTTAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAg                           >  2:375125/1‑93 (MQ=255)

                                                               ttAGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGc                          <  2:5722732/93‑1 (MQ=255)

                                                                 aGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGa                        >  2:381151/1‑93 (MQ=255)

                                                                 aGAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGa                        <  1:7440352/93‑1 (MQ=255)

                                                                  gAACAGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGaa                       >  1:4220668/1‑93 (MQ=255)

                                                                      aGGATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAgg                   >  1:583991/1‑93 (MQ=255)

                                                                       ggATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGt                  <  2:2326863/93‑1 (MQ=255)

                                                                       ggATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGt                  <  2:1830258/93‑1 (MQ=255)

                                                                       ggATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGt                  <  2:1023991/93‑1 (MQ=255)

                                                                        gATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGtt                 <  1:530537/93‑1 (MQ=255)

                                                                        gATGGTTTTCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGtt                 >  2:1201262/1‑93 (MQ=255)

                                                                         aTGGTTTTCTTAACCGTGTCGCGTATCCGGGGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGGtg                >  2:1579795/1‑93 (MQ=255)

                                                                             ttttCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtgctgc            <  1:2839108/93‑1 (MQ=255)

                                                                             ttttCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtgctgc            >  1:2560262/1‑93 (MQ=255)

                                                                             ttttCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTtgctgc            >  2:7720959/1‑93 (MQ=255)

                                                                               ttCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCAc          >  2:2941361/1‑93 (MQ=255)

                                                                                tCTTATCCGTGTCGCGTATCGGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACt         <  2:2804618/93‑1 (MQ=255)

                                                                                tCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACt         >  2:7925105/1‑93 (MQ=255)

                                                                                tCTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACt         >  1:2476640/1‑93 (MQ=255)

                                                                                 cTTAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTg        <  1:6052467/93‑1 (MQ=255)

                                                                                  ttAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTgg       >  1:3524971/1‑93 (MQ=255)

                                                                                  ttAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTgg       >  1:2078731/1‑93 (MQ=255)

                                                                                  ttAACCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTgg       >  1:1319960/1‑93 (MQ=255)

                                                                                    aaCCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCg     >  2:425631/1‑93 (MQ=255)

                                                                                    aaCCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACAGGCg     <  1:3639763/93‑1 (MQ=255)

                                                                                     aCCGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCgg    >  2:7476219/1‑93 (MQ=255)

                                                                                      ccGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGa   <  2:5278896/93‑1 (MQ=255)

                                                                                       cGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAc  >  1:7803841/1‑93 (MQ=255)

                                                                                       cGTGTCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAc  <  1:2279224/93‑1 (MQ=255)

AGCGACCTGCGGCGCAAAATCGGTGGGGTGAGTGAAAAGATGCTTGCGCAGTCGTTACAGGCGTTAGAACAGGATGGTTTTCTTAACCGTATCGCGTATCCGGTGGTGCCGCCGCATGTGGAATATAGCCTCACGCCGCTGGGCGAGCAGGTGAGCGAAAAGGTTGCTGCACTGGCGGAC  >  NC_000913/4431835‑4432014

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