BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	716,235	A→C	pseudogene (363/651 nt)	ybfG ←	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	716,235	0	A	C	96.3%	438.8 / 2.5	134	pseudogene (363/651 nt)	ybfG	pseudogene
Reads supporting (aligned to +/- strand): ref base A (3/2); new base C (54/75); total (57/77)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

AATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGT  >  NC_000913/716138‑716333

aaaaaCCGGATAATGCCCACTGAAGAAACTCAGAGGATCCAACTTCAAGCCTGTCCCAAGTCAACGGTTCTGGCGCAGAATCATAAAGGCAGCCAGGATc                                                                                                  >  2:1780361/5‑100 (MQ=255)

aaTCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGtctc                                                                                                  <  1:2351863/100‑1 (MQ=255)

    aCCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTcc                                                                                              >  1:1233199/1‑100 (MQ=255)

    aCCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTcc                                                                                              >  2:2892411/1‑100 (MQ=255)

     ccGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTccc                                                                                             <  2:1787560/100‑1 (MQ=255)

       ggATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCa                                                                                           <  1:776320/100‑1 (MQ=255)

         aTAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAg                                                                                         >  2:2267374/1‑100 (MQ=255)

         aTAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAACAAACGCAGCCAGTCTCTTCCCCAAg                                                                                         <  2:3227413/100‑1 (MQ=255)

           aaTGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGAt                                                                                       <  1:2939638/100‑1 (MQ=255)

           aaTGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGAt                                                                                       <  2:2242169/100‑1 (MQ=255)

            aTGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATc                                                                                      <  1:3089727/100‑1 (MQ=255)

            aTGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATc                                                                                      >  2:2590448/1‑100 (MQ=255)

            aTGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATc                                                                                      <  2:2928816/100‑1 (MQ=255)

             tGCCCACTCAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTCTCTGGCGCAAACCACTAAACGCAGCCAGTCTCTTCCCCAAGATCg                                                                                     >  2:818536/1‑100 (MQ=255)

               cccACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGcc                                                                                   <  2:3973009/100‑1 (MQ=255)

               cccACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGcc                                                                                   <  2:1979723/100‑1 (MQ=255)

                ccACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGcca                                                                                  >  1:1956854/1‑100 (MQ=255)

                 cACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGccac                                                                                 <  1:3131503/100‑1 (MQ=255)

                    tGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCAt                                                                              <  1:1461086/100‑1 (MQ=255)

                    tGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCAt                                                                              <  2:1321017/100‑1 (MQ=255)

                     gaagaaACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCCAGATCGCCACCAtt                                                                             >  2:2071198/1‑100 (MQ=255)

                     gaagaaACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCAtt                                                                             >  2:2849327/1‑100 (MQ=255)

                      aagaaACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTa                                                                            <  2:2685079/100‑1 (MQ=255)

                        gaaACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAag                                                                          <  2:2088763/100‑1 (MQ=255)

                        gaaACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAag                                                                          <  1:2646795/100‑1 (MQ=255)

                        gaaACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAag                                                                          <  1:2813353/100‑1 (MQ=255)

                          aaCTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAagag                                                                        <  1:3428375/100‑1 (MQ=255)

                            cTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCa                                                                      >  2:1462723/1‑100 (MQ=255)

                            cTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCa                                                                      >  2:3518572/1‑100 (MQ=255)

                             tCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACTCAACCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCaa                                                                     >  2:1972735/1‑100 (MQ=255)

                               aGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCaaaa                                                                   >  2:1603260/1‑100 (MQ=255)

                               aGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCaaaa                                                                   >  1:3169437/1‑100 (MQ=255)

                                gAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCaaaaa                                                                  <  1:3081284/100‑1 (MQ=255)

                                gAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCaaaaa                                                                  >  1:782601/1‑100 (MQ=255)

                                gAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCaaaaa                                                                  <  1:2274092/100‑1 (MQ=255)

                                 aaTATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAAt                                                                 <  1:2139065/100‑1 (MQ=255)

                                 aaTATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAAt                                                                 <  1:1085926/100‑1 (MQ=255)

                                  atatCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAAtt                                                                >  1:2429641/1‑100 (MQ=255)

                                   tatCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTa                                                               <  2:1702184/100‑1 (MQ=255)

                                   cttCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTa                                                               <  1:158195/98‑1 (MQ=255)

                                    atCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATACACGCAGCCATTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTAc                                                              >  1:2783436/1‑100 (MQ=255)

                                    atCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTAc                                                              <  1:268739/100‑1 (MQ=255)

                                    atCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTAc                                                              <  1:2655651/100‑1 (MQ=255)

                                    atCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTAc                                                              <  2:612350/100‑1 (MQ=255)

                                     tCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTaag                                                             >  2:155575/1‑98 (MQ=255)

                                     tCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTAcc                                                             <  1:2046540/100‑1 (MQ=255)

                                      cccaCTTCCAGAATTTCCGAATACAGCGTATCTGGCGCAAAATAATAAACGCAGACAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAGAATTAccc                                                            <  2:1298632/97‑1 (MQ=255)

                                      ccAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTAccc                                                            <  2:3597091/100‑1 (MQ=255)

                                       cAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTAcccc                                                           <  1:2207596/100‑1 (MQ=255)

                                        aaCTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCt                                                          >  1:3552237/1‑100 (MQ=255)

                                          cTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGc                                                        >  1:468029/1‑100 (MQ=255)

                                          cTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGc                                                        <  1:2892822/100‑1 (MQ=255)

                                           ttCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCg                                                       <  2:1750121/100‑1 (MQ=255)

                                              cAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACa                                                    <  1:2511559/100‑1 (MQ=255)

                                              cAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACa                                                    >  1:1301977/1‑100 (MQ=255)

                                                gACTTGCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATc                                                  <  1:916043/100‑1 (MQ=255)

                                                  cTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGt                                                >  1:2999355/1‑100 (MQ=255)

                                                     tCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAg                                             <  1:954596/100‑1 (MQ=255)

                                                     tCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAg                                             >  1:781803/1‑100 (MQ=255)

                                                      cccAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGc                                            >  2:1496569/1‑100 (MQ=255)

                                                      cccAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGc                                            <  1:3167095/100‑1 (MQ=255)

                                                       ccAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCg                                           >  1:3057685/1‑100 (MQ=255)

                                                         aaTACAGCGTATCTGGCGCAAAATAATAAAAGCAGCAAGTCTCTTCCCCTTGATCGCAAAAAAAAAAAAAAAAAATTACCCCTGCGACATCGTCAGCGAt                                         <  2:2712120/100‑1 (MQ=255)

                                                          aTTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATg                                        >  1:1893628/1‑100 (MQ=255)

                                                          aTTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATg                                        >  2:2813014/1‑100 (MQ=255)

                                                           ttCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATgg                                       >  2:2408708/1‑100 (MQ=255)

                                                           ttCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATgg                                       <  1:3932969/100‑1 (MQ=255)

                                                             cAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGt                                     >  1:2427674/1‑100 (MQ=255)

                                                               gCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTaa                                   >  1:923892/1‑100 (MQ=255)

                                                               gCGTATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTaa                                   <  2:1890276/100‑1 (MQ=255)

                                                                  tATCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCAACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGc                                <  2:3689903/100‑1 (MQ=255)

                                                                   aTCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGcc                               >  1:1533405/1‑100 (MQ=255)

                                                                   aTCTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGcc                               <  2:2581895/100‑1 (MQ=255)

                                                                     cTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCg                             >  2:3732077/1‑100 (MQ=255)

                                                                     cTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCg                             >  1:636650/1‑100 (MQ=255)

                                                                     cTGGCGCAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCg                             >  1:1678393/1‑100 (MQ=255)

                                                                        gcgcAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCAc                          <  2:1022904/100‑1 (MQ=255)

                                                                         cgcAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACg                         >  1:1669814/1‑100 (MQ=255)

                                                                          gcAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGCTGGGTAAAGCCCGCACGc                        >  2:1186604/1‑100 (MQ=255)

                                                                           cAAAATAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGcc                       >  1:2060947/1‑100 (MQ=255)

                                                                             aaaTAATAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTg                     >  2:606372/1‑100 (MQ=255)

                                                                              aataatAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGc                    >  1:2067347/1‑100 (MQ=255)

                                                                              aataatAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGc                    <  2:1590841/100‑1 (MQ=255)

                                                                              aataatAAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGc                    <  2:1590904/100‑1 (MQ=255)

                                                                               ataataAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCTAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCg                   <  1:4052456/100‑1 (MQ=255)

                                                                               ataataAACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCg                   >  2:593047/1‑100 (MQ=255)

                                                                               ataataAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGTAAAAATTACCCCTGCGACAGCGTCAGCGATGGTTAAAGCCCGCACGCCTGCg                   <  2:1701216/100‑1 (MQ=255)

                                                                                   taaACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAAtt               >  1:1711521/1‑100 (MQ=255)

                                                                                   taaACGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAAtt               >  2:886943/1‑100 (MQ=255)

                                                                                    aaaCGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTc              <  2:2577209/100‑1 (MQ=255)

                                                                                     aaCGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCa             >  1:2135723/1‑100 (MQ=255)

                                                                                     aaCGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCa             <  1:3442866/100‑1 (MQ=255)

                                                                                      gatctgCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAg            <  2:3081405/95‑1 (MQ=255)

                                                                                      aCGCAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAg            <  1:3841707/100‑1 (MQ=255)

                                                                                         cAGCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGc         >  1:3490016/1‑100 (MQ=255)

                                                                                           gCCAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGt       >  2:1237203/1‑100 (MQ=255)

                                                                                            ccAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGtt      >  1:637285/1‑100 (MQ=255)

                                                                                             cAGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTtg     <  2:2922118/100‑1 (MQ=255)

                                                                                              aGTCTCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTtgt    >  2:3351445/1‑100 (MQ=255)

                                                                                                tATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTgtgt  <  2:2647985/100‑1 (MQ=255)

AATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGT  >  NC_000913/716138‑716333

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