BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,892,505	T→C	G75G (GGT→GGC)	yieF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,892,505	0	T	C	99.6%	785.7 / ‑6.1	252	G75G (GGT→GGC)	yieF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (121/130); minor base G (0/1); total (121/131)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 271 total aligned reads displayed.

CATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACC  >  NC_000913/3892415‑3892594

cATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGgcgt                                                                                         >  1:2631336/1‑93 (MQ=255)

 aTTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGcgtg                                                                                        <  2:2300635/93‑1 (MQ=255)

  ttCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCgtgg                                                                                       <  2:2542501/93‑1 (MQ=255)

   tCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCgtggt                                                                                      <  2:2669521/93‑1 (MQ=255)

   tCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCgtggt                                                                                      >  2:2774460/1‑93 (MQ=255)

   tCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCgtggt                                                                                      <  2:3099420/93‑1 (MQ=255)

     ccTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGa                                                                                    <  1:246269/93‑1 (MQ=255)

     ccTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGa                                                                                    >  2:3691090/1‑93 (MQ=255)

     ccTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGa                                                                                    <  2:4446219/93‑1 (MQ=255)

     ccTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGa                                                                                    >  2:854390/1‑93 (MQ=255)

      cTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGAt                                                                                   >  2:3274059/1‑93 (MQ=255)

      cTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGAt                                                                                   >  1:3487746/1‑93 (MQ=255)

       ttGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATc                                                                                  >  1:1558329/1‑93 (MQ=255)

       ttGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATc                                                                                  >  2:3606994/1‑93 (MQ=255)

       ttGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATc                                                                                  <  2:2212734/93‑1 (MQ=255)

       ttGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATc                                                                                  >  2:2091793/1‑93 (MQ=255)

         gTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGt                                                                                <  2:1254912/93‑1 (MQ=255)

          tATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTc                                                                               <  2:1767195/93‑1 (MQ=255)

          tATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTc                                                                               >  2:1524902/1‑93 (MQ=255)

          tATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGAGTGGTGATCGTc                                                                               >  1:3254050/1‑93 (MQ=255)

            tGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCAc                                                                             <  1:1892233/93‑1 (MQ=255)

            tGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCAc                                                                             >  2:4287526/1‑93 (MQ=255)

             gACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACg                                                                            >  1:780928/1‑93 (MQ=255)

              aCGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGc                                                                           >  1:407415/1‑93 (MQ=255)

              aCGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGc                                                                           >  2:1359274/1‑93 (MQ=255)

                gatGACGTACACCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCg                                                                         <  2:3314655/91‑1 (MQ=255)

                 cTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCgg                                                                        <  1:1758143/93‑1 (MQ=255)

                  tGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGa                                                                       >  1:3671612/1‑93 (MQ=255)

                  tGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGa                                                                       <  1:3723868/93‑1 (MQ=255)

                   gACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGaa                                                                      <  1:3319315/93‑1 (MQ=255)

                   gACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGaa                                                                      <  1:1448711/93‑1 (MQ=255)

                   gACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGaa                                                                      >  1:4586961/1‑93 (MQ=255)

                    aCGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAat                                                                     <  1:975883/93‑1 (MQ=255)

                    aCGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAat                                                                     <  2:4329217/93‑1 (MQ=255)

                     cGTACAGCAGGAAGAAGGTTTTCCAGCAACGTTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAata                                                                    <  1:312568/93‑1 (MQ=255)

                      gTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAatat                                                                   <  1:2139406/93‑1 (MQ=255)

                        aCAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATaa                                                                 >  1:2583380/1‑93 (MQ=255)

                           gcagGAAGAAGGCTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGCGATCGGCATGCCGGACTATAactc                                                              >  2:2388703/1‑92 (MQ=255)

                            gagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAactac                                                             <  2:3403343/92‑1 (MQ=255)

                            cagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAactac                                                             <  2:2718050/93‑1 (MQ=255)

                             agGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAactact                                                            >  1:137786/1‑93 (MQ=255)

                              ggAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTc                                                           >  2:1315003/1‑93 (MQ=255)

                                aagaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCgg                                                         <  2:1086691/93‑1 (MQ=255)

                                 agaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGt                                                        >  2:665727/1‑93 (MQ=255)

                                  gaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTa                                                       >  2:3742229/1‑93 (MQ=255)

                                  gaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTa                                                       >  1:1104050/1‑93 (MQ=255)

                                    agGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTAcc                                                     <  1:13179/93‑1 (MQ=255)

                                     ggTTTTCAAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCg                                                    >  1:3977008/1‑93 (MQ=255)

                                      gTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCgg                                                   >  2:3836768/1‑93 (MQ=255)

                                      gTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCgg                                                   >  1:2972453/1‑93 (MQ=255)

                                         ttCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTg                                                <  2:2882449/93‑1 (MQ=255)

                                          tCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTgg                                               >  2:4191293/1‑93 (MQ=255)

                                           ccAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTggg                                              >  2:1046563/1‑93 (MQ=255)

                                           ccAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTggg                                              >  2:460792/1‑93 (MQ=255)

                                             aGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCt                                            >  1:1634319/1‑93 (MQ=255)

                                              gCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTg                                           >  2:1699844/1‑93 (MQ=255)

                                                aaCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATTGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaa                                         <  1:4312940/93‑1 (MQ=255)

                                                aaCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGGGGGCTGaa                                         >  2:3188727/1‑93 (MQ=255)

                                                 aCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaa                                        >  2:812926/1‑93 (MQ=255)

                                                 aCGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaa                                        >  1:1520923/1‑93 (MQ=255)

                                                  cGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaa                                       <  2:4496648/93‑1 (MQ=255)

                                                  cGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGAGGTGATCGGCACGCCGGCATATAACTACTCGGTACCGGGTGGCCTGaaaa                                       >  1:3073311/1‑93 (MQ=255)

                                                   gggtGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaaa                                      <  2:1337771/90‑1 (MQ=255)

                                                   ggTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaaa                                      <  2:1430506/93‑1 (MQ=255)

                                                   ggTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGaaaaa                                      <  2:1620616/93‑1 (MQ=255)

                                                    gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATTGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAAt                                     >  2:2971807/1‑93 (MQ=255)

                                                     ttGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATg                                    >  1:1419506/1‑93 (MQ=255)

                                                     ttGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATg                                    <  2:1942507/93‑1 (MQ=255)

                                                       gAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGcc                                  >  1:1748237/1‑93 (MQ=255)

                                                       gAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGcc                                  >  1:1693692/1‑93 (MQ=255)

                                                       gAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGcc                                  >  2:3024895/1‑93 (MQ=255)

                                                          gCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATc                               >  2:1804048/1‑93 (MQ=255)

                                                          gCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATc                               >  2:3563508/1‑93 (MQ=255)

                                                          gCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATc                               >  2:4137483/1‑93 (MQ=255)

                                                          gCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATc                               <  2:486425/93‑1 (MQ=255)

                                                           ctctGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCg                              >  1:4249598/1‑93 (MQ=255)

                                                            tctGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGa                             <  1:2159632/93‑1 (MQ=255)

                                                            tctGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGa                             >  2:1993436/1‑93 (MQ=255)

                                                               ggCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTg                          >  1:4624849/1‑93 (MQ=255)

                                                                gCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTgg                         >  2:3054339/1‑93 (MQ=255)

                                                                 cGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGc                        >  2:4040687/1‑93 (MQ=255)

                                                                   gAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCtt                      >  1:2345956/1‑93 (MQ=255)

                                                                    aaCAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCttt                     <  2:3625196/93‑1 (MQ=255)

                                                                    aaCAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTCGCttt                     >  1:334835/1‑93 (MQ=255)

                                                                     aCAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTc                    >  1:684330/1‑93 (MQ=255)

                                                                     aCAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTc                    >  2:2913251/1‑93 (MQ=255)

                                                                      cAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTcc                   >  1:615984/1‑93 (MQ=255)

                                                                         aTCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGc                >  1:2331819/1‑93 (MQ=255)

                                                                          tCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGcc               >  1:4155631/1‑93 (MQ=255)

                                                                           ccGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCt              <  1:3574969/93‑1 (MQ=255)

                                                                             gTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGc            >  1:2875841/1‑93 (MQ=255)

                                                                              tCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGcc           <  2:61037/93‑1 (MQ=255)

                                                                                  gCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGAt       <  1:2061349/93‑1 (MQ=255)

                                                                                   cTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATc      >  2:4008073/1‑93 (MQ=255)

                                                                                    tGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCa     >  1:1040200/1‑93 (MQ=255)

                                                                                    tGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCa     >  2:183500/1‑93 (MQ=255)

                                                                                     gACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCaa    <  2:3297366/93‑1 (MQ=255)

                                                                                      gcGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc   <  1:4432614/92‑1 (MQ=255)

                                                                                      aCGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc   <  2:4530551/93‑1 (MQ=255)

                                                                                       cGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAcc  <  1:1990629/93‑1 (MQ=255)

CATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACC  >  NC_000913/3892415‑3892594

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