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.2%	1126.5 / 1.7	246	coding (504/672 nt)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base G (1/0); major base . (126/116); minor base T (2/0); total (130/116)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 274 total aligned reads displayed.

AAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTA  >  NC_000913/3873533‑3873713

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

  aaaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCgga                                                                                        >  2:836176/1‑92 (MQ=255)

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

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

   aaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggat                                                                                       <  2:1945755/93‑3 (MQ=255)

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

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

    aTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatg                                                                                      <  1:5013771/93‑4 (MQ=255)

     tAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatgg                                                                                     >  1:1009164/1‑89 (MQ=255)

     tAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatgg                                                                                     >  1:4410053/1‑89 (MQ=255)

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

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

       aaaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgtg                                                                                  <  1:722214/93‑1 (MQ=255)

       aaaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgtg                                                                                  <  1:3082605/93‑1 (MQ=255)

        aaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgtgt                                                                                 >  2:1583029/1‑93 (MQ=255)

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

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

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

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

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

              ggATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgc                                                                           <  1:4734526/93‑1 (MQ=255)

                 ttACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgcagc                                                                        >  1:5758038/1‑93 (MQ=255)

                  tACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgcagca                                                                       >  1:2787846/1‑93 (MQ=255)

                  tACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgcagca                                                                       >  2:4469225/1‑93 (MQ=255)

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

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

                      cAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTAc                                                                   >  2:1733201/1‑93 (MQ=255)

                        aaGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATGCTATTCACCCAGAGGCATTCC‑GAATGGTGTTTCCGCAGCATTACTg                                                                 <  2:6061392/93‑1 (MQ=255)

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

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

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

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

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

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

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

                           gTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCa                                                              >  1:1736386/1‑93 (MQ=255)

                             gATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGt                                                            >  2:353042/1‑93 (MQ=255)

                               tGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCCCCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGa                                                          <  2:2439501/93‑1 (MQ=255)

                                ggCGTTCTCCCACCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGAc                                                         <  2:2219567/93‑1 (MQ=255)

                                ggCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGAc                                                         >  1:3044295/1‑93 (MQ=255)

                                   gTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGt                                                      <  1:5348899/93‑1 (MQ=255)

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

                                    ttCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGtt                                                     <  2:4717403/93‑1 (MQ=255)

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

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

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

                                        aCCATCCATTTAAAGTACGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGc                                                 <  1:2476684/93‑1 (MQ=255)

                                         ccatccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCt                                                <  1:5212678/93‑1 (MQ=255)

                                         ccatccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCt                                                >  2:2337204/1‑93 (MQ=255)

                                         ccatccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGAATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCt                                                >  1:443126/1‑93 (MQ=255)

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

                                           atccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGGTTCCGCAGCATTACTGCCAGTGACTGTTACGCTgg                                              >  1:4694298/1‑93 (MQ=255)

                                             ccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAc                                            <  1:4310003/93‑1 (MQ=255)

                                             ccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAc                                            >  2:5012802/1‑93 (MQ=255)

                                                 ttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatca                                        >  2:2751148/1‑93 (MQ=255)

                                                 ttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatca                                        >  2:1348971/1‑93 (MQ=255)

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

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

                                                   aaaGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGGTTCCGCAGCATTACTGCCAGTGATTGTGACGCTGGATATCAtt                                      >  1:5487667/1‑93 (MQ=255)

                                                     aGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTAc                                    >  2:1882521/1‑93 (MQ=255)

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

                                                        gCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGc                                 >  1:1134049/1‑93 (MQ=255)

                                                          gATTTTATGAATACTATTCACCCAGAGGCATTCC‑CGATGGGGTTTCCACAGCATTACTGCCAAGGACTGTTACGGTGGACATCATTACTGCAc                               >  2:3084087/1‑93 (MQ=255)

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

                                                           aTTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCAcc                              <  2:5528456/93‑1 (MQ=255)

                                                           aTTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCAcc                              <  1:3791129/93‑1 (MQ=255)

                                                            ttttATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCg                             >  1:2855193/1‑93 (MQ=255)

                                                             tttATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGc                            <  2:265448/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                  gAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCaa                       <  1:1837843/93‑1 (MQ=255)

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

                                                                   aaTACTATTCACCCAGAGGCATTCC‑GGATGGTGATTCCGCCGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAt                      >  1:410251/1‑93 (MQ=255)

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

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

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

                                                                        tATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTg                 <  1:2829441/93‑1 (MQ=255)

                                                                        tATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTg                 <  1:3802621/93‑1 (MQ=255)

                                                                          ttCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGt               <  2:4305861/93‑1 (MQ=255)

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

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

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

                                                                               ccAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGtat          <  1:2228647/93‑1 (MQ=255)

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

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

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

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

                                                                                   agGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGc      <  1:4833889/93‑1 (MQ=255)

                                                                                   agGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGc      <  2:4158705/93‑1 (MQ=255)

                                                                                   agGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGc      <  1:1382042/93‑1 (MQ=255)

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

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

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

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

                                                                                       aTTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTa  >  1:4523624/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: `‑`