Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	788,389	A→C	7.1%	L150R (CTG→CGG)	galM ←	galactose‑1‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	788,389	0	A	C	7.1%	256.5 / 10.6	85	L150R (CTG→CGG)	galM	galactose‑1‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (44/35);  new base C (3/3);  total (47/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

GATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTT  >  NC_000913/788240‑788537                                                                                                                                                      |                                                                                                                                                     gATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCa                                                                                                                                                      >  1:1471547/1‑150 (MQ=255)   tatTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGa                                                                                                                                                    >  2:86486/1‑150 (MQ=255)   tatTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGa                                                                                                                                                    >  1:716940/1‑150 (MQ=255)           cGCCAGAATCTGCTACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAATATTGAAATAGACGTGATTAGTCATATTCACCGGGCTAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATAGGTTTATCGTCGGTCAGACGATATTg                                                                                                                                            <  2:1024781/150‑1 (MQ=255)           cGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGACCGACGataacc                                                                                                                                            <  2:1499385/150‑4 (MQ=255)           cGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGGCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGACCGACGataacc                                                                                                                                            <  2:1431671/150‑4 (MQ=255)             ccAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCa                                                                                                                                          <  1:562527/150‑1 (MQ=255)              cAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCTGACCGACGataaccgta                                                                                                                                         <  1:888378/150‑7 (MQ=255)                taaTCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACTATATTGCACCg                                                                                                                                       <  2:786694/149‑1 (MQ=255)                  aTCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTc                                                                                                                                     <  1:286105/150‑1 (MQ=255)                   tCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcg                                                                                                                                    >  1:407268/1‑150 (MQ=255)                   tCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTcg                                                                                                                                    <  1:399636/150‑1 (MQ=255)                          cTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGGTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCgag                                                                                                                             <  1:142233/150‑1 (MQ=255)                           ttGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCgaga                                                                                                                            >  1:591766/1‑150 (MQ=255)                              tgATTGCGCACGTCAGACTGCTAGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTa                                                                                                                         <  2:1457212/150‑1 (MQ=255)                                  tGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCg                                                                                                                     <  2:672360/150‑1 (MQ=255)                                   gcgcACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCgg                                                                                                                    >  1:665017/1‑150 (MQ=255)                                    cgcACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCggg                                                                                                                   <  2:201027/150‑1 (MQ=255)                                    cgcACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCggg                                                                                                                   <  1:793748/150‑1 (MQ=255)                                      cACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCGCCGTCGCGCCGAGATTACCCGGGaa                                                                                                                 >  1:1315045/1‑150 (MQ=255)                                          tCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGccc                                                                                                             >  1:156552/1‑150 (MQ=255)                                          tCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGccc                                                                                                             <  1:13640/150‑1 (MQ=255)                                          tCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGccc                                                                                                             >  1:708240/1‑150 (MQ=255)                                           cAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCt                                                                                                            <  1:86486/150‑1 (MQ=255)                                             gACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGa                                                                                                          >  2:25624/1‑150 (MQ=255)                                               cTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATc                                                                                                        <  2:318415/150‑1 (MQ=255)                                                tGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCa                                                                                                       <  1:999826/150‑1 (MQ=255)                                                    cGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcat                                                                                                   >  2:1438380/1‑150 (MQ=255)                                                       cGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATGACcatcat                                                                                                >  1:1433674/1‑150 (MQ=255)                                                       cGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACcatcat                                                                                                <  1:570969/150‑1 (MQ=255)                                                          cAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTg                                                                                             >  1:1394144/1‑150 (MQ=255)                                                            aGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGaa                                                                                           >  2:588898/1‑150 (MQ=255)                                                             gATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAAc                                                                                          >  2:148357/1‑150 (MQ=255)                                                                tGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCa                                                                                       >  2:952184/1‑150 (MQ=255)                                                                 gAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAg                                                                                      >  2:473765/1‑150 (MQ=255)                                                                 gAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAg                                                                                      >  2:1190975/1‑150 (MQ=255)                                                                 gAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAg                                                                                      >  2:589025/1‑150 (MQ=255)                                                                 gAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAg                                                                                      >  1:423601/1‑150 (MQ=255)                                                                  aaaTAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATAAACTGTGTCGCGATAAGTAATGTAGATACTGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAgg                                                                                     <  1:1008952/150‑1 (MQ=255)                                                                     tAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCa                                                                                  <  2:851340/150‑1 (MQ=255)                                                                       gACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCaaa                                                                                >  2:1413309/1‑150 (MQ=255)                                                                           tGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGc                                                                            >  2:617050/1‑150 (MQ=255)                                                                           tGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGc                                                                            >  1:1184371/1‑150 (MQ=255)                                                                            gATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCa                                                                           >  2:5414/1‑150 (MQ=255)                                                                                 gTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGa                                                                      >  2:755063/1‑150 (MQ=255)                                                                                 gTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGa                                                                      >  2:1184225/1‑150 (MQ=255)                                                                                 gTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGACAAACAGCACCTGa                                                                      >  1:1128324/1‑150 (MQ=255)                                                                                 ctCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGa                                                                      >  1:1419685/2‑150 (MQ=255)                                                                                     tatTCACCGGGCATGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGAt                                                                  <  2:1538486/150‑1 (MQ=255)                                                                                          aCCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCTCCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTc                                                             <  2:1315045/150‑1 (MQ=255)                                                                                          aCCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACATGACGATCGTTc                                                             <  1:639488/150‑1 (MQ=255)                                                                                          aCCGGGCAAGGTTTAGCATCTGTGGCGCGATTAGTTATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGATCTCAGGTCAAACAGCACCTGACGATCGTTc                                                             <  2:1180441/150‑1 (MQ=255)                                                                                           ccGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCt                                                            <  1:952184/150‑1 (MQ=255)                                                                                             gggCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTgg                                                          >  2:1083001/1‑150 (MQ=255)                                                                                              ggCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGt                                                         <  1:1303492/150‑1 (MQ=255)                                                                                              ggCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGt                                                         <  1:1413309/150‑1 (MQ=255)                                                                                               gCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGtt                                                        >  2:897314/1‑150 (MQ=255)                                                                                                 aaGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTca                                                      <  2:31727/150‑1 (MQ=255)                                                                                                        aTCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTg                                               <  1:1310063/150‑1 (MQ=255)                                                                                                            aCTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAg                                           <  2:1046375/150‑1 (MQ=255)                                                                                                             cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAAt                                                  <  1:1574155/142‑1 (MQ=255)                                                                                                             cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAAt                                                  >  2:1574155/1‑142 (MQ=255)                                                                                                             cTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGc                                          >  1:401732/1‑150 (MQ=255)                                                                                                               gtgGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGa                                        >  2:1030681/1‑150 (MQ=255)                                                                                                                   cgcgATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGCTTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCCCAATCTGCCAGCGACGtt                                    >  2:1404217/1‑150 (MQ=255)                                                                                                                   cgcgATAAGAAATGAAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCAACATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGtt                                    >  2:588716/1‑150 (MQ=255)                                                                                                                      gATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGt                                 >  1:730735/1‑150 (MQ=255)                                                                                                                         aaGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGa                              >  1:1504959/1‑150 (MQ=255)                                                                                                                          aGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGaa                             >  1:1054007/1‑150 (MQ=255)                                                                                                                           gTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAAc                            <  1:5414/150‑1 (MQ=255)                                                                                                                            tAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAAcc                           >  1:1516841/1‑150 (MQ=255)                                                                                                                               tGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCtt                        <  1:1438380/150‑1 (MQ=255)                                                                                                                                ggAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTc                       >  1:1544408/1‑150 (MQ=255)                                                                                                                                 gagaTACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTcc                      >  1:1529869/1‑150 (MQ=255)                                                                                                                                  agaAACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATGACCCGGGAAGCCCTGATCACCATCCTCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCg                     >  2:1537779/1‑150 (MQ=255)                                                                                                                                   taTACGGTTATCGTCGGTCAGACGATATTGCACCTTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGTGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCgg                    <  2:716940/149‑1 (MQ=255)                                                                                                                                    aTACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGc                   <  2:156552/150‑1 (MQ=255)                                                                                                                                     tACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGcc                  >  1:925820/1‑150 (MQ=255)                                                                                                                                     gtgcaatATCGTCTGACCGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGcc                  >  2:1515268/7‑150 (MQ=255)                                                                                                                                      aCGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGccc                 >  1:852157/1‑150 (MQ=255)                                                                                                                                       gcaatATCGTCTGACCGACGATATTGCACCGTCGCGCCGAGCTTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCccg                >  1:1334592/5‑150 (MQ=255)                                                                                                                                        caatATCGTCTGACCGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAACCAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCccgc               >  1:406941/4‑150 (MQ=255)                                                                                                                                           tATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGt            >  2:409827/1‑150 (MQ=255)                                                                                                                                                  ggTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATTGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTg     <  1:473765/150‑1 (MQ=255)                                                                                                                                                    tCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGt   <  1:755063/150‑1 (MQ=255)                                                                                                                                                     cAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGtt  <  1:25624/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GATATTCGTCCGCCAGAATCTGCAACTTGTGATTGCGCACGTCAGACTGCTCGCCGTCAAGATTGAAATAGACGTGATTAGTCATATTCACCGGGCAAGGTTTATCAACTGTGGCGCGATAAGTAATGGAGATACGGTTATCGTCGGTCAGACGATATTGCACCGTCGCGCCGAGATTACCCGGGAAGCCCTGATCACCATCATCTGAACTCAGGGCAAACAGCACCTGACGATCGTTCTGGTTCACAATCTGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTT  >  NC_000913/788240‑788537

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