BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	673,104	T→C	K301K (AAA→AAG)	leuS ←	leucine tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	673,104	0	T	C	99.7%	940.3 / NA	296	K301K (AAA→AAG)	leuS	leucine tRNA synthetase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (149/146); total (149/147)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 321 total aligned reads displayed.

GCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCTTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTCCGC  >  NC_000913/673014‑673193

gCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCctt                                                                                         <  2:1291380/93‑1 (MQ=255)

gCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCctt                                                                                         >  1:2728751/1‑93 (MQ=255)

  cGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCctttt                                                                                       <  1:926649/93‑1 (MQ=255)

  cGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCctttt                                                                                       <  1:1880749/93‑1 (MQ=255)

  cGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCctttt                                                                                       <  2:6660128/93‑1 (MQ=255)

   gTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCcttttt                                                                                      >  2:3357019/1‑93 (MQ=255)

    tACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTtc                                                                                     >  1:8465333/1‑93 (MQ=255)

    tACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTtc                                                                                     >  2:3986337/1‑93 (MQ=255)

      cTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTtctc                                                                                   >  2:6436573/1‑93 (MQ=255)

        ccATCAATACGAAGTTTGCTGCCCAAGCGGGAATTTCTTCGCCAGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCa                                                                                 <  1:5729488/93‑1 (MQ=255)

         cATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCAt                                                                                >  2:848009/1‑93 (MQ=255)

         cATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCAt                                                                                <  2:2975616/93‑1 (MQ=255)

         cATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCAt                                                                                <  2:5168568/93‑1 (MQ=255)

              atgcGAAGTTTGCCGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCAtcgtcg                                                                           <  2:3785089/90‑1 (MQ=255)

               tACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGc                                                                          >  2:3188904/1‑93 (MQ=255)

                 cGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCa                                                                        >  1:2416478/1‑93 (MQ=255)

                  gAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCAt                                                                       >  1:7219570/1‑93 (MQ=255)

                  gAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCAt                                                                       >  1:6015553/1‑93 (MQ=255)

                  gAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCAt                                                                       >  2:7473230/1‑93 (MQ=255)

                   aaGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCAtt                                                                      >  2:6533785/1‑93 (MQ=255)

                   aaGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCAtt                                                                      >  1:9754816/1‑93 (MQ=255)

                    aGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCAttt                                                                     <  2:2444628/93‑1 (MQ=255)

                    aGTTTGCTGCCCAAACGGGAATTTCCTCGCCCGTTAATTGGTGAACCGCTTTAAAGCCAGTCTCAAC‑GCCCTTTTTCTCCATCGTCGCCAttt                                                                     <  2:3475340/93‑1 (MQ=255)

                      tttGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCa                                                                   <  1:664734/93‑1 (MQ=255)

                      tttGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCa                                                                   >  1:7109143/1‑93 (MQ=255)

                        tgctgcCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGc                                                                 <  1:153529/93‑1 (MQ=255)

                         gctgcCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGCCGCCATTTCAGCt                                                                <  2:7259968/93‑1 (MQ=255)

                         gctgcCCAAAAGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCt                                                                <  1:8591947/93‑1 (MQ=255)

                          ctgcCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCtt                                                               >  2:5277224/1‑93 (MQ=255)

                            gcCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCg                                                             >  2:742927/1‑93 (MQ=255)

                                 aaCGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAAc                                                        >  1:1101163/1‑93 (MQ=255)

                                   cGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACtt                                                      >  2:4621083/1‑93 (MQ=255)

                                    gggAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACttt                                                     >  1:5212884/1‑93 (MQ=255)

                                     ggAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTg                                                    >  2:9598491/1‑93 (MQ=255)

                                     ggAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTg                                                    >  1:6213332/1‑93 (MQ=255)

                                      gAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTgg                                                   <  2:2106345/93‑1 (MQ=255)

                                       aaTTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGgt                                                  >  2:9408735/1‑93 (MQ=255)

                                        aTTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGgtg                                                 >  2:5911315/1‑93 (MQ=255)

                                        aTTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGgtg                                                 >  2:9206877/1‑93 (MQ=255)

                                        aTTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGgtg                                                 >  1:6491065/1‑93 (MQ=255)

                                           tcttcGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTa                                              >  2:9153328/1‑93 (MQ=255)

                                            cttcGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTAc                                             <  2:760975/93‑1 (MQ=255)

                                            cttcGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTAc                                             >  2:1158307/1‑93 (MQ=255)

                                            cttcGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTAc                                             <  1:8555922/93‑1 (MQ=255)

                                             ttcGCCCGTTAATGGGTGAACCGCTTTAAAGCGAGTATCGCC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACg                                            >  1:2172931/1‑93 (MQ=255)

                                             ttcGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACg                                            <  1:3492225/93‑1 (MQ=255)

                                              tcGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACgg                                           >  1:2848651/1‑93 (MQ=255)

                                                gCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCa                                         <  1:5078507/93‑1 (MQ=255)

                                                  ccGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCAtt                                       >  2:4888175/1‑93 (MQ=255)

                                                  ccGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCAtt                                       <  2:1378405/93‑1 (MQ=255)

                                                  ccGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCAtt                                       >  2:15424/1‑93 (MQ=255)

                                                  ccGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCAtt                                       >  1:3251347/1‑93 (MQ=255)

                                                  ccGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTGTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCAtt                                       <  1:6005047/93‑1 (MQ=255)

                                                   cGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTc                                      >  2:4500186/1‑93 (MQ=255)

                                                    gTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCg                                     <  2:4591911/93‑1 (MQ=255)

                                                    gTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCg                                     >  2:5554313/1‑93 (MQ=255)

                                                     ttAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGt                                    >  2:8790661/1‑93 (MQ=255)

                                                     ttAATGGGGGAACCGCTTTAAAGGCAGTTTTGGCGGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCg                                     <  2:361856/93‑1 (MQ=255)

                                                      tAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTc                                   <  1:3555793/93‑1 (MQ=255)

                                                       aaTGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCa                                  >  1:4426041/1‑93 (MQ=255)

                                                           ggTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATaa                              <  1:4104251/93‑1 (MQ=255)

                                                           ggTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATaa                              >  1:3060563/1‑93 (MQ=255)

                                                            gTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATaaa                             >  2:7640993/1‑93 (MQ=255)

                                                            gTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATaaa                             <  1:9906225/93‑1 (MQ=255)

                                                            gTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATaaa                             <  2:8332868/93‑1 (MQ=255)

                                                            gTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATaaa                             <  1:936115/93‑1 (MQ=255)

                                                            gTGAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATaaa                             <  2:4006642/93‑1 (MQ=255)

                                                              gAACCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAgg                           >  2:2853763/1‑93 (MQ=255)

                                                               aaCCGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgc                          <  2:6135240/93‑1 (MQ=255)

                                                                  cGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgccgc                       <  1:3945662/93‑1 (MQ=255)

                                                                  cGCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgccgc                       <  1:9534617/93‑1 (MQ=255)

                                                                   gCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgccgcc                      <  2:6839317/93‑1 (MQ=255)

                                                                   gCTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgccgcc                      <  2:8088943/93‑1 (MQ=255)

                                                                    cTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCa                     >  1:4775168/1‑93 (MQ=255)

                                                                    cTTTAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCa                     >  1:7590251/1‑93 (MQ=255)

                                                                     tttAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGACATTCGTCAATAAAGGCCGCCAg                    >  1:1073203/1‑93 (MQ=255)

                                                                      ttAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGt                   >  2:5043927/1‑93 (MQ=255)

                                                                      ttAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGt                   >  2:6264188/1‑93 (MQ=255)

                                                                      ttAAAGCCAGGATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGt                   <  2:3039291/93‑1 (MQ=255)

                                                                       tAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGtt                  >  1:7471493/1‑93 (MQ=255)

                                                                       tAAAGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGtt                  >  2:1902979/1‑93 (MQ=255)

                                                                        aaaGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTc                 <  1:7342761/93‑1 (MQ=255)

                                                                         aaGCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCa                <  1:1609207/93‑1 (MQ=255)

                                                                           gCCAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAgg              >  2:5528418/1‑93 (MQ=255)

                                                                            ccAGTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGa             <  1:5624213/93‑1 (MQ=255)

                                                                             cAGTATCGAC‑GCCCTTTTTCTCGATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGat            >  1:6947666/1‑93 (MQ=255)

                                                                               gTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGatta          <  2:4308190/93‑1 (MQ=255)

                                                                               gTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGatta          <  2:1818609/93‑1 (MQ=255)

                                                                               gTATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACCTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGatta          >  1:4211186/1‑93 (MQ=255)

                                                                                tATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGattat         <  2:7845537/93‑1 (MQ=255)

                                                                                tATCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGattat         >  2:6019239/1‑93 (MQ=255)

                                                                                 aTCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGattatt        <  2:227326/93‑1 (MQ=255)

                                                                                 aTCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGattatt        <  1:7430809/93‑1 (MQ=255)

                                                                                  tCGAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTAttt       <  1:7683505/93‑1 (MQ=255)

                                                                                  tCGAC‑GCACTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTAttt       >  1:8033607/1‑93 (MQ=255)

                                                                                    gAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTc     >  2:2029222/1‑93 (MQ=255)

                                                                                    gAC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTc     <  1:4661951/93‑1 (MQ=255)

                                                                                     aC‑GCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTcc    <  1:6360312/93‑1 (MQ=255)

                                                                                        gCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTccgc  <  1:1470425/93‑1 (MQ=255)

                                                                                        gCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTccgc  <  2:8563400/93‑1 (MQ=255)

GCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGAC‑GCCTTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTCCGC  >  NC_000913/673014‑673193

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