BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,486,394	T→C	A278A (GCA→GCG)	yfdE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,486,394	0	T	C	99.7%	1020.8 / ‑5.1	298	A278A (GCA→GCG)	yfdE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (150/147); total (150/148)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 310 total aligned reads displayed.

CTCAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCAT  >  NC_000913/2486304‑2486483

 tCAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGcc                                                                                        >  2:7103726/1‑93 (MQ=255)

 tCAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGcc                                                                                        >  2:1804466/1‑93 (MQ=255)

 tCAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGcc                                                                                        <  1:1926353/93‑1 (MQ=255)

 tCAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGcc                                                                                        <  1:8237189/93‑1 (MQ=255)

   aaTATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTg                                                                                      >  2:3924697/1‑93 (MQ=255)

    atatatTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTgg                                                                                     >  2:6008838/1‑93 (MQ=255)

      atatTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCa                                                                                   >  1:8580148/1‑93 (MQ=255)

      atatTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCa                                                                                   >  1:7427875/1‑93 (MQ=255)

       tatTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCAt                                                                                  >  1:383477/1‑93 (MQ=255)

       tatTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCAt                                                                                  >  2:4747529/1‑93 (MQ=255)

            tttAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGc                                                                             <  2:2892207/93‑1 (MQ=255)

              tAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAg                                                                           <  2:5748039/93‑1 (MQ=255)

               aaGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGa                                                                          >  2:2659504/1‑93 (MQ=255)

               aaGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGa                                                                          >  1:7087241/1‑93 (MQ=255)

                   aTAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTACCCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGaaaaa                                                                      <  2:4311116/93‑1 (MQ=255)

                   aTAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGaaaaa                                                                      <  2:7811714/93‑1 (MQ=255)

                    tAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGaaaaaa                                                                     <  2:1515879/93‑1 (MQ=255)

                     aGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAg                                                                    <  2:5288229/93‑1 (MQ=255)

                      gCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGc                                                                   >  1:1066333/1‑93 (MQ=255)

                      gCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGc                                                                   >  1:3319665/1‑93 (MQ=255)

                       ccTGGTTTTGTACGCGTAAAATATTGCTGCTAAGTCGGGGATCATTACCCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCt                                                                  <  1:5172152/93‑1 (MQ=255)

                       ccTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCt                                                                  >  2:2309475/1‑93 (MQ=255)

                       ccTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCt                                                                  >  1:5666762/1‑93 (MQ=255)

                       ccTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCt                                                                  <  2:8141052/93‑1 (MQ=255)

                        cTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCtt                                                                 >  2:4356503/1‑93 (MQ=255)

                        cTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCtt                                                                 >  1:2690536/1‑93 (MQ=255)

                        cTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCtt                                                                 >  2:1043066/1‑93 (MQ=255)

                         tGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAACAGCTTg                                                                >  1:6843886/1‑93 (MQ=255)

                         tGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTg                                                                >  2:663760/1‑93 (MQ=255)

                          ggTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAATCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGt                                                               >  1:6689455/1‑93 (MQ=255)

                          ggTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGt                                                               >  2:1107628/1‑93 (MQ=255)

                          ggTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGt                                                               >  1:8613698/1‑93 (MQ=255)

                          ggTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGt                                                               >  2:2139011/1‑93 (MQ=255)

                           gTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTc                                                              >  2:2768290/1‑93 (MQ=255)

                           gTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTc                                                              >  2:6402136/1‑93 (MQ=255)

                           gTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTc                                                              <  2:4183852/93‑1 (MQ=255)

                              ttGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCAtt                                                           >  1:6123227/1‑93 (MQ=255)

                              ttGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAAAGCAGAAAAAAGCTTGTCAtt                                                           <  2:2870797/93‑1 (MQ=255)

                                gTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTAc                                                         <  1:376797/93‑1 (MQ=255)

                                 tACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTAcc                                                        <  2:3046155/93‑1 (MQ=255)

                                  aCGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACca                                                       >  1:6224728/1‑93 (MQ=255)

                                    gcgTAAAATATTGCTGATAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACcaca                                                     >  2:8366121/1‑93 (MQ=255)

                                     cgTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCacaa                                                    >  2:3262987/1‑93 (MQ=255)

                                     cgTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCacaa                                                    >  2:3439497/1‑93 (MQ=255)

                                        aaaaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCAcaacaa                                                 <  1:4389442/93‑1 (MQ=255)

                                         aaaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACaaa                                                >  1:1305599/1‑93 (MQ=255)

                                         aaaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACaaa                                                >  1:2569440/1‑93 (MQ=255)

                                         aaaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACaaa                                                <  1:4741828/93‑1 (MQ=255)

                                          aaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAAt                                               >  2:1201832/1‑93 (MQ=255)

                                          aaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAAt                                               >  1:2826918/1‑93 (MQ=255)

                                          aaTATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAAt                                               >  2:3733011/1‑93 (MQ=255)

                                              ttGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTa                                           >  1:6282352/1‑93 (MQ=255)

                                              ttGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTa                                           >  2:7794140/1‑93 (MQ=255)

                                                gctgctAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAAt                                         <  1:7648601/93‑1 (MQ=255)

                                                gctgctAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAAt                                         <  1:5205818/93‑1 (MQ=255)

                                                 ctgctAAATCGGGGATCATTACCCAGTTCCGTAAGCTCCAGCGCCTGCCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATc                                        <  1:6486799/93‑1 (MQ=255)

                                                 ctgctAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCTCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATc                                        <  2:5508009/93‑1 (MQ=255)

                                                 ctgctAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATc                                        <  2:4057132/93‑1 (MQ=255)

                                                 ctgctAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATc                                        <  2:1644775/93‑1 (MQ=255)

                                                  tgctAAATCGGGGATCATTAACCATTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCg                                       >  1:7458685/1‑93 (MQ=255)

                                                  tgctAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCg                                       >  2:2510776/1‑93 (MQ=255)

                                                       aaTCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTa                                  <  2:7326282/93‑1 (MQ=255)

                                                        aTCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTAt                                 >  2:4869236/1‑93 (MQ=255)

                                                        aTCGGGGATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTAt                                 <  2:3298470/93‑1 (MQ=255)

                                                           ggggATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCt                              >  2:4443929/1‑93 (MQ=255)

                                                            gggATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAACAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTCACCTg                             >  2:7434916/1‑93 (MQ=255)

                                                            gggATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTg                             >  1:697962/1‑93 (MQ=255)

                                                            gggATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTg                             >  1:4576094/1‑93 (MQ=255)

                                                             ggATCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGa                            >  1:8788528/1‑93 (MQ=255)

                                                                tCATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTa                         >  2:4954233/1‑93 (MQ=255)

                                                                 cATTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTAt                        >  1:896841/1‑93 (MQ=255)

                                                                  aTTAACCAGTTCCGTAAGCTCCAGCGCCTGGCATGACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTAtt                       <  2:135753/93‑1 (MQ=255)

                                                                    tAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGa                     <  2:591360/93‑1 (MQ=255)

                                                                    tAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGa                     >  2:5067316/1‑93 (MQ=255)

                                                                    tAACCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGa                     <  1:3670321/93‑1 (MQ=255)

                                                                     aaCCAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGaa                    <  2:7922874/93‑1 (MQ=255)

                                                                       ccAGTTCCGTAAGGTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGaaaa                  <  1:7592659/93‑1 (MQ=255)

                                                                       ccAGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGaaaa                  >  2:4803294/1‑93 (MQ=255)

                                                                         aGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACa                <  2:6952561/93‑1 (MQ=255)

                                                                         aGTTCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACa                >  1:3086233/1‑93 (MQ=255)

                                                                           ttCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATc              <  2:8629823/93‑1 (MQ=255)

                                                                            tCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCa             >  1:4151657/1‑93 (MQ=255)

                                                                            tCCGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCa             <  1:5386836/93‑1 (MQ=255)

                                                                             ccGTAATCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCaa            <  2:625890/93‑1 (MQ=255)

                                                                             ccGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATcga            >  2:7890833/1‑91 (MQ=255)

                                                                             ccGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCaa            <  2:4652734/93‑1 (MQ=255)

                                                                              cGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCaaa           <  2:4541978/93‑1 (MQ=255)

                                                                              cGTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCaaa           >  2:5463488/1‑93 (MQ=255)

                                                                               gTAGGCTCCAACGCGTGGCATAACGGAGAAAAAAGCTTGTCATTACCACAACAAATGGTAGTCGGCGTATCCTGAGTATTGAAAACATCaatg          >  2:8762987/1‑91 (MQ=255)

                                                                               gTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCaaaa          <  1:8510055/93‑1 (MQ=255)

                                                                               gTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCaaaa          <  2:5886746/93‑1 (MQ=255)

                                                                               gTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCaaaa          <  2:2071396/93‑1 (MQ=255)

                                                                               gTAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCaaaa          <  1:6008838/93‑1 (MQ=255)

                                                                                tAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAg         <  1:1735242/93‑1 (MQ=255)

                                                                                tAAGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACAACAACAAATCGTAATCGGTTTATCCTGAGTATTGAAAACATCAACAg         >  1:4488444/1‑93 (MQ=255)

                                                                                 aaGCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAgg        <  1:2456760/93‑1 (MQ=255)

                                                                                   gCTCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTg      <  2:1489782/93‑1 (MQ=255)

                                                                                     tCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGcc    <  2:3824339/93‑1 (MQ=255)

                                                                                     tCCAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGcc    >  2:8516764/1‑93 (MQ=255)

                                                                                       cAGCGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCAt  <  1:5867365/93‑1 (MQ=255)

CTCAATATATTGTTTAAGAATAGCCTGGTTTTGTACGCGTAAAATATTGCTGCTAAATCGGGGATCATTAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCAT  >  NC_000913/2486304‑2486483

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