BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,724,929	T→C	E235G (GAG→GGG)	ydhF ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,724,929	0	T	C	89.8%	863.4 / 72.3	264	E235G (GAG→GGG)	ydhF	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (9/18); new base C (135/102); total (144/120)
Fisher's exact test for biased strand distribution p-value = 2.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 278 total aligned reads displayed.

GCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCATCGGACGAACACG  >  NC_000913/1724831‑1725026

gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCctct                                                                                                 >  1:3116461/1‑101 (MQ=255)

gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCccct                                                                                                 >  2:3116457/1‑101 (MQ=255)

gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCccct                                                                                                 <  1:3115977/101‑1 (MQ=255)

gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCccct                                                                                                 >  1:3116455/1‑101 (MQ=255)

 cTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTg                                                                                                <  2:3116191/101‑1 (MQ=255)

 cTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTg                                                                                                >  2:3116463/1‑101 (MQ=255)

 cTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTg                                                                                                <  1:3116029/101‑1 (MQ=255)

   cAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGcc                                                                                              >  1:3116466/1‑101 (MQ=255)

    aaTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGcca                                                                                             >  1:3116468/1‑101 (MQ=255)

     aTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCCCCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGccac                                                                                            <  2:3116130/101‑1 (MQ=255)

     aTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCCCCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGccac                                                                                            <  2:3115991/101‑1 (MQ=255)

     aTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGccac                                                                                            >  2:3116470/1‑101 (MQ=255)

      ttttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGccacc                                                                                           <  1:3116041/101‑1 (MQ=255)

        ttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCcaccac                                                                                         >  1:3116475/1‑101 (MQ=255)

          aCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAg                                                                                       >  1:3116478/1‑101 (MQ=255)

           ccTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGc                                                                                      >  1:3116484/1‑101 (MQ=255)

           ccTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGc                                                                                      >  1:3116482/1‑101 (MQ=255)

            cTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGcc                                                                                     >  2:3116485/1‑101 (MQ=255)

            cTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGcc                                                                                     <  2:3115938/101‑1 (MQ=255)

             tGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCa                                                                                    >  2:3116488/1‑101 (MQ=255)

             tGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCa                                                                                    <  1:3116095/101‑1 (MQ=255)

               aaCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGt                                                                                  <  1:3116102/101‑1 (MQ=255)

               aaCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGt                                                                                  <  2:3115649/101‑1 (MQ=255)

                aCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGtt                                                                                 >  1:3116490/1‑101 (MQ=255)

                 ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTtc                                                                                <  1:3116024/101‑1 (MQ=255)

                   gATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTtcat                                                                              >  2:3116495/1‑101 (MQ=255)

                   gATAATTGGCAGCAGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTtcat                                                                              >  1:3116496/1‑101 (MQ=255)

                     tAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTtcatca                                                                            >  2:3116499/1‑101 (MQ=255)

                     tAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTtcatca                                                                            >  2:3116501/1‑101 (MQ=255)

                      aaTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCAc                                                                           <  2:3116251/101‑1 (MQ=255)

                      aaTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCAc                                                                           <  1:3116228/101‑1 (MQ=255)

                        ttGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGc                                                                         <  1:3116106/101‑1 (MQ=255)

                        ttGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGc                                                                         >  1:3116506/1‑101 (MQ=255)

                         tGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCa                                                                        >  2:3116509/1‑101 (MQ=255)

                          ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAg                                                                       >  1:3116510/1‑101 (MQ=255)

                           gCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGc                                                                      <  2:3116046/101‑1 (MQ=255)

                           gCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGc                                                                      <  1:3116116/101‑1 (MQ=255)

                            cAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCg                                                                     >  1:3116512/1‑101 (MQ=255)

                               cGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCt                                                                  >  2:3116517/1‑101 (MQ=255)

                                ggCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                                                                 <  2:3116109/101‑1 (MQ=255)

                                ggCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                                                                 >  1:3116522/1‑101 (MQ=255)

                                ggCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCAAAGCCAGTTCATCACGCAGCGGCTg                                                                 >  2:3116520/1‑101 (MQ=255)

                                   tGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaa                                                              >  2:3116526/1‑101 (MQ=255)

                                    gCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGaaa                                                             >  1:3116528/1‑101 (MQ=255)

                                     cGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg                                                                 <  1:3115761/96‑1 (MQ=255)

                                     cGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaat                                                            >  2:3116530/1‑101 (MQ=255)

                                      gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaata                                                           >  1:3116535/1‑101 (MQ=255)

                                       aTGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGc                                                                   <  1:3116069/92‑1 (MQ=255)

                                       aTGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAaataa                                                          >  2:3116537/1‑101 (MQ=255)

                                          gTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatc                                                        <  2:3116089/100‑1 (MQ=255)

                                           tAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcat                                                      >  2:3116540/1‑101 (MQ=255)

                                            aaaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAatcatc                                                     >  1:3116541/1‑101 (MQ=255)

                                               cGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCAtt                                                  <  1:3116354/101‑1 (MQ=255)

                                                gTAATCCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTa                                                 <  2:3115961/101‑1 (MQ=255)

                                                gTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTa                                                 <  1:3116087/101‑1 (MQ=255)

                                                  aaTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa                                               <  1:3116135/101‑1 (MQ=255)

                                                  aaTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa                                               >  1:3116549/1‑101 (MQ=255)

                                                  aaTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa                                               >  1:3116546/1‑101 (MQ=255)

                                                  aaTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa                                               <  2:3116030/101‑1 (MQ=255)

                                                    tACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACa                                             >  2:3116552/1‑101 (MQ=255)

                                                      cccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTaaa                                               <  1:3115968/97‑1 (MQ=255)

                                                      cccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAga                                           >  1:3116554/1‑101 (MQ=255)

                                                       ccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAgac                                          >  2:3116555/1‑101 (MQ=255)

                                                       ccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAgac                                          <  2:3116304/101‑1 (MQ=255)

                                                         aGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAga                                           >  2:3116558/1‑98 (MQ=255)

                                                          ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAgacgac                                       <  2:3116002/101‑1 (MQ=255)

                                                          ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAgacgac                                       <  2:3116162/101‑1 (MQ=255)

                                                             gTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACa                                             >  2:3116565/1‑92 (MQ=255)

                                                             gTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGaccac                                    >  2:3116562/1‑101 (MQ=255)

                                                               aaaCCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGaccacca                                  >  1:3116567/1‑101 (MQ=255)

                                                                aaCCACCTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGaccaccac                                 >  1:3116570/1‑101 (MQ=255)

                                                                 accaccTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGaccaccacc                                <  1:3116049/101‑1 (MQ=255)

                                                                 accaccTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGaccaccacc                                <  1:3116311/101‑1 (MQ=255)

                                                                  ccaccTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATTATCATTAAACAGACGAccaccacca                               >  2:3116571/1‑101 (MQ=255)

                                                                      cTGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGc                           >  1:3116575/1‑101 (MQ=255)

                                                                       tGTTCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCaa                              >  2:3116578/1‑97 (MQ=255)

                                                                         ttCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAgg                        <  1:3116257/101‑1 (MQ=255)

                                                                          tCAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGa                       >  2:3116582/1‑101 (MQ=255)

                                                                           cAATCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGAc                      >  1:3116584/1‑101 (MQ=255)

                                                                            aaTCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAg                             >  2:3116589/1‑93 (MQ=255)

                                                                            aaTCGAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGAcc                     >  1:3116588/1‑101 (MQ=255)

                                                                                gAGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGc                 >  2:3116595/1‑101 (MQ=255)

                                                                                gAGCCCGCGTTTAACTCCCATGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGc                 >  2:3116593/1‑101 (MQ=255)

                                                                                 aGCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGcc                <  2:3116018/101‑1 (MQ=255)

                                                                                  gCCCGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAgacgac                                       <  1:3116149/77‑1 (MQ=255)

                                                                                   cccGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCAt              <  2:3116237/101‑1 (MQ=255)

                                                                                    ccGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCa                    <  1:3116154/94‑1 (MQ=255)

                                                                                    ccGCGTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCATc             <  2:3116078/101‑1 (MQ=255)

                                                                                     cgcgTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCATCg            >  2:3116599/1‑101 (MQ=255)

                                                                                     cgcgTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCATCg            >  1:3116600/1‑101 (MQ=255)

                                                                                      gcgTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAg                             <  1:3115995/83‑1 (MQ=255)

                                                                                      gcgTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCa                                                                        <  2:3115955/40‑1 (MQ=38)

                                                                                      gcgTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCATCgg           <  1:3116200/101‑1 (MQ=255)

                                                                                      gcgTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCATCgg           <  1:3115894/101‑1 (MQ=255)

                                                                                        gTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCATCGGAc         <  1:3115940/101‑1 (MQ=255)

                                                                                        gTTTAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCa               >  1:3116604/1‑95 (MQ=255)

                                                                                         tttAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCATCGGACg        <  2:3116120/101‑1 (MQ=255)

                                                                                          ttAACTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCATCGGACGa       >  2:3116607/1‑101 (MQ=255)

                                                                                             aCTCCCCTGCCACCACAGCCAGTTCATCACGCAGCGGCTgg                                                                <  1:3115678/41‑1 (MQ=255)

                                                                                              ctccccTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGc                 >  2:3116610/1‑87 (MQ=255)

GCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGGAAATAATCATCATTAAACAGACGACCACCACCAAGGCAGGACCACGCCATCGGACGAACACG  >  NC_000913/1724831‑1725026

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