BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,713,565	G→A	H117Y (CAT→TAT)	pdxY ←	pyridoxal kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,713,565	0	G	A	99.5%	683.1 / ‑6.6	188	H117Y (CAT→TAT)	pdxY	pyridoxal kinase 2
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (91/96); total (92/96)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 192 total aligned reads displayed.

GATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCG  >  NC_000913/1713475‑1713654

  tttGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACTCCCGGTGCAACGATACAGCCTTTTTCCGGATAAccc                                                                                       <  1:3201603/93‑1 (MQ=255)

    tGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCAGGATAACCAAt                                                                                     >  1:4405179/1‑93 (MQ=255)

     ggCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCAtt                                                                                    >  2:2147955/1‑93 (MQ=255)

       cgcAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTAc                                                                                  <  1:3247345/93‑1 (MQ=255)

           aTGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACCCCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGa                                                                              <  2:3476056/93‑1 (MQ=255)

           aTGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCTGa                                                                              >  1:5746474/1‑93 (MQ=255)

           aTGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGa                                                                              <  1:3162449/93‑1 (MQ=255)

           aTGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGa                                                                              <  2:804064/93‑1 (MQ=255)

           aTGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGa                                                                              >  2:3100284/1‑93 (MQ=255)

           aTGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGa                                                                              >  1:5510736/1‑93 (MQ=255)

            tGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAGCTCTGCGACCCCCGGTGCAACGATACAGCCTTTTTCCGGATAGCCCATTACCGGAt                                                                             <  1:4968981/93‑1 (MQ=255)

            tGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGAt                                                                             <  2:2250537/93‑1 (MQ=255)

            tGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGAt                                                                             <  2:4614266/93‑1 (MQ=255)

             gATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATc                                                                            <  1:1851428/93‑1 (MQ=255)

               tatCGCTGGCAGGCAAACCGTGCCGCACATGAAGCTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGc                                                                          <  1:5864102/93‑1 (MQ=255)

                atCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCa                                                                         >  1:3452209/1‑93 (MQ=255)

                  cGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTGTTTCCGGATAACCCATTACCGGATCGCaaa                                                                       <  2:1137498/93‑1 (MQ=255)

                    cTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCaaaaa                                                                     >  2:6146469/1‑93 (MQ=255)

                    cTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCaaaaa                                                                     >  2:5736523/1‑93 (MQ=255)

                     tGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAat                                                                    >  1:6237572/1‑93 (MQ=255)

                     tGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAat                                                                    >  2:5207869/1‑93 (MQ=255)

                          ggcaCACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTc                                                               <  2:4656260/93‑1 (MQ=255)

                          ggcaAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTc                                                               <  2:1764208/93‑1 (MQ=255)

                          ggcaAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTc                                                               <  2:4796776/93‑1 (MQ=255)

                             aaaCCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGcc                                                            <  1:44719/93‑1 (MQ=255)

                              aaCCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCt                                                           >  1:361624/1‑93 (MQ=255)

                              aaCCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCt                                                           <  1:6607137/93‑1 (MQ=255)

                              aaCCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCt                                                           >  1:1509167/1‑93 (MQ=255)

                              aaCCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCt                                                           <  2:4150246/93‑1 (MQ=255)

                               aCCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTg                                                          >  2:2465983/1‑93 (MQ=255)

                               aCCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTg                                                          >  1:423440/1‑93 (MQ=255)

                                  gTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCgg                                                       >  2:3019958/1‑93 (MQ=255)

                                   tGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACTGGATCGCAAAAATATTTCGCCTGCGGa                                                      <  1:4160192/93‑1 (MQ=255)

                                   tGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGa                                                      <  1:6699348/93‑1 (MQ=255)

                                     ccGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGAtt                                                    <  1:4616272/93‑1 (MQ=255)

                                        cacaTGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTcgc                                                 >  1:1198862/1‑93 (MQ=255)

                                         acaTGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTcgcg                                                >  2:701484/1‑93 (MQ=255)

                                          caTGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCgg                                               >  1:4571990/1‑93 (MQ=255)

                                           aTGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGc                                              >  2:442565/1‑93 (MQ=255)

                                           aTGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGc                                              <  1:6361768/93‑1 (MQ=255)

                                           aTGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGc                                              <  1:2923875/93‑1 (MQ=255)

                                           aTGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGc                                              <  1:6647284/93‑1 (MQ=255)

                                             gAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCtt                                            <  1:2778826/93‑1 (MQ=255)

                                                 ctctGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCAc                                        <  2:552808/93‑1 (MQ=255)

                                                  tctGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCAcc                                       >  1:3763731/1‑93 (MQ=255)

                                                    tGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTg                                     <  1:1299747/93‑1 (MQ=255)

                                                      cGACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTgac                                   <  2:5347866/93‑1 (MQ=255)

                                                       gACACCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTgacg                                  >  1:5372012/1‑93 (MQ=255)

                                                         cacCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTgacgga                                <  1:2969172/93‑1 (MQ=255)

                                                         cacCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTgacgga                                <  1:2733397/93‑1 (MQ=255)

                                                          acCCGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTgacggac                               >  1:3593902/1‑93 (MQ=255)

                                                            ccGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGa                             <  2:1506798/93‑1 (MQ=255)

                                                            ccGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGa                             <  1:1682548/93‑1 (MQ=255)

                                                             cGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGAt                            <  2:989285/93‑1 (MQ=255)

                                                             cGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGAt                            <  2:2354014/93‑1 (MQ=255)

                                                             cGGTGCAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGAt                            >  2:3622382/1‑93 (MQ=255)

                                                                  cAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGa                       <  2:3002155/93‑1 (MQ=255)

                                                                  cAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGa                       <  1:1968085/93‑1 (MQ=255)

                                                                  cAACGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGa                       <  2:5139812/93‑1 (MQ=255)

                                                                   aaCGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAg                      >  2:6660638/1‑93 (MQ=255)

                                                                    aCGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAgg                     >  2:327293/1‑93 (MQ=255)

                                                                    aCGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAgg                     >  1:5963313/1‑93 (MQ=255)

                                                                    aCGATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAgg                     >  2:6454978/1‑93 (MQ=255)

                                                                      gATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGat                   <  1:1074061/93‑1 (MQ=255)

                                                                      gATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGat                   >  1:484079/1‑93 (MQ=255)

                                                                      gATACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGat                   >  2:3777695/1‑93 (MQ=255)

                                                                       aTACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGata                  >  2:1884627/1‑93 (MQ=255)

                                                                       aTACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGata                  >  1:145049/1‑93 (MQ=255)

                                                                       aTACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGata                  <  2:2432269/93‑1 (MQ=255)

                                                                       aTACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGata                  <  1:1666826/93‑1 (MQ=255)

                                                                        tACAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGatat                 <  1:2689461/93‑1 (MQ=255)

                                                                          cAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGt               <  1:5103348/93‑1 (MQ=255)

                                                                          cAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGt               <  2:4771587/93‑1 (MQ=255)

                                                                          cAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGt               >  2:3521645/1‑93 (MQ=255)

                                                                          cAGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGt               <  2:1210432/93‑1 (MQ=255)

                                                                           aGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGtt              >  1:964177/1‑93 (MQ=255)

                                                                           aGCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGtt              >  2:2260993/1‑93 (MQ=255)

                                                                            gCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTc             >  2:3392746/1‑93 (MQ=255)

                                                                            gCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTc             >  1:753836/1‑93 (MQ=255)

                                                                            gCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTc             <  1:3388799/93‑1 (MQ=255)

                                                                            gCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTc             >  1:2128393/1‑93 (MQ=255)

                                                                            gCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTc             >  2:4841550/1‑93 (MQ=255)

                                                                            gCCTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTc             <  2:5115616/93‑1 (MQ=255)

                                                                              cTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGTGGCTTTCACCTGACGGACGATAACGAGGATATGTTCAc           >  1:1749900/1‑93 (MQ=255)

                                                                              cTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCAc           >  1:3675670/1‑93 (MQ=255)

                                                                              cTTTTTCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCAc           >  2:2878966/1‑93 (MQ=255)

                                                                               tttttCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCAcc          >  2:1438522/1‑93 (MQ=255)

                                                                               tttttCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCAcc          >  2:5576138/1‑93 (MQ=255)

                                                                                ttttCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCAccc         <  1:1075539/93‑1 (MQ=255)

                                                                                ttttCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCAccc         <  2:128019/93‑1 (MQ=255)

                                                                                ttttCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCAccc         <  1:515585/93‑1 (MQ=255)

                                                                                 tttCCGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCt        >  1:2511094/1‑93 (MQ=255)

                                                                                     cGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTc    <  2:4465031/93‑1 (MQ=255)

                                                                                     cGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTc    <  2:1572011/93‑1 (MQ=255)

                                                                                     cGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTc    >  1:702405/1‑93 (MQ=255)

                                                                                     cGGATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTc    <  2:620976/93‑1 (MQ=255)

                                                                                      ggATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTcc   >  2:5967062/1‑93 (MQ=255)

                                                                                      ggATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTcc   >  1:2019882/1‑93 (MQ=255)

                                                                                       gATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTccg  >  1:5690849/1‑93 (MQ=255)

                                                                                       gATAACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTccg  >  1:1334587/1‑93 (MQ=255)

GATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCG  >  NC_000913/1713475‑1713654

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