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	100.0%	925.3 / NA	283	Y171H (TAC→CAC)	spoU	putative RNA methylase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (130/153); total (130/153)

Only 100 of 285 total aligned reads displayed.

TGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAA  >  NC_000913/3822563‑3822741

tgatCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCAcc                                                                                        <  1:4845985/93‑1 (MQ=255)

  atCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTg                                                                                      >  1:167825/1‑93 (MQ=255)

      gCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGGCCAGCGTCAGCGGCAAAATACAGGCATGCAACTGCGTg                                                                                  >  2:3583197/1‑93 (MQ=255)

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

        aTGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGaa                                                                                >  1:4551244/1‑93 (MQ=255)

         tGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGaaa                                                                               >  1:5347782/1‑93 (MQ=255)

          ggTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGaaaa                                                                              <  1:4927861/93‑1 (MQ=255)

           gTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCCCCTGCGTGAAAAc                                                                             <  2:5784698/93‑1 (MQ=255)

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

            tGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACa                                                                            <  2:4092154/93‑1 (MQ=255)

            tGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACa                                                                            <  1:5766682/93‑1 (MQ=255)

            tGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACa                                                                            >  1:4863224/1‑93 (MQ=255)

             gCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAg                                                                           >  1:3893148/1‑93 (MQ=255)

             gCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAg                                                                           <  1:2527559/93‑1 (MQ=255)

               aGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGACAACAGCa                                                                         >  2:2763854/1‑93 (MQ=255)

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

                gTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCAt                                                                        >  1:3338975/1‑93 (MQ=255)

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

                  cGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGt                                                                      <  1:5149170/93‑1 (MQ=255)

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

                   gCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGtt                                                                     <  1:4077475/93‑1 (MQ=255)

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

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

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

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

                        aTGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCCGCGTCAGCGGCACAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCgg                                                                >  1:3947741/1‑93 (MQ=255)

                        aTGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCgg                                                                >  2:435389/1‑93 (MQ=255)

                         tGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGa                                                               >  1:4233313/1‑93 (MQ=255)

                            ttCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGc                                                            >  2:1771505/1‑93 (MQ=255)

                             tCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCa                                                           >  2:4200543/1‑93 (MQ=255)

                              ccGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGACAACAGCATGTTGCCGGAAGGag                                                          >  2:885797/1‑93 (MQ=255)

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

                               cGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGAGAAAACAGCATGTTGCCGGAAGCaga                                                         <  1:3576268/93‑1 (MQ=255)

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

                                 ttGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGc                                                       >  1:3230616/1‑93 (MQ=255)

                                 ttGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGc                                                       <  1:1884389/93‑1 (MQ=255)

                                 ttGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGc                                                       >  2:2226967/1‑93 (MQ=255)

                                  tGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGca                                                      >  2:3710253/1‑93 (MQ=255)

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

                                    ccTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaac                                                    >  1:6015409/1‑93 (MQ=255)

                                    ccTCACCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGATcaac                                                    <  2:2970464/93‑1 (MQ=255)

                                      tCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                  <  2:1652196/93‑1 (MQ=255)

                                      tCAGCCCTCATTCTTTACGAAGCACAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                  >  2:4610294/1‑93 (MQ=255)

                                        aGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACg                                                >  2:4773552/1‑93 (MQ=255)

                                         gCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGc                                               <  1:2860994/93‑1 (MQ=255)

                                           ccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCt                                             >  1:836/1‑93 (MQ=255)

                                           ccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCt                                             >  2:1479000/1‑93 (MQ=255)

                                           ccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCt                                             >  2:5720829/1‑93 (MQ=255)

                                           ccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCt                                             <  2:5600693/93‑1 (MQ=255)

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

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

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

                                              cATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCAACTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgtt                                          <  1:6248710/93‑1 (MQ=255)

                                               aTTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttg                                         >  1:2729295/1‑93 (MQ=255)

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

                                                 tCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTgttgtt                                       >  1:17434/1‑93 (MQ=255)

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

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

                                                   tttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCGCCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTg                                     <  2:4967639/93‑1 (MQ=255)

                                                   tttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTg                                     >  1:3123446/1‑93 (MQ=255)

                                                    ttACGAAGCCCAGCGTCAGCGGCAAAATGCAGGAATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAAAAACGCCTGTTGTTTGa                                    >  2:532030/1‑93 (MQ=255)

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

                                                     tACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGaa                                   >  2:5014938/1‑93 (MQ=255)

                                                        gAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggc                                >  2:794477/1‑93 (MQ=255)

                                                         aaGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcg                               <  2:6194478/93‑1 (MQ=255)

                                                         aaGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAggcg                               >  1:994578/1‑93 (MQ=255)

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

                                                             ccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTa                           <  1:2914456/93‑1 (MQ=255)

                                                             ccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTa                           >  1:2223757/1‑93 (MQ=255)

                                                              cAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTAt                          <  1:2444084/93‑1 (MQ=255)

                                                               aGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATc                         >  1:4171702/1‑93 (MQ=255)

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

                                                                gCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATcc                        >  2:2392468/1‑93 (MQ=255)

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

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

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

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

                                                                   tCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCTGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGt                     <  1:104268/93‑1 (MQ=255)

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

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

                                                                    cAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTg                    >  2:1020722/1‑93 (MQ=255)

                                                                     aGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGc                   <  1:1141066/93‑1 (MQ=255)

                                                                     aGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGc                   <  2:692169/93‑1 (MQ=255)

                                                                     aGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGc                   <  2:2056674/93‑1 (MQ=255)

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

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

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

                                                                             aaTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGaaa           >  1:6275965/1‑93 (MQ=255)

                                                                             aaTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGaaa           >  1:5161873/1‑93 (MQ=255)

                                                                              aTGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAg          >  1:3271992/1‑93 (MQ=255)

                                                                               tGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGt         <  1:6103791/93‑1 (MQ=255)

                                                                               tGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGt         <  1:2965726/93‑1 (MQ=255)

                                                                                gCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTc        >  2:305468/1‑93 (MQ=255)

                                                                                gCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTc        <  2:4486260/93‑1 (MQ=255)

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

                                                                                   ggCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCa     >  2:3462021/1‑93 (MQ=255)

                                                                                      aTGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGGGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaaa  >  1:141587/1‑93 (MQ=255)

                                                                                      aTGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaaa  <  2:606631/93‑1 (MQ=255)

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

                                                                                      aTGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaaa  <  1:3793449/93‑1 (MQ=255)

TGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAA  >  NC_000913/3822563‑3822741

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