BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,725,180	G→A	I70I (ATC→ATT)	yiaB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,725,180	0	G	A	99.5%	1573.8 / ‑4.8	436	I70I (ATC→ATT)	yiaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (190/244); total (191/245)
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.76e-01

Only 100 of 444 total aligned reads displayed.

GGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAG  >  NC_000913/3725090‑3725269

  gATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAatat                                                                                        >  1:4073401/1‑93 (MQ=255)

    tAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATgg                                                                                      <  1:1709191/93‑1 (MQ=255)

    tAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATgg                                                                                      >  2:585147/1‑93 (MQ=255)

    tAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATgg                                                                                      >  1:2568317/1‑93 (MQ=255)

      gATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATggga                                                                                    >  2:5645485/1‑92 (MQ=255)

      gATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGt                                                                                    >  2:1091557/1‑93 (MQ=255)

        tCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaa                                                                                  >  2:2456175/1‑93 (MQ=255)

        tCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaa                                                                                  >  1:3612418/1‑93 (MQ=255)

         cACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaa                                                                                 <  2:4157334/93‑1 (MQ=255)

          aCCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaa                                                                                >  2:2827040/1‑93 (MQ=255)

          aCCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaa                                                                                >  2:6902103/1‑93 (MQ=255)

            cATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaa                                                                              >  2:4887976/1‑93 (MQ=255)

            cATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaa                                                                              >  2:4737396/1‑93 (MQ=255)

             aTTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaaa                                                                             <  1:6140932/93‑1 (MQ=255)

             aTTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaaa                                                                             <  2:6842330/93‑1 (MQ=255)

             aTTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaaa                                                                             <  2:4016621/93‑1 (MQ=255)

             aTTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTaaaaaaa                                                                             <  1:4325921/93‑1 (MQ=255)

              tttCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAat                                                                            <  1:7720305/93‑1 (MQ=255)

              tttCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAat                                                                            >  1:4889103/1‑93 (MQ=255)

                tCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatca                                                                          <  2:4643409/93‑1 (MQ=255)

                tCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatca                                                                          <  2:760658/93‑1 (MQ=255)

                tCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatca                                                                          <  1:7510642/93‑1 (MQ=255)

                  aTAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAatcatc                                                                        >  1:5072483/1‑93 (MQ=255)

                   tAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                       >  1:7068532/1‑93 (MQ=255)

                   tAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATcc                                                                       >  2:32979/1‑93 (MQ=255)

                    aaGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCa                                                                      <  1:6266183/93‑1 (MQ=255)

                    aaGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCa                                                                      >  1:3757057/1‑93 (MQ=255)

                    aaGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCa                                                                      <  2:1167185/93‑1 (MQ=255)

                      gTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGt                                                                    <  2:770472/93‑1 (MQ=255)

                         ttGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGc                                                                 >  2:3566715/1‑93 (MQ=255)

                            aTAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCg                                                              >  1:943738/1‑93 (MQ=255)

                            aTAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCg                                                              <  2:6614990/93‑1 (MQ=255)

                             tAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGa                                                             >  2:1035377/1‑93 (MQ=255)

                               ggTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑t                                                          >  2:6129635/1‑93 (MQ=255)

                               ggTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑t                                                          >  1:6356939/1‑93 (MQ=255)

                               ggTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑t                                                          >  1:1333002/1‑93 (MQ=255)

                                 tGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑ttt                                                        >  2:1831805/1‑93 (MQ=255)

                                 tGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑ttt                                                        >  1:6688865/1‑93 (MQ=255)

                                  gCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑tttt                                                       <  2:7434253/93‑1 (MQ=255)

                                  gCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑tttt                                                       <  1:7416270/93‑1 (MQ=255)

                                   cGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑ttttt                                                      <  2:451568/93‑1 (MQ=255)

                                    gTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTc                                                     >  2:2090368/1‑93 (MQ=255)

                                    gTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTc                                                     <  2:5181500/93‑1 (MQ=255)

                                    gTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTc                                                     >  1:83202/1‑93 (MQ=255)

                                    gTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTc                                                     <  1:1255012/93‑1 (MQ=255)

                                    gTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTc                                                     <  1:5684017/93‑1 (MQ=255)

                                     tttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCt                                                    >  1:3796494/1‑93 (MQ=255)

                                     tttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCt                                                    >  2:1258697/1‑93 (MQ=255)

                                     tttAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCt                                                    <  1:2346580/93‑1 (MQ=255)

                                         aaaCACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTgcgc                                                >  2:4091177/1‑93 (MQ=255)

                                            cacCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGaa                                             >  2:3144206/1‑93 (MQ=255)

                                                gATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGata                                         <  1:1667936/93‑1 (MQ=255)

                                                gATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGata                                         <  1:5912988/93‑1 (MQ=255)

                                                gATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGata                                         <  1:3994201/93‑1 (MQ=255)

                                                 aTAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGatat                                        <  2:7297294/93‑1 (MQ=255)

                                                  tAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATg                                       >  1:284689/1‑93 (MQ=255)

                                                  tAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATg                                       <  2:1542781/93‑1 (MQ=255)

                                                  tAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATg                                       <  1:4758604/93‑1 (MQ=255)

                                                   aaaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGCG‑TTTTTCTGCGCGAGGCTATga                                      >  1:6528661/1‑92 (MQ=255)

                                                   aaaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGc                                      <  2:4203695/93‑1 (MQ=255)

                                                    aaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCa                                     <  1:7120199/93‑1 (MQ=255)

                                                    aaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCa                                     <  2:7795507/93‑1 (MQ=255)

                                                    aaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCa                                     >  1:7818390/1‑93 (MQ=255)

                                                    aaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCa                                     <  1:538815/93‑1 (MQ=255)

                                                    aaCAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCa                                     >  2:1362572/1‑93 (MQ=255)

                                                      cAAGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATa                                   <  2:5062258/93‑1 (MQ=255)

                                                       aaGAGACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAg                                  >  1:1971452/1‑93 (MQ=255)

                                                        agagACCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGt                                 <  2:427768/93‑1 (MQ=255)

                                                            aCCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGcc                             >  2:7061263/1‑93 (MQ=255)

                                                            aCCGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGcc                             <  1:3198100/93‑1 (MQ=255)

                                                              cGATAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCaa                           >  1:1496193/1‑93 (MQ=255)

                                                                aTAGTGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCaaaa                         >  1:3071743/1‑93 (MQ=255)

                                                                    tGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTg                     >  1:5244354/1‑93 (MQ=255)

                                                                    tGATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTg                     <  2:2249304/93‑1 (MQ=255)

                                                                     gATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGc                    <  1:519567/93‑1 (MQ=255)

                                                                     gATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGc                    <  2:6376823/93‑1 (MQ=255)

                                                                     gATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGc                    >  1:2802324/1‑93 (MQ=255)

                                                                     gATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGc                    <  1:4605258/93‑1 (MQ=255)

                                                                     gATTAACGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGc                    <  2:1956845/93‑1 (MQ=255)

                                                                        tAACGCAACTAACTGGCAAATATGGATAAAAAAAACATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGt                 <  1:3913853/93‑1 (MQ=255)

                                                                          aCGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCa               <  1:1070879/93‑1 (MQ=255)

                                                                           cGCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAt              <  2:3082897/93‑1 (MQ=255)

                                                                           agCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAt              <  2:4584678/92‑1 (MQ=255)

                                                                            gCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAGTTTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAt              <  1:5591203/93‑1 (MQ=255)

                                                                            gCAACTAACTGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCAtt             >  2:5455065/1‑93 (MQ=255)

                                                                             cAACTAACTGGCAAAGATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTa            <  2:2758322/93‑1 (MQ=255)

                                                                              aactaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTaa           >  2:433198/1‑93 (MQ=255)

                                                                              aactaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTaa           >  1:1569851/1‑93 (MQ=255)

                                                                               actaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAac          <  2:6286208/93‑1 (MQ=255)

                                                                               actaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAac          >  1:208410/1‑93 (MQ=255)

                                                                                ctaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAaca         <  1:5485014/93‑1 (MQ=255)

                                                                                ctaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAaca         >  2:7225489/1‑93 (MQ=255)

                                                                                ctaactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAaca         <  2:2990989/93‑1 (MQ=255)

                                                                                 taactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAacac        <  1:2243819/93‑1 (MQ=255)

                                                                                 taactGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAacac        >  2:2597262/1‑93 (MQ=255)

                                                                                   actGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGc      <  2:5731265/93‑1 (MQ=255)

                                                                                   actGGCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGc      <  1:5800038/93‑1 (MQ=255)

                                                                                      ggCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCaa   <  2:6689669/93‑1 (MQ=255)

                                                                                      ggCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCaa   <  1:3397194/93‑1 (MQ=255)

                                                                                      ggCAAATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCaa   <  1:7057827/93‑1 (MQ=255)

GGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACCGATAAACAAGAGACCGATAGTGATTAACGCAACTAACTGGCAGATATGGGTAAAAAAATCATCCAGTTGCCCGAG‑TTTTTCTGCGCGAAGATATGCATAGTTGCCAAAAGTTGCTGTCATTAACACGCCAAG  >  NC_000913/3725090‑3725269

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