BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,725,180	G→A	I70I (ATC→ATT)	yiaB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,725,180	0	G	A	99.6%	919.8 / ‑5.7	252	I70I (ATC→ATT)	yiaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (105/146); total (105/147)
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.09e-01

Only 100 of 256 total aligned reads displayed.

GGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAG  >  NC_000913/3725090‑3725269

gggATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAat                                                                                         <  1:855966/93‑1 (MQ=255)

gggATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAat                                                                                         >  2:4391581/1‑93 (MQ=255)

 ggATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAata                                                                                        >  2:4339281/1‑93 (MQ=255)

 ggATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAata                                                                                        <  1:2502070/93‑1 (MQ=255)

 ggATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAata                                                                                        <  1:2295780/93‑1 (MQ=255)

  gATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAatat                                                                                       <  1:647953/93‑1 (MQ=255)

  gATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAatat                                                                                       <  1:1555067/93‑1 (MQ=255)

     aGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATggg                                                                                    <  1:175999/93‑1 (MQ=255)

        tCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaa                                                                                 >  1:3588483/1‑93 (MQ=255)

            cATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaa                                                                             >  1:743824/1‑93 (MQ=255)

             aTTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaaa                                                                            >  2:733769/1‑93 (MQ=255)

             aTTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaaa                                                                            <  1:299049/93‑1 (MQ=255)

              tttCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAat                                                                           <  1:4211588/93‑1 (MQ=255)

              tttCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAat                                                                           <  1:3990481/93‑1 (MQ=255)

                tCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatca                                                                         <  2:2474141/93‑1 (MQ=255)

                 cATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatcat                                                                        <  2:419139/93‑1 (MQ=255)

                  aTAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatcatc                                                                       <  1:4285032/93‑1 (MQ=255)

                  aTAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatcatc                                                                       <  1:3626075/93‑1 (MQ=255)

                  aTAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatcatc                                                                       >  1:2525795/1‑93 (MQ=255)

                   tAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                      >  2:3817762/1‑93 (MQ=255)

                   tAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                      <  1:3781080/93‑1 (MQ=255)

                    aaGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCa                                                                     <  2:1064911/93‑1 (MQ=255)

                    aaGTATTCACAGTTGCGTTAACAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACGATCAAAAATGGGTAAAAAAATCATCCa                                                                     <  2:2522858/93‑1 (MQ=255)

                     aGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAg                                                                    >  1:593114/1‑93 (MQ=255)

                      gTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGt                                                                   <  2:2764710/93‑1 (MQ=255)

                        aTTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTg                                                                 >  2:3277689/1‑93 (MQ=255)

                         ttGATAGGTGCGTTTAAAACACCGTTAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGc                                                                <  2:4598097/93‑1 (MQ=255)

                         ttGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGc                                                                <  2:3224327/93‑1 (MQ=255)

                         ttGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGc                                                                >  2:3196394/1‑93 (MQ=255)

                           gATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGccc                                                              >  1:2708094/1‑93 (MQ=255)

                            aTAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCg                                                             >  2:1153032/1‑93 (MQ=255)

                            aTAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCg                                                             <  1:1346223/93‑1 (MQ=255)

                             tAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCgg                                                            >  2:1458164/1‑92 (MQ=255)

                             tAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGa                                                            <  1:1871968/93‑1 (MQ=255)

                               ggTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGt                                                          <  1:34074/93‑1 (MQ=255)

                               ggTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGt                                                          >  2:4140648/1‑93 (MQ=255)

                               ggTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGt                                                          >  2:4636419/1‑93 (MQ=255)

                                   cGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGttttt                                                      >  1:3395946/1‑93 (MQ=255)

                                    gTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTc                                                     <  2:1819136/93‑1 (MQ=255)

                                    gTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTc                                                     <  1:3920819/93‑1 (MQ=255)

                                    gTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTc                                                     <  1:4406705/93‑1 (MQ=255)

                                     tttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCt                                                    <  1:470713/93‑1 (MQ=255)

                                     tttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCt                                                    >  1:3026110/1‑93 (MQ=255)

                                     tttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCt                                                    <  1:2789542/93‑1 (MQ=255)

                                      ttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTg                                                   <  2:1128074/93‑1 (MQ=255)

                                       tAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTgc                                                  <  2:3413500/93‑1 (MQ=255)

                                         aaaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTgcgc                                                <  1:1090459/93‑1 (MQ=255)

                                         aaaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTgcgc                                                >  1:1984329/1‑93 (MQ=255)

                                           acacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGa                                              >  1:3252621/1‑93 (MQ=255)

                                           acacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGa                                              >  2:1625946/1‑93 (MQ=255)

                                            cacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGaa                                             <  1:4136764/93‑1 (MQ=255)

                                            cacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGaa                                             <  2:3923526/93‑1 (MQ=255)

                                             acCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAg                                            >  2:3121963/1‑93 (MQ=255)

                                             acCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAg                                            >  1:3206697/1‑93 (MQ=255)

                                             acCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAg                                            <  2:2309756/93‑1 (MQ=255)

                                              ccGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGa                                           <  2:4094639/93‑1 (MQ=255)

                                               cGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGat                                          >  1:2189485/1‑93 (MQ=255)

                                               cGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGat                                          <  1:2604490/93‑1 (MQ=255)

                                                gATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGata                                         >  2:4292601/1‑93 (MQ=255)

                                                gATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGata                                         <  2:191902/93‑1 (MQ=255)

                                                gATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGata                                         <  2:1758120/93‑1 (MQ=255)

                                                gATAAACAAGAGACCGATAGTGAATAACGCAACGAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGata                                         <  1:39653/93‑1 (MQ=255)

                                                    aaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCa                                     >  1:4657024/1‑93 (MQ=255)

                                                    aaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCa                                     <  2:2994509/93‑1 (MQ=255)

                                                     aCAAGAGCGCGATAGTGATGAACGCAAGTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCAt                                    <  2:1992910/93‑1 (MQ=255)

                                                     aCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCAt                                    <  1:3842845/93‑1 (MQ=255)

                                                     aCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCAt                                    >  1:1323037/1‑93 (MQ=255)

                                                      cAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATa                                   <  1:1392316/93‑1 (MQ=255)

                                                      cAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATa                                   <  1:408118/93‑1 (MQ=255)

                                                      cAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATa                                   <  2:1020890/93‑1 (MQ=255)

                                                      cAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATa                                   <  2:3729383/93‑1 (MQ=255)

                                                       aaGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAg                                  >  2:4233878/1‑93 (MQ=255)

                                                       aaGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAg                                  <  2:754709/93‑1 (MQ=255)

                                                        agagACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGt                                 >  1:901341/1‑93 (MQ=255)

                                                         gagaCCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGtt                                >  2:1805300/1‑93 (MQ=255)

                                                         gagaCCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGtt                                <  1:208188/93‑1 (MQ=255)

                                                         gagaCCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGtt                                <  1:345558/93‑1 (MQ=255)

                                                           gaCCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGc                              >  2:2216719/1‑93 (MQ=255)

                                                                aTAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCaaaa                         <  1:1521090/93‑1 (MQ=255)

                                                                aTAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCaaaa                         >  2:881441/1‑93 (MQ=255)

                                                                 tAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAg                        <  1:2698612/93‑1 (MQ=255)

                                                                  aGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGt                       <  2:348309/93‑1 (MQ=255)

                                                                  aGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGt                       <  1:933507/93‑1 (MQ=255)

                                                                  aGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGt                       <  2:2740840/93‑1 (MQ=255)

                                                                     gATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGc                    >  2:2245971/1‑93 (MQ=255)

                                                                     gATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGc                    <  1:1424009/93‑1 (MQ=255)

                                                                       ttAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTg                  <  1:3655225/93‑1 (MQ=255)

                                                                        tAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGt                 <  1:4112428/93‑1 (MQ=255)

                                                                        tAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGt                 <  2:3645958/93‑1 (MQ=255)

                                                                         aaCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTc                <  2:1354721/93‑1 (MQ=255)

                                                                          aCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCa               <  2:2914598/93‑1 (MQ=255)

                                                                          aCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCa               >  1:3433387/1‑93 (MQ=255)

                                                                            gCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAtt             <  2:2653407/93‑1 (MQ=255)

                                                                            gCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAtt             <  1:1007242/93‑1 (MQ=255)

                                                                             cAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTa            <  1:2439602/93‑1 (MQ=255)

                                                                              aactaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTaa           >  2:1241940/1‑93 (MQ=255)

                                                                                ctaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAaca         <  1:2369814/93‑1 (MQ=255)

                                                                                    ctGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGcc     >  1:96557/1‑93 (MQ=255)

                                                                                    ctGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGcc     <  1:1641147/93‑1 (MQ=255)

                                                                                      ggCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCaa   <  2:3574669/93‑1 (MQ=255)

GGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAG  >  NC_000913/3725090‑3725269

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