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	.	98.8%	1111.9 / 4.2	246	coding (504/672 nt)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base . (129/114); minor base A (2/1); total (131/115)
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 269 total aligned reads displayed.

AAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTA  >  NC_000913/3873533‑3873713

aaaaaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCg                                                                                          <  2:3712574/93‑1 (MQ=255)

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

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

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

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

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

     tAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatgg                                                                                     <  2:3903108/93‑5 (MQ=255)

      aaaaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgt                                                                                   >  2:2672050/1‑93 (MQ=255)

      aaaaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgt                                                                                   >  1:1693605/1‑93 (MQ=255)

        aaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCATCCAGAGGCATTCC‑GGATGgtgt                                                                                 <  2:2990649/93‑1 (MQ=255)

        aaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgtgt                                                                                 <  1:290804/93‑1 (MQ=255)

        aaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATT‑CAGGATGgtgt                                                                                 >  2:808423/1‑93 (MQ=255)

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

          aCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGttt                                                                               <  2:709166/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                       aaaGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACt                                                                  <  2:263845/93‑1 (MQ=255)

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

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

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

                             gATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCTTTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGt                                                            >  2:1444542/1‑93 (MQ=255)

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

                               tGGCGTACTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCGCCCAGAGGCATTCC‑GTATGGTGTTTCCGCAGCATTACTGCCAGTGa                                                          <  2:3155799/93‑1 (MQ=255)

                                ggCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGAc                                                         <  1:260671/93‑1 (MQ=255)

                                ggCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGAc                                                         <  2:2923862/93‑1 (MQ=255)

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

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

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

                                       tACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACg                                                  <  1:3120467/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                               atTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCGGGAcgt                                          >  1:1092360/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                      gtccGATTTTATGAAGACTATTCACCCAGAGGCATTCC‑GGATGGTGTTCCCGTAGCATTACTGCCAGTGGCTGTTACTATAGTCATCATTACt                                   <  2:1896038/90‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                              cccAGAGGCATTCC‑GGATGGTGTTGCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGta           <  2:3480891/93‑1 (MQ=255)

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

                                                                               ccAGAGGCATTCC‑GGATGATGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGtat          >  2:3807552/1‑93 (MQ=255)

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

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

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

                                                                                  gagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATg       >  2:2538974/1‑93 (MQ=255)

                                                                                  gagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATg       >  2:1157783/1‑93 (MQ=255)

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

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

                                                                                       aTTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTa  <  1:896294/93‑1 (MQ=255)

AAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTA  >  NC_000913/3873533‑3873713

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