Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,174,410	A→T	7.7%	L53M (TTG→ATG) ‡	gatB ←	galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,174,410	0	A	T	7.7%	232.1 / 4.8	76	L53M (TTG→ATG) ‡	gatB	galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base A (37/33);  new base T (4/2);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCAAGCGACAATAATCTT  >  NC_000913/2174261‑2174557                                                                                                                                                      |                                                                                                                                                    cATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCaa                                                                                                                                                     >  1:768344/1‑150 (MQ=255)      cTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGca                                                                                                                                                >  1:1363831/1‑150 (MQ=255)         cTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCacac                                                                                                                                             >  2:757164/1‑150 (MQ=255)         cTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCacac                                                                                                                                             >  2:658265/1‑150 (MQ=255)          tGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCacaca                                                                                                                                            >  1:1147242/1‑149 (MQ=255)          tGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAcc                                                                                                                                            >  1:934656/1‑150 (MQ=255)           gAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACAcca                                                                                                                                           >  2:1005978/1‑150 (MQ=255)            aaaaCATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATTTGATATGcaccactg                                                                                                                                          <  1:535344/150‑6 (MQ=255)             aaaCATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATTTGATATGcaccactgc                                                                                                                                         <  2:1270117/150‑7 (MQ=255)                   aGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAAAGGAATATCGCCAAAACTACGATCCACTTTCGCAGTGGTGCATATCAAATGCACACCATCCAtata                                                                                                                                   >  2:497371/1‑150 (MQ=255)                    ggTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAg                                                                                                                                  >  1:1096214/1‑150 (MQ=255)                    ggTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAg                                                                                                                                  >  1:185722/1‑150 (MQ=255)                     gTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAgg                                                                                                                                 <  1:495797/150‑1 (MQ=255)                        aCCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGc                                                                                                                                                 >  1:830544/1‑131 (MQ=255)                        aCCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGc                                                                                                                                                 <  2:830544/131‑1 (MQ=255)                         ccccTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTc                                                                                                                             >  1:238566/1‑150 (MQ=255)                         ccccTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTc                                                                                                                             >  2:412287/1‑150 (MQ=255)                          cccTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCt                                                                                                                            <  2:1520965/150‑1 (MQ=255)                             tGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGTAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTAtt                                                                                                                         <  1:1484196/150‑1 (MQ=255)                             tGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTAtt                                                                                                                         <  1:1494867/150‑1 (MQ=255)                              gTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTAttt                                                                                                                        <  1:878407/150‑1 (MQ=255)                              gTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTAttt                                                                                                                        <  2:64120/150‑1 (MQ=255)                               tAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTc                                                                                                                       <  1:222001/150‑1 (MQ=255)                                  gATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCAtt                                                                                                                    <  1:566381/150‑1 (MQ=255)                                  gATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCAtt                                                                                                                    >  2:1353739/1‑150 (MQ=255)                                    tAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTaa                                                                                                                  >  2:1139805/1‑150 (MQ=255)                                    tAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTaa                                                                                                                  >  1:155242/1‑150 (MQ=255)                                      gTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAAcc                                                                                                                <  2:1363831/150‑1 (MQ=255)                                        cAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCg                                                                                                              >  1:1157846/1‑150 (MQ=255)                                             ttttattttGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACt                                                                                                         <  1:1534052/150‑1 (MQ=255)                                                tattttGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGa                                                                                                      >  2:1493356/1‑150 (MQ=255)                                                 attttGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGAt                                                                                                     <  1:93679/150‑1 (MQ=255)                                                      gTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTCACGGAATATCGCCAAAACTACGATCCACTTTGGTTGTTTATTTTTTCTAATTCACACGATCGATATAGGGTTCTATTTGATTTACCCGACAGTGTATTAAtt                                                                                                <  2:715195/150‑1 (MQ=255)                                                       tAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTc                                                                                               <  1:757164/150‑1 (MQ=255)                                                        aaTGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCa                                                                                              <  1:229214/150‑1 (MQ=255)                                                        aaTGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCa                                                                                              >  1:868770/1‑150 (MQ=255)                                                          tGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAAc                                                                                            <  2:956597/150‑1 (MQ=255)                                                           gCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACa                                                                                           <  1:1139805/150‑1 (MQ=255)                                                                gataccgataccAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGTCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAat                                                                                      <  2:698616/150‑1 (MQ=255)                                                                gataccgataccAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAat                                                                                      >  1:50697/1‑150 (MQ=255)                                                                  taccgataccAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAatat                                                                                    >  2:830866/1‑150 (MQ=255)                                                                      gataccAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATg                                                                                <  1:310178/150‑1 (MQ=255)                                                                            aGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTg                                                                          <  1:520169/150‑1 (MQ=255)                                                                                 tAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGAcaca                                                                     >  1:810636/1‑150 (MQ=255)                                                                                        cATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCttt                                                              <  2:155242/150‑1 (MQ=255)                                                                                           gCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAAt                                                           >  1:324796/1‑150 (MQ=255)                                                                                           gCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAAt                                                           >  2:796547/1‑150 (MQ=255)                                                                                            ccGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAAtt                                                          >  2:687458/1‑150 (MQ=255)                                                                                              gTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATttc                                                        <  1:317668/150‑1 (MQ=255)                                                                                                    tAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAg                                                                                          >  2:1193622/1‑110 (MQ=255)                                                                                                    tAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAg                                                                                          <  1:1193622/110‑1 (MQ=255)                                                                                                    tAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTccgc                                                  <  1:658265/150‑1 (MQ=255)                                                                                                            aatCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCg                                          <  2:197242/149‑1 (MQ=255)                                                                                                              tCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTa                                        >  1:1293660/1‑150 (MQ=255)                                                                                                               cGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTag                                       <  2:934656/150‑1 (MQ=255)                                                                                                               cGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTag                                       <  1:697567/150‑1 (MQ=255)                                                                                                                gCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGCTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCAGCCGCCATCGTaga                                      >  1:1530272/1‑150 (MQ=255)                                                                                                                           gATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATGCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACcg                           <  2:633301/150‑1 (MQ=255)                                                                                                                           gATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCCGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACcg                           >  2:566094/1‑150 (MQ=255)                                                                                                                             tCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACcgcg                         <  2:238566/150‑1 (MQ=255)                                                                                                                             tCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACcgcg                         <  1:1353739/150‑1 (MQ=255)                                                                                                                               cACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGcc                       <  1:1005978/150‑1 (MQ=255)                                                                                                                               cACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGcc                       >  2:851448/1‑150 (MQ=255)                                                                                                                               cACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGcc                       <  1:412287/150‑1 (MQ=255)                                                                                                                                aCTTTGGCAGTGGTGCATCTCAACTGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTGGAGGTCGCAACCGCGCCt                      >  2:1596445/1‑150 (MQ=255)                                                                                                                                    tGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGc                  <  2:768344/150‑1 (MQ=255)                                                                                                                                     ggCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCa                 >  1:103721/1‑150 (MQ=255)                                                                                                                                      gatggtgtgCATTTGATATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCaa                >  2:356093/7‑150 (MQ=255)                                                                                                                                       atggtgtgCATTTGATATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCAAg               >  2:1482732/6‑150 (MQ=255)                                                                                                                                       atggtgtgCATATGATATGCACACCATCCATATACGATTCTATTTCATTAACCCGACACTGGATTAATTCAACACGAATATTATGATTCTCACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCAAg               >  2:1488664/6‑150 (MQ=255)                                                                                                                                        aGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCAAGc              >  1:1445357/1‑150 (MQ=255)                                                                                                                                           tgtgCATTTGATATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCAAGCGAc           >  1:1165241/2‑150 (MQ=255)                                                                                                                                           ggtgCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCAAGCGAc           >  2:1598322/1‑150 (MQ=255)                                                                                                                                             tgCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCAAGCGACaa         >  1:1034623/1‑150 (MQ=255)                                                                                                                                             tgCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCAAGCGACaa         >  1:903892/1‑150 (MQ=255)                                                                                                                                                  atCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCAAGCGACAATAATc    >  1:844519/1‑150 (MQ=255)                                                                                                                                                    cAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCAAGCGACAATAATCtt  >  2:89656/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CATGACTTCTGAAAACATAGGTCACCCCTGTAAGATAGTCAGAATTTTATTTTGTAATGCTTCGATACCGATACCAGAAATAAAAGGCATGCCGTGAACTAACGGAATATCGCCAAAACTACGATCCACTTTGGCAGTGGTGCATATCAAATGCACACCATCCATATAGGTTTCTATTTCATTAACCCGACACTGGATTAATTCAACAGGAATATTATGATTCTGACACAACTCTTTAATTTCTTCCGCCGCCATCGTAGAGGTCGCAACCGCGCCTCCGCAAGCGACAATAATCTT  >  NC_000913/2174261‑2174557

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