BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,043,403	A→G	H260H (CAT→CAC)	yihG ←	putative endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,043,403	0	A	G	99.5%	1190.9 / ‑4.4	375	H260H (CAT→CAC)	yihG	putative endonuclease
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (203/169); total (206/169)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 395 total aligned reads displayed.

TGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAGTT  >  NC_000913/4043313‑4043494

tGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTg                                                                                           >  1:3724031/1‑93 (MQ=255)

tGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTg                                                                                           <  2:7606609/93‑1 (MQ=255)

 gAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACgtgt                                                                                          <  2:2173458/93‑1 (MQ=255)

 gAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACgtgt                                                                                          >  2:3330202/1‑93 (MQ=255)

 gAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACgtgt                                                                                          <  1:7312497/93‑1 (MQ=255)

  aGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTa                                                                                         >  2:3961502/1‑93 (MQ=255)

  aGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTa                                                                                         <  1:233132/93‑1 (MQ=255)

   gTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTAc                                                                                        <  1:2898963/93‑1 (MQ=255)

      cAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGAc                                                                                     <  1:770288/93‑1 (MQ=255)

       aGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACa                                                                                    <  1:2489083/93‑1 (MQ=255)

         ccATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAAt                                                                                  <  2:7377446/93‑1 (MQ=255)

         ccATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAAt                                                                                  <  2:4602484/93‑1 (MQ=255)

          cATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATc                                                                                 >  1:7849077/1‑93 (MQ=255)

           aTTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATcc                                                                                <  1:4881316/93‑1 (MQ=255)

            ttGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATcag                                                                               >  2:2687386/1‑91 (MQ=255)

            ttGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCcg                                                                               >  2:6602698/1‑93 (MQ=255)

            ttGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCcg                                                                               >  2:6787675/1‑93 (MQ=255)

              gCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCcgcg                                                                             >  1:7099439/1‑93 (MQ=255)

              gCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCcgcg                                                                             <  2:5143328/93‑1 (MQ=255)

                tGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTc                                                                           >  1:3614385/1‑93 (MQ=255)

                tGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTc                                                                           >  1:6814121/1‑93 (MQ=255)

                 gAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCa                                                                          >  2:5735712/1‑93 (MQ=255)

                  aaaaTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAg                                                                         >  2:3544024/1‑93 (MQ=255)

                  aaaaTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAg                                                                         >  2:5389635/1‑93 (MQ=255)

                       ggCGCTTAAAGCTTGTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTAc                                                                    <  2:7915950/93‑1 (MQ=255)

                          gcTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGc                                                                 >  1:5410180/1‑93 (MQ=255)

                          gcTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGc                                                                 >  2:723388/1‑93 (MQ=255)

                           cTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCt                                                                >  1:5609141/1‑93 (MQ=255)

                           cTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCt                                                                >  2:4099512/1‑93 (MQ=255)

                            ttAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGCCAATCCGCGTCAGTTTACCGCtt                                                               >  2:2017111/1‑93 (MQ=255)

                            ttAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCtt                                                               <  1:6892901/93‑1 (MQ=255)

                            ttAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCtt                                                               >  1:2151180/1‑93 (MQ=255)

                            ttAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCtt                                                               >  1:7926004/1‑93 (MQ=255)

                             tAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTa                                                              <  2:1702530/93‑1 (MQ=255)

                             tAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTa                                                              >  1:1594956/1‑93 (MQ=255)

                             tAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTa                                                              >  2:4471396/1‑93 (MQ=255)

                                aGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACa                                                           <  2:104469/93‑1 (MQ=255)

                                aGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACa                                                           >  1:3809304/1‑93 (MQ=255)

                                  cTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACata                                                         >  2:2601583/1‑93 (MQ=255)

                                  cTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACata                                                         <  1:1428083/93‑1 (MQ=255)

                                  cTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACata                                                         <  1:1968261/93‑1 (MQ=255)

                                  cTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATACGCTGTAAATCCACGTGTCCGACAATCCGCGTCAGTTTACCGCTTAACata                                                         >  1:5917512/1‑93 (MQ=255)

                                   ttttGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACatat                                                        >  1:7402307/1‑93 (MQ=255)

                                    tttGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATc                                                       <  2:1861077/93‑1 (MQ=255)

                                    tttGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATc                                                       <  2:3272003/93‑1 (MQ=255)

                                    tttGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATc                                                       >  1:3854496/1‑93 (MQ=255)

                                      tGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCTACGTGTACGACAATCCGGGTAAGTTTACCGCTAAACCAATCga                                                     >  1:6318605/1‑93 (MQ=255)

                                      tGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCga                                                     <  1:7034301/93‑1 (MQ=255)

                                          aTTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaagaa                                                 >  1:1293677/1‑93 (MQ=255)

                                          aTTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaagaa                                                 >  2:682809/1‑93 (MQ=255)

                                           ttGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAAc                                                >  1:3098465/1‑93 (MQ=255)

                                            tGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACg                                               <  1:4483854/93‑1 (MQ=255)

                                             gATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACgg                                              <  1:7253949/93‑1 (MQ=255)

                                             gATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACgg                                              >  2:4848342/1‑93 (MQ=255)

                                              aTGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGc                                             <  1:1771576/93‑1 (MQ=255)

                                              aTGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGc                                             <  2:6917802/93‑1 (MQ=255)

                                                 tAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGa                                          <  2:914810/93‑1 (MQ=255)

                                                  aGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGAc                                         <  1:4582483/93‑1 (MQ=255)

                                                      gCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGtt                                     <  2:2364335/93‑1 (MQ=255)

                                                      gCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGtt                                     >  1:323228/1‑93 (MQ=255)

                                                         aTGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTAtt                                  >  1:5565352/1‑93 (MQ=255)

                                                         aTGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTAtt                                  >  1:7755424/1‑93 (MQ=255)

                                                          tGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTAttt                                 >  2:6113597/1‑92 (MQ=255)

                                                          tGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTg                                 <  1:6555482/93‑1 (MQ=255)

                                                          tGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTg                                 >  1:1121867/1‑93 (MQ=255)

                                                          tGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTg                                 >  1:5795193/1‑93 (MQ=255)

                                                          tGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTg                                 >  2:50692/1‑93 (MQ=255)

                                                          tGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGCTATTg                                 >  2:3664/1‑93 (MQ=255)

                                                            tAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTc                               <  1:5251777/93‑1 (MQ=255)

                                                             aaCTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTcc                              <  1:3011600/93‑1 (MQ=255)

                                                                tCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGa                           <  1:6715172/93‑1 (MQ=255)

                                                                tCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGa                           <  2:7756695/93‑1 (MQ=255)

                                                                  aTCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATa                         >  1:7621914/1‑93 (MQ=255)

                                                                   tCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGa                                          >  1:239020/1‑75 (MQ=255)

                                                                   tCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATaa                        >  2:5478025/1‑93 (MQ=255)

                                                                   tCGGCAATAGGCTGTAAATCCACCTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATaa                        >  1:1004406/1‑93 (MQ=255)

                                                                    cGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAac                       <  2:582759/93‑1 (MQ=255)

                                                                     ggCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAaca                      >  1:4370739/1‑93 (MQ=255)

                                                                       cAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAAcaca                    <  2:4386467/93‑1 (MQ=255)

                                                                       cAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAAcaca                    <  1:2737148/93‑1 (MQ=255)

                                                                       cAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAAcaca                    >  2:3748364/1‑93 (MQ=255)

                                                                        aaTAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAg                   >  2:4339435/1‑93 (MQ=255)

                                                                         aTAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGc                  <  2:5303602/93‑1 (MQ=255)

                                                                         aTAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGc                  >  2:3188557/1‑93 (MQ=255)

                                                                          tAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCg                 >  1:5666615/1‑93 (MQ=255)

                                                                          tAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCg                 >  2:626428/1‑93 (MQ=255)

                                                                            ggCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTa               >  1:3328758/1‑93 (MQ=255)

                                                                            ggCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTa               >  2:142448/1‑93 (MQ=255)

                                                                            ggCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTa               <  1:6153923/93‑1 (MQ=255)

                                                                               tGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACa            >  2:6483768/1‑93 (MQ=255)

                                                                               tGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACa            <  1:426383/93‑1 (MQ=255)

                                                                               tGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACa            >  2:1567026/1‑93 (MQ=255)

                                                                                 tAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACAtt          >  1:4776687/1‑93 (MQ=255)

                                                                                  aaaTCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTa         >  1:4732420/1‑93 (MQ=255)

                                                                                     tCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGc      >  1:3938443/1‑93 (MQ=255)

                                                                                     tCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGc      <  2:3870958/93‑1 (MQ=255)

                                                                                     tCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGc      >  2:2825556/1‑93 (MQ=255)

                                                                                       cacgTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAg    <  2:4798559/93‑1 (MQ=255)

                                                                                       cacgTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAg    <  2:421411/93‑1 (MQ=255)

                                                                                         caTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAGtt  <  2:4979115/93‑1 (MQ=255)

TGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAGTT  >  NC_000913/4043313‑4043494

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