BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,683,566	(C)_6→7	coding (1047/1404 nt)	fumC ←	fumarase C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,683,560	1	.	C	97.4%	639.4 / 8.0	193	coding (1053/1404 nt)	fumC	fumarase C
Reads supporting (aligned to +/- strand): ref base . (4/1); new base C (95/93); total (99/94)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 230 total aligned reads displayed.

TTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑C‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCC  >  NC_000913/1683469‑1683647

tttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAgcg‑c                                                                                          >  2:554716/1‑93 (MQ=255)

tttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAgcg‑c                                                                                          >  2:2173892/1‑93 (MQ=255)

ttcccATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGTAAGACGTTCAGTTCAAAGTTACCGGAAgcg‑c                                                                                          <  1:2588959/90‑1 (MQ=255)

 ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑c‑c                                                                                        <  1:826452/93‑1 (MQ=255)

   ccATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑c‑c‑cc                                                                                     <  1:4880374/93‑1 (MQ=255)

    cATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑c‑c‑ccc                                                                                    <  1:1343429/93‑1 (MQ=255)

    cATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑c‑c‑ccc                                                                                    >  1:5340915/1‑93 (MQ=255)

     aTGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑c‑c‑cccc                                                                                   >  2:5028658/1‑93 (MQ=255)

     aTGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑c‑c‑cccc                                                                                   >  1:3641328/1‑93 (MQ=255)

     aTGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑c‑c‑cccc                                                                                   >  1:6359782/1‑93 (MQ=255)

     aTGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑c‑c‑cccc                                                                                   >  1:921576/1‑93 (MQ=255)

     aTGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑c‑c‑cccc                                                                                   >  2:3248907/1‑93 (MQ=255)

     aTGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑c‑c‑cccc                                                                                   >  1:1229619/1‑93 (MQ=255)

      tGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑c‑c‑ccccc                                                                                  >  1:1053338/1‑92 (MQ=255)

           tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGtt                                                                              >  2:1405986/1‑93 (MQ=255)

           tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑C‑CACCCCATGtt                                                                              >  1:6389547/1‑93 (MQ=255)

            cTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTg                                                                             <  2:3104999/93‑1 (MQ=255)

              gCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGAt                                                                           <  1:473561/93‑1 (MQ=255)

              gCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGAt                                                                           <  1:5696853/93‑1 (MQ=255)

              gCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGAt                                                                           >  2:5805181/1‑93 (MQ=255)

               ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATc                                                                          >  2:840473/1‑93 (MQ=255)

               ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATc                                                                          <  2:2736767/93‑1 (MQ=255)

                cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCg                                                                         >  2:6052198/1‑93 (MQ=255)

                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGcc                                                                       <  2:4228993/93‑1 (MQ=255)

                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGcc                                                                       >  1:2808859/1‑93 (MQ=255)

                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGcc                                                                       <  2:2233528/93‑1 (MQ=255)

                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGcc                                                                       >  1:1881912/1‑93 (MQ=255)

                     aGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCAcg                                                                    >  1:3510777/1‑93 (MQ=255)

                      gcgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCAcgt                                                                   <  2:661157/93‑1 (MQ=255)

                       cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCAcgtc                                                                  <  2:4181238/93‑1 (MQ=255)

                       cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCAcgtc                                                                  <  2:5903100/93‑1 (MQ=255)

                        gcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCAcgtcg                                                                 >  2:2710194/1‑93 (MQ=255)

                        gcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAAACGTTCAGTTCAAAGTTACCGGAAGCG‑CGC‑CCCCATGTTGATCGCCAcgtcg                                                                 >  1:5249863/1‑93 (MQ=255)

                         cACCGATTGCTGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCAcgtcgt                                                                >  1:194312/1‑93 (MQ=255)

                         cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑C‑CACCCCATGTTGATAGCTAcgtcgt                                                                >  1:3751061/1‑93 (MQ=255)

                          aCCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGtt                                                               >  1:665459/1‑93 (MQ=255)

                             gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTccc                                                            <  2:6195813/93‑1 (MQ=255)

                               ttGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCa                                                          >  2:3387584/1‑93 (MQ=255)

                                   aGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCAc                                                      <  2:2508201/93‑1 (MQ=255)

                                   aGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCAc                                                      >  1:3197871/1‑93 (MQ=255)

                                   aGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCAc                                                      >  1:2971656/1‑93 (MQ=255)

                                     gAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                    >  2:3036886/1‑93 (MQ=255)

                                     gAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                    <  1:2369943/93‑1 (MQ=255)

                                      aaaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTg                                                   >  2:2927445/1‑93 (MQ=255)

                                      aaaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTg                                                   <  2:120539/93‑1 (MQ=255)

                                      aaaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTg                                                   >  2:3454853/1‑93 (MQ=255)

                                        aTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAc                                                 <  2:2389244/93‑1 (MQ=255)

                                        aTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAc                                                 <  2:3793212/93‑1 (MQ=255)

                                        aTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAc                                                 >  1:1538463/1‑93 (MQ=255)

                                        aTTGTGGATCACCAGTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAc                                                 <  2:790371/93‑1 (MQ=255)

                                            tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagca                                             <  1:1161315/93‑1 (MQ=255)

                                             ggATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagcag                                            <  2:4575450/93‑1 (MQ=255)

                                              gATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTCCCGGACCCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCaga                                           <  1:2522780/93‑1 (MQ=255)

                                              gATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCaga                                           >  1:2722870/1‑93 (MQ=255)

                                               ttCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCCCCTGACAGCagag                                          <  1:4506223/92‑1 (MQ=255)

                                                 cACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                        <  2:5250060/93‑1 (MQ=255)

                                                 cACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                        >  2:3346588/1‑93 (MQ=255)

                                                  aCCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                       >  2:2532833/1‑93 (MQ=255)

                                                   ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATg                                      >  2:3947229/1‑93 (MQ=255)

                                                   ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATg                                      >  1:2029213/1‑93 (MQ=255)

                                                   ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATg                                      <  2:6344575/93‑1 (MQ=255)

                                                    cATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACCGCAGAGCATgg                                     >  1:4442414/1‑93 (MQ=255)

                                                    cATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATgg                                     <  1:397692/93‑1 (MQ=255)

                                                    cATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATgg                                     >  1:4284966/1‑93 (MQ=255)

                                                     aTTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                    <  2:5238431/93‑1 (MQ=255)

                                                      ttGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGtt                                   <  1:5813557/93‑1 (MQ=255)

                                                       tGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                  <  1:3582129/93‑1 (MQ=255)

                                                        ggACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTaa                                 <  1:4514061/93‑1 (MQ=255)

                                                        ggACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTaa                                 <  2:4351069/93‑1 (MQ=255)

                                                         gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCGTCACCTGACAGCAGAGCATGGTTAAt                                >  2:3970563/1‑93 (MQ=255)

                                                         gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAAt                                >  2:897346/1‑93 (MQ=255)

                                                         gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAAt                                <  2:4865428/93‑1 (MQ=255)

                                                         gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAAt                                <  1:4102841/93‑1 (MQ=255)

                                                            ggAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                             >  2:4095276/1‑93 (MQ=255)

                                                            ggAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                             >  2:5008316/1‑93 (MQ=255)

                                                             gAAGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                            <  2:1725424/93‑1 (MQ=255)

                                                              aaGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                           >  2:137642/1‑93 (MQ=255)

                                                              aaGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                           <  1:80621/93‑1 (MQ=255)

                                                               aGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTca                          >  1:3091886/1‑93 (MQ=255)

                                                               aGACGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTca                          >  2:5619903/1‑93 (MQ=255)

                                                                  cGTTCAGTTCAAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCacac                       <  1:6216369/93‑1 (MQ=255)

                                                                   gttcagttcaAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACt                      >  2:153393/1‑93 (MQ=255)

                                                                   gttcagttcaAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACt                      >  2:5445136/1‑93 (MQ=255)

                                                                    ttcagttcaAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACtg                     <  1:5084411/93‑1 (MQ=255)

                                                                     tcagttcaAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGACTCACACtgt                    >  2:3768408/1‑93 (MQ=255)

                                                                         ttcaAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTg                >  2:5113919/1‑93 (MQ=255)

                                                                         ttcaAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTg                <  1:2209788/93‑1 (MQ=255)

                                                                          tcaAAGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg               >  1:5965763/1‑93 (MQ=255)

                                                                            aaaGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGt             <  2:1228129/93‑1 (MQ=255)

                                                                            aaaGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGt             >  1:5166684/1‑93 (MQ=255)

                                                                             aaGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGtt            <  2:2591279/93‑1 (MQ=255)

                                                                              aGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTc           >  2:1074718/1‑93 (MQ=255)

                                                                              aGTTACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTc           >  2:2077639/1‑93 (MQ=255)

                                                                                ttACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCAc         <  1:1192247/93‑1 (MQ=255)

                                                                                 tACCGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACt        >  1:5924553/1‑93 (MQ=255)

                                                                                   ccGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACttt      >  2:4669047/1‑93 (MQ=255)

                                                                                    cGGAAGCTCC‑C‑CCCCATGTTGATCGCCACCTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc     <  1:1779859/93‑1 (MQ=255)

                                                                                    cGGAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc     >  2:5302329/1‑93 (MQ=255)

                                                                                     ggAAGCGCC‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTcc    <  2:5710159/93‑1 (MQ=255)

                                                                                      gAAGCG‑C‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTcccc  >  1:3791973/1‑93 (MQ=255)

TTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑C‑C‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCC  >  NC_000913/1683469‑1683647

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