BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,822,652	T→C	Y171H (TAC→CAC)	spoU →	putative RNA methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,822,652	0	T	C	99.5%	1428.9 / ‑4.2	437	Y171H (TAC→CAC)	spoU	putative RNA methylase
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (212/223); total (212/225)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 445 total aligned reads displayed.

ATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAA  >  NC_000913/3822562‑3822741

  ggtCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCt                                                                                       <  1:1580982/91‑1 (MQ=255)

  gatCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCt                                                                                       >  1:1687580/1‑93 (MQ=255)

    tCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGc                                                                                     <  2:74591/93‑1 (MQ=255)

     cGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCg                                                                                    <  1:8222252/93‑1 (MQ=255)

     cGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCg                                                                                    <  1:6554543/93‑1 (MQ=255)

      ggCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGt                                                                                   >  2:10047549/1‑93 (MQ=255)

      ggCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGt                                                                                   <  2:3436434/93‑1 (MQ=255)

       gCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTg                                                                                  >  2:7918587/1‑93 (MQ=255)

        cATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCCGGCATGCACCTGCGTGa                                                                                 >  1:4171187/1‑93 (MQ=255)

        cATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGa                                                                                 <  2:96086/93‑1 (MQ=255)

        cATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGa                                                                                 >  2:5002671/1‑93 (MQ=255)

         aTGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCATAATGCAGGCATGCACCTGCGTGaa                                                                                >  1:7721923/1‑93 (MQ=255)

          tGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACTAAGCCCAGCGTCAGCGGCAAAATGCAGGCATCCACCTGCGTGaaa                                                                               >  1:67666/1‑93 (MQ=255)

            gTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAAc                                                                             <  2:9123811/93‑1 (MQ=255)

                aGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCATCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCa                                                                         <  2:4019942/93‑1 (MQ=255)

                aGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCa                                                                         >  2:7145250/1‑93 (MQ=255)

                 gTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCAt                                                                        <  1:5173794/93‑1 (MQ=255)

                  tCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATg                                                                       <  1:8283205/93‑1 (MQ=255)

                  tCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATg                                                                       <  1:1000046/93‑1 (MQ=255)

                  tCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATg                                                                       >  2:3789376/1‑93 (MQ=255)

                   cGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGt                                                                      >  1:8464512/1‑93 (MQ=255)

                    gCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGtt                                                                     >  2:5808711/1‑93 (MQ=255)

                    gCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGtt                                                                     >  2:2249468/1‑93 (MQ=255)

                      tGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGc                                                                   >  1:9058444/1‑93 (MQ=255)

                       gAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGcc                                                                  >  1:5067619/1‑93 (MQ=255)

                        aaTGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCg                                                                 <  2:3509750/93‑1 (MQ=255)

                          tGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGa                                                               <  2:8321997/93‑1 (MQ=255)

                            tttCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAg                                                             >  1:458963/1‑93 (MQ=255)

                             ttCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGc                                                            <  1:3799842/93‑1 (MQ=255)

                              tCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCa                                                           <  2:8517236/93‑1 (MQ=255)

                              tCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCa                                                           <  1:2360345/93‑1 (MQ=255)

                               ccGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCag                                                          <  1:1509724/93‑1 (MQ=255)

                               ccGTCGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCag                                                          <  2:3895993/93‑1 (MQ=255)

                                cGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCaga                                                         <  2:4270570/93‑1 (MQ=255)

                                cGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCaga                                                         <  2:1581079/93‑1 (MQ=255)

                                cGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCaga                                                         <  1:6976364/93‑1 (MQ=255)

                                cGTTGCCTCAGCACTCATTCTTTCGGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAGCAGCATGTTGCCGGAAGCaga                                                         <  2:6842365/93‑1 (MQ=255)

                                 gTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCagag                                                        >  1:3893297/1‑93 (MQ=255)

                                 gTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCagag                                                        <  2:1290024/93‑1 (MQ=255)

                                 gTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCagag                                                        <  1:4431882/93‑1 (MQ=255)

                                   tGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGca                                                      <  2:8376907/93‑1 (MQ=255)

                                    gCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaa                                                     <  1:9306644/93‑1 (MQ=255)

                                       tCAGCCCTCATTCTTTGCGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                  <  1:5450419/93‑1 (MQ=255)

                                       tCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                  >  1:8616258/1‑93 (MQ=255)

                                        cAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAACCAGCATGTTGCCGGAAGCAGAGCaacaac                                                 <  2:4577167/93‑1 (MQ=255)

                                        cAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCaacaac                                                 >  2:8802800/1‑93 (MQ=255)

                                         aGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACg                                                <  2:2473757/93‑1 (MQ=255)

                                          gCCCTCATTCTTTACGAAGCCCAGCGTCAGCGTCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGc                                               >  2:7750065/1‑93 (MQ=255)

                                           ncnTCATNCTTTACGAAGCCCNGNGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGcc                                              <  1:8800732/92‑1 (MQ=255)

                                             cTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtg                                            >  1:1850282/1‑93 (MQ=255)

                                              tCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgt                                           <  1:8657281/93‑1 (MQ=255)

                                              tCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgt                                           >  1:5914624/1‑93 (MQ=255)

                                               cATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgtt                                          <  2:666056/93‑1 (MQ=255)

                                                aTTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGAAACAACGCCtgttg                                         >  1:2635457/1‑93 (MQ=255)

                                                 ttCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttgt                                        >  2:3067144/1‑93 (MQ=255)

                                                 ttCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttgt                                        >  1:7513256/1‑93 (MQ=255)

                                                 ttCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAAAGCAACAAAGCCtgttgt                                        >  1:2026802/1‑93 (MQ=255)

                                                  tCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTgttgtt                                       <  2:6142922/93‑1 (MQ=255)

                                                  tCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTgttgtt                                       <  2:8009163/93‑1 (MQ=255)

                                                  tCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTgttgtt                                       <  2:9000853/93‑1 (MQ=255)

                                                   cTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGttt                                      <  2:5697949/93‑1 (MQ=255)

                                                    tttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTg                                     <  2:1754168/93‑1 (MQ=255)

                                                     ttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGa                                    >  2:8760261/1‑93 (MQ=255)

                                                      tACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGaa                                   <  2:4952606/93‑1 (MQ=255)

                                                       aCGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAg                                  <  1:3024995/93‑1 (MQ=255)

                                                        cGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAgg                                 <  2:1968463/93‑1 (MQ=255)

                                                        cGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAgg                                 >  2:4218459/1‑93 (MQ=255)

                                                         gAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggc                                <  1:7389953/93‑1 (MQ=255)

                                                         gAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggc                                <  1:2801660/93‑1 (MQ=255)

                                                           aGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcgg                              >  1:2053165/1‑93 (MQ=255)

                                                            gCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcggc                             >  2:5307422/1‑93 (MQ=255)

                                                             cccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCt                            <  2:4406758/93‑1 (MQ=255)

                                                              ccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGGTGTTTGAAGGCGGCTa                           >  2:9361168/1‑93 (MQ=255)

                                                               cAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTac                          >  1:9624119/1‑92 (MQ=255)

                                                                aGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATc                         <  2:6727510/93‑1 (MQ=255)

                                                                 gCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATcc                        <  2:2743263/93‑1 (MQ=255)

                                                                  cGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCg                       <  1:6084531/93‑1 (MQ=255)

                                                                   gTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCgg                      >  2:3340307/1‑93 (MQ=255)

                                                                   gTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCgg                      <  2:2985007/93‑1 (MQ=255)

                                                                    tCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                     >  1:4820248/1‑93 (MQ=255)

                                                                    tCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                     >  1:1030613/1‑93 (MQ=255)

                                                                    tCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                     <  1:7891859/93‑1 (MQ=255)

                                                                       gCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCt                  >  2:4302081/1‑93 (MQ=255)

                                                                        cGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTg                 <  2:2817656/93‑1 (MQ=255)

                                                                        cGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTg                 >  2:7616393/1‑93 (MQ=255)

                                                                         ggCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTgg                >  1:9017137/1‑93 (MQ=255)

                                                                          gCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGc               >  1:6203819/1‑93 (MQ=255)

                                                                           cAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGACTGTTGTTTGAAGGGGGATACCCGGggcgggcg              >  1:1260623/1‑93 (MQ=255)

                                                                           cAAAATGCAGGCATGCACCTGCGCGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCg              <  2:5526001/93‑1 (MQ=255)

                                                                            aaaaTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGa             >  1:8082470/1‑93 (MQ=255)

                                                                            aaaaTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGa             >  2:6598621/1‑93 (MQ=255)

                                                                             aaaTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGaa            >  2:26052/1‑93 (MQ=255)

                                                                              aaTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGaaa           <  1:7255114/93‑1 (MQ=255)

                                                                              aaTGAAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGaaa           <  2:4650862/93‑1 (MQ=255)

                                                                               aTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAg          <  1:3419315/93‑1 (MQ=255)

                                                                                    ggCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCa     <  1:5718677/93‑1 (MQ=255)

                                                                                     gCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaa    >  2:8239831/1‑93 (MQ=255)

                                                                                      cATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTAGCCGGTGCTGGCGAAAGTCTCaaa   <  2:5146906/93‑1 (MQ=255)

                                                                                       aTGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaaa  >  1:1145784/1‑93 (MQ=255)

                                                                                       aTGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaaa  >  1:9949095/1‑93 (MQ=255)

ATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAA  >  NC_000913/3822562‑3822741

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