BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,953,902	T→C	Q377R (CAG→CGG)	b1872 ←	biotin sulfoxide reductase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,953,902	0	T	C	100.0%	615.8 / NA	190	Q377R (CAG→CGG)	b1872	biotin sulfoxide reductase 2
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (106/84); total (106/84)

Only 100 of 195 total aligned reads displayed.

ATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATT  >  NC_000913/1953812‑1953991

aTTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGc                                                                                         >  1:629955/1‑93 (MQ=255)

aTTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGc                                                                                         >  1:1931240/1‑93 (MQ=255)

      ggCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAAc                                                                                   >  2:1677815/1‑93 (MQ=255)

      ggCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAAc                                                                                   >  1:305240/1‑93 (MQ=255)

      ggCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAAc                                                                                   <  1:5474878/93‑1 (MQ=255)

       gcagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACa                                                                                  <  1:5232859/93‑1 (MQ=255)

       gcagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACa                                                                                  >  1:2144368/1‑93 (MQ=255)

       gcagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACa                                                                                  <  2:2495657/93‑1 (MQ=255)

        cagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACAt                                                                                 <  2:606781/93‑1 (MQ=255)

        cagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACAt                                                                                 >  1:4399953/1‑93 (MQ=255)

        cagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACAt                                                                                 <  2:131662/93‑1 (MQ=255)

        cagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACAt                                                                                 <  1:5650143/93‑1 (MQ=255)

         agcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                                <  1:5343490/93‑1 (MQ=255)

         agcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                                <  1:3213658/93‑1 (MQ=255)

          gcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcg                                                                               <  2:4389639/93‑1 (MQ=255)

          gcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcg                                                                               <  1:3358095/93‑1 (MQ=255)

           caCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcgc                                                                              <  1:2351090/93‑1 (MQ=255)

               ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCgccgcc                                                                          <  2:1009927/93‑1 (MQ=255)

               ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCgccgcc                                                                          <  2:4402002/93‑1 (MQ=255)

               ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCgccgcc                                                                          >  2:3655640/1‑93 (MQ=255)

               ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCgccgcc                                                                          <  2:5958244/93‑1 (MQ=255)

                cgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCa                                                                         >  2:2814066/1‑93 (MQ=255)

                cgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCa                                                                         >  2:3990538/1‑93 (MQ=255)

                  ccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGc                                                                       >  2:2661276/1‑93 (MQ=255)

                    aaCACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGt                                                                     >  2:4270112/1‑93 (MQ=255)

                     acacGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTc                                                                    <  2:4965977/93‑1 (MQ=255)

                     acacGCGTCGGGTTGACACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTc                                                                    >  2:3087317/1‑93 (MQ=255)

                      cacGCGTCTGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCa                                                                   >  2:6148334/1‑93 (MQ=255)

                      cacGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCa                                                                   >  2:2686994/1‑93 (MQ=255)

                       acGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAc                                                                  >  1:742188/1‑93 (MQ=255)

                        cgcgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAcc                                                                 <  2:5069249/93‑1 (MQ=255)

                           gTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGc                                                              >  1:2011134/1‑93 (MQ=255)

                            tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCa                                                             >  1:2233225/1‑93 (MQ=255)

                             cGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCAt                                                            <  2:3802010/93‑1 (MQ=255)

                              gggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAACCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATc                                                           <  2:4307228/93‑1 (MQ=255)

                              gggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATc                                                           <  1:228730/93‑1 (MQ=255)

                              gggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATc                                                           >  1:2870560/1‑93 (MQ=255)

                               ggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATcc                                                          >  1:5178900/1‑93 (MQ=255)

                               ggTGGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATcc                                                          <  2:2382454/93‑1 (MQ=255)

                                  tGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                                                       >  2:4313493/1‑93 (MQ=255)

                                  tGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                                                       <  2:5028984/93‑1 (MQ=255)

                                  tGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                                                       <  2:2074889/93‑1 (MQ=255)

                                   gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtg                                                      >  2:3481347/1‑93 (MQ=255)

                                    ccACCGTTGGAGTAATGATATGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                     >  2:630741/1‑93 (MQ=255)

                                    ccACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                     <  2:5226626/93‑1 (MQ=255)

                                    ccACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                     <  2:1037808/93‑1 (MQ=255)

                                     cACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                                                    <  1:578458/93‑1 (MQ=255)

                                     cACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                                                    >  1:140404/1‑93 (MQ=255)

                                           tGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                              >  1:4328645/1‑93 (MQ=255)

                                             gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa                                            >  1:908329/1‑93 (MQ=255)

                                               gTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAcc                                          >  1:2947281/1‑93 (MQ=255)

                                               gTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAcc                                          >  1:6201896/1‑93 (MQ=255)

                                                tAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                         >  2:2128655/1‑93 (MQ=255)

                                                tAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                         >  1:2315411/1‑93 (MQ=255)

                                                   tGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCatat                                      >  1:5872458/1‑93 (MQ=255)

                                                     aTAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                                    >  1:2929277/1‑93 (MQ=255)

                                                       aaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                  <  2:3731228/93‑1 (MQ=255)

                                                       aaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                  <  2:1611846/93‑1 (MQ=255)

                                                        agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                 <  1:1048452/93‑1 (MQ=255)

                                                        agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                 >  1:5839830/1‑93 (MQ=255)

                                                         gagaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTgg                                >  2:5437222/1‑93 (MQ=255)

                                                         gagaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTgg                                >  2:4654949/1‑93 (MQ=255)

                                                           gaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcg                              >  2:5846549/1‑93 (MQ=255)

                                                             aaCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCt                            <  1:3510984/93‑1 (MQ=255)

                                                              aCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                           >  1:1228476/1‑93 (MQ=255)

                                                                cAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGaa                         >  2:2226542/1‑93 (MQ=255)

                                                                cAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGaa                         <  1:4764462/93‑1 (MQ=255)

                                                                 aaaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAt                        >  2:1544928/1‑93 (MQ=255)

                                                                 aaaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAt                        >  1:5794958/1‑93 (MQ=255)

                                                                  aaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAtt                       >  1:1197761/1‑93 (MQ=255)

                                                                   aaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc                      <  1:1489929/93‑1 (MQ=255)

                                                                    aCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc                     >  2:4523123/1‑93 (MQ=255)

                                                                      cgccgcCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                   <  2:5659712/93‑1 (MQ=255)

                                                                      cgccgcCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                   <  1:5926588/93‑1 (MQ=255)

                                                                       gctgccCGGTGTACCAATCCGCCCTACCAGCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                  <  2:124216/90‑1 (MQ=255)

                                                                       gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                  <  2:3400764/93‑1 (MQ=255)

                                                                       gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                  <  2:1129420/93‑1 (MQ=255)

                                                                       gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                  <  2:4088465/93‑1 (MQ=255)

                                                                       cccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                  <  2:3871487/92‑1 (MQ=255)

                                                                        ccgccCGGTGTCCCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                 <  1:5142738/93‑1 (MQ=255)

                                                                        ccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                 >  1:2097968/1‑93 (MQ=255)

                                                                         cgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTTTTCACCATATTGCTGGCGCTGAATTCCCCAGc                >  1:498985/1‑93 (MQ=255)

                                                                         cgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGc                >  2:360161/1‑93 (MQ=255)

                                                                           cccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCt              <  2:4123256/93‑1 (MQ=255)

                                                                            ccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCGGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg             >  2:5903218/1‑93 (MQ=255)

                                                                            ccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg             >  1:4049822/1‑93 (MQ=255)

                                                                            ccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg             >  1:2814792/1‑93 (MQ=255)

                                                                             cGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc            >  1:1705514/1‑93 (MQ=255)

                                                                             cGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc            >  1:3929656/1‑93 (MQ=255)

                                                                             cGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc            >  2:1801769/1‑93 (MQ=255)

                                                                               gtgtACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCa          <  1:3452241/93‑1 (MQ=255)

                                                                               gtgtACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCa          >  1:2475627/1‑93 (MQ=255)

                                                                                 gtACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATc        <  2:2188788/93‑1 (MQ=255)

                                                                                  tACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCAAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCa       <  1:3850930/93‑1 (MQ=255)

                                                                                   aCCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGCCGCTGAATTCCCCAGCCTGCCCTCAg      >  1:3053267/1‑93 (MQ=255)

                                                                                    ccAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc     <  2:5493398/93‑1 (MQ=255)

                                                                                     cAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa    <  2:3914111/93‑1 (MQ=255)

                                                                                      aaTCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt   >  1:2994782/1‑93 (MQ=255)

                                                                                      aaTCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt   >  2:709470/1‑93 (MQ=255)

                                                                                       aTCCGCCCTAACAGCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAtt  <  1:1858279/93‑1 (MQ=255)

ATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATT  >  NC_000913/1953812‑1953991

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