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.6%	1169.1 / ‑3.3	256	coding (504/672 nt)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base . (154/101); minor base A (1/0); total (155/101)
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 283 total aligned reads displayed.

AAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTA  >  NC_000913/3873533‑3873713

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

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

  aaaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCgga                                                                                        <  2:3802632/93‑2 (MQ=255)

  aaaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTANNAATACTATTCACCCAGANGNATTCCnga                                                                                        >  2:2446801/1‑92 (MQ=255)

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

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

    aTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatg                                                                                      <  2:1067641/93‑4 (MQ=255)

    aTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatg                                                                                      >  2:1561877/1‑90 (MQ=255)

     tAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatgg                                                                                     <  1:2969517/93‑5 (MQ=255)

      aaaaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTATAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgt                                                                                   >  2:4168087/1‑93 (MQ=255)

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

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

         aaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGtt                                                                                <  2:3565914/93‑1 (MQ=255)

         aaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGtt                                                                                <  2:1683751/93‑1 (MQ=255)

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

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

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

             tGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCg                                                                            >  1:739238/1‑93 (MQ=255)

              ggATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgc                                                                           >  2:2237083/1‑93 (MQ=255)

                 ttACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTTACCCAGAGGCATTCC‑GGATGGTGTTTCCgcagc                                                                        >  2:1761674/1‑93 (MQ=255)

                  tACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCtcagca                                                                       >  2:1629993/1‑93 (MQ=255)

                   aCTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCAt                                                                      >  1:4010509/1‑93 (MQ=255)

                      cAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTAc                                                                   <  2:143168/93‑1 (MQ=255)

                      cAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTAc                                                                   <  1:2826089/93‑1 (MQ=255)

                      cAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTAc                                                                   <  2:1889176/93‑1 (MQ=255)

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

                        aaGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTg                                                                 <  1:317026/93‑1 (MQ=255)

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

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

                         aGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGc                                                                >  1:1767641/1‑93 (MQ=255)

                         aGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGc                                                                <  2:3025086/93‑1 (MQ=255)

                          ggTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGcc                                                               >  2:608744/1‑93 (MQ=255)

                              aTGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTg                                                           <  1:2027281/93‑1 (MQ=255)

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

                               tGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGa                                                          >  2:4049817/1‑93 (MQ=255)

                                 gCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGCCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACt                                                        <  1:1702720/93‑1 (MQ=255)

                                    ttCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGtt                                                     >  1:2682154/1‑93 (MQ=255)

                                     tCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTa                                                    >  2:3175027/1‑93 (MQ=255)

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

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

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

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

                                        aCCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGc                                                 <  2:460598/93‑1 (MQ=255)

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

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

                                           atccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAAGCATTCC‑GGCTGGTGTTTCGGCAGCATTACTGCCAGTGACAGCTACGCCgg                                              >  1:1500408/1‑93 (MQ=255)

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

                                            tccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATACC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGa                                             >  1:3041000/1‑93 (MQ=255)

                                              catTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAca                                           >  1:24499/1‑93 (MQ=255)

                                              catTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAca                                           >  1:4551874/1‑93 (MQ=255)

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

                                                tttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatc                                         <  1:913990/93‑1 (MQ=255)

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

                                                  tAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatcat                                       >  2:2358634/1‑93 (MQ=255)

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

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

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

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

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

                                                         cGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCa                                <  2:3662810/93‑1 (MQ=255)

                                                           aTTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCAcc                              >  2:4344497/1‑93 (MQ=255)

                                                           aTTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCAcc                              >  2:490295/1‑93 (MQ=255)

                                                              ttATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCt                           <  1:3963575/93‑1 (MQ=255)

                                                               tATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTg                          >  1:3561681/1‑93 (MQ=255)

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

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

                                                                  gAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCaa                       >  2:3261808/1‑93 (MQ=255)

                                                                   aaTACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAt                      >  1:2304703/1‑93 (MQ=255)

                                                                    aTACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAtt                     <  2:2006091/93‑1 (MQ=255)

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

                                                                     tACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAttt                    >  2:1363075/1‑93 (MQ=255)

                                                                      aCTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTc                   <  1:117194/93‑1 (MQ=255)

                                                                      aCTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTc                   <  1:4440401/93‑1 (MQ=255)

                                                                      aCTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATAACTGCACCGCTGCAATTTc                   <  2:4566731/93‑1 (MQ=255)

                                                                       cTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCt                  <  2:1569771/93‑1 (MQ=255)

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

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

                                                                         aTTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTgg                >  2:2730575/1‑93 (MQ=255)

                                                                         aTTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTgg                >  2:3787694/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                                                 agagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTatat        >  2:4593447/1‑93 (MQ=255)

                                                                                 agagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTatat        >  2:1123338/1‑93 (MQ=255)

                                                                                 agagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTatat        >  1:1001392/1‑93 (MQ=255)

                                                                                  gagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACAGCTGCAATTTCTGGTTGTATATg       >  2:3969271/1‑93 (MQ=255)

                                                                                   agGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGc      >  2:1953572/1‑93 (MQ=255)

                                                                                    ggCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCa     <  1:3335119/93‑1 (MQ=255)

                                                                                    ggCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCa     >  1:2116330/1‑93 (MQ=255)

                                                                                     gCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAg    <  1:3790483/93‑1 (MQ=255)

                                                                                      cATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGt   <  1:2174606/93‑1 (MQ=255)

                                                                                      cATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGt   >  2:68495/1‑93 (MQ=255)

AAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTA  >  NC_000913/3873533‑3873713

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