BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,098,623	T→C	R174R (CGT→CGC)	ycdX →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,098,623	0	T	C	99.5%	666.0 / ‑4.6	205	R174R (CGT→CGC)	ycdX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (95/109); minor base A (0/1); total (95/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

Only 100 of 209 total aligned reads displayed.

GCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGTGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAATTTGAAGA  >  NC_000913/1098534‑1098712

gCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAg                                                                                        >  2:1663739/1‑93 (MQ=255)

 cTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAgg                                                                                       >  1:1766905/1‑93 (MQ=255)

  tGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGt                                                                                      >  2:4253914/1‑93 (MQ=255)

   gAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAgcga                                                                                     >  1:3233723/1‑90 (MQ=255)

   gAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTa                                                                                     <  2:4387429/93‑1 (MQ=255)

   gAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTa                                                                                     >  1:5451353/1‑93 (MQ=255)

   gAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTa                                                                                     >  2:5072200/1‑93 (MQ=255)

   gAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTa                                                                                     >  2:5743234/1‑93 (MQ=255)

   gAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTa                                                                                     >  1:4336929/1‑93 (MQ=255)

     agcagcCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGc                                                                                   <  2:3443138/93‑1 (MQ=255)

     agcagcCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGc                                                                                   >  1:3460134/1‑93 (MQ=255)

      gcagcCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGcc                                                                                  <  1:1271746/93‑1 (MQ=255)

        agcCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGc                                                                                >  1:6109984/1‑93 (MQ=255)

         gcCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCa                                                                               >  2:1463387/1‑93 (MQ=255)

           cgcgAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGc                                                                             >  1:4548476/1‑93 (MQ=255)

           cgcgAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGc                                                                             <  1:4354090/93‑1 (MQ=255)

            gcgAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCg                                                                            <  1:1173882/93‑1 (MQ=255)

              gAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGt                                                                          >  1:4694424/1‑93 (MQ=255)

                aaCATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTAc                                                                        >  1:2387588/1‑93 (MQ=255)

                aaCATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTAc                                                                        >  1:78612/1‑93 (MQ=255)

                 aCATCAGGTGGCGCTGGAAAGCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACg                                                                       <  1:2860229/93‑1 (MQ=255)

                  cATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGt                                                                      >  1:6644305/1‑93 (MQ=255)

                   aTCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTg                                                                     >  1:38794/1‑93 (MQ=255)

                    tCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGa                                                                    >  1:4599721/1‑93 (MQ=255)

                    tCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGa                                                                    >  1:2160944/1‑93 (MQ=255)

                      aGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATg                                                                  >  1:2695249/1‑93 (MQ=255)

                      aGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATg                                                                  >  1:4413057/1‑93 (MQ=255)

                      aGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATg                                                                  >  1:4300003/1‑93 (MQ=255)

                      aGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATg                                                                  >  2:3500288/1‑93 (MQ=255)

                      aGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATg                                                                  <  1:320269/93‑1 (MQ=255)

                       ggTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGc                                                                 >  1:4164503/1‑93 (MQ=255)

                        gTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCt                                                                <  2:405950/93‑1 (MQ=255)

                        gTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCt                                                                <  2:4761061/93‑1 (MQ=255)

                             gcTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCGGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAgg                                                           <  2:6565600/93‑1 (MQ=255)

                             gcTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAgg                                                           <  1:3552675/93‑1 (MQ=255)

                               tGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGtt                                                         >  1:4426893/1‑93 (MQ=255)

                                ggAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTg                                                        >  1:2841801/1‑93 (MQ=255)

                                  aaaTCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTggg                                                      >  1:2437851/1‑93 (MQ=255)

                                    aTCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTg                                                    <  2:2448947/93‑1 (MQ=255)

                                     tCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTgg                                                   >  2:3082289/1‑93 (MQ=255)

                                       aataatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCa                                                 >  2:5029392/1‑93 (MQ=255)

                                        ataatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCAt                                                <  2:3666800/93‑1 (MQ=255)

                                        ataatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCAt                                                >  2:1025616/1‑93 (MQ=255)

                                         taatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCAtt                                               >  1:3754991/1‑93 (MQ=255)

                                         taatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCAtt                                               >  1:6440162/1‑93 (MQ=255)

                                          aatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTa                                              <  2:530660/93‑1 (MQ=255)

                                          aatTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTa                                              >  1:5971736/1‑93 (MQ=255)

                                           atTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAg                                             <  2:1735461/93‑1 (MQ=255)

                                            ttCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAgg                                            <  1:5927811/93‑1 (MQ=255)

                                             tCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGc                                           >  1:5241864/1‑93 (MQ=255)

                                             tCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGc                                           <  2:2811742/93‑1 (MQ=255)

                                             tCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGc                                           >  1:1466744/1‑93 (MQ=255)

                                             tCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGc                                           <  2:4592297/93‑1 (MQ=255)

                                             tCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTACGc                                           >  1:169226/1‑93 (MQ=255)

                                              ccTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCt                                          >  2:6175078/1‑93 (MQ=255)

                                              ccTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCt                                          <  2:1701354/93‑1 (MQ=255)

                                                tCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCg                                        >  1:403108/1‑93 (MQ=255)

                                                tCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCg                                        <  2:6464968/93‑1 (MQ=255)

                                                  aTTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGa                                      >  2:4002551/1‑93 (MQ=255)

                                                  aTTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGa                                      <  1:5797382/93‑1 (MQ=255)

                                                   tttttACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGAt                                     <  2:3775801/93‑1 (MQ=255)

                                                     tttACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATtc                                   <  2:1296027/93‑1 (MQ=255)

                                                     tttACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTCGAGCATGGGTGGCATTAGGCTCGGATtc                                   >  2:3365549/1‑93 (MQ=255)

                                                      ttACACTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATtct                                  >  2:2425277/1‑93 (MQ=255)

                                                        acacTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTca                                <  1:1535420/93‑1 (MQ=255)

                                                        acacTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTca                                <  2:2033820/93‑1 (MQ=255)

                                                         cacTCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTcac                               >  2:121102/1‑93 (MQ=255)

                                                            tCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACAcc                            <  1:6064753/93‑1 (MQ=255)

                                                            tCACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACAcc                            <  1:3806189/93‑1 (MQ=255)

                                                             cACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACcg                           >  2:1858130/1‑93 (MQ=255)

                                                             cACGTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACcg                           >  1:2754611/1‑93 (MQ=255)

                                                                gTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGt                        >  2:4142632/1‑93 (MQ=255)

                                                                gTAAGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGt                        >  2:141999/1‑93 (MQ=255)

                                                                  aaGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGttt                      <  2:2486446/93‑1 (MQ=255)

                                                                  aaGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGttt                      <  1:188586/93‑1 (MQ=255)

                                                                   aGGGCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTa                     <  2:904723/93‑1 (MQ=255)

                                                                    gggCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTAc                    >  2:382850/1‑93 (MQ=255)

                                                                     ggCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTAcc                   <  2:230395/93‑1 (MQ=255)

                                                                      gCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCa                  <  2:5984483/93‑1 (MQ=255)

                                                                      gCAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCa                  <  2:1759640/93‑1 (MQ=255)

                                                                      gCAGTGAAGACAACTGGCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCa                  <  1:6318712/93‑1 (MQ=255)

                                                                       cAGTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCAt                 <  1:3463145/93‑1 (MQ=255)

                                                                         gTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATgg               <  1:4785563/93‑1 (MQ=255)

                                                                         gTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATgg               <  1:1120412/93‑1 (MQ=255)

                                                                         gTGAAGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATgg               >  1:5572426/1‑93 (MQ=255)

                                                                            aaGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATggggg            <  1:103572/93‑1 (MQ=255)

                                                                            aaGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATggggg            <  1:4509643/93‑1 (MQ=255)

                                                                             aGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGa           <  1:6492593/93‑1 (MQ=255)

                                                                             aGACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGa           >  1:2169669/1‑93 (MQ=255)

                                                                              gACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGaa          <  1:725322/93‑1 (MQ=255)

                                                                              gACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGaa          <  1:4961123/93‑1 (MQ=255)

                                                                              gACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGaa          >  2:5530641/1‑93 (MQ=255)

                                                                              gACAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGaa          <  2:631502/93‑1 (MQ=255)

                                                                                cAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAAtt        <  1:5545162/93‑1 (MQ=255)

                                                                                cAACTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAAtt        <  1:4476549/93‑1 (MQ=255)

                                                                                 aaCTGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAAttt       <  2:2912471/93‑1 (MQ=255)

                                                                                    tGTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAATTTGaa    >  2:2573807/1‑93 (MQ=255)

                                                                                     gTCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAATTTGAag   <  1:214669/93‑1 (MQ=255)

                                                                                      tCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAATTTGAaga  >  1:367568/1‑93 (MQ=255)

                                                                                      tCGCGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAATTTGAaga  <  1:930179/93‑1 (MQ=255)

GCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTCACGTAAGGGCAGTGAAGACAACTGTCGTGAGGTAGCCGCAGCGGTACGTGATGCTGGAGGTTGGGTGGCATTAGGCTCGGATTCTCACACCGCGTTTACCATGGGGGAATTTGAAGA  >  NC_000913/1098534‑1098712

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