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.7%	1354.6 / ‑3.5	287	coding (504/672 nt)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base G (1/0); new base . (176/110); total (177/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 300 total aligned reads displayed.

AAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTA  >  NC_000913/3873534‑3873713

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

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

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

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

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

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

        aaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGtt                                                                                >  2:18832/1‑93 (MQ=255)

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

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

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

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

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

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

               aTTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgcag                                                                         <  1:3694986/93‑1 (MQ=255)

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

               aTTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgcag                                                                         <  2:5075306/93‑1 (MQ=255)

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

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

                   cTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCAtt                                                                     <  2:5395354/93‑1 (MQ=255)

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

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

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

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

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

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

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

                          gTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGCCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCa                                                              <  2:3013520/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                      tACCATCCATTTAAAGTGCGATATTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACg                                                  >  1:752246/1‑93 (MQ=255)

                                       aCCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGc                                                 >  1:2485343/1‑93 (MQ=255)

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

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

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

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

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

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

                                             catTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAca                                           <  2:440568/93‑1 (MQ=255)

                                             catTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAca                                           <  1:2016366/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                   aaGTGCGATTTTATGAATACCATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTa                                     <  1:6102767/93‑1 (MQ=255)

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

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

                                                      tGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTg                                  <  2:2525378/93‑1 (MQ=255)

                                                      tGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTg                                  <  2:2827320/93‑1 (MQ=255)

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

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

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

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

                                                           ttttATGAATACTATTCACTCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCg                             >  1:6502328/1‑93 (MQ=255)

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

                                                             ttATGAATACTATTCTCCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCt                           <  2:2906475/93‑1 (MQ=255)

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

                                                              tATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTg                          <  1:1599509/93‑1 (MQ=255)

                                                              tATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTCCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTg                          <  2:4788978/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                           cACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTg             >  1:6602588/1‑93 (MQ=255)

                                                                            aCCCAGAGGCATTCC‑GTATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACTCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGt            <  2:6104250/93‑1 (MQ=255)

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

                                                                            aCCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGt            <  2:458964/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                    gCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGAc                                                         <  2:1182245/40‑1 (MQ=255)

                                                                                      aTTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTa  >  1:2671763/1‑93 (MQ=255)

                                                                                      aTTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTa  >  2:3572809/1‑93 (MQ=255)

AAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTA  >  NC_000913/3873534‑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: `‑`