BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,146,393	A→G	R484R (CGT→CGC)	yijP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,146,393	0	A	G	92.1%	1019.3 / 59.4	379	R484R (CGT→CGC)	yijP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (20/10); new base G (168/181); total (188/191)
Fisher's exact test for biased strand distribution p-value = 5.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-01

Only 100 of 399 total aligned reads displayed.

TATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAA  >  NC_000913/4146301‑4146482

  tttACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGcg                                                                                         <  1:7249722/93‑1 (MQ=255)

   ttACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGgcgc                                                                                        >  2:9444804/1‑93 (MQ=255)

   ttACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGgcgc                                                                                        >  2:1346241/1‑93 (MQ=255)

   ttACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGgcgc                                                                                        >  2:3963500/1‑93 (MQ=255)

    tACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGcgcg                                                                                       <  2:7672069/93‑1 (MQ=255)

      cGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTc                                                                                     <  1:9831572/93‑1 (MQ=255)

      cGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTc                                                                                     <  2:8652682/93‑1 (MQ=255)

      cGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTc                                                                                     <  1:7837499/93‑1 (MQ=255)

      cGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTc                                                                                     <  1:7572779/93‑1 (MQ=255)

      cGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTc                                                                                     <  2:7271876/93‑1 (MQ=255)

        aTCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCgg                                                                                   >  2:6717064/1‑93 (MQ=255)

        aTCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCgg                                                                                   >  1:3685074/1‑93 (MQ=255)

        aTCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCgg                                                                                   <  2:6188221/93‑1 (MQ=255)

         tCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCggg                                                                                  >  2:841708/1‑93 (MQ=255)

          cAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGt                                                                                 <  2:2428240/93‑1 (MQ=255)

          cAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTAAATATGGCGCGTCGGGt                                                                                 >  2:2822474/1‑93 (MQ=255)

           aaCATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGtt                                                                                >  1:939682/1‑93 (MQ=255)

           aaCATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGtt                                                                                <  2:6118501/93‑1 (MQ=255)

           aaCATCCTGCGAGAAATAACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGtt                                                                                <  1:3166004/93‑1 (MQ=255)

            aCATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTg                                                                               <  2:5898245/93‑1 (MQ=255)

            aCATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTg                                                                               >  2:8585461/1‑93 (MQ=255)

              aTCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTc                                                                             >  1:1537934/1‑93 (MQ=255)

               tCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTcc                                                                            >  1:4690420/1‑93 (MQ=255)

               tCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTcc                                                                            <  1:2128295/93‑1 (MQ=255)

               tCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTcc                                                                            <  2:3802036/93‑1 (MQ=255)

                ccTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCt                                                                           <  1:6118467/93‑1 (MQ=255)

                 cTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGAGCGTCGGGTTGTCCTc                                                                          >  2:4433359/1‑93 (MQ=255)

                  tGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCa                                                                         <  2:1783400/93‑1 (MQ=255)

                  cgCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCa                                                                         <  1:7476063/92‑1 (MQ=255)

                     gagaAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTa                                                                      >  1:1901328/1‑93 (MQ=255)

                     gagaAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTa                                                                      >  2:2549919/1‑93 (MQ=255)

                          aTCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGccc                                                                 >  2:9391277/1‑93 (MQ=255)

                           tCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGcccc                                                                >  1:5161167/1‑93 (MQ=255)

                             acagggATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGCTGTCCTCATTATGCCCCTg                                                              >  2:797825/4‑93 (MQ=255)

                              cGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGa                                                             >  1:7097431/1‑93 (MQ=255)

                               ggggATGAGTCGCTTGCCATTGTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAg                                                            <  2:7043409/93‑1 (MQ=255)

                                 ggATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTATCCTCATTACGCCCCTGAggt                                                          >  1:3299957/1‑91 (MQ=255)

                                  gATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGttt                                                         <  1:7399193/93‑1 (MQ=255)

                                  gATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGttt                                                         <  2:2043224/93‑1 (MQ=255)

                                   aTGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGtttt                                                        <  2:3016471/93‑1 (MQ=255)

                                    tGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTtg                                                       >  1:5917667/1‑93 (MQ=255)

                                     gAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTtgt                                                      <  2:4380186/93‑1 (MQ=255)

                                      aGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTtgtg                                                     >  1:2418257/1‑93 (MQ=255)

                                      aGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTtgtg                                                     <  2:6825460/93‑1 (MQ=255)

                                         cGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTgcgg                                                  <  2:9632820/93‑1 (MQ=255)

                                          gCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTgcggc                                                 >  2:4803980/1‑93 (MQ=255)

                                          gCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTgcggc                                                 >  2:875371/1‑93 (MQ=255)

                                          gCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGGTTTGTgcggg                                                 >  2:262512/1‑92 (MQ=255)

                                            ttGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTgcggcgg                                               >  2:9173568/1‑93 (MQ=255)

                                              gCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTgcggcggcg                                             <  2:5320127/93‑1 (MQ=255)

                                              gCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTgcggcggcg                                             >  1:8564238/1‑93 (MQ=255)

                                               ccATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCgt                                            <  2:4516398/93‑1 (MQ=255)

                                               ccATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCgt                                            <  1:658207/93‑1 (MQ=255)

                                                cATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCgtg                                           <  1:1157740/93‑1 (MQ=255)

                                                   ttttCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCg                                        >  1:10017483/1‑93 (MQ=255)

                                                     ttCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTa                                      >  1:573067/1‑93 (MQ=255)

                                                      tCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTaa                                     >  2:7435973/1‑93 (MQ=255)

                                                      tCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTaa                                     <  1:6976649/93‑1 (MQ=255)

                                                      tCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTaa                                     >  1:2010668/1‑93 (MQ=255)

                                                      tCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTaa                                     >  1:2339125/1‑93 (MQ=255)

                                                       cTGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTaaa                                    <  1:8136601/93‑1 (MQ=255)

                                                        tGACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAAc                                   <  2:1891301/93‑1 (MQ=255)

                                                         gACGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAAcc                                  >  1:3099580/1‑93 (MQ=255)

                                                          aCGTCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACct                                 >  2:10906/1‑93 (MQ=255)

                                                             tCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTACACCTCtt                              >  1:8842052/1‑93 (MQ=255)

                                                             tCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCtt                              <  2:3271496/93‑1 (MQ=255)

                                                             tCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCtt                              <  1:6836708/93‑1 (MQ=255)

                                                             tCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCtt                              >  2:1520529/1‑93 (MQ=255)

                                                             tCCACAGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCtt                              >  2:7222432/1‑93 (MQ=255)

                                                               cacaGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCa                            <  1:2198708/93‑1 (MQ=255)

                                                                ccaggagGAACGGGATGGTGGACATATGGCGCGTGCGGTGGTCCTCAGTACGCCCCTGAGTTTTGTGCGGCGGCTTGTCGTAAACCTCTTCAc                           <  1:9679577/92‑1 (MQ=255)

                                                                acaGCAGGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCAc                           >  2:4151311/1‑93 (MQ=255)

                                                                   gcagGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGt                        <  2:4872587/93‑1 (MQ=255)

                                                                    cagGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTg                       <  2:8806170/93‑1 (MQ=255)

                                                                     agGAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTgg                      >  2:6633245/1‑93 (MQ=255)

                                                                      ggAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGt                     >  2:8933463/1‑93 (MQ=255)

                                                                      ggAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGt                     >  2:6592395/1‑93 (MQ=255)

                                                                      ggAACGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGt                     <  2:5668901/93‑1 (MQ=255)

                                                                          cGGAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCaga                 >  1:7725198/1‑93 (MQ=255)

                                                                            gAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAgaga               <  2:1179294/93‑1 (MQ=255)

                                                                            gAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAgaga               <  1:130489/93‑1 (MQ=255)

                                                                            gAATGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAgaga               >  2:8288014/1‑93 (MQ=255)

                                                                             aaTGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGaa              <  1:7963926/93‑1 (MQ=255)

                                                                              aTGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGaaa             <  1:451358/93‑1 (MQ=255)

                                                                              aTGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGaaa             <  2:215318/93‑1 (MQ=255)

                                                                               tGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAAt            >  2:353112/1‑93 (MQ=255)

                                                                               tGGTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTAACCGGGGTCAGCGaaaa            >  1:2718706/1‑92 (MQ=255)

                                                                                ggTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATa           <  1:4052394/93‑1 (MQ=255)

                                                                                ggTGTACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAAACTCTTCACCGTGGTCAGAGAAATa           >  1:4959273/1‑93 (MQ=255)

                                                                                 gtgtACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATaa          <  2:5061855/93‑1 (MQ=255)

                                                                                  tgtACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATaaa         >  1:172126/1‑93 (MQ=255)

                                                                                  tgtACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAAAAATaaa         >  1:4218964/1‑93 (MQ=255)

                                                                                   gtACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAAc        <  2:2257852/93‑1 (MQ=255)

                                                                                   gtACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAAc        <  1:8321387/93‑1 (MQ=255)

                                                                                    tACATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAAcc       <  1:5386907/93‑1 (MQ=255)

                                                                                    tACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAAcc       <  1:6396856/93‑1 (MQ=255)

                                                                                     aCATATGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCa      >  1:3836783/1‑93 (MQ=255)

                                                                                     aCATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCa      >  2:6420812/1‑93 (MQ=255)

                                                                                       atatGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAgg    <  1:9094862/93‑1 (MQ=255)

                                                                                        tatGGCGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGa   >  1:8996910/1‑93 (MQ=255)

TATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTGTACATATGACGCGTCGGGTTGTCCTCATTACGCCCCTGAGTTTTGTGCGGCGGCGTGTCGTAAACCTCTTCACCGTGGTCAGAGAAATAAACCAGGAA  >  NC_000913/4146301‑4146482

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