Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	769,691	C→T	P572S (CCG→TCG)	mngB →	alpha‑mannosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	769,691	0	C	T	98.5%	246.4 / ‑5.2	67	P572S (CCG→TCG)	mngB	alpha‑mannosidase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (33/33);  total (33/34)
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.71e-01

TGAAGAAAGCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGG  >  NC_000913/769543‑769838                                                                                                                                                     |                                                                                                                                                    tGAAGAAAGCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCg                                                                                                                                                   >  2:396387/1‑151 (MQ=255)  gaagaaAGCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGc                                                                                                                                                  >  2:924632/1‑151 (MQ=255)      aaaGCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCaa                                                                                                                                              <  1:788586/151‑1 (MQ=255)         gCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCGGTGTCGCTCAAtt                                                                                                                                            >  1:696093/1‑150 (MQ=255)             tgatgaTGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTc                                                                                                                                        >  1:349501/1‑150 (MQ=255)               atgatgGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGa                                                                                                                                     <  2:850864/151‑1 (MQ=255)                  atgGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACg                                                                                                                                  <  1:841725/151‑1 (MQ=255)                  atgGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACg                                                                                                                                  >  1:144772/1‑151 (MQ=255)                    ggTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCa                                                                                                                                <  1:777718/151‑1 (MQ=255)                       gatgaATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAgcg                                                                                                                             <  1:161750/151‑1 (MQ=255)                       gatgaATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAgcg                                                                                                                             <  1:4744/151‑1 (MQ=255)                              atGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGc                                                                                                                       >  1:288142/1‑150 (MQ=255)                                    aTTCACCCGCAAATGAAGAGTGGGCAATTACCGCAGCGCACGCGAAACCGCAATGCGATATTATTCGTGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACACGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCAGGCAATCCAc                                                                                                                >  2:726508/1‑151 (MQ=255)                                    aTTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATcc                                                                                                                  <  1:629955/149‑1 (MQ=255)                                      tCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACt                                                                                                               <  2:308936/150‑1 (MQ=255)                                        aCCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGc                                                                                                            >  1:1017685/1‑151 (MQ=255)                                              aaaaGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCag                                                                                                          >  2:419687/1‑147 (MQ=255)                                              aaaaGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTa                                                                                                      <  1:48800/151‑1 (MQ=255)                                               aaaGCATAGTGGGTAATTACCGACGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGACCAGGGCTGTTATCCGCCATGACATCTCAGTGTCGCTCAATTTGTCAGAACGCACCGCCCGGCAATCCACTgcccgcg                                                                                                        >  2:764711/1‑142 (MQ=255)                                                aagaagAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAgg                                                                                                    >  1:647263/1‑151 (MQ=255)                                                   aagaGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAgggg                                                                                                  >  1:829900/1‑150 (MQ=255)                                                          ggTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTgttg                                                                                          >  1:446875/1‑151 (MQ=255)                                                             aaTTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTgttggt                                                                                        <  1:1063036/150‑1 (MQ=255)                                                              aTTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTgttggttg                                                                                      <  1:232913/151‑1 (MQ=255)                                                                tACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTc                                                                                    >  2:177684/1‑151 (MQ=255)                                                                 aCCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTc                                                                                    >  1:632342/1‑150 (MQ=255)                                                                   cGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCAct                                                                                 <  1:481410/151‑1 (MQ=255)                                                                   cGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCAct                                                                                 >  1:1017741/1‑151 (MQ=255)                                                                        cGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCtt                                                                              >  2:358830/1‑149 (MQ=255)                                                                        cGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAg                                                                            <  2:981066/151‑1 (MQ=255)                                                                         gAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAg                                                                            >  1:765709/1‑150 (MQ=255)                                                                         gAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGt                                                                           >  1:556168/1‑151 (MQ=255)                                                                         gAACGCGAAACCGCAATGCGAAATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGt                                                                           <  1:755853/151‑1 (MQ=255)                                                                          aaCGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGt                                                                           <  2:365244/150‑1 (MQ=255)                                                                           aagcgAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTATCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGt                                                                           >  1:333794/3‑149 (MQ=255)                                                                            cgcgAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCAt                                                                        <  1:238197/151‑1 (MQ=255)                                                                            cgcgAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCAt                                                                        <  1:376655/151‑1 (MQ=255)                                                                             gcgAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATa                                                                       <  2:589050/151‑1 (MQ=255)                                                                                aaaCCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAAc                                                                     >  2:576667/1‑150 (MQ=255)                                                                                    cGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAAcagcag                                                                <  2:427660/151‑1 (MQ=255)                                                                                    cGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAAcagc                                                                  >  2:756110/1‑149 (MQ=255)                                                                                     gCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAAcagcag                                                                <  2:954178/150‑1 (MQ=255)                                                                                       aaTGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTg                                                                                                              <  1:966938/102‑1 (MQ=255)                                                                                       aaTGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTg                                                                                                              >  2:966889/1‑102 (MQ=255)                                                                                       aaTGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCg                                                             >  1:51169/1‑151 (MQ=255)                                                                                       aaTGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCg                                                             <  1:945963/151‑1 (MQ=255)                                                                                       aaTGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCg                                                             <  1:27191/151‑1 (MQ=255)                                                                                          gCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTa                                                           <  2:275798/150‑1 (MQ=255)                                                                                          gCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTAt                                                          >  1:100233/1‑151 (MQ=255)                                                                                          gCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTAt                                                          <  1:609011/151‑1 (MQ=255)                                                                                             atatTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGGTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTg                                                        >  2:465212/1‑150 (MQ=255)                                                                                             atatTATTCATGAAGCCTGGCAGAGCAGGGCTGTTAGCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTa                                                           <  2:22465/147‑1 (MQ=255)                                                                                                    tCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGa                                                 <  2:502848/150‑1 (MQ=255)                                                                                                    tCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGa                                                 >  2:659541/1‑150 (MQ=255)                                                                                                          aGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAAt                                          >  1:154811/1‑151 (MQ=255)                                                                                                           gCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATc                                         >  2:387314/1‑151 (MQ=255)                                                                                                                 cAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGAt                                    >  1:928059/1‑150 (MQ=255)                                                                                                                      cAGGGCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAAcc                                >  1:1048882/1‑149 (MQ=255)                                                                                                                          gCTGTTATCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTg                          <  1:662107/151‑1 (MQ=255)                                                                                                                                aTCCGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGatc                    >  2:582994/1‑151 (MQ=255)                                                                                                                                  ccGCTATGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGatca                   <  1:262406/150‑1 (MQ=255)                                                                                                                                       aTGACATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCt             <  1:1017417/151‑1 (MQ=255)                                                                                                                                        tGACATGGCAGTGTCGCGCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGGTGACGATCATCGcc              <  2:972399/149‑1 (MQ=255)                                                                                                                                           cATGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCgt         <  1:516044/151‑1 (MQ=255)                                                                                                                                            aTGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCTCTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCgt         <  2:542775/150‑1 (MQ=255)                                                                                                                                            aTGGCAGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCgtg        <  2:136886/151‑1 (MQ=255)                                                                                                                                                 aGTGTCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTg   >  2:662280/1‑151 (MQ=255)                                                                                                                                                  gtgtCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTgg  <  2:749710/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGAAGAAAGCTCTGATGATGGTGATGAATATGACTATTCACCCGCAAAAGAAGAGTGGGTAATTACCGCAGCGAACGCGAAACCGCAATGCGATATTATTCATGAAGCCTGGCAGAGCAGGGCTGTTATCCGCTATGACATGGCAGTGCCGCTCAATTTGTCAGAACGCAGCGCCCGGCAATCCACTGGCAGAGTAGGGGTGGTGTTGGTTGTCACTCTTAGTCATAACAGCAGGCGTATTGATGTGGATATCAATCTTGATAACCAGGCTGACGATCATCGCCTTCGTGTCCTGG  >  NC_000913/769543‑769838

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