BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,858,325	A→G	D344D (GAT→GAC)	b1776 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,858,325	0	A	G	98.7%	507.5 / ‑3.4	159	D344D (GAT→GAC)	b1776	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (71/86); total (73/86)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 166 total aligned reads displayed.

GTTTTGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTATCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGCCGCTGGCTAA  >  NC_000913/1858235‑1858414

gTTTTGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTc                                                                                         <  2:1768972/93‑1 (MQ=255)

gTTTTGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTc                                                                                         <  1:781197/93‑1 (MQ=255)

gTTTTGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTc                                                                                         <  1:3508808/93‑1 (MQ=255)

  tttGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGa                                                                                       <  2:1412496/93‑1 (MQ=255)

   ttGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGAc                                                                                      >  2:1503523/1‑93 (MQ=255)

         tCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGc                                                                                >  2:1253344/1‑93 (MQ=255)

          cTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGct                                                                               >  2:1653594/1‑92 (MQ=255)

           tGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCAt                                                                              >  2:3477292/1‑93 (MQ=255)

           tGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCAt                                                                              >  2:4047257/1‑93 (MQ=255)

           tGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCAt                                                                              >  1:1956076/1‑93 (MQ=255)

             ttCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCa                                                                            >  2:380439/1‑93 (MQ=255)

               ccATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCaaa                                                                          >  2:3018004/1‑93 (MQ=255)

               ccATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCaaa                                                                          >  2:3018003/1‑93 (MQ=255)

                cATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAAc                                                                         >  1:3527785/1‑93 (MQ=255)

                cATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAAc                                                                         >  1:1690253/1‑93 (MQ=255)

                cATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAAc                                                                         <  2:2552251/93‑1 (MQ=255)

                 aTGCGATAACGTTCCTTCAGGAGACTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAAcc                                                                        <  1:3381446/93‑1 (MQ=255)

                  tGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAAccc                                                                       <  2:4592778/93‑1 (MQ=255)

                   tcGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCt                                                                      <  2:1375943/92‑1 (MQ=255)

                    cGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCtt                                                                     <  1:561645/93‑1 (MQ=255)

                    cGATAACGGTTCTTCAGGAGAGTGGTTATTTAATAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCtt                                                                     <  2:3166040/93‑1 (MQ=255)

                     gATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTc                                                                    <  2:3496681/93‑1 (MQ=255)

                     gATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTc                                                                    <  2:3681355/93‑1 (MQ=255)

                       tAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTcgc                                                                  >  1:4586588/1‑93 (MQ=255)

                        aaCGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTcgcg                                                                 >  1:2850278/1‑93 (MQ=255)

                           gTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCa                                                              >  2:2837985/1‑93 (MQ=255)

                            ttCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCAt                                                             <  1:3989396/93‑1 (MQ=255)

                            ttCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCAt                                                             <  1:1802236/93‑1 (MQ=255)

                            ttCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCAt                                                             <  1:3660636/93‑1 (MQ=255)

                             tCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCAtt                                                            >  1:3986766/1‑93 (MQ=255)

                               cTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGc                                                          >  2:4094639/1‑93 (MQ=255)

                                 tCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGa                                                        <  1:1433938/93‑1 (MQ=255)

                                 tCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGa                                                        >  1:1402844/1‑93 (MQ=255)

                                  cAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGAc                                                       <  1:3675349/93‑1 (MQ=255)

                                  cAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGAc                                                       <  1:4594595/93‑1 (MQ=255)

                                  cAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGAc                                                       <  1:3324373/93‑1 (MQ=255)

                                   aGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACa                                                      <  2:2484352/93‑1 (MQ=255)

                                   aGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACa                                                      <  2:1082813/93‑1 (MQ=255)

                                   aGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACa                                                      >  1:4415373/1‑93 (MQ=255)

                                    ggAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTTCGGTTTTATCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAg                                                     >  2:4108843/1‑93 (MQ=255)

                                    ggAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAg                                                     >  1:4051809/1‑93 (MQ=255)

                                     gagaGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAgg                                                    <  1:334943/93‑1 (MQ=255)

                                      ggagTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGc                                                   <  1:1186207/92‑1 (MQ=255)

                                      agagTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGc                                                   <  2:4433701/93‑1 (MQ=255)

                                         gTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGa                                                >  1:3380597/1‑93 (MQ=255)

                                         gTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGa                                                >  2:3594499/1‑93 (MQ=255)

                                          tGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGAt                                               <  2:3105372/93‑1 (MQ=255)

                                          tGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGAt                                               <  1:890086/93‑1 (MQ=255)

                                          tGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGAt                                               <  1:4670021/93‑1 (MQ=255)

                                          tGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGAt                                               <  1:505966/93‑1 (MQ=255)

                                             ttATTCATCAAAGTCGTAAGTCATGATAACTTTGATTGCGGTTTTGTCAACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACg                                            >  1:4336407/1‑93 (MQ=255)

                                                ttCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATg                                         >  2:1999048/1‑93 (MQ=255)

                                                 tcatcaAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGc                                        <  2:3283525/93‑1 (MQ=255)

                                                  catcaAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCg                                       <  1:1908954/93‑1 (MQ=255)

                                                   atcaAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGt                                      >  2:2112595/1‑93 (MQ=255)

                                                   atcaAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGt                                      <  2:4313314/93‑1 (MQ=255)

                                                   atcaAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGt                                      >  2:2028284/1‑93 (MQ=255)

                                                   atcaAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGt                                      >  1:3859944/1‑93 (MQ=255)

                                                    tgaaaGGCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTg                                     <  1:949047/91‑1 (MQ=255)

                                                     caAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGa                                    <  2:4162535/93‑1 (MQ=255)

                                                     caAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGa                                    >  1:2459603/1‑93 (MQ=255)

                                                     caAAATCGTAAGTCAAGATCTCTTTGATTGCGGTTTTGGCGACCATCGCATCACACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGa                                    <  2:1672347/93‑1 (MQ=255)

                                                      aaaGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGat                                   >  2:4067957/1‑93 (MQ=255)

                                                       aaGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTCCGACAGGCCGATACGATGCGTGatc                                  <  2:1607377/93‑1 (MQ=255)

                                                         gTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGatcat                                <  1:3946714/93‑1 (MQ=255)

                                                          tCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGatcatc                               <  1:1062342/93‑1 (MQ=255)

                                                          tCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGatcatc                               <  1:1595923/93‑1 (MQ=255)

                                                           cGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCg                              <  1:3923178/93‑1 (MQ=255)

                                                           cGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCg                              <  1:3923179/93‑1 (MQ=255)

                                                               aGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGttt                          >  1:2177930/1‑93 (MQ=255)

                                                                 tCATGATCACTTTGATTGCGGTTTTGTCGCCCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGa                        <  1:2139355/93‑1 (MQ=255)

                                                                 tCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGa                        <  1:2200884/93‑1 (MQ=255)

                                                                 tCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCATTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGa                        >  1:1770964/1‑93 (MQ=255)

                                                                  cATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGAc                       >  2:2007450/1‑93 (MQ=255)

                                                                   aTGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACt                      >  2:72061/1‑93 (MQ=255)

                                                                    tGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACtt                     <  2:4267762/93‑1 (MQ=255)

                                                                    tGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACtt                     <  1:4620591/93‑1 (MQ=255)

                                                                    tGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACtt                     >  1:1049984/1‑93 (MQ=255)

                                                                    tGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACtt                     <  1:2546223/93‑1 (MQ=255)

                                                                    agATCACTGTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACtt                     <  1:3359008/92‑1 (MQ=255)

                                                                     gATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACttt                    <  2:615790/93‑1 (MQ=255)

                                                                     gATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACttt                    <  1:2165507/93‑1 (MQ=255)

                                                                       tCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGGGTGATCATCGGTTTGACTTTGa                  >  1:365914/1‑93 (MQ=255)

                                                                       tCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGa                  <  2:851951/93‑1 (MQ=255)

                                                                        cACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGAt                 <  1:574665/93‑1 (MQ=255)

                                                                         ccTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATa                <  1:3098293/92‑1 (MQ=255)

                                                                         aCTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATa                >  2:2647361/1‑93 (MQ=255)

                                                                          cTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAg               >  2:3789340/1‑93 (MQ=255)

                                                                          cTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAg               <  1:843601/93‑1 (MQ=255)

                                                                           tttGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAgc              <  1:2673297/93‑1 (MQ=255)

                                                                           tttGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAgc              >  2:390655/1‑93 (MQ=255)

                                                                           tttGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAgc              <  2:2389827/93‑1 (MQ=255)

                                                                             tGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAgcgc            >  1:4530826/1‑93 (MQ=255)

                                                                               aTTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGcgccg          >  1:2857367/1‑93 (MQ=255)

                                                                               aTTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGcgccg          >  1:1490222/1‑93 (MQ=255)

                                                                               aTTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGcgccg          >  2:1753968/1‑93 (MQ=255)

                                                                                ttGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGcgccgc         >  2:2686156/1‑93 (MQ=255)

                                                                                ttGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCAGATACGATGCGTGATCATCGGTTTGACTTTGATAGcgccgc         >  1:1236145/1‑93 (MQ=255)

                                                                                       tttGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGCCGCTGGCTaa  >  1:3559340/1‑93 (MQ=255)

                                                                                       tttGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGCCGCTGGCTaa  >  1:425912/1‑93 (MQ=255)

GTTTTGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTATCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGCCGCTGGCTAA  >  NC_000913/1858235‑1858414

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