BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,429,526	C→T	A12A (GCG→GCA)	ytfE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,429,526	0	C	T	98.1%	964.7 / 1.0	268	A12A (GCG→GCA)	ytfE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (4/1); new base T (113/150); total (117/151)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 279 total aligned reads displayed.

TGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGG  >  NC_000913/4429436‑4429615

tGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAgtgc                                                                                         >  1:647886/1‑93 (MQ=255)

  ccgccgCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCa                                                                                       <  2:3537499/93‑1 (MQ=255)

    gccgcGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGt                                                                                     >  1:5680779/1‑93 (MQ=255)

     ccgcGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGtt                                                                                    <  2:5504700/93‑1 (MQ=255)

       gcgcgcgcCAGCTTCTGCTTACCGCCGCAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCa                                                                                  <  1:3118844/93‑1 (MQ=255)

         gcgcgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCAcc                                                                                <  1:5060970/93‑1 (MQ=255)

          cgcgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCt                                                                               >  2:2411418/1‑93 (MQ=255)

          cgcgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCt                                                                               >  2:2882628/1‑93 (MQ=255)

           gcgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTa                                                                              >  2:2235944/1‑93 (MQ=255)

            cgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTaa                                                                             <  2:2367911/93‑1 (MQ=255)

            cgcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTaa                                                                             >  1:3864333/1‑93 (MQ=255)

             gcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTaaa                                                                            >  2:1287606/1‑93 (MQ=255)

             gcCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTaaa                                                                            <  1:1058840/93‑1 (MQ=255)

              ccAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAg                                                                           >  2:5348586/1‑93 (MQ=255)

                aGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGt                                                                         <  2:1882301/93‑1 (MQ=255)

                 gCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGtt                                                                        <  2:4780104/93‑1 (MQ=255)

                   gTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTgg                                                                      >  2:2442321/1‑93 (MQ=255)

                   gTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTgg                                                                      >  2:4627436/1‑93 (MQ=255)

                    tCTGCTTACCGCCACAGCAGTAATCCATATCATGTTGACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGt                                                                     <  2:2643002/93‑1 (MQ=255)

                    tCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGt                                                                     >  1:6114371/1‑93 (MQ=255)

                    tCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGt                                                                     >  1:448740/1‑93 (MQ=255)

                    tCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGt                                                                     >  1:3037970/1‑93 (MQ=255)

                     cTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTc                                                                    <  1:812169/93‑1 (MQ=255)

                         ttACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCTCGa                                                                >  1:5835346/1‑93 (MQ=255)

                         ttACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGa                                                                <  1:5988474/93‑1 (MQ=255)

                         ttACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGa                                                                <  1:5387108/93‑1 (MQ=255)

                         ttACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGa                                                                >  1:4382045/1‑93 (MQ=255)

                          tACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGAt                                                               >  1:159493/1‑93 (MQ=255)

                          tACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGAt                                                               >  2:74633/1‑93 (MQ=255)

                              gCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGc                                                           >  1:3593697/1‑93 (MQ=255)

                              gCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGc                                                           <  2:3819437/93‑1 (MQ=255)

                                  cagcagTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATa                                                       <  2:5637199/93‑1 (MQ=255)

                                   agtagTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAg                                                      <  2:3637137/90‑1 (MQ=255)

                                       gTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGa                                                  <  1:4570217/93‑1 (MQ=255)

                                       gTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGa                                                  <  1:4682468/93‑1 (MQ=255)

                                       gTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGa                                                  >  1:4841790/1‑93 (MQ=255)

                                        tAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGAt                                                 <  2:4389299/93‑1 (MQ=255)

                                        tAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGAt                                                 <  2:1348106/93‑1 (MQ=255)

                                          aTCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATAc                                               <  1:1899892/93‑1 (MQ=255)

                                          aTCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATAc                                               <  1:735704/93‑1 (MQ=255)

                                           tCCATATCATATTTACGAAACAGAGCTGAAGAGCGCGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATAcc                                              >  1:2235280/1‑93 (MQ=255)

                                             catatcatatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTc                                            >  1:1992644/1‑93 (MQ=255)

                                              atatcatatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCa                                           <  2:4479364/93‑1 (MQ=255)

                                               tatcatatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCAt                                          >  1:2687075/1‑93 (MQ=255)

                                               tatcatatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCAt                                          <  1:259554/93‑1 (MQ=255)

                                               tatcatatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCAt                                          <  1:420491/93‑1 (MQ=255)

                                               tatcatatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCAt                                          >  1:1449796/1‑93 (MQ=255)

                                               tatcatatTTACGAAACAGAGCGGAAGCGCGAGGAAAAGAGAGTCCCAGTTGACCTAAAGGTTGGTCGCGATAAACCATAGCTGAAACCTCAt                                          >  2:1189062/1‑93 (MQ=255)

                                                atcatatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCAtt                                         >  1:4776221/1‑93 (MQ=255)

                                                   atatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCtt                                      <  1:1184328/93‑1 (MQ=255)

                                                    tatTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCATCTAAAGGTTGGTGGCGATAAGCCATAGCTGATATCTCATTCTta                                     >  1:5108884/1‑93 (MQ=255)

                                                     atTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaa                                    <  1:2762179/93‑1 (MQ=255)

                                                      tttACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaat                                   >  1:3957113/1‑93 (MQ=255)

                                                      tttACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaat                                   <  1:4602171/93‑1 (MQ=255)

                                                       ttACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaata                                  >  2:3255282/1‑93 (MQ=255)

                                                       ttACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCCtaata                                  <  1:3438234/93‑1 (MQ=255)

                                                        tACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaataa                                 <  2:3689310/93‑1 (MQ=255)

                                                        tACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTtaataa                                 <  1:5322496/93‑1 (MQ=255)

                                                         aCGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataag                                <  1:2068178/93‑1 (MQ=255)

                                                         aCGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataag                                <  1:667853/93‑1 (MQ=255)

                                                          cGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataaga                               >  1:3322923/1‑93 (MQ=255)

                                                          cGAAACAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataaga                               <  1:1471408/93‑1 (MQ=255)

                                                            aaaCAGGGCGAAAGCGAGAGGAATAGAGAGGGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataagata                             <  1:244279/93‑1 (MQ=255)

                                                            aaaCAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataagata                             <  1:1881172/93‑1 (MQ=255)

                                                              aCAGAGCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAataagataag                           <  2:514254/93‑1 (MQ=255)

                                                                agagCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATaagataagat                         >  2:48656/1‑93 (MQ=255)

                                                                agagCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATaagataagat                         >  2:906853/1‑93 (MQ=255)

                                                                agagCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATaagataagat                         >  2:3045316/1‑93 (MQ=255)

                                                                 gagCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATg                        >  2:85038/1‑93 (MQ=255)

                                                                  agCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGt                       <  1:2911631/93‑1 (MQ=255)

                                                                  agCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGt                       <  1:1888775/93‑1 (MQ=255)

                                                                  agCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGt                       >  2:1959485/1‑93 (MQ=255)

                                                                   gCTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATCTa                      >  2:6005946/1‑93 (MQ=255)

                                                                    cTGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTAt                     >  2:4253829/1‑93 (MQ=255)

                                                                     tGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTAtt                    >  2:6108621/1‑93 (MQ=255)

                                                                     tGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTAtt                    >  2:5210757/1‑93 (MQ=255)

                                                                     tGAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTAtt                    >  2:258806/1‑93 (MQ=255)

                                                                      gAAGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTAttt                   >  1:250661/1‑93 (MQ=255)

                                                                        aGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTa                 <  1:5792561/93‑1 (MQ=255)

                                                                        aGCGCGAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTa                 >  2:5092137/1‑93 (MQ=255)

                                                                          cgcgAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTaaa               <  2:2930613/93‑1 (MQ=255)

                                                                          cgcgAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTaaa               <  2:251379/93‑1 (MQ=255)

                                                                          cgcgAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTaaa               <  1:5830206/93‑1 (MQ=255)

                                                                           gcgAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAAt              <  2:3896670/93‑1 (MQ=255)

                                                                           gcgAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAAt              <  1:956044/93‑1 (MQ=255)

                                                                             gAGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGc            >  1:5225741/1‑93 (MQ=255)

                                                                              aGGAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCa           >  1:5884394/1‑93 (MQ=255)

                                                                               ggAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCAt          <  2:4002627/93‑1 (MQ=255)

                                                                                gAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATc         <  2:1912140/93‑1 (MQ=255)

                                                                                gAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATc         <  1:5727405/93‑1 (MQ=255)

                                                                                gAATAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATc         >  1:281298/1‑93 (MQ=255)

                                                                                 aaTAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCt        <  2:1624351/93‑1 (MQ=255)

                                                                                   tAGAGAGTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCttt      >  2:1032065/1‑93 (MQ=255)

                                                                                   tAGAGAGTGCCAGTTCACCTAAAGGTAGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCttt      <  1:5481737/93‑1 (MQ=255)

                                                                                     gagagTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTaa    <  1:5643203/93‑1 (MQ=255)

                                                                                     gagagTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTaa    <  2:5859176/93‑1 (MQ=255)

                                                                                      agagTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAg   <  2:4717696/93‑1 (MQ=255)

                                                                                      agagTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAg   <  1:1393247/93‑1 (MQ=255)

                                                                                      agagTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAg   <  1:1036950/93‑1 (MQ=255)

                                                                                       gagTGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAgg  >  1:1612076/1‑93 (MQ=255)

TGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCCATATCATATTTACGAAACAGAGCTGAAGCGCGAGGAATAGAGAGCGCCAGTTCACCTAAAGGTTGGTCGCGATAAGCCATAGCTGATACCTCATTCTTAATAAGATAAGATGTATTTTAAATGCATCTTTAAGG  >  NC_000913/4429436‑4429615

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