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.8%	1528.1 / NA	473	Y171H (TAC→CAC)	spoU	putative RNA methylase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (252/220); total (253/220)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 487 total aligned reads displayed.

ATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCAAAA  >  NC_000913/3822562‑3822741

atgatCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCAc                                                                                         <  2:5395592/93‑1 (MQ=255)

  gatCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCt                                                                                       <  2:6877443/93‑1 (MQ=255)

   atCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTg                                                                                      <  1:7547247/93‑1 (MQ=255)

    tCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGc                                                                                     >  2:7526515/1‑93 (MQ=255)

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

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

          tGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGaga                                                                               >  1:6052887/1‑91 (MQ=255)

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

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

           ggTGCAGTCGCTGATTGTTTCCGTTGCCTCAGCCCTCATTCTTTCCGACGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGaaaa                                                                              <  1:6996708/93‑1 (MQ=255)

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

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

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

                aGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCa                                                                         <  2:2194511/93‑1 (MQ=255)

                aGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCa                                                                         <  2:1632884/93‑1 (MQ=255)

                aGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCa                                                                         <  2:1352258/93‑1 (MQ=255)

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

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

                    gCTGAATGTTTCCGTCGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGtt                                                                     <  1:3883407/93‑1 (MQ=255)

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

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

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

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

                       gAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGcc                                                                  <  2:4008498/93‑1 (MQ=255)

                        aaTGTTTTCGTTGTCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCg                                                                 >  2:6995624/1‑93 (MQ=255)

                        aaTGTTTTCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCg                                                                 >  2:6198054/1‑93 (MQ=255)

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

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

                           gTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGaa                                                              <  2:5635101/93‑1 (MQ=255)

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

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

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

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

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

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

                               ccGTTGCCTCAGCCCTCATGCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCag                                                          <  1:1867185/93‑1 (MQ=255)

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

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

                                cGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCaga                                                         >  2:3601495/1‑93 (MQ=255)

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

                                 cttGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCagag                                                        <  2:5902413/92‑1 (MQ=255)

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

                                    gCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaa                                                     >  1:3590273/1‑93 (MQ=255)

                                      cTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaaca                                                   >  2:2648170/1‑93 (MQ=255)

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

                                       tCAGCCCTCATTCTTTACGAAGCACAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGcaacaa                                                  <  2:4368652/93‑1 (MQ=255)

                                        cAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCaaccac                                                 >  1:5424103/1‑93 (MQ=255)

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

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

                                          gCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGc                                               >  1:5289805/1‑93 (MQ=255)

                                           cccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGcc                                              >  2:4191489/1‑93 (MQ=255)

                                           cccTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGcc                                              <  2:3119875/93‑1 (MQ=255)

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

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

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

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

                                               cATTCTTTACGAAGCCCAGCGTCAGCGTCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgtt                                          <  1:4607458/93‑1 (MQ=255)

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

                                                aTTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttg                                         <  2:57745/93‑1 (MQ=255)

                                                 ttCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGc                                                                   >  1:1044827/1‑66 (MQ=255)

                                                 ttCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGc                                                                   <  2:1044827/66‑1 (MQ=255)

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

                                                 ttCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCtgttgt                                        <  1:1555662/93‑1 (MQ=255)

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

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

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

                                                        tgAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAgg                                 >  1:6616652/2‑93 (MQ=255)

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

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

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

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

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

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

                                                             cccAGCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCt                            >  2:4759053/1‑93 (MQ=255)

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

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

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

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

                                                                 gCGTCAGCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATcg                        >  2:6706267/1‑92 (MQ=255)

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

                                                                 acctCAGCGGCAAAATGCAGGCATGCACCTGCGTTAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATcc                        <  2:1848236/90‑1 (MQ=255)

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

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

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

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

                                                                       gCGGCAAAATGCAGGCATGCACGTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCt                  <  2:7778526/93‑1 (MQ=255)

                                                                       gCGGCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACCACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCt                  >  1:7610747/1‑93 (MQ=255)

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

                                                                         ggCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTgg                <  2:3289322/93‑1 (MQ=255)

                                                                          gCAAAATGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTggg               >  2:2780252/1‑92 (MQ=255)

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

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

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

                                                                                tGCAGGCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGt         >  1:4293936/1‑93 (MQ=255)

                                                                                    ggCATGCACCTGCGTGAAAACAGCGTGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCa     <  2:1144969/93‑1 (MQ=255)

                                                                                     gCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaa    <  1:4004713/93‑1 (MQ=255)

                                                                                     gCATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaa    <  2:5329771/93‑1 (MQ=255)

                                                                                      cATGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaa   >  1:249613/1‑93 (MQ=255)

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

                                                                                       aTGCACCTGCGTGAAAACAGCATGTTGCCGGAAGCAGAGCAACAACGCCTGTTGTTTGAAGGCGGCTATCCGGTGCTGGCGAAAGTCGCaaaa  <  2:2099544/93‑1 (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` `< 30 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`