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.8%	1295.0 / ‑6.3	415	G75G (GGT→GGC)	yieF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (193/221); total (194/221)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 447 total aligned reads displayed.

CATTCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGTGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACC  >  NC_000913/3892415‑3892594

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

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

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

  ttCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCgtgg                                                                                       >  2:4735894/1‑93 (MQ=255)

  ttCCCTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGAGGAACAGATCCGTCATGCTGACGGTgtgg                                                                                       >  1:481725/1‑93 (MQ=255)

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

      cTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGAt                                                                                   <  1:3826042/93‑1 (MQ=255)

      cTTGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGAt                                                                                   <  1:2353408/93‑1 (MQ=255)

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

        tGTATGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCg                                                                                 <  2:1172765/93‑1 (MQ=255)

           aTGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCa                                                                              >  2:7931503/1‑93 (MQ=255)

           aTGACGCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCa                                                                              <  2:2586740/93‑1 (MQ=255)

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

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

                gCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCg                                                                         <  2:2509378/93‑1 (MQ=255)

                gCTGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCg                                                                         <  2:4348110/93‑1 (MQ=255)

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

                  tGACGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTCGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGa                                                                       >  1:5974195/1‑93 (MQ=255)

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

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

                     cGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAata                                                                    <  1:7743220/93‑1 (MQ=255)

                     cGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAata                                                                    <  2:4141046/93‑1 (MQ=255)

                     cGTACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAata                                                                    <  1:1069696/93‑1 (MQ=255)

                       tACAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAAtata                                                                  >  2:6941860/1‑93 (MQ=255)

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

                        aCAGCAGGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATaa                                                                 <  2:1853634/93‑1 (MQ=255)

                         cagcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGACTATAac                                                                >  1:685485/1‑93 (MQ=255)

                         cagcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAac                                                                >  2:2686904/1‑93 (MQ=255)

                          agcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAact                                                               >  2:2453617/1‑93 (MQ=255)

                           gcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAactc                                                              >  1:6282625/1‑92 (MQ=255)

                           gcagGAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAacta                                                              >  2:3398488/1‑93 (MQ=255)

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

                              ggAAGAAGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTc                                                           <  2:1498072/93‑1 (MQ=255)

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

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

                               gaagaaGGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAAATACTCg                                                          >  2:1377021/1‑93 (MQ=255)

                                 agaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGCCGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGt                                                        >  2:2090192/1‑93 (MQ=255)

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

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

                                  gaagGTTTTCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTa                                                       <  1:24195/93‑1 (MQ=255)

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

                                       ttttCCAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGCACCggg                                                  >  1:477939/1‑93 (MQ=255)

                                        ttgccAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGt                                                 <  1:6463405/90‑1 (MQ=255)

                                           ccAGCAACGGTTGAAGCTCTGGCTGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTggg                                              <  2:4290018/93‑1 (MQ=255)

                                            cAGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGc                                             <  2:3121438/93‑1 (MQ=255)

                                            cAGCAACGGTTGAAGCTCTGGCGGAACAGATCCATCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGc                                             <  1:7160176/93‑1 (MQ=255)

                                             aGCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGCGATCGTCACGCCGGAATATAACTACTTGGTACCGGGTGGGCt                                            >  2:171362/1‑93 (MQ=255)

                                              gCAACGGTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTg                                           <  1:5121022/93‑1 (MQ=255)

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

                                                    gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAAt                                     <  1:209607/93‑1 (MQ=255)

                                                    gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAAt                                     >  2:3888593/1‑93 (MQ=255)

                                                    gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAAt                                     >  1:2888429/1‑93 (MQ=255)

                                                    gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAAt                                     >  2:7481249/1‑93 (MQ=255)

                                                    gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAAt                                     >  1:4093586/1‑93 (MQ=255)

                                                    gTTGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAAt                                     >  1:3935906/1‑93 (MQ=255)

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

                                                      tGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGc                                   <  1:4186864/93‑1 (MQ=255)

                                                      tGAAGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGc                                   >  2:7339541/1‑93 (MQ=255)

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

                                                        aaGCTCTGTGGGAACAGATTCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCa                                 <  1:5886893/93‑1 (MQ=255)

                                                        aaGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCa                                 <  2:4008/93‑1 (MQ=255)

                                                        aaGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCa                                 <  1:2524829/93‑1 (MQ=255)

                                                        aaGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCa                                 >  2:2627252/1‑93 (MQ=255)

                                                         aGCTCTGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCAt                                <  2:6113429/93‑1 (MQ=255)

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

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

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

                                                          gCTCTGGCAGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGAAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATc                               <  2:6674747/93‑1 (MQ=255)

                                                           ctctGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCg                              <  1:7128378/93‑1 (MQ=255)

                                                           ctctGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCg                              <  1:837875/93‑1 (MQ=255)

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

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

                                                            cctGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGa                             <  1:2621849/92‑1 (MQ=255)

                                                             ctGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGAc                            <  2:7831187/93‑1 (MQ=255)

                                                              tGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACt                           <  2:3258859/93‑1 (MQ=255)

                                                              tGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACt                           <  1:1926824/93‑1 (MQ=255)

                                                              tGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACt                           >  1:7528466/1‑93 (MQ=255)

                                                              tGGCGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACt                           <  1:4608984/93‑1 (MQ=255)

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

                                                                 cGGAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGc                        <  1:1407000/93‑1 (MQ=255)

                                                                  ggAACAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCt                       <  2:5490482/93‑1 (MQ=255)

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

                                                                    aaCAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCttt                     >  2:7604939/1‑93 (MQ=255)

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

                                                                      cAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTcc                   <  1:6359045/93‑1 (MQ=255)

                                                                      cAGATCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTcc                   <  2:6277773/93‑1 (MQ=255)

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

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

                                                                          tCCGTCAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGcc               <  2:349822/93‑1 (MQ=255)

                                                                               cAGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCg          <  1:1632358/93‑1 (MQ=255)

                                                                                aGGCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCgg         <  2:22836/93‑1 (MQ=255)

                                                                                 ggCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGa        >  2:373896/1‑93 (MQ=255)

                                                                                  gCTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGAt       >  1:4918605/1‑93 (MQ=255)

                                                                                   cTGACGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATc      <  2:5771394/93‑1 (MQ=255)

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

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

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

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

                                                                                      aCGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAc   >  1:5486523/1‑93 (MQ=255)

                                                                                       cGGCGTGGTGATCGTCACGCCGGAATATAACTACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAAcc  >  1:1529792/1‑93 (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` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`