Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,793,850	(C)6→7	coding (260/552 nt)	btuE ←	glutathione peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,793,844	1	.	C	98.9%	342.9 / ‑3.3	90	coding (266/552 nt)	btuE	glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (37/52);  total (37/53)
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.99e-01

GCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑A‑CCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATAT  >  NC_000913/1793698‑1793984                                                                                                                                                     |                                                                                                                                             gCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑A‑cccc                                                                                                                                          >  4:446460/1‑151 (MQ=255)  cGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑A‑ccccc                                                                                                                                         >  2:807850/1‑151 (MQ=255)   ggTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑A‑cccccc                                                                                                                                        >  2:1126452/1‑151 (MQ=255)          ggCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATgtgg                                                                                                                                  <  2:255843/151‑1 (MQ=255)           gCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATgtggt                                                                                                                                 <  1:251256/151‑1 (MQ=255)               ttGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGtggtgg                                                                                                                               <  2:580820/149‑1 (MQ=255)               ttGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTggtggt                                                                                                                              >  2:835162/1‑150 (MQ=255)                tGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTa                                                                                                                             <  1:501494/150‑1 (MQ=255)                 gCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTAc                                                                                                                            <  1:206101/150‑1 (MQ=255)                    gACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGt                                                                                                                         <  2:476762/150‑1 (MQ=255)                    gACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTa                                                                                                                        <  2:303458/151‑1 (MQ=255)                        aTACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGtt                                                                                                                    <  1:1131823/151‑1 (MQ=255)                          aCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGtttt                                                                                                                  <  3:117056/151‑1 (MQ=255)                           cGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGtttt                                                                                                                  >  1:854032/1‑150 (MQ=255)                            gggCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTaa                                                                                                                <  2:842519/151‑1 (MQ=255)                            gggCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTaa                                                                                                                <  2:653234/151‑1 (MQ=255)                              gCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAAtc                                                                                                              <  2:454465/151‑1 (MQ=255)                               cATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAAtct                                                                                                             <  1:1195563/151‑1 (MQ=255)                               cATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAAtct                                                                                                             <  1:940636/151‑1 (MQ=255)                                  aGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATctct                                                                                                           <  3:254962/150‑1 (MQ=255)                                  aGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATctct                                                                                                           <  3:50614/150‑1 (MQ=255)                                    aaTCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCa                                                                                                        <  4:365955/151‑1 (MQ=255)                                     aTCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCAt                                                                                                       <  1:399697/151‑1 (MQ=255)                                       ccGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCg                                                                                                     <  2:199996/151‑1 (MQ=255)                                       ccGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCg                                                                                                     <  3:161208/151‑1 (MQ=255)                                       ccGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCg                                                                                                     <  4:433893/151‑1 (MQ=255)                                        cGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGc                                                                                                    <  1:955516/151‑1 (MQ=255)                                             tctTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGccc                                                                                               >  1:707068/1‑151 (MQ=255)                                             tctTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGccc                                                                                               >  1:167325/1‑151 (MQ=255)                                              ctTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGcc                                                                                                <  1:376427/149‑1 (MQ=255)                                              ctTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGc                                                                                                 <  1:389263/148‑1 (MQ=255)                                                tCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCg                                                                                              <  4:24151/149‑1 (MQ=255)                                                 ccGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGt                                                                                            <  4:260508/150‑1 (MQ=255)                                                      gcgACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGt                                                                                      >  2:426450/1‑151 (MQ=255)                                                         aCTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTcc                                                                                   >  1:1065834/1‑151 (MQ=255)                                                                 cGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACt                                                                            <  1:359090/150‑1 (MQ=255)                                                                  ggCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTgg                                                                          >  1:584782/1‑151 (MQ=255)                                                                    cgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGtt                                                                        >  3:143486/1‑151 (MQ=255)                                                                    cgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGtt                                                                        <  2:160225/151‑1 (MQ=255)                                                                     gcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTg                                                                       >  2:362669/1‑151 (MQ=255)                                                                       tGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGc                                                                      >  1:880524/1‑150 (MQ=255)                                                                       tGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGc                                                                      >  4:270581/1‑150 (MQ=255)                                                                        gCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGc                                                                      <  1:877161/149‑1 (MQ=255)                                                                         cGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACg                                                                   <  1:444907/151‑1 (MQ=255)                                                                         cGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACg                                                                   >  4:6358/1‑151 (MQ=255)                                                                            caatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGa                                                                <  1:447368/151‑1 (MQ=255)                                                                             aatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGaa                                                               <  3:248184/151‑1 (MQ=255)                                                                             aatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGaa                                                               <  1:924724/151‑1 (MQ=255)                                                                              atcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAAt                                                              >  4:196903/1‑151 (MQ=255)                                                                               tcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATc                                                             >  2:1267020/1‑151 (MQ=255)                                                                                caatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATcc                                                            <  3:69832/151‑1 (MQ=255)                                                                                     tttGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGc                                                        >  2:495724/1‑150 (MQ=255)                                                                                       tGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCAc                                                      <  1:847849/150‑1 (MQ=255)                                                                                        gATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCAcc                                                     >  2:578082/1‑150 (MQ=255)                                                                                        gATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCAcc                                                     >  4:354154/1‑150 (MQ=255)                                                                                            cAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaa                                                <  1:616654/151‑1 (MQ=255)                                                                                            cAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaa                                                 >  2:642859/1‑150 (MQ=255)                                                                                             aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                               <  2:1265109/151‑1 (MQ=255)                                                                                             aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                               <  2:1102173/151‑1 (MQ=255)                                                                                             aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                               <  2:579892/151‑1 (MQ=255)                                                                                               cGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAAcc                                             <  2:685185/151‑1 (MQ=255)                                                                                                ggATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAAcc                                             >  2:232312/1‑150 (MQ=255)                                                                                                ggATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCt                                            >  2:623162/1‑151 (MQ=255)                                                                                                 gATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtc                                           <  2:711853/151‑1 (MQ=255)                                                                                                 gATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtc                                           <  2:796664/151‑1 (MQ=255)                                                                                                  aTGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcg                                          <  1:621830/151‑1 (MQ=255)                                                                                                  aTGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcg                                          >  2:113215/1‑151 (MQ=255)                                                                                                    ggCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                        >  1:522011/1‑151 (MQ=255)                                                                                                    ggCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcga                                         <  1:1270703/150‑1 (MQ=255)                                                                                                     gCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatc                                       <  3:143066/151‑1 (MQ=255)                                                                                                      cGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatcg                                      <  3:114592/151‑1 (MQ=255)                                                                                                        tCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGAc                                    >  3:411314/1‑151 (MQ=255)                                                                                                             cGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAgg                               >  2:578219/1‑151 (MQ=255)                                                                                                             cGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAg                                <  2:706923/150‑1 (MQ=255)                                                                                                               ccATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCt                             <  1:143864/151‑1 (MQ=255)                                                                                                               ccATTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCt                             >  1:986912/1‑151 (MQ=255)                                                                                                               caaTTAACTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGAGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCt                             >  3:352321/3‑151 (MQ=255)                                                                                                                cATTAACTTCAATCTTACTGAACATCGGGAACGTCTC‑CCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCCCCATAAAACCTCGATCGACCCAGGCtt                            <  2:1051525/151‑1 (MQ=255)                                                                                                                    aaCTTCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTg                        >  1:248208/1‑151 (MQ=255)                                                                                                                       ttCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                     >  3:163519/1‑151 (MQ=255)                                                                                                                        tCAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                      <  2:133424/149‑1 (MQ=255)                                                                                                                         cAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAatat                   >  1:1101736/1‑151 (MQ=255)                                                                                                                         cAATCTTACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAata                    <  4:405147/150‑1 (MQ=255)                                                                                                                              ttACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTcc              >  4:121007/1‑151 (MQ=255)                                                                                                                              ttACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTcc              >  2:769059/1‑151 (MQ=255)                                                                                                                               tACTGAACATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTcc              <  2:974532/150‑1 (MQ=255)                                                                                                                                     aCATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCt       >  3:321255/1‑151 (MQ=255)                                                                                                                                     aCATCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCt       >  2:243196/1‑151 (MQ=255)                                                                                                                                       aTCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTc      >  2:147353/1‑150 (MQ=255)                                                                                                                                       aTCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCa     <  3:272367/151‑1 (MQ=255)                                                                                                                                        tCGGGAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCat    >  2:1280246/1‑151 (MQ=255)                                                                                                                                          gggAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCatat  >  2:311291/1‑151 (MQ=255)                                                                                                                                          gggAACGTC‑ACCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCatat  >  2:1048307/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                             GCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTC‑A‑CCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATAT  >  NC_000913/1793698‑1793984

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