Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,054,193:1	+T	6.5%	coding (767/1320 nt)	PP_3573 →	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,054,193	1	.	T	6.5%	249.6 / 9.8	62	coding (767/1320 nt)	PP_3573	monooxygenase
Reads supporting (aligned to +/- strand):  ref base . (28/30);  new base T (2/2);  total (30/32)
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.23e-01

ACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCT  >  NC_002947/4054046‑4054337                                                                                                                                                      |                                                                                                                                                 aCCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAggg                                                                                                                                                   >  4:337144/1‑148 (MQ=255)  ccGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAggg                                                                                                                                                   <  4:243165/147‑1 (MQ=255)  ccGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑c                                                                                                                                                 <  6:30498/148‑1 (MQ=255)        ccGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGAccc                                                                                                                                          <  3:337144/149‑1 (MQ=255)         cgcgACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCt                                                                                                                                         >  8:87948/1‑149 (MQ=255)          gcgTCTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTg                                                                                                                                        <  6:357319/149‑1 (MQ=255)              cTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCg                                                                                                                                    <  7:22455/149‑1 (MQ=255)                tCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGGTCATG‑CCCTGGtgg                                                                                                                                    <  4:92052/147‑3 (MQ=255)                tCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGGTCATG‑CCCTGGtgg                                                                                                                                    >  3:92052/1‑145 (MQ=255)                   gctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCt                                                                                                                                >  8:276342/1‑148 (MQ=255)                          gCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGGTCATG‑CCCTGGTGGGc                                                                                                                                  >  1:147015/1‑139 (MQ=255)                          gCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGGTCATG‑CCCTGGTGGGc                                                                                                                                  <  2:147015/139‑1 (MQ=255)                           cTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTc                                                                                                                       <  7:276342/149‑1 (MQ=255)                            tGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCg                                                                                                                      >  6:358703/1‑149 (MQ=255)                              ggCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGc                                                                                                                                  >  2:247815/1‑135 (MQ=255)                              ggCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGtt                                                                                                                    <  7:149210/149‑1 (MQ=255)                              ggCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGt                                                                                                                     <  6:183043/148‑1 (MQ=255)                              ggCGTGCGGGGCGAGTGGGAAGCGGAAATTGCCAAGCCCGGTCGCGATCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGCTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGc                                                                                                                                  <  1:247815/136‑1 (MQ=255)                                    cTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGaaaa                                                                                                              >  3:320000/1‑149 (MQ=255)                                     tGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTg                                                                                                                                        >  2:191578/1‑122 (MQ=255)                                     tGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTg                                                                                                                                        <  1:191578/122‑1 (MQ=255)                                     tGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAAc                                                                                                             >  7:238367/1‑149 (MQ=255)                                       ggCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACg                                                                                                            <  7:87948/148‑1 (MQ=255)                                            gTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGA‑CATCGCGGTCAGAGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTcgc                                                                                                     >  1:92807/1‑149 (MQ=255)                                                gACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAAc                                                                                                             >  1:143717/1‑138 (MQ=255)                                                gACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAAc                                                                                                             <  2:143717/138‑1 (MQ=255)                                                     ggAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTcc                                                                                             <  4:320000/149‑1 (MQ=255)                                                     ggAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTc                                                                                              >  4:281619/1‑148 (MQ=255)                                                      gAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCa                                                                                            <  2:152103/149‑1 (MQ=255)                                                       aaaTTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAg                                                                                           <  5:358703/149‑1 (MQ=255)                                                              cAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAAcc                                                                                    >  3:187633/1‑149 (MQ=255)                                                               aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTc                                                                                              <  1:179493/138‑1 (MQ=255)                                                               aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTc                                                                                              >  2:179493/1‑138 (MQ=255)                                                                aGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTg                                                                                  >  8:56691/1‑149 (MQ=255)                                                                 gCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTgg                                                                                 <  2:183011/149‑1 (MQ=255)                                                                   ccGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TTGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTc                                                                               <  6:3565/149‑1 (MQ=255)                                                                     ggTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCACGTGGTGGCCACAACG‑TGGACTTACGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAAACAGTCGTTCGAAAAc                                                                                                             >  5:136101/1‑117 (MQ=255)                                                                      gTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCgag                                                                            <  4:187633/149‑1 (MQ=255)                                                                           gAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACAt                                                                       <  3:34764/149‑1 (MQ=255)                                                                            aGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACAt                                                                       >  5:322696/1‑148 (MQ=255)                                                                                     ccATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                             <  3:281619/149‑1 (MQ=255)                                                                                        tCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcgac                                                          >  8:235554/1‑149 (MQ=255)                                                                                                  tGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAcc                                                >  5:92662/1‑149 (MQ=255)                                                                                                     cgcCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTgg                                             <  2:92807/149‑1 (MQ=255)                                                                                                          gtggtgGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGc                                            <  3:326169/145‑1 (MQ=255)                                                                                                          gtggtgGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGc                                            >  4:326169/1‑145 (MQ=255)                                                                                                                  cacaCCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcga                                                           <  4:137338/122‑1 (MQ=255)                                                                                                                  cacaCCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGc                                >  3:308003/1‑149 (MQ=255)                                                                                                                    cacCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATc                                                                      >  6:2051/1‑109 (MQ=255)                                                                                                                    cacCG‑TGGACTTCCGTGCCCTAGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATc                                                                      <  5:2051/109‑1 (MQ=255)                                                                                                                          tGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACg                                                                                                            >  4:323424/1‑66 (MQ=255)                                                                                                                          tGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGAt                         <  6:92662/149‑1 (MQ=255)                                                                                                                              cTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTcc                                                                                             >  6:53656/1‑77 (MQ=255)                                                                                                                                tCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACGTGGACGAGACCAGCAAACCGGGCGAGGACACGAACCTGGCCCTCGGGGTTGCGCCCGAGGCGTa                    >  4:210381/1‑148 (MQ=255)                                                                                                                                 ccGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGa                          >  6:254994/1‑141 (MQ=255)                                                                                                                                 ccGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGa                          <  5:254994/141‑1 (MQ=255)                                                                                                                                   gTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATc                <  6:322696/149‑1 (MQ=255)                                                                                                                                         tCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGAt                                  <  8:332112/125‑1 (MQ=255)                                                                                                                                         tCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGAt                                  >  7:332112/1‑125 (MQ=255)                                                                                                                                         tCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGaa             <  5:251698/146‑1 (MQ=255)                                                                                                                                         tCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGaa             >  6:251698/1‑146 (MQ=255)                                                                                                                                         tCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACg           >  5:157006/1‑148 (MQ=255)                                                                                                                                              caccaGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACgg     <  7:235554/149‑1 (MQ=255)                                                                                                                                                 caGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGcc                                           <  6:364281/108‑1 (MQ=255)                                                                                                                                                 caGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCt  >  3:141094/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                 ACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCG‑TGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCT  >  NC_002947/4054046‑4054337

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: ‑