BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,648,042	T→C	L59P (CTG→CCG)	arsC →	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,648,042	0	T	C	99.8%	1428.4 / ‑5.4	455	L59P (CTG→CCG)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (226/228); total (226/229)
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.79e-01

Only 100 of 491 total aligned reads displayed.

CGACTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCTGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCACC  >  NC_000913/3647952‑3648131

 gACTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCcgcg                                                                                        >  2:1277541/1‑93 (MQ=255)

 gACTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCcgcg                                                                                        >  1:5433389/1‑93 (MQ=255)

 gACTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCcgcg                                                                                        >  1:2353282/1‑93 (MQ=255)

 gACTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCcgcg                                                                                        >  2:7341926/1‑93 (MQ=255)

 gACTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCcgcg                                                                                        >  1:2032855/1‑93 (MQ=255)

  aCTATTATCCATTATCTGGAAACTACGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTATGCGCGCTGCCGCGt                                                                                       >  1:1039514/1‑93 (MQ=255)

   cTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTa                                                                                      >  2:3915666/1‑93 (MQ=255)

    tattatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaa                                                                                     >  1:4158555/1‑93 (MQ=255)

    tattatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaa                                                                                     >  1:3992820/1‑93 (MQ=255)

     attatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaaa                                                                                    <  2:1405309/93‑1 (MQ=255)

     attatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaaa                                                                                    <  2:1494449/93‑1 (MQ=255)

      ttatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaaaa                                                                                   >  1:2491232/1‑93 (MQ=255)

       tatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaaaaa                                                                                  <  2:2195415/93‑1 (MQ=255)

       tatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaaaaa                                                                                  <  2:2275948/93‑1 (MQ=255)

       tatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaaaaa                                                                                  <  2:5527770/93‑1 (MQ=255)

       tatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaaaaa                                                                                  >  1:7294788/1‑93 (MQ=255)

       tatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCCTaaaaa                                                                                  <  2:476084/93‑1 (MQ=255)

         tCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACg                                                                                >  2:5697554/1‑93 (MQ=255)

          ccATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCTGCGTAAAAACGt                                                                               <  2:2801013/93‑1 (MQ=255)

          ccATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGt                                                                               >  1:2717287/1‑93 (MQ=255)

            aTTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTca                                                                             >  1:6857614/1‑92 (MQ=255)

              tATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGaa                                                                           <  1:5000767/93‑1 (MQ=255)

                tCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAAcc                                                                         >  1:2582315/1‑93 (MQ=255)

                  tGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGt                                                                       >  1:7777474/1‑93 (MQ=255)

                    gAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTAt                                                                     <  2:2450093/93‑1 (MQ=255)

                    gAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTAt                                                                     <  1:6442810/93‑1 (MQ=255)

                    gAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTAt                                                                     <  1:5462582/93‑1 (MQ=255)

                      aaCTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATga                                                                   <  2:3365829/93‑1 (MQ=255)

                        cTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATgagg                                                                 <  1:3292468/93‑1 (MQ=255)

                            gCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCt                                                             <  1:7140710/93‑1 (MQ=255)

                             ccAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGACCCGTATGAGGAGCTg                                                            >  2:2333884/1‑93 (MQ=255)

                             ccAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTg                                                            <  2:1776072/93‑1 (MQ=255)

                             ccAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTg                                                            >  1:4499229/1‑93 (MQ=255)

                               aaCGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTggg                                                          >  1:6284526/1‑93 (MQ=255)

                               aaCGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTggg                                                          >  2:7458182/1‑93 (MQ=255)

                                aCGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGc                                                         >  2:302428/1‑93 (MQ=255)

                                  gcgcgATGAACTGGTCTAACTCATTGCCGATATGGGGATTTCCGTCCGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCt                                                       <  1:6994815/93‑1 (MQ=255)

                                  gcgcgATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCt                                                       <  2:3443642/93‑1 (MQ=255)

                                  gcgcgATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCt                                                       <  2:7686634/93‑1 (MQ=255)

                                    gcgATGAACGGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTg                                                     <  1:5270576/93‑1 (MQ=255)

                                     cgATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGc                                                    <  1:1802261/93‑1 (MQ=255)

                                      gATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCg                                                   >  2:2538284/1‑93 (MQ=255)

                                       aTGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCgg                                                  >  2:6464084/1‑93 (MQ=255)

                                       aTGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCgg                                                  <  2:6706561/93‑1 (MQ=255)

                                         gAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGaa                                                <  2:434201/93‑1 (MQ=255)

                                          aaCTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAg                                               >  1:3963099/1‑93 (MQ=255)

                                          aaCTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAg                                               >  2:199053/1‑93 (MQ=255)

                                           aCTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGa                                              <  1:3678282/93‑1 (MQ=255)

                                           aCTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGa                                              <  2:7104620/93‑1 (MQ=255)

                                            cTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGAt                                             >  2:580940/1‑93 (MQ=255)

                                            cTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGAt                                             >  1:920048/1‑93 (MQ=255)

                                             tGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATa                                            <  2:1842483/93‑1 (MQ=255)

                                             tGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATa                                            >  1:558166/1‑93 (MQ=255)

                                             tGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATa                                            >  1:4606721/1‑93 (MQ=255)

                                              ggTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATaa                                           >  2:6184141/1‑93 (MQ=255)

                                               gTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATaaa                                          <  2:5415148/93‑1 (MQ=255)

                                               gTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATaaa                                          <  1:2769739/93‑1 (MQ=255)

                                                 cAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAAtt                                        >  1:4849051/1‑93 (MQ=255)

                                                 cAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAAtt                                        >  2:2913890/1‑93 (MQ=255)

                                                 cAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAAtt                                        <  2:3237487/93‑1 (MQ=255)

                                                  aaaCTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAAttt                                       >  1:2240140/1‑93 (MQ=255)

                                                     cTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACt                                    >  2:6279490/1‑93 (MQ=255)

                                                      tCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTg                                   <  2:4058328/93‑1 (MQ=255)

                                                       cATTGCCGATATGTGGATTGCCGTACGCGCGCTGCCGCGTAAAAGCGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGa                                  <  1:7730326/93‑1 (MQ=255)

                                                       cATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGa                                  >  2:6012329/1‑93 (MQ=255)

                                                         ttGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACg                                <  1:6765139/93‑1 (MQ=255)

                                                         ttGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACg                                <  1:6652650/93‑1 (MQ=255)

                                                           gCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGAt                              <  2:5037935/93‑1 (MQ=255)

                                                            ccGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATc                             <  2:798741/93‑1 (MQ=255)

                                                            ccGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATc                             <  2:1037103/93‑1 (MQ=255)

                                                            ccGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATc                             >  1:6040263/1‑93 (MQ=255)

                                                                tatGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCAGAAGATAAATTTACTGACGATCGGtt                         <  2:726181/93‑1 (MQ=255)

                                                                tatGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGACCTTGCGGAAGATAAATTTACTGACGATCGGtt                         >  2:2052038/1‑93 (MQ=255)

                                                                 atGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTa                        <  2:6622094/93‑1 (MQ=255)

                                                                  tGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTaa                       >  1:1405497/1‑93 (MQ=255)

                                                                   ggggATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAAt                      <  2:6982263/93‑1 (MQ=255)

                                                                   ggggATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAAt                      <  1:4717414/93‑1 (MQ=255)

                                                                     ggATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTACTca                    >  2:3135751/1‑92 (MQ=255)

                                                                      gATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGa                   <  1:3336915/93‑1 (MQ=255)

                                                                        tttCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACt                 >  2:4579266/1‑93 (MQ=255)

                                                                        tttCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACt                 <  1:2112117/93‑1 (MQ=255)

                                                                           ccGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGGGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTa              <  1:3187134/93‑1 (MQ=255)

                                                                           ccGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTa              <  1:1879693/93‑1 (MQ=255)

                                                                           ccGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTa              <  1:5848905/93‑1 (MQ=255)

                                                                           ccGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGAGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTa              <  2:2690010/93‑1 (MQ=255)

                                                                            cGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTAt             >  1:4256876/1‑93 (MQ=255)

                                                                            cGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTAt             <  2:469452/93‑1 (MQ=255)

                                                                            cGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTAt             <  1:7883088/93‑1 (MQ=255)

                                                                             gTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATg            >  1:1284548/1‑93 (MQ=255)

                                                                             gTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATg            >  2:678899/1‑93 (MQ=255)

                                                                              tACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATACATTTACTGACGATCGGTTAATCGACTTTATGc           >  1:7409043/1‑93 (MQ=255)

                                                                               aCGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCt          <  1:408504/93‑1 (MQ=255)

                                                                                cgcgcgCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCtt         >  2:7837596/1‑93 (MQ=255)

                                                                                 gcgcgcTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTc        <  1:3125427/93‑1 (MQ=255)

                                                                                 gcgcgcTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTc        <  1:387089/93‑1 (MQ=255)

                                                                                  cgcgcTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCa       <  1:3473643/93‑1 (MQ=255)

                                                                                  cgcgcTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCa       <  1:2999725/93‑1 (MQ=255)

                                                                                  cgcgcTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCa       <  1:5745885/93‑1 (MQ=255)

                                                                                     gctgccGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCa    <  2:3671127/93‑1 (MQ=255)

                                                                                       tgccgcGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCAcc  <  2:4886411/93‑1 (MQ=255)

CGACTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCTGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCACC  >  NC_000913/3647952‑3648131

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