BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,444,741	A→G	E1018G (GAG→GGG)	ytfN →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,444,741	0	A	G	99.7%	1179.3 / ‑5.6	347	E1018G (GAG→GGG)	ytfN	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (170/176); minor base T (1/0); total (171/176)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

Only 100 of 356 total aligned reads displayed.

GGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGAAAGCGC  >  NC_000913/4444651‑4444830

 ggCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGcc                                                                                        >  2:4135198/1‑93 (MQ=255)

 ggCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGcc                                                                                        >  2:2038397/1‑93 (MQ=255)

 ggCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGcc                                                                                        >  1:1008740/1‑93 (MQ=255)

    gcgcATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGaa                                                                                     >  1:5977005/1‑93 (MQ=255)

      gcATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGaaaa                                                                                   <  2:7403784/93‑1 (MQ=255)

      gcATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGaaaa                                                                                   <  2:7513453/93‑1 (MQ=255)

       cATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAAc                                                                                  <  1:661283/93‑1 (MQ=255)

        aTCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACg                                                                                 <  2:2832103/93‑1 (MQ=255)

          cGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGc                                                                               >  2:5104121/1‑93 (MQ=255)

           gtggtgCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCg                                                                              <  2:8783483/93‑1 (MQ=255)

           gtggtgCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCg                                                                              >  2:8556107/1‑93 (MQ=255)

               tgCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGa                                                                          >  1:3158396/1‑93 (MQ=255)

                  aCGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTc                                                                       >  2:669075/1‑93 (MQ=255)

                   cGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTcc                                                                      <  2:6079448/93‑1 (MQ=255)

                   cGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTcc                                                                      <  2:6398089/93‑1 (MQ=255)

                   cGATCTGCCAGAAAGCGCAGTAGGCGGCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTcc                                                                      <  1:2359327/93‑1 (MQ=255)

                    gATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCg                                                                     <  1:2587716/93‑1 (MQ=255)

                     aTCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGa                                                                    <  1:2049501/93‑1 (MQ=255)

                       cTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGAtt                                                                  >  1:2470336/1‑93 (MQ=255)

                           cGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGCAGGGCCGAAAACGGCGTCGATTCCGATTAACa                                                              >  1:3630803/1‑93 (MQ=255)

                           cGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCAGAAAACGGCGTCGATTCCGGTTAACa                                                              >  1:6720923/1‑93 (MQ=255)

                             gAAAGCGCAGTAGGCGCCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGt                                                            >  2:1554478/1‑93 (MQ=255)

                               aaGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTaa                                                          <  1:7591147/93‑1 (MQ=255)

                                 gcgcAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAAcc                                                        <  2:6323885/93‑1 (MQ=255)

                                  cgcAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCt                                                       >  1:2225175/1‑93 (MQ=255)

                                      gTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGAtt                                                   >  1:1647212/1‑93 (MQ=255)

                                       tAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTg                                                  <  1:4817985/93‑1 (MQ=255)

                                        aGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGt                                                 >  2:353694/1‑93 (MQ=255)

                                         ggCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTc                                                >  1:938640/1‑93 (MQ=255)

                                         ggCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCATAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTc                                                <  1:5731765/93‑1 (MQ=255)

                                         ggCGTCTCCAGCGATGTGGTGATGCTTAACGATAAACTGCAACCGGAAGGGCCGAAAACAGCGTCGATTCCGATTAACAGTAACCTGATTGTc                                                >  1:3792746/1‑93 (MQ=255)

                                           cGTCCCCAGCCATGTGGTGATGCTTAACGATAACCTGCAGCCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCa                                              <  2:229436/93‑1 (MQ=255)

                                            gTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCAc                                             >  1:306648/1‑93 (MQ=255)

                                            gTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCAc                                             <  1:5464261/93‑1 (MQ=255)

                                            gTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCAc                                             <  1:7064287/93‑1 (MQ=255)

                                            gTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCAc                                             >  1:6141225/1‑93 (MQ=255)

                                             tctcCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACg                                            >  1:4478325/1‑93 (MQ=255)

                                              ctcCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGt                                           <  1:6274899/93‑1 (MQ=255)

                                               tcCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCCCGtt                                          >  2:5979707/1‑93 (MQ=255)

                                               tcCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGtt                                          <  2:2475826/93‑1 (MQ=255)

                                                ccAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTg                                         >  2:4859217/1‑93 (MQ=255)

                                                ccAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTg                                         <  2:6983742/93‑1 (MQ=255)

                                                ccAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTg                                         <  1:2855163/93‑1 (MQ=255)

                                                  aGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGc                                       >  1:3225023/1‑93 (MQ=255)

                                                   gCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGca                                      <  2:5352417/93‑1 (MQ=255)

                                                     gATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGcaac                                    >  2:3345847/1‑93 (MQ=255)

                                                      aTGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGcaaca                                   >  2:3134158/1‑93 (MQ=255)

                                                      aTGTGGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGcaaca                                   >  1:4149306/1‑93 (MQ=255)

                                                       tgtgGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGcaacaa                                  >  1:1193446/1‑93 (MQ=255)

                                                       tgtgGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGcaacaa                                  >  1:3476711/1‑93 (MQ=255)

                                                       tgtgGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGcaacaa                                  <  1:7699314/93‑1 (MQ=255)

                                                       tgtgGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGcaacaa                                  >  2:2176758/1‑93 (MQ=255)

                                                       tgtgGTGATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGcaacaa                                  <  2:3699646/93‑1 (MQ=255)

                                                       ggtggtgATGCTTAACGATAGCCTGCAACCGGAAGGGCCGAAAACCGCGTCGATTCCGATTAACAGTAACCTGATTGACCACGTTGGcaacaa                                  <  1:1165345/92‑1 (MQ=255)

                                                        gtggtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAt                                 >  1:8304027/1‑93 (MQ=255)

                                                        gtggtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGACGATTCCGATTAACAGAAACCTGATTGTCCACGTTGGCAACAAt                                 >  1:1861793/1‑93 (MQ=255)

                                                         tggtgATGGTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTACCCTGATTGTCCACGTTGGCAACAAtg                                <  1:1974478/93‑1 (MQ=255)

                                                          ggtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAtgt                               <  2:6735992/93‑1 (MQ=255)

                                                          ggtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAtgt                               >  1:6867258/1‑93 (MQ=255)

                                                           gtgATGCTTAACGATAACCTGCAGCCGGAAGGGCCGAAAACGGAGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAtgtg                              <  1:6519793/93‑1 (MQ=255)

                                                           gtgATGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAAtgtg                              >  1:6444351/1‑93 (MQ=255)

                                                               tGCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCa                          >  1:2653863/1‑93 (MQ=255)

                                                                gCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCAt                         <  2:5777568/93‑1 (MQ=255)

                                                                gCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCAt                         <  1:2956795/93‑1 (MQ=255)

                                                                gCTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCAt                         <  1:3850821/93‑1 (MQ=255)

                                                                 cTTAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCAtt                        <  1:5018421/93‑1 (MQ=255)

                                                                  ttAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTg                       >  2:5880195/1‑93 (MQ=255)

                                                                  ttAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTg                       <  1:5192008/93‑1 (MQ=255)

                                                                   tAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCCCGTTGGCAACAATGTGGGCATTGa                      <  1:7213720/93‑1 (MQ=255)

                                                                   tAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCAttta                      >  1:5263263/1‑91 (MQ=255)

                                                                   tAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGa                      >  2:5702634/1‑93 (MQ=255)

                                                                   tAACGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGa                      >  1:7941412/1‑93 (MQ=255)

                                                                    aaCGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGAc                     <  2:6574576/93‑1 (MQ=255)

                                                                     aCGATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACg                    <  1:1081519/93‑1 (MQ=255)

                                                                       gATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGcc                  >  2:4091746/1‑93 (MQ=255)

                                                                       gATAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGcc                  >  2:1320238/1‑93 (MQ=255)

                                                                        aTAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCt                 <  1:4534811/93‑1 (MQ=255)

                                                                         tAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCtt                >  1:3887738/1‑93 (MQ=255)

                                                                         tAACCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCtt                <  2:8452691/93‑1 (MQ=255)

                                                                          aaCCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCttt               <  1:8530893/93‑1 (MQ=255)

                                                                           aCCTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTg              <  1:212381/93‑1 (MQ=255)

                                                                            ccTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTgg             <  2:2393385/93‑1 (MQ=255)

                                                                            ccTGCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTgg             <  2:3252215/93‑1 (MQ=255)

                                                                               gCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCCCGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCt          <  2:2613141/93‑1 (MQ=255)

                                                                               gCAACCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCt          >  2:1699265/1‑93 (MQ=255)

                                                                                 aaCCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGa        >  2:5201391/1‑93 (MQ=255)

                                                                                 aaCCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGa        >  2:1046906/1‑93 (MQ=255)

                                                                                  aCCGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGaa       <  2:5383600/93‑1 (MQ=255)

                                                                                  aCCGGAAGGGCCGAAAACGGAGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGaa       >  1:4779279/1‑93 (MQ=255)

                                                                                   ccGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGaag      >  1:7491277/1‑92 (MQ=255)

                                                                                   ccGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGaaa      <  1:3263223/93‑1 (MQ=255)

                                                                                    cGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGAAAg     <  2:7932447/93‑1 (MQ=255)

                                                                                    cGGAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGAAAg     >  1:8437975/1‑93 (MQ=255)

                                                                                    cGGAAGGGCCAAAACCGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGAAAg     <  2:3008408/93‑1 (MQ=255)

                                                                                     ggAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGAAAgc    >  1:6699905/1‑93 (MQ=255)

                                                                                     ggAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGAAAgc    >  1:3947401/1‑93 (MQ=255)

                                                                                      gAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTACCCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTCGCCTCACAgcg   >  2:4176154/1‑93 (MQ=255)

                                                                                      gAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGACAgcg   >  2:3442957/1‑93 (MQ=255)

                                                                                      gAAGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGAGCGCATTGACGCCTTTGCCCTGAAAgcg   >  2:6291052/1‑93 (MQ=255)

                                                                                       aaGGGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGAAAgcgc  <  2:8301178/93‑1 (MQ=255)

GGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTCCACGTTGGCAACAATGTGCGCATTGACGCCTTTGGCCTGAAAGCGC  >  NC_000913/4444651‑4444830

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