BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,627	Δ1 bp	coding (506/672 nt)	yidX →	putative replicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,625	0	G	.	99.8%	2032.1 / ‑3.7	445	coding (504/672 nt)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base . (250/194); minor base A (1/0); total (251/194)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 485 total aligned reads displayed.

AAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTA  >  NC_000913/3873533‑3873713

 aaaaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCgg                                                                                         >  1:5865458/1‑93 (MQ=255)

 aaaaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCgg                                                                                         >  1:7857797/1‑93 (MQ=255)

   aaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggat                                                                                       <  1:3968259/93‑3 (MQ=255)

   aaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggat                                                                                       >  1:883175/1‑91 (MQ=255)

   aaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggat                                                                                       >  1:4125851/1‑91 (MQ=255)

       aaaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgtg                                                                                  >  2:3417293/1‑93 (MQ=255)

       aaaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgtg                                                                                  >  2:1453460/1‑93 (MQ=255)

          aCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGttt                                                                               >  1:7014702/1‑93 (MQ=255)

          aCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAATGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGttt                                                                               <  2:2147324/93‑1 (MQ=255)

           cATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTc                                                                              >  2:5880792/1‑93 (MQ=255)

           cATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTc                                                                              <  1:174604/93‑1 (MQ=255)

           cATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTc                                                                              <  2:6776541/93‑1 (MQ=255)

            aTGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTcc                                                                             <  1:6274321/93‑1 (MQ=255)

            aTGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTcc                                                                             >  1:6734314/1‑93 (MQ=255)

            aTGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTcc                                                                             <  2:5935012/93‑1 (MQ=255)

               gATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgca                                                                          >  2:457795/1‑93 (MQ=255)

                aTTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgcag                                                                         >  2:1910863/1‑93 (MQ=255)

                 ttACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgcagc                                                                        <  1:6124767/93‑1 (MQ=255)

                   aCTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCAt                                                                      <  2:6446840/93‑1 (MQ=255)

                   aCTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCAt                                                                      <  1:4012368/93‑1 (MQ=255)

                   aCTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCAt                                                                      <  1:7672486/93‑1 (MQ=255)

                    cTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCAtt                                                                     >  1:3238986/1‑93 (MQ=255)

                    cTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCAtt                                                                     >  1:2020247/1‑93 (MQ=255)

                    cTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCAtt                                                                     >  2:6136284/1‑93 (MQ=255)

                     tCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTa                                                                    >  2:1766292/1‑93 (MQ=255)

                     tCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTa                                                                    >  2:923715/1‑93 (MQ=255)

                     tCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTa                                                                    >  2:233869/1‑93 (MQ=255)

                     tCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTa                                                                    >  2:5298447/1‑93 (MQ=255)

                     tCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTa                                                                    >  2:7648752/1‑93 (MQ=255)

                     tCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTa                                                                    <  1:4479280/93‑1 (MQ=255)

                      cAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATCCTACTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTAc                                                                   <  2:795136/93‑1 (MQ=255)

                       aaaGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACt                                                                  >  1:7932891/1‑93 (MQ=255)

                        aaGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTg                                                                 >  2:1565401/1‑93 (MQ=255)

                          ggTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGcc                                                               <  1:1045154/93‑1 (MQ=255)

                           gTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCa                                                              <  1:7543833/93‑1 (MQ=255)

                            tGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAg                                                             <  1:1556740/93‑1 (MQ=255)

                             gATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGt                                                            <  1:1485890/93‑1 (MQ=255)

                             gATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGt                                                            <  1:2231518/93‑1 (MQ=255)

                              aTGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTg                                                           >  2:4828331/1‑93 (MQ=255)

                               tGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGa                                                          <  2:7315874/93‑1 (MQ=255)

                                     tCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTa                                                    <  1:574445/93‑1 (MQ=255)

                                      cTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTAc                                                   >  2:2297694/1‑93 (MQ=255)

                                      cTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTAc                                                   <  1:1837754/93‑1 (MQ=255)

                                       tACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACg                                                  >  1:7149675/1‑93 (MQ=255)

                                          catccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTg                                               <  1:7443585/93‑1 (MQ=255)

                                          catccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTg                                               >  2:2470920/1‑93 (MQ=255)

                                           atccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTgg                                              >  2:4783453/1‑93 (MQ=255)

                                           atccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTgg                                              >  1:3870475/1‑93 (MQ=255)

                                           atccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTgg                                              >  1:115884/1‑93 (MQ=255)

                                            tccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGa                                             >  2:4315601/1‑93 (MQ=255)

                                            tccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGa                                             >  2:3510033/1‑93 (MQ=255)

                                            tccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGa                                             >  1:3040793/1‑93 (MQ=255)

                                              catTTAAAGGGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGCGACTGTTACGCCGGAca                                           >  2:7044652/1‑93 (MQ=255)

                                               atTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcat                                          >  2:5677894/1‑93 (MQ=255)

                                               atTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcat                                          >  1:5438384/1‑93 (MQ=255)

                                                tttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatc                                         >  2:2692893/1‑93 (MQ=255)

                                                tttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatc                                         >  1:2543144/1‑93 (MQ=255)

                                                tttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatc                                         >  2:3610140/1‑93 (MQ=255)

                                                 ttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatca                                        <  1:3719015/93‑1 (MQ=255)

                                                 ttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatca                                        <  1:6668106/93‑1 (MQ=255)

                                                 ttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatca                                        <  2:2805519/93‑1 (MQ=255)

                                                  tAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatcat                                       <  2:4177400/93‑1 (MQ=255)

                                                   aaaGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCAtt                                      >  2:5509299/1‑93 (MQ=255)

                                                   aaaGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCAtt                                      <  2:114480/93‑1 (MQ=255)

                                                   aaaGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCAtt                                      <  2:2979812/93‑1 (MQ=255)

                                                     aGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTAc                                    <  1:5953839/93‑1 (MQ=255)

                                                     aGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTAc                                    <  2:6232098/93‑1 (MQ=255)

                                                      atGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACt                                   <  2:2577306/92‑1 (MQ=255)

                                                        gCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGc                                 <  2:5210062/93‑1 (MQ=255)

                                                         cGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCa                                >  2:3825233/1‑93 (MQ=255)

                                                          gATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCAc                               >  2:1386650/1‑93 (MQ=255)

                                                          gATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCAc                               >  1:500695/1‑93 (MQ=255)

                                                           aTTTTATAAATACTATGCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCAcc                              <  1:5317715/93‑1 (MQ=255)

                                                             tttATTAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGc                            <  2:1202952/93‑1 (MQ=255)

                                                                aTGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGc                         <  1:5697408/93‑1 (MQ=255)

                                                                aTGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGc                         >  1:406217/1‑93 (MQ=255)

                                                                 tGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCa                        >  2:29290/1‑93 (MQ=255)

                                                                   aaTACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAt                      <  1:4237975/93‑1 (MQ=255)

                                                                    aTACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAtt                     >  1:2942720/1‑93 (MQ=255)

                                                                      aCTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTc                   >  1:3424263/1‑93 (MQ=255)

                                                                       cTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCt                  >  1:779612/1‑93 (MQ=255)

                                                                        tATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTg                 >  2:401103/1‑93 (MQ=255)

                                                                          ttCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGt               >  2:6694693/1‑93 (MQ=255)

                                                                           tCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGtt              <  2:1035437/93‑1 (MQ=255)

                                                                            cACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTg             <  1:1952483/93‑1 (MQ=255)

                                                                             aCCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGt            >  2:3853554/1‑93 (MQ=255)

                                                                             aCCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGt            >  1:6484470/1‑93 (MQ=255)

                                                                              cccAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGta           >  2:7199596/1‑93 (MQ=255)

                                                                              cccAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGta           <  1:5549478/93‑1 (MQ=255)

                                                                               ccAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGtat          >  2:7866133/1‑93 (MQ=255)

                                                                                cAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGtata         >  1:4353392/1‑93 (MQ=255)

                                                                                cAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGtata         <  1:6772939/93‑1 (MQ=255)

                                                                                cAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGtata         <  1:6812206/93‑1 (MQ=255)

                                                                                cAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGtata         <  1:2379681/93‑1 (MQ=255)

                                                                                 agagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTatat        <  2:3246612/93‑1 (MQ=255)

                                                                                 agagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTatat        <  1:3541371/93‑1 (MQ=255)

                                                                                  gagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATg       <  2:7510911/93‑1 (MQ=255)

                                                                                  gagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACAGTTACGCTGGACATCATTACTGGACCGCTGCAATTTCTGGTTGTATATg       >  1:50989/1‑93 (MQ=255)

                                                                                    ggCATTCAGGG‑TGGTGTTTCCGCAGCGTTACAGACAGTGAATGTTACGCTGGACATCATTACTGCAACGCTGCAATTTCTGGTTGTATATGCa     >  2:49876/1‑93 (MQ=255)

                                                                                     gCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAg    >  1:2798678/1‑93 (MQ=255)

AAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTA  >  NC_000913/3873533‑3873713

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