BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	409,218	A→G	V9A (GTT→GCT)	yaiD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	409,218	0	A	G	99.4%	550.8 / ‑4.8	174	V9A (GTT→GCT)	yaiD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (78/95); total (79/95)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 175 total aligned reads displayed.

TCCTGGCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAACCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCTTGC  >  NC_000913/409128‑409307

tCCTGGCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGcc                                                                                         <  2:42888/93‑1 (MQ=255)

  cTGGCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCAt                                                                                       >  2:287584/1‑93 (MQ=255)

    ggCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTa                                                                                     <  1:4367756/93‑1 (MQ=255)

    ggCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTa                                                                                     >  1:431270/1‑93 (MQ=255)

      cTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTaaa                                                                                   <  2:1317058/93‑1 (MQ=255)

      cTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTaaa                                                                                   <  1:524294/93‑1 (MQ=255)

       tGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAAt                                                                                  >  2:930793/1‑93 (MQ=255)

       tGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAAt                                                                                  <  2:1984711/93‑1 (MQ=255)

       tGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAAt                                                                                  >  2:4604609/1‑93 (MQ=255)

        gCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAAtt                                                                                 <  2:2198085/93‑1 (MQ=255)

        gCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAAtt                                                                                 >  2:2384873/1‑93 (MQ=255)

        gCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAAtt                                                                                 >  1:2556528/1‑93 (MQ=255)

        gCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAAtt                                                                                 <  2:2548152/93‑1 (MQ=255)

        gCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAAtt                                                                                 <  1:850730/93‑1 (MQ=255)

           gCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAAttttt                                                                              >  1:2810028/1‑93 (MQ=255)

             aTGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGa                                                                            >  1:2753408/1‑93 (MQ=255)

               ggggTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAAc                                                                          <  1:191828/93‑1 (MQ=255)

                gggTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAAcc                                                                         >  2:3549402/1‑93 (MQ=255)

                 ggTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAACCATTAAATTTTTGAACca                                                                        >  1:377355/1‑93 (MQ=255)

                   tAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACcaca                                                                      >  1:52579/1‑93 (MQ=255)

                     aaTGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGc                                                                    >  2:2913267/1‑93 (MQ=255)

                      aTGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCa                                                                   >  1:1148008/1‑93 (MQ=255)

                        gCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATg                                                                 <  1:3701144/93‑1 (MQ=255)

                          cATCTAGTCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGTCGGTAAGCCATTAAATTTTTGAACCACAGgatgat                                                               <  2:2304498/93‑1 (MQ=255)

                          cATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGAt                                                               >  2:2939370/1‑93 (MQ=255)

                          cATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGAt                                                               >  1:3839590/1‑93 (MQ=255)

                          cATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGAt                                                               >  2:4418455/1‑93 (MQ=255)

                           aTCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATa                                                              <  1:4556608/93‑1 (MQ=255)

                            tCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATaa                                                             >  2:3928888/1‑93 (MQ=255)

                              gAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGTTAAGCCATTAAATTTTTGAACCACAGCATGATAAtt                                                           <  2:2805054/93‑1 (MQ=255)

                              gAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAAtt                                                           >  1:2919435/1‑93 (MQ=255)

                                ggCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTc                                                         >  1:2700330/1‑93 (MQ=255)

                                ggCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTc                                                         >  2:3006384/1‑93 (MQ=255)

                                ggCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTc                                                         <  1:4462945/93‑1 (MQ=255)

                                  cTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCa                                                       <  2:900101/93‑1 (MQ=255)

                                  cTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCa                                                       >  2:1120394/1‑93 (MQ=255)

                                  cTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCa                                                       >  2:1254818/1‑93 (MQ=255)

                                   tAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCAc                                                      <  1:2169995/93‑1 (MQ=255)

                                    aGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACg                                                     >  2:1370418/1‑93 (MQ=255)

                                    aGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACg                                                     >  2:1513372/1‑93 (MQ=255)

                                    aGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACg                                                     <  1:1855879/93‑1 (MQ=255)

                                    aGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACg                                                     <  1:4121562/93‑1 (MQ=255)

                                     gCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTACATTTTTGAACCACAGCATGATAATTTCCACgg                                                    >  2:1129918/1‑93 (MQ=255)

                                      cTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGc                                                   <  2:3621409/93‑1 (MQ=255)

                                        gTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCt                                                 <  1:1959190/93‑1 (MQ=255)

                                        gTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCt                                                 <  1:395780/93‑1 (MQ=255)

                                           tttCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGt                                              <  1:1629938/93‑1 (MQ=255)

                                            ttCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTc                                             <  1:19882/93‑1 (MQ=255)

                                             tCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCg                                            >  2:3249450/1‑93 (MQ=255)

                                               cATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGtt                                          >  1:1820532/1‑93 (MQ=255)

                                               cATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGtt                                          >  1:762728/1‑93 (MQ=255)

                                               cATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGtt                                          >  2:3793046/1‑93 (MQ=255)

                                               cATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGtt                                          <  1:2547843/93‑1 (MQ=255)

                                                aTCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTa                                         >  2:1652347/1‑93 (MQ=255)

                                                aTCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTa                                         <  1:2969183/93‑1 (MQ=255)

                                                 tctcTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTaa                                        >  1:4649800/1‑93 (MQ=255)

                                                 tctcTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTaa                                        <  1:1271244/93‑1 (MQ=255)

                                                 tctcTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTaa                                        >  1:3487567/1‑93 (MQ=255)

                                                 tctcTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTaa                                        <  2:509306/93‑1 (MQ=255)

                                                  ctctTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTaaa                                       <  2:4164821/93‑1 (MQ=255)

                                                  ctctTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTaaa                                       >  2:3681913/1‑93 (MQ=255)

                                                  ctctTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTaaa                                       >  2:3325803/1‑93 (MQ=255)

                                                   tcttctGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAAt                                      <  2:4294245/93‑1 (MQ=255)

                                                   tcttctGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAAt                                      >  2:1589517/1‑93 (MQ=255)

                                                   tcttctGCACGCAGCGAAATCTCGCGGCTAAGACGGAAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAAt                                      >  1:1380792/1‑93 (MQ=255)

                                                    cttctGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAAtt                                     <  1:92714/93‑1 (MQ=255)

                                                      tctGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTa                                   <  1:3921523/93‑1 (MQ=255)

                                                       ctGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAg                                  <  2:2669087/93‑1 (MQ=255)

                                                        tGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGc                                 >  2:1885530/1‑93 (MQ=255)

                                                          cACGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCgg                               <  1:244280/93‑1 (MQ=255)

                                                            cGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGc                             <  2:2717260/93‑1 (MQ=255)

                                                            cGCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGc                             <  1:4594723/93‑1 (MQ=255)

                                                             gCAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCa                            >  2:1523868/1‑93 (MQ=255)

                                                              cAGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCAt                           <  2:2961433/93‑1 (MQ=255)

                                                               aGCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATg                          <  1:1895335/93‑1 (MQ=255)

                                                                gCGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGa                         >  1:4175165/1‑93 (MQ=255)

                                                                 cGAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGAt                        >  1:2995470/1‑93 (MQ=255)

                                                                  gAAATCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATa                       >  2:1197843/1‑93 (MQ=255)

                                                                   aaaTCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATaa                      <  1:954279/93‑1 (MQ=255)

                                                                   aaaTCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATaa                      <  1:1367063/93‑1 (MQ=255)

                                                                   aaaTCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATaa                      <  2:3259192/93‑1 (MQ=255)

                                                                   aaaTCTCGCGGCTAAGAC‑GTAAGCCCTTAAATTTTTGAACCACAGCATGATAATTTTCCCCGCCTTGTCGTTAAATTTCGCGGGCATGATAAc                     >  2:4046245/1‑93 (MQ=255)

                                                                     aTCTCGCGGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACg                    >  1:3267581/1‑93 (MQ=255)

                                                                         cgcgGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAAtt                <  1:4295839/93‑1 (MQ=255)

                                                                         cgcgGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAAtt                <  1:132388/93‑1 (MQ=255)

                                                                         cgcgGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAAtt                <  2:1306622/93‑1 (MQ=255)

                                                                           cgGCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGt              <  1:1180060/93‑1 (MQ=255)

                                                                            ggCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTc             <  1:1565707/93‑1 (MQ=255)

                                                                            ggCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTc             <  2:2273614/93‑1 (MQ=255)

                                                                            ggCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTc             >  1:866575/1‑93 (MQ=255)

                                                                             gCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTcg            >  2:3399113/1‑93 (MQ=255)

                                                                             gCTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTcg            >  1:4667172/1‑93 (MQ=255)

                                                                              cTAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTcgg           <  2:590186/93‑1 (MQ=255)

                                                                               tAAGACGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTcggc          >  2:2458704/1‑93 (MQ=255)

                                                                                   aCGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGcc      >  2:1165250/1‑93 (MQ=255)

                                                                                    cGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCt     <  1:3979579/93‑1 (MQ=255)

                                                                                    cGGTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCt     >  1:704531/1‑93 (MQ=255)

                                                                                     ggTAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCtt    >  1:3264076/1‑93 (MQ=255)

                                                                                       tAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCTTGc  >  1:333297/1‑93 (MQ=255)

                                                                                       tAAGCCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCTTGc  >  2:3586627/1‑93 (MQ=255)

TCCTGGCTGCCGCATGGGGTAAATGCCATCGAGGCTAGCTGTTTTTCCATCTCTTCTGCACGCAGCGAAATCTCGCGGCTAAGACGGTAAACCATTAAATTTTTGAACCACAGCATGATAATTTCCACGGCCTTGTCGTTAAATTTAGCGGGCATGATAACGAATTGTCGGCGGCCTTGC  >  NC_000913/409128‑409307

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