BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	649,964:1	+G	100%	coding (527/909 nt)	citE ←	citrate lyase, citryl‑ACP lyase (beta) subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	649,960	1	.	G	89.6%	287.4 / 30.9	96	coding (531/909 nt)	citE	citrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand): ref base . (6/4); new base G (51/34); total (58/38)
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 108 total aligned reads displayed.

GTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑G‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  NC_000913/649861‑650054

gttgttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGttcc‑ttc‑c                                                                                                    <  1:583835/100‑1 (MQ=255)

   gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑g‑gg                                                                                               <  1:1249179/100‑1 (MQ=255)

   gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑g‑gg                                                                                               <  2:40550/100‑1 (MQ=255)

     tAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑g‑gg‑gg                                                                                            >  2:2900953/1‑99 (MQ=255)

     tAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑g‑gg‑gg                                                                                            <  1:2744847/100‑2 (MQ=255)

      aGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑g‑gg‑gga                                                                                           >  1:2842646/1‑98 (MQ=255)

       gCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACATCGTGCGAAAAGCAGTTCAGTTTC‑TTT‑T‑G‑GG‑G‑AGc                                                                                         >  1:2846397/1‑100 (MQ=255)

       gCGTCGGAATAGACGGTATCGAACGACTTAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGTACAGTAGTTCAGTTCC‑TTC‑C‑g‑gg‑ggag                                                                                          <  2:388026/100‑4 (MQ=255)

         gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑g‑gg‑ggagcg                                                                                        <  1:1523511/100‑6 (MQ=255)

            ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTc                                                                                     >  1:390492/1‑100 (MQ=255)

            ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTc                                                                                     >  1:2217696/1‑100 (MQ=255)

            ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTc                                                                                     <  1:2080509/100‑1 (MQ=255)

               aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCAAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGt                                                                                  >  1:2132552/1‑100 (MQ=255)

                 aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTAc                                                                                >  2:445443/1‑100 (MQ=255)

                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACg                                                                               >  1:70469/1‑100 (MQ=255)

                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCATGCGAACAGCAGCTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACg                                                                               >  2:920459/1‑100 (MQ=255)

                   aCGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGc                                                                              >  2:2141800/1‑100 (MQ=255)

                   aCGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAACAGTTCAGTTCA‑TTC‑CAG‑GA‑G‑AGCGTTCTGTACGc                                                                              >  1:647384/1‑100 (MQ=255)

                    cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGc                                                                              <  2:1867814/99‑1 (MQ=255)

                    cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCa                                                                             >  1:1289279/1‑100 (MQ=255)

                    cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCa                                                                             >  2:2903585/1‑100 (MQ=255)

                    cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCCAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCa                                                                             >  2:593695/1‑100 (MQ=255)

                     ggTATCGAACGTCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAg                                                                            <  1:1939187/100‑1 (MQ=255)

                       tATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGt                                                                          >  2:2848943/1‑100 (MQ=255)

                         tCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTg                                                                        >  1:1374365/1‑100 (MQ=255)

                          cGAACGCCTGAATACCCGCAAAGCGAGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGttac                                                                       >  2:1140062/1‑98 (MQ=255)

                           ggaCGCCTGAATACCCGCAGAGCGCGCGGCCTGTAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTgcg                                                                      >  1:2058214/3‑100 (MQ=255)

                           gAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAATAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTgcg                                                                      <  2:895651/100‑1 (MQ=255)

                           gAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTgcg                                                                      >  1:2229706/1‑100 (MQ=255)

                            aaCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTgcgc                                                                     >  2:2519442/1‑100 (MQ=255)

                             aCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGca                                                                    <  2:1401823/100‑1 (MQ=255)

                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑GTGG‑G‑AGCGTTCTGTACGCAGGTTGCGca                                                                    <  2:1012304/99‑1 (MQ=255)

                                ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACAt                                                                 >  2:766267/1‑100 (MQ=255)

                                ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGCAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑CGCGTTCTGTACGCAGGTTGCGCACAt                                                                 >  2:652401/1‑100 (MQ=255)

                                catGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAAGAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCTTTCTGTACGCAGGTTGCGCACAt                                                                 <  2:931591/98‑1 (MQ=255)

                                 cTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACATTGTGTTAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACAt                                                                 <  1:805333/99‑1 (MQ=255)

                                 cTGAATACCCGCAGAGCGCGCGGCCTGAAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATa                                                                >  2:2091350/1‑100 (MQ=255)

                                  tGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAg                                                               <  2:2229683/100‑1 (MQ=255)

                                   gAATACCCGTTGTGCGCGAGGCCTGCAAAAGGGAACAGCTTGCGAACAGCAGTTCAGTTCC‑GTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGca                                                                    <  1:1238241/94‑1 (MQ=255)

                                      tACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGtct                                                            <  2:2821898/99‑1 (MQ=255)

                                       aCCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          >  2:1468215/1‑100 (MQ=255)

                                       aCCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          <  1:1016854/100‑1 (MQ=255)

                                        cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTCAG‑G‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  1:757457/1‑100 (MQ=255)

                                        cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  2:1206117/1‑100 (MQ=255)

                                        cccGCAGAGCGCGCGACCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTGC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          >  2:1604471/1‑99 (MQ=255)

                                         ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          >  1:2030853/1‑98 (MQ=255)

                                         ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTAAGCAGGTTGCGCACATAGTCTTCTg                                                        >  1:2540435/1‑100 (MQ=255)

                                         ccGCAGAGCGCGCGGACTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCATTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  2:1096838/1‑99 (MQ=255)

                                          cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc                                                       >  2:1450459/1‑100 (MQ=255)

                                          cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc                                                       <  1:1450566/100‑1 (MQ=255)

                                           gCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCa                                                      >  2:2113546/1‑100 (MQ=255)

                                              gagCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCATGTTGTGCACATAGTCTTCAGCACCg                                                   <  2:1205683/100‑1 (MQ=255)

                                              gagCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑G‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGa                                                  >  2:740197/1‑100 (MQ=255)

                                                 cgcgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGGGCACATAGTCTTCTGCACCGAgg                                                <  2:2163862/100‑1 (MQ=255)

                                                    gcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCg                                             <  2:1657768/100‑1 (MQ=255)

                                                    gcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCg                                             >  2:1269168/1‑100 (MQ=255)

                                                     cgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGAACATAGTCTTCTGCACCGAGGGCGa                                            >  1:2598806/1‑100 (MQ=255)

                                                      ggCCTGCAAAATGTAACAGCGTGGGAACAGTAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCATGTTGCGCACATAGTCTTCTGCACCGAGGGc                                              <  2:2060758/97‑1 (MQ=255)

                                                         cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        <  2:1347036/100‑1 (MQ=255)

                                                         cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCACACATAGTCTTCTGCACCGAGGGCGATAcc                                        <  1:2422445/100‑1 (MQ=255)

                                                         cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑G‑GG‑GTAGCGTTCTGTACGAAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        >  1:103630/1‑100 (MQ=255)

                                                          tGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCg                                       <  2:2522851/100‑1 (MQ=255)

                                                           gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                      >  1:2580583/1‑100 (MQ=255)

                                                           gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                      <  1:127130/100‑1 (MQ=255)

                                                           gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGACTTCTGCACCGAGGGCGATACCGa                                      >  1:1659898/1‑100 (MQ=255)

                                                            cAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGATTGCAAACATAGTCTTCTGCACCGATGGCGATACCGAt                                     >  1:889355/1‑100 (MQ=255)

                                                             aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    <  2:298572/100‑1 (MQ=255)

                                                              aaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCa                                   >  2:2021197/1‑100 (MQ=255)

                                                              aaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCa                                   <  2:70469/100‑1 (MQ=255)

                                                              aaaTGGAACAGCGTGCGAACAGCAGTTAAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCa                                   >  2:2325844/1‑100 (MQ=255)

                                                                 tGGAGCAGCGTGCGAACGGCAGTTCAGTTCC‑GTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaaa                                 <  2:2491958/99‑1 (MQ=255)

                                                                  tgAACAGCGTGCGAACAGCAGTTAAGTTCT‑TGCAG‑G‑GG‑G‑AGCGTTCTGTGTGCAGGTTGCGTTCATAGTCTTCTTCATCGAGGGCGTTACCGATCAAAc                                >  2:1988280/2‑99 (MQ=255)

                                                                  ggAACATCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑GTGG‑G‑AGCGTTCTGTACGCAGGTTGCGAAGATAGTCTTCTGCACCGAGGGCGATACCGATCAAAc                                <  1:19692/99‑1 (MQ=255)

                                                                  ggAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACg                               >  2:1831480/1‑100 (MQ=255)

                                                                  ggAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCAATACCGATCAAACg                               <  2:1996979/100‑1 (MQ=255)

                                                                        gCGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGa                         >  1:957461/1‑100 (MQ=255)

                                                                         cGTGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGaa                        >  1:478750/1‑100 (MQ=255)

                                                                           tGCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGc                      >  1:2472760/1‑100 (MQ=255)

                                                                            gCGAACAGCAGTTCAGTTCCTTTC‑G‑G‑GG‑G‑AGTTTTTTGTACGCAGGTTGTGCACATAGTCTGCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     <  1:427463/100‑1 (MQ=255)

                                                                            gCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     >  2:16781/1‑100 (MQ=255)

                                                                            gCGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     >  1:2214643/1‑100 (MQ=255)

                                                                             cGAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGt                    >  1:1618089/1‑100 (MQ=255)

                                                                             cGAACAGCAGTTCAGTCGC‑TTC‑C‑G‑GGAG‑AGCGTTCTGTACGCAGATTGCGTACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGTAAGCGt                    <  2:836258/100‑1 (MQ=255)

                                                                              gAACAGCAGTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTg                   <  1:2938920/100‑1 (MQ=255)

                                                                                 cagcagTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTGCGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATTCTGATCAATCGCTCGGAAGCGTGAg                 <  2:1412706/99‑1 (MQ=255)

                                                                                  agcagTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCg               <  1:2894858/100‑1 (MQ=255)

                                                                                   gcagTTCAGTTCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGa              >  1:1936857/1‑100 (MQ=255)

                                                                                    cagttcagttCC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAt             >  1:1831194/1‑100 (MQ=255)

                                                                                     agttcagttcC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAtt            <  2:1877143/100‑1 (MQ=255)

                                                                                        tcagttcC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTcc         <  2:195879/100‑1 (MQ=255)

                                                                                         cagttcC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCGCATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCa        >  1:2908768/1‑100 (MQ=255)

                                                                                         cagttcC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCa        <  2:2977929/100‑1 (MQ=255)

                                                                                          agttcC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCAc       >  1:1698515/1‑100 (MQ=255)

                                                                                          agttcC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCAc       <  1:993858/100‑1 (MQ=255)

                                                                                          agttcC‑TTC‑CGG‑GG‑G‑AGCGTTCTGGACGCAGGTTGCGCACATAGTGTTCTGCACAGAGGGGAATACCGATCGAAAGCTCGGAAGCGGGAGCGACTTcaa        >  2:2487105/1‑97 (MQ=255)

                                                                                           gttcC‑TTC‑CGG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACt      >  2:2407259/1‑100 (MQ=255)

                                                                                              cc‑ttc‑cgG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcg   >  2:2887314/1‑100 (MQ=255)

                                                                                              cc‑ttc‑cgG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcg   <  1:2783311/100‑1 (MQ=255)

                                                                                               c‑ttc‑cgG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc  >  1:773137/1‑100 (MQ=255)

                                                                                               c‑ttc‑cgG‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGCAAGCGTGAGCGATTTCCACTgcgc  >  2:1304927/1‑100 (MQ=255)

GTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCC‑TTC‑C‑G‑GG‑G‑AGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  NC_000913/649861‑650054

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