Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,198,900	G→A	6.0%	G37G (GGG→GGA)	nfi →	endonuclease V; deoxyinosine 3' endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,198,900	0	G	A	6.0%	214.4 / 2.7	67	G37G (GGG→GGA)	nfi	endonuclease V; deoxyinosine 3' endonuclease
Reads supporting (aligned to +/- strand):  ref base G (36/27);  new base A (2/2);  total (38/29)
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.11e-01

ATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGC  >  NC_000913/4198752‑4199048                                                                                                                                                      |                                                                                                                                                     aTCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTAAGCGCTCAACAAA‑TCGAACTGGCTTCTTATGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCggg                                                                                                                                                      >  3:5708/1‑149 (MQ=255)    gACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGCTCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGttt                                                                                                                                                   <  1:173048/149‑1 (MQ=255)     aCTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTg                                                                                                                                                  <  2:356773/149‑1 (MQ=255)      cTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGttt                                                                                                                                                   >  7:323259/1‑147 (MQ=255)      cTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 <  5:307827/149‑1 (MQ=255)      cTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 <  8:323259/149‑1 (MQ=255)         tCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAg                                                                                                                                                >  1:130173/1‑147 (MQ=255)         tCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGATTTGAg                                                                                                                                                >  7:275981/1‑147 (MQ=255)         tCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGATTTGAg                                                                                                                                                >  7:281163/1‑147 (MQ=255)           tGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAgg                                                                                                                                            >  7:376499/1‑149 (MQ=255)           tGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAgg                                                                                                                                            >  2:47739/1‑149 (MQ=255)                      aCCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGt                                                                                                                                  >  5:329586/1‑148 (MQ=255)                      aCCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGt                                                                                                                                  >  8:306058/1‑148 (MQ=255)                      aCCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTg                                                                                                                                 <  4:117874/149‑1 (MQ=255)                             aGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGATTTGAGCAGGGCGGAGAAGTGAcgcg                                                                                                                            <  8:275981/147‑1 (MQ=255)                             aGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGATTTGAGCAGGGCGGAGAAGTGAcgcg                                                                                                                            <  8:281163/147‑1 (MQ=255)                              ggAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAgc                                                                                                                         >  1:289550/1‑149 (MQ=255)                                    aTTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGAt                                                                                                                   <  8:150091/149‑1 (MQ=255)                                     ttATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                                                                    <  8:109146/131‑1 (MQ=255)                                     ttATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                                                                    >  7:109146/1‑131 (MQ=255)                                     ttATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATg                                                                                                                  <  2:289550/149‑1 (MQ=255)                                        tGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGttt                                                                                                                                                   <  2:311306/113‑1 (MQ=255)                                        tGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGttt                                                                                                                                                   >  1:311306/1‑113 (MQ=255)                                        tGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                                                                    >  2:180543/1‑128 (MQ=255)                                        tGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                                                                    <  1:180543/128‑1 (MQ=255)                                         ggATCTCGCGTCATTACGCGCTCACAAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGtgc                                                                                                              >  5:28516/1‑149 (MQ=255)                                         ggATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 <  4:303690/114‑1 (MQ=255)                                         ggATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 >  3:303690/1‑114 (MQ=255)                                           aTCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGtgctg                                                                                                            >  7:275073/1‑149 (MQ=255)                                              tcGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 >  8:103714/1‑109 (MQ=255)                                              tcGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 <  7:103714/109‑1 (MQ=255)                                              tcGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTgctgct                                                                                                          >  6:24257/1‑148 (MQ=255)                                                 cgTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATg                                                                                                                  >  5:263057/1‑137 (MQ=255)                                                 cgTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATg                                                                                                                  <  6:263057/137‑1 (MQ=255)                                                        aCGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATccc                                                                                                <  7:34008/148‑1 (MQ=255)                                                         cgcgCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCt                                                                                               <  7:306058/148‑1 (MQ=255)                                                                aacaaA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAg                                                                                       >  1:235392/1‑149 (MQ=255)                                                                aacaaA‑TCGAACTGGCTTCTTCTGTGATCAGCGAGGATCGACTCGATAAAGATACACCGGATCTGATCTCCGGAGCCGATGTCGGGTTTGAGCAGGGCGCAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGc                                                                                            >  5:89990/1‑144 (MQ=255)                                                                 acaaATTCGAACTGGCTTGTTCTGTGCTGCGCGAGGATCGACTCGATCTAGTTCGACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGGTGAAATATCCCTCGc                                                                                            <  6:89990/144‑1 (MQ=255)                                                                  caaA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCt                                                                                     >  7:51517/1‑149 (MQ=255)                                                                         gAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGa                                                                               <  6:28516/149‑1 (MQ=255)                                                                           aCTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGt                                                                             <  6:329586/149‑1 (MQ=255)                                                                            cTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCAACCTGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGt                                                                             <  4:5708/148‑1 (MQ=255)                                                                             tGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGta                                                                            >  4:59716/1‑148 (MQ=255)                                                                               gCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGCAATATCCCTCGCTTGAGCTGGTc                                                                                 >  7:180151/1‑141 (MQ=255)                                                                               gCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGCAATATCCCTCGCTTGAGCTGGTc                                                                                 <  8:180151/141‑1 (MQ=255)                                                                               gCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGtata                                                                          <  5:24257/148‑1 (MQ=255)                                                                                      ctGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGtt                                                                     <  8:275073/146‑1 (MQ=255)                                                                                      ctGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGcc                                                                  >  5:92687/1‑149 (MQ=255)                                                                                         tgATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATaaa                                                                        >  4:328362/1‑140 (MQ=255)                                                                                         tgATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATaaa                                                                        <  3:328362/140‑1 (MQ=255)                                                                                              cgcgAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGc                                                          >  5:136102/1‑149 (MQ=255)                                                                                               gcgAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGcc                                                         >  2:130790/1‑149 (MQ=255)                                                                                                   ggATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccacca                                                     >  2:321565/1‑149 (MQ=255)                                                                                                                 aaGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACAtt                                       >  7:122097/1‑149 (MQ=255)                                                                                                                  aGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTc                                      >  8:124310/1‑149 (MQ=255)                                                                                                                   gATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTc                                      >  6:209267/1‑148 (MQ=255)                                                                                                                    aTCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCa                                    <  3:59716/149‑1 (MQ=255)                                                                                                                    aTCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCa                                    <  8:122097/149‑1 (MQ=255)                                                                                                                     tCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCt                                                                                               >  8:68831/1‑89 (MQ=255)                                                                                                                      ccACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAgg                                  >  2:127981/1‑149 (MQ=255)                                                                                                                       cACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGt                                                                             <  7:18634/105‑1 (MQ=255)                                                                                                                       cACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGt                                                                             >  8:18634/1‑105 (MQ=255)                                                                                                                       cACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGt                                 <  7:124310/149‑1 (MQ=255)                                                                                                                        aCCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGtt                                >  1:349554/1‑149 (MQ=255)                                                                                                                                 gATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTttcct                       >  4:339845/1‑149 (MQ=255)                                                                                                                                              gATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATcc          >  2:314709/1‑149 (MQ=255)                                                                                                                                              gATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATcc          >  2:263757/1‑149 (MQ=255)                                                                                                                                                    gggTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCgct    >  2:75616/1‑149 (MQ=255)                                                                                                                                                      gTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCgctgc  <  5:209267/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     ATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAA‑TCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGC  >  NC_000913/4198752‑4199048

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑