BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,137,208	A→G	T27T (ACT→ACC)	talC ←	putative transaldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,137,208	0	A	G	99.4%	1119.8 / ‑5.3	361	T27T (ACT→ACC)	talC	putative transaldolase
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (173/185); total (176/185)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 395 total aligned reads displayed.

AAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC  >  NC_000913/4137118‑4137297

 aaCAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGggttggtt                                                                                        <  1:4233315/93‑1 (MQ=255)

 aaCAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGggttggtt                                                                                        >  1:6151642/1‑93 (MQ=255)

  aCAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGgttggttg                                                                                       >  1:1947647/1‑93 (MQ=255)

   cAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGt                                                                                      >  1:6390905/1‑93 (MQ=255)

    aGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTc                                                                                     >  1:2022175/1‑93 (MQ=255)

     gAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGtta                                                                                    >  1:74242/1‑91 (MQ=255)

     gAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGtga                                                                                    >  2:149804/1‑91 (MQ=255)

     gAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCa                                                                                    >  2:7378570/1‑93 (MQ=255)

     gAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCa                                                                                    >  2:5253737/1‑93 (MQ=255)

      aaTGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCAc                                                                                   >  2:1430627/1‑93 (MQ=255)

              catccCCAAGCCTTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCaa                                                                           <  1:895635/93‑1 (MQ=255)

              catcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCaa                                                                           >  2:5145702/1‑93 (MQ=255)

               atcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAAt                                                                          >  2:4400394/1‑93 (MQ=255)

                 caCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATgg                                                                        <  1:2130944/93‑1 (MQ=255)

                  cccAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATggg                                                                       <  1:3075490/92‑1 (MQ=255)

                  aCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATggg                                                                       >  1:3485646/1‑93 (MQ=255)

                   ccAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATgggg                                                                      <  2:3496584/93‑1 (MQ=255)

                      aTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCTGCGATAATGCTCGGGTTGGTTGTCTCCCTGGCAATGGGGAat                                                                   >  2:7048383/1‑93 (MQ=255)

                      aTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAat                                                                   <  1:6899444/93‑1 (MQ=255)

                      aTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAat                                                                   >  2:2639719/1‑93 (MQ=255)

                       tCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                  >  1:508158/1‑93 (MQ=255)

                       tCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTAGCAATGGGGAata                                                                  >  1:1325168/1‑93 (MQ=255)

                        cGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAatat                                                                 <  1:3718160/93‑1 (MQ=255)

                        cGCTTTTTGCAGACGCGGCAGCACCTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAatat                                                                 <  2:6365026/93‑1 (MQ=255)

                         gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATg                                                                <  2:7500534/93‑1 (MQ=255)

                          cTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGc                                                               >  2:5072729/1‑93 (MQ=255)

                          cTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGc                                                               >  2:3789957/1‑93 (MQ=255)

                           tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCg                                                              >  1:6782327/1‑93 (MQ=255)

                           tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCg                                                              >  2:4449003/1‑93 (MQ=255)

                           tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCg                                                              >  1:1140641/1‑93 (MQ=255)

                            ttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGt                                                             >  1:5569406/1‑93 (MQ=255)

                            ttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGt                                                             <  1:3140458/93‑1 (MQ=255)

                             tttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTg                                                            <  1:1388284/93‑1 (MQ=255)

                              ttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGc                                                           <  1:3015445/93‑1 (MQ=255)

                               tGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGcc                                                          <  2:6785410/93‑1 (MQ=255)

                                gCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCa                                                         <  2:3155302/93‑1 (MQ=255)

                                 cAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAg                                                        <  1:7671499/93‑1 (MQ=255)

                                    aCGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACg                                                     <  1:4696175/93‑1 (MQ=255)

                                    aCGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACg                                                     >  1:3247632/1‑93 (MQ=255)

                                     cgcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGt                                                    >  1:5589134/1‑93 (MQ=255)

                                      gcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGtt                                                   >  2:3616017/1‑93 (MQ=255)

                                      gcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGtt                                                   >  2:1651058/1‑93 (MQ=255)

                                      gcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGtt                                                   <  1:2316960/93‑1 (MQ=255)

                                        ggCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCg                                                 >  2:4214955/1‑93 (MQ=255)

                                         gcagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGa                                                <  1:798769/93‑1 (MQ=255)

                                         gcagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGa                                                >  1:477417/1‑93 (MQ=255)

                                         gcagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCAACCCGGCAATGGGGCATATGCGTGCCAGACGTTCgt                                                >  2:2480858/1‑92 (MQ=255)

                                          cagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGAc                                               <  1:4318011/93‑1 (MQ=255)

                                            gcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACtt                                             >  1:4424094/1‑93 (MQ=255)

                                             caCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTc                                            <  1:2710606/93‑1 (MQ=255)

                                             caCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTc                                            <  2:4067876/93‑1 (MQ=255)

                                             caCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTc                                            <  1:6259206/93‑1 (MQ=255)

                                               cTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTg                                          <  2:853016/93‑1 (MQ=255)

                                               cTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTg                                          <  2:5596410/93‑1 (MQ=255)

                                               cTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTg                                          >  2:298290/1‑93 (MQ=255)

                                                  cccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGa                                       <  2:657310/93‑1 (MQ=255)

                                                  cccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGa                                       <  1:5679585/93‑1 (MQ=255)

                                                  cccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGa                                       >  2:5730337/1‑93 (MQ=255)

                                                  cccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGa                                       >  2:3347639/1‑93 (MQ=255)

                                                   ccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGAc                                      >  2:7674322/1‑93 (MQ=255)

                                                    cATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACg                                     >  1:3949384/1‑93 (MQ=255)

                                                        tGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAg                                 >  2:5970098/1‑93 (MQ=255)

                                                         ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGc                                >  2:7574844/1‑93 (MQ=255)

                                                         ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCCGCGACGTTAGc                                >  2:4712212/1‑93 (MQ=255)

                                                          gACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCg                               >  1:4969547/1‑93 (MQ=255)

                                                           aCTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCATGCCAGACGTTCGACTTCTGCGACGTTAGCgg                              <  2:157255/93‑1 (MQ=255)

                                                            cTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgt                             >  2:1912244/1‑93 (MQ=255)

                                                            cTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgt                             <  1:1636019/93‑1 (MQ=255)

                                                             tCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgtg                            >  1:3842939/1‑93 (MQ=255)

                                                             tCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgtg                            >  2:2370451/1‑93 (MQ=255)

                                                              ccTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgtgt                           >  2:1254439/1‑93 (MQ=255)

                                                               cTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTc                          >  1:1781412/1‑93 (MQ=255)

                                                               cTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTc                          >  2:7152641/1‑93 (MQ=255)

                                                                 tGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCa                        <  1:1426360/93‑1 (MQ=255)

                                                                 tGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCa                        <  2:7769817/93‑1 (MQ=255)

                                                                 tGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCa                        >  2:3267518/1‑93 (MQ=255)

                                                                   cTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGa                      <  1:3977063/93‑1 (MQ=255)

                                                                   cTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGa                      <  2:2851234/93‑1 (MQ=255)

                                                                   cTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGa                      <  1:169350/93‑1 (MQ=255)

                                                                    tGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGAt                     <  2:1053476/93‑1 (MQ=255)

                                                                      gCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATAc                   <  1:6043294/93‑1 (MQ=255)

                                                                       cAGCTATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACa                  <  1:7594491/93‑1 (MQ=255)

                                                                       cAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGGGTCCAGATACa                  >  1:1024102/1‑93 (MQ=255)

                                                                         gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                <  1:588381/93‑1 (MQ=255)

                                                                         gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                >  2:6736576/1‑93 (MQ=255)

                                                                           gATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTc              <  2:6098218/93‑1 (MQ=255)

                                                                           gATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTc              >  2:622699/1‑93 (MQ=255)

                                                                           gATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTc              <  1:2491229/93‑1 (MQ=255)

                                                                            aTAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTcc             >  1:6727652/1‑93 (MQ=255)

                                                                             tAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCa            >  1:6517953/1‑93 (MQ=255)

                                                                                tGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatga         <  1:7048384/93‑1 (MQ=255)

                                                                                tGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatga         <  1:7311425/93‑1 (MQ=255)

                                                                                  cTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatgatg       <  2:2800750/93‑1 (MQ=255)

                                                                                   tCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGt      <  1:5370561/93‑1 (MQ=255)

                                                                                   tCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGt      >  1:661992/1‑93 (MQ=255)

                                                                                    cGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGtt     <  2:6940007/93‑1 (MQ=255)

                                                                                      ggttggttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTc   <  2:2187027/93‑1 (MQ=255)

                                                                                      ggttggttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTc   >  1:2872021/1‑93 (MQ=255)

                                                                                       gttggttgTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTcc  <  2:5810108/93‑1 (MQ=255)

                                                                                       gttggttgTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTcc  >  2:5515692/1‑93 (MQ=255)

AAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCC  >  NC_000913/4137118‑4137297

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