BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,848,059	T→C	V400A (GTT→GCT)	sppA →	protease IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,848,059	0	T	C	83.5%	539.0 / 96.1	243	V400A (GTT→GCT)	sppA	protease IV
Reads supporting (aligned to +/- strand): ref base T (18/22); new base C (88/115); total (106/137)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.95e-01

Only 100 of 246 total aligned reads displayed.

CTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCG  >  NC_000913/1847968‑1848148

 tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTg                                                                                         <  1:2082940/93‑1 (MQ=255)

 tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTg                                                                                         >  2:7007797/1‑93 (MQ=255)

 tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTg                                                                                         >  2:3907998/1‑93 (MQ=255)

   cGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTa                                                                                       >  2:5730301/1‑93 (MQ=255)

     tgtTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATc                                                                                     >  2:7470087/1‑93 (MQ=255)

      gtTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCg                                                                                    >  2:3809122/1‑93 (MQ=255)

       ttAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGa                                                                                   <  1:2321975/93‑1 (MQ=255)

       ttAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGa                                                                                   <  1:949300/93‑1 (MQ=255)

         aaTAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATg                                                                                 >  1:2858643/1‑93 (MQ=255)

          aTAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATgg                                                                                <  2:2730719/93‑1 (MQ=255)

          aTAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATgg                                                                                >  1:2483518/1‑93 (MQ=255)

           tAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATggg                                                                               <  1:7012602/93‑1 (MQ=255)

           tAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATggg                                                                               <  1:1679412/93‑1 (MQ=255)

             gCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcg                                                                             <  1:5316321/93‑1 (MQ=255)

               ccAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggc                                                                           >  1:1270919/1‑93 (MQ=255)

               ccAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggc                                                                           >  2:7285259/1‑93 (MQ=255)

               ccAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggc                                                                           >  1:6879647/1‑93 (MQ=255)

                cAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggca                                                                          <  2:4053681/93‑1 (MQ=255)

                cAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggca                                                                          >  1:520471/1‑93 (MQ=255)

                cAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggca                                                                          <  2:3532856/93‑1 (MQ=255)

                  ggcggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatg                                                                        >  2:6055819/1‑93 (MQ=255)

                   gcggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatgg                                                                       <  1:4017867/93‑1 (MQ=255)

                   gcggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatgg                                                                       <  2:4720148/93‑1 (MQ=255)

                    cggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggc                                                                      <  2:5247265/93‑1 (MQ=255)

                     ggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcg                                                                     <  2:3729021/93‑1 (MQ=255)

                      gcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcgg                                                                    <  2:1458859/93‑1 (MQ=255)

                      gcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcgg                                                                    >  2:6363970/1‑93 (MQ=255)

                         gCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggcggcatggcggcat                                                                 <  2:1606084/93‑1 (MQ=255)

                         gCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcggcat                                                                 >  1:5918514/1‑93 (MQ=255)

                         gCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcggcat                                                                 >  1:346218/1‑93 (MQ=255)

                          cGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATc                                                                <  2:231525/93‑1 (MQ=255)

                          cGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATc                                                                <  2:2774184/93‑1 (MQ=255)

                            ttACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCtg                                                              >  2:1766510/1‑93 (MQ=255)

                            ttACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCtg                                                              <  1:302190/93‑1 (MQ=255)

                            ttACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCtg                                                              <  1:1608807/93‑1 (MQ=255)

                            ttACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCtg                                                              >  2:647910/1‑93 (MQ=255)

                               ccGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCtggtg                                                           <  1:4231307/93‑1 (MQ=255)

                               ccGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCtggtg                                                           >  1:2947662/1‑93 (MQ=255)

                                cgcgTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCtggtgg                                                          <  1:1747460/93‑1 (MQ=255)

                                 gcgTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTggtggt                                                         <  2:6119049/93‑1 (MQ=255)

                                 gcgTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTggtggt                                                         >  2:6778102/1‑93 (MQ=255)

                                 gcgTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTggtggt                                                         <  2:6539533/93‑1 (MQ=255)

                                 gcgTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTggtggt                                                         <  2:5690948/93‑1 (MQ=255)

                                  cgTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGtt                                                        >  2:2635952/1‑93 (MQ=255)

                                    tCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTAc                                                      <  2:1669689/93‑1 (MQ=255)

                                     cTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACt                                                     >  2:4918087/1‑93 (MQ=255)

                                     cTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACt                                                     >  2:2534092/1‑93 (MQ=255)

                                      tGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTg                                                    <  1:18570/93‑1 (MQ=255)

                                       gAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTgg                                                   >  1:7290782/1‑93 (MQ=255)

                                         aGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGAt                                                 >  2:5355443/1‑93 (MQ=255)

                                            gATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTc                                              <  1:3171955/93‑1 (MQ=255)

                                            gATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTc                                              >  1:2806204/1‑93 (MQ=255)

                                            gATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTc                                              <  1:1998844/93‑1 (MQ=255)

                                            gATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTc                                              <  1:3599956/93‑1 (MQ=255)

                                            gATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTc                                              <  1:6316171/93‑1 (MQ=255)

                                               tCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCAc                                           <  2:7193934/93‑1 (MQ=255)

                                                    cTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAg                                      <  2:1253443/93‑1 (MQ=255)

                                                     tGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc                                     <  1:2437393/93‑1 (MQ=255)

                                                     tGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc                                     <  2:5138366/93‑1 (MQ=255)

                                                     tGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc                                     >  2:4511462/1‑93 (MQ=255)

                                                     tGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc                                     <  1:6533988/93‑1 (MQ=255)

                                                     tGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc                                     <  2:1642086/93‑1 (MQ=255)

                                                        aCTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTaa                                  <  2:2461572/93‑1 (MQ=255)

                                                        aCTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTaa                                  <  2:5389130/93‑1 (MQ=255)

                                                         cTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAAt                                 <  2:3623722/93‑1 (MQ=255)

                                                         cTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAAt                                 >  1:6695197/1‑93 (MQ=255)

                                                         cTGGCAGCAGCCCGGGCAGCGGGTAAGACTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGTTAAt                                 >  1:5742882/1‑93 (MQ=255)

                                                          tGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAAtt                                <  2:2893717/93‑1 (MQ=255)

                                                          tGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAAtt                                <  1:1516155/93‑1 (MQ=255)

                                                           ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTa                               <  2:5122128/93‑1 (MQ=255)

                                                            gcagcaGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTAc                              <  2:3095896/93‑1 (MQ=255)

                                                            gcagcaGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTAc                              <  2:1341143/93‑1 (MQ=255)

                                                              agcagcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAt                            <  1:617601/93‑1 (MQ=255)

                                                              agcagcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAt                            <  1:5133971/93‑1 (MQ=255)

                                                               gcagcaCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCAGCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAtt                           <  1:4987023/93‑1 (MQ=255)

                                                                cagcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtg                          <  2:3504656/93‑1 (MQ=255)

                                                                cagcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtg                          <  2:4462483/93‑1 (MQ=255)

                                                                 agcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgt                         <  1:6817894/93‑1 (MQ=255)

                                                                 agcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgt                         <  1:3824039/93‑1 (MQ=255)

                                                                 agcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgt                         <  1:5522420/93‑1 (MQ=255)

                                                                  gcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgtg                        >  2:3964458/1‑93 (MQ=255)

                                                                   cccGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTgg                       <  2:854700/93‑1 (MQ=255)

                                                                   cccGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTgg                       <  1:803411/93‑1 (MQ=255)

                                                                   cccGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTgg                       <  1:1146878/93‑1 (MQ=255)

                                                                    ccGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGc                      >  2:2171164/1‑93 (MQ=255)

                                                                       ggCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGACATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTaa                   >  2:5183973/1‑93 (MQ=255)

                                                                         cAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAcc                 <  1:1092115/93‑1 (MQ=255)

                                                                           gCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAcccc               >  1:1354915/1‑93 (MQ=255)

                                                                            cGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCa              <  2:3459126/93‑1 (MQ=255)

                                                                            cGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCa              <  2:5948209/93‑1 (MQ=255)

                                                                             gggTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAg             >  1:4568167/1‑93 (MQ=255)

                                                                              ggTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGc            >  1:1037369/1‑93 (MQ=255)

                                                                              ggTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGCATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGc            >  1:3353334/1‑93 (MQ=255)

                                                                               gTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCa           <  1:6083904/93‑1 (MQ=255)

                                                                                tAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCAc          <  1:4593926/93‑1 (MQ=255)

                                                                                  aGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCAccc        <  2:4640226/93‑1 (MQ=255)

                                                                                    ccTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTg      >  1:2241456/1‑93 (MQ=255)

                                                                                      tgtgGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGAc    <  2:1857971/93‑1 (MQ=255)

                                                                                      tgtgGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGAc    <  2:6907252/93‑1 (MQ=255)

                                                                                      tgtgGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGAc    >  1:5484835/1‑93 (MQ=255)

CTGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACCG  >  NC_000913/1847968‑1848148

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