BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,212,282	T→C	S392S (AGT→AGC)	dnaG →	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,212,282	0	T	C	99.7%	1549.6 / ‑4.9	365	S392S (AGT→AGC)	dnaG	DNA primase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (177/187); minor base A (0/1); total (177/188)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01

Only 100 of 365 total aligned reads displayed.

GCGCGGATGGAGCAGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACA  >  NC_000913/3212184‑3212379

 cgcgGATGGAGCAGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACg                                                                                                >  2:5968023/1‑101 (MQ=255)

  gcgGATGGAGCAGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGc                                                                                               <  1:5967468/101‑1 (MQ=255)

    ggATGGAGCAGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTg                                                                                             >  1:5968025/1‑101 (MQ=255)

      aTGGAGCAGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGc                                                                                           >  1:5968030/1‑101 (MQ=255)

        ggAGCAGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCAc                                                                                         >  2:5968034/1‑101 (MQ=255)

        ggAGCAGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCAc                                                                                         <  2:5967334/101‑1 (MQ=255)

         gAGCAGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACt                                                                                        <  2:5967504/101‑1 (MQ=255)

          aGCAGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTa                                                                                       <  2:5967705/101‑1 (MQ=255)

           gCAGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTAc                                                                                      >  1:5968037/1‑101 (MQ=255)

             aGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCa                                                                                    >  2:5968042/1‑101 (MQ=255)

              ggCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCAt                                                                                   <  2:5967541/101‑1 (MQ=255)

              ggCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCAt                                                                                   >  1:5968044/1‑101 (MQ=255)

                  aTGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGat                                                                               <  2:5967417/101‑1 (MQ=255)

                  aTGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGat                                                                               >  2:5968048/1‑101 (MQ=255)

                     ccACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATc                                                                            <  2:5967567/101‑1 (MQ=255)

                     ccACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATc                                                                            >  1:5968053/1‑101 (MQ=255)

                     ccACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATc                                                                            <  1:5967463/101‑1 (MQ=255)

                      cACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCg                                                                           >  2:5968057/1‑101 (MQ=255)

                       aCTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGc                                                                          >  1:5968060/1‑101 (MQ=255)

                        ctctCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCa                                                                         <  1:5967067/101‑1 (MQ=255)

                         tctcCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCaa                                                                        >  1:5968066/1‑101 (MQ=255)

                          ctcCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGAACGCTGGCACTACCATTGATATCGCAAg                                                                       <  1:5967694/101‑1 (MQ=255)

                           tcCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGt                                                                      >  2:5968067/1‑101 (MQ=255)

                            ccGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTg                                                                     <  2:5967236/101‑1 (MQ=255)

                            ccGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTg                                                                     <  1:5967556/101‑1 (MQ=255)

                             cGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGc                                                                    <  1:5967579/101‑1 (MQ=255)

                             cGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGc                                                                    >  2:5968072/1‑101 (MQ=255)

                                aTTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCgg                                                                 >  1:5968076/1‑101 (MQ=255)

                                 tttCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCggg                                                                <  2:5967675/101‑1 (MQ=255)

                                 tttCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCggg                                                                <  2:5967863/101‑1 (MQ=255)

                                 tttCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCggg                                                                >  2:5968079/1‑101 (MQ=255)

                                   tCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCg                                                              <  2:5967596/101‑1 (MQ=255)

                                   tCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCg                                                              <  2:5967529/101‑1 (MQ=255)

                                   tCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCg                                                              <  1:5967986/101‑1 (MQ=255)

                                    cTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGa                                                             >  1:5968085/1‑101 (MQ=255)

                                    cTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGa                                                             >  2:5968086/1‑101 (MQ=255)

                                      gTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGaaa                                                           >  2:5968089/1‑101 (MQ=255)

                                          aaCAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTg                                                                                             >  2:5968095/1‑63 (MQ=255)

                                          aaCAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCt                                                       <  2:5967303/101‑1 (MQ=255)

                                           aCAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTg                                                      >  1:5968098/1‑101 (MQ=255)

                                           aCAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTg                                                      <  1:5967381/101‑1 (MQ=255)

                                             aGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCg                                                    >  2:5968102/1‑101 (MQ=255)

                                              gTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGa                                                   >  1:5968106/1‑101 (MQ=255)

                                                cTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAat                                                 <  1:5967482/101‑1 (MQ=255)

                                                  gATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAatat                                               <  2:5967472/101‑1 (MQ=255)

                                                  gATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCCGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAatat                                               >  1:5968111/1‑101 (MQ=255)

                                                    tGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTATCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAatatat                                             <  1:5967563/101‑1 (MQ=255)

                                                    tGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAatatat                                             >  2:5968113/1‑101 (MQ=255)

                                                      ccGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCt                                           >  1:5968116/1‑101 (MQ=255)

                                                       cGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCtt                                          >  2:5968119/1‑101 (MQ=255)

                                                        gCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTc                                         >  2:5968124/1‑101 (MQ=255)

                                                         cAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCg                                        <  2:5967811/101‑1 (MQ=255)

                                                          aaGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGt                                       >  2:5968129/1‑101 (MQ=255)

                                                          aaGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGt                                       <  1:5967614/101‑1 (MQ=255)

                                                           aGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTc                                      <  1:5967314/101‑1 (MQ=255)

                                                           aGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTc                                      <  1:5967805/101‑1 (MQ=255)

                                                           aGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTc                                      >  1:5968130/1‑101 (MQ=255)

                                                             ttGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAg                                    <  2:5967494/101‑1 (MQ=255)

                                                             ttGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAg                                    <  1:5967641/101‑1 (MQ=255)

                                                             ttGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAg                                    <  1:5967633/101‑1 (MQ=255)

                                                             ttGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAg                                    >  2:5968136/1‑101 (MQ=255)

                                                              tGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAgg                                   <  1:5967660/101‑1 (MQ=255)

                                                              tGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAgg                                   <  1:5967599/101‑1 (MQ=255)

                                                               gATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGa                                  >  1:5968144/1‑101 (MQ=255)

                                                               gATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGa                                  >  1:5968141/1‑101 (MQ=255)

                                                                aTCTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGaa                                 <  1:5967754/101‑1 (MQ=255)

                                                                  cTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAAtt                               <  2:5967429/101‑1 (MQ=255)

                                                                  cTGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAAtt                               >  1:5968148/1‑101 (MQ=255)

                                                                   tGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTa                              <  2:5967555/101‑1 (MQ=255)

                                                                   tGAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTa                              <  2:5967449/101‑1 (MQ=255)

                                                                    gAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAg                             <  1:5967540/101‑1 (MQ=255)

                                                                    gAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAg                             <  1:5967514/101‑1 (MQ=255)

                                                                    gAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAg                             <  2:5967516/101‑1 (MQ=255)

                                                                    gAGTACCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAg                             >  2:5968150/1‑101 (MQ=255)

                                                                        aCCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGcaa                         <  1:5967257/101‑1 (MQ=255)

                                                                        aCCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGcaa                         <  2:5967631/101‑1 (MQ=255)

                                                                        aCCCCTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGcaa                         >  1:5968153/1‑101 (MQ=255)

                                                                          cccTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGcaaca                       <  2:5967387/101‑1 (MQ=255)

                                                                          cccTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGcaaca                       >  2:5968157/1‑101 (MQ=255)

                                                                           ccTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGcaacaa                      >  2:5968162/1‑101 (MQ=255)

                                                                           ccTGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGcaacaa                      >  1:5968163/1‑101 (MQ=255)

                                                                             tGACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAAt                    >  2:5968167/1‑101 (MQ=255)

                                                                              gACGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAAtt                   >  2:5968172/1‑101 (MQ=255)

                                                                               aCGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAAt                    >  1:5968175/1‑99 (MQ=255)

                                                                               aCGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTa                  <  2:5967663/101‑1 (MQ=255)

                                                                               aCGGGCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTa                  <  1:5967408/101‑1 (MQ=255)

                                                                                 gggCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAgg                >  2:5968180/1‑101 (MQ=255)

                                                                                 gggCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAgg                >  1:5968181/1‑101 (MQ=255)

                                                                                 gggCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAgg                <  2:5967785/101‑1 (MQ=255)

                                                                                  ggCGCGCACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAAtt                   <  1:5967844/97‑1 (MQ=255)

                                                                                      cgcACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATAc           <  1:5967522/101‑1 (MQ=255)

                                                                                       gcACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACt          <  1:5967717/101‑1 (MQ=255)

                                                                                       gcACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACt          <  2:5967607/101‑1 (MQ=255)

                                                                                       gcACGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACt          >  2:5968188/1‑101 (MQ=255)

                                                                                         aCGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAg                             <  2:5967587/80‑1 (MQ=255)

                                                                                         aCGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtg        >  1:5968191/1‑101 (MQ=255)

                                                                                         aCGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtg        >  2:5968192/1‑101 (MQ=255)

                                                                                          cGTTTGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtga       <  1:5967646/101‑1 (MQ=255)

                                                                                            tttGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTtgatg     >  2:5968197/1‑101 (MQ=255)

                                                                                              tGAGCACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACaaa                     <  2:5967366/83‑1 (MQ=255)

GCGCGGATGGAGCAGGCGATGCCACTCTCCGCATTTCTGTTTAACAGTCTGATGCCGCAAGTTGATCTGAGTACCCCTGACGGGCGCGCACGTTTGAGTACGCTGGCACTACCATTGATATCGCAAGTGCCGGGCGAAACGCTGCGAATATATCTTCGTCAGGAATTAGGCAACAAATTAGGCATACTTGATGACA  >  NC_000913/3212184‑3212379

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