Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,855,020	C→T	6.9%	A347V (GCC→GTC)	tcyP →	sodium/cystine symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,855,020	0	C	T	6.9%	169.1 / 2.8	58	A347V (GCC→GTC)	tcyP	sodium/cystine symporter
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base T (2/2);  total (29/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.12e-01

CTGACCTTCGCCTTCACCAGCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTCACGCTGAGCTCTGCAGGGGTTGCCGGCGTGGGCGGCGGTGCCACCTT  >  NC_002947/5854878‑5855164                                                                                                                                               |                                                                                                                                                 cTGACCTTCGCCTTCACCAGCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGct                                                                                                                                            >  6:39718/1‑149 (MQ=255)  tGACCTTCGCCTTCACCAGCAGATCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGTCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCCGCAACACCATAGGCCAGAACTGATGCGTCGGGatc                                                                                                                                           >  2:152209/1‑149 (MQ=255)  tGACATTCGCATTCACCAGCCGCTCCAGCGCAGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTTCCGCAATCGATCGCCAGCTTCACCGCATCGTTCGGCAACACCATTGGCCAGAACGGAAGCGTCGGtctc                                                                                                                                           >  2:148108/1‑149 (MQ=255)   gACCTTCGCCTTCACCAGCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGctc                                                                                                                                           >  6:72551/1‑148 (MQ=255)    aCCTTCGCCTTCACCAGCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGcc                                                                                                                                                 <  5:104927/141‑1 (MQ=255)    aCCTTCGCCTTCACCAGCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGcc                                                                                                                                                 >  6:104927/1‑141 (MQ=255)      cTTCGCCTTCACCAGCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATc                                                                                                                                       <  8:224893/149‑1 (MQ=255)      cTTCGCCTTCACCAGCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATc                                                                                                                                       >  6:169599/1‑149 (MQ=255)       ttCGCCTTCACCAGCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATcc                                                                                                                                      >  2:14054/1‑149 (MQ=255)            cTTCACCAGCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCa                                                                                                                                 <  4:82582/149‑1 (MQ=255)                  cAGCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTgg                                                                                                                           >  3:243321/1‑149 (MQ=255)                    gCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCg                                                                                                                         >  4:163649/1‑149 (MQ=255)                     ccGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACCAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCCGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCgg                                                                                                                        >  3:154740/1‑149 (MQ=255)                      ggCTCCAGCGCCGACAGCTTTCCGCTGAACAGCGAGGCAGAATCCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCAGCGTTCGGCACCACCATTGGCGAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGt                                                                                                                       <  3:177513/148‑1 (MQ=255)                       gCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGt                                                                                                                       >  2:205071/1‑148 (MQ=255)                            aGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGt                                                                                                                 <  1:205071/149‑1 (MQ=255)                             gcgcCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTc                                                                                                                <  3:33977/149‑1 (MQ=255)                               gccgccAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTcgc                                                                                                              <  4:16032/149‑1 (MQ=255)                               gccgccAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTcgc                                                                                                              <  4:19596/149‑1 (MQ=255)                                 cgccAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCgcgc                                                                                                            <  1:150406/149‑1 (MQ=255)                                   ccAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCTGGCTCTATCCGGCCAt                                                                                                                                <  4:26657/127‑1 (MQ=255)                                   ccAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCAt                                                                                                                                >  3:26657/1‑127 (MQ=255)                                   ccAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCgcgc                                                                                                            >  3:41836/1‑147 (MQ=255)                                    cAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCcag                                                                                                         <  5:39718/149‑1 (MQ=255)                                    cAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCcag                                                                                                         <  1:100253/149‑1 (MQ=255)                                         ttCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTa                                                                                                    <  5:169599/149‑1 (MQ=255)                                          tCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCgcgc                                                                                                            >  5:174249/1‑140 (MQ=255)                                          tCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCgcgc                                                                                                            <  6:174249/140‑1 (MQ=255)                                            cGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAgg                                                                                                  >  5:3919/1‑148 (MQ=255)                                             gCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCa                                                                                                >  1:167381/1‑149 (MQ=255)                                                    aTCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGacac                                                                                         <  4:41836/149‑1 (MQ=255)                                                        aGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTc                                                                                     >  5:93134/1‑149 (MQ=255)                                                        aGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTc                                                                                     <  1:14054/149‑1 (MQ=255)                                                          gCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGa                                                                                   <  4:238781/149‑1 (MQ=255)                                                          gCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACAGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGa                                                                                   >  5:193642/1‑149 (MQ=255)                                                           cacaAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGt                                                                                                     >  1:156438/1‑130 (MQ=255)                                                           cacaAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGt                                                                                                     <  2:156438/130‑1 (MQ=255)                                                             caAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGAccc                                                                                >  5:159349/1‑149 (MQ=255)                                                                aCCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCt                                                                             >  2:219116/1‑149 (MQ=255)                                                                               gTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCt                                                              <  3:163649/149‑1 (MQ=255)                                                                                tGCCGCAATCGATGGCCAGCTTCAGCTCATCGTTCGGCAACACGATTGGCCAGAACGGCTGCGTCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCTCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCtg                                                             <  1:148108/149‑1 (MQ=255)                                                                                tGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCAACACCATTGTCCAGAACGGCTGCGTCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCtg                                                             <  1:152209/149‑1 (MQ=255)                                                                                    gCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCtggtgg                                                         <  5:72551/149‑1 (MQ=255)                                                                                     cAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCa                                                                  <  8:102826/139‑1 (MQ=255)                                                                                     cAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCa                                                                  >  7:102826/1‑139 (MQ=255)                                                                                       atcgatcgCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCa                                                      >  3:201053/1‑149 (MQ=255)                                                                                       atcgatcgCCAGCTCCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGc                                                        <  5:11855/147‑1 (MQ=255)                                                                                       atcgatcgCCAGCTCCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGc                                                        >  6:11855/1‑147 (MQ=255)                                                                                        tcgatcgCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCAt                                                     >  7:107519/1‑149 (MQ=255)                                                                                          gatcgCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATc                                                    <  2:167381/148‑1 (MQ=255)                                                                                           atcgCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGt                                                  >  7:35952/1‑149 (MQ=255)                                                                                            tcgCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTc                                                 <  4:243321/149‑1 (MQ=255)                                                                                            tcgCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTc                                                 <  2:213941/149‑1 (MQ=255)                                                                                              gCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTCa                                                <  4:154740/148‑1 (MQ=255)                                                                                                 aGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTCACGCt                                            <  8:35952/149‑1 (MQ=255)                                                                                                                ttCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTCACGCTGAGCTCTGCAGGGGt                             >  1:234667/1‑149 (MQ=255)                                                                                                                      accaccATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTCACGCTGAGCTCTGCAGGGGTTGCCgg                       >  3:105477/1‑149 (MQ=255)                                                                                                                          ccaTTGGCCAGAACGGCTGCGCCGGGCTCTATACGGCGATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTTGATCGCCACCCTGGTGGCCTTCGTCACGCTGAGCTCTGCAGGGGTTGCCGGCGTg                   <  6:193642/149‑1 (MQ=255)                                                                                                                                           tGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTCACGCTGAGCTCTGCAGGGGTTGCCGGCGTGGGCGGCGGTGCCACCtt  >  6:4361/1‑149 (MQ=255)                                                                                                                                           tGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTCACGCTGAGCTCTGCAGGGGTTGCCGGCGTGGGCGGCGGTGCCACCtt  >  3:96533/1‑149 (MQ=255)                                                                                                                                               ccGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTCACGCTGAGCTCTGCAGGGGTTGCCGGCGTGGGCGGCGGt          <  1:203298/137‑1 (MQ=255)                                                                                                                                               ccGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTCACGCTGAGCTCTGCAGGGGTTGCCGGCGTGGGCGGCGGt          >  2:203298/1‑137 (MQ=255)                                                                                                                                               |                                                                                                                                                 CTGACCTTCGCCTTCACCAGCCGCTCCAGCGCCGCCAGCATTCCGCTGAACATCGAGGCACAAACCCGCCGTCTGGGCGTGCCGCAATCGATCGCCAGCTTCAGCGCATCGTTCGGCACCACCATTGGCCAGAACGGCTGCGCCGGGCTCTATCCGGCCATGCTGGCGGTGATGGTCGCGCCAGCAGTAGGCATCGACACCTTCGACCCGCTGTGGATCGCCACCCTGGTGGCCATCGTCACGCTGAGCTCTGCAGGGGTTGCCGGCGTGGGCGGCGGTGCCACCTT  >  NC_002947/5854878‑5855164

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