Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,068,830	A→G	F337S (TTT→TCT)	rutG ←	pyrimidine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,068,830	0	A	G	97.9%	335.4 / ‑1.5	93	F337S (TTT→TCT)	rutG	pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (37/54);  total (37/56)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

ATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAAAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCG  >  NC_000913/1068683‑1068976                                                                                                                                                    |                                                                                                                                                   aTCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGc                                                                                                                                                 <  2:60806/151‑1 (MQ=255)   cAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGccc                                                                                                                                               <  1:1084257/151‑1 (MQ=255)    aaaTTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGccc                                                                                                                                               >  2:347950/1‑150 (MQ=255)    aaaTTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCa                                                                                                                                              <  1:686766/151‑1 (MQ=255)     aaTTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCaa                                                                                                                                             <  1:518032/151‑1 (MQ=255)        tACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACa                                                                                                                                           <  2:1519973/150‑1 (MQ=255)          ccGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCa                                                                                                                                        <  1:240992/151‑1 (MQ=255)           cGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCAt                                                                                                                                       <  1:615718/151‑1 (MQ=255)            gTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATg                                                                                                                                      <  1:887100/151‑1 (MQ=255)            gTTCTGGCTTAAATCGACACGGTTTTGTACCCATAACCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCTCCAATAACCGCGGCCGGAATGGTATGAATCAGCGCCCCAAAGTTCGGTGAAGAGCCCAGCAGCATg                                                                                                                                      <  1:20299/151‑1 (MQ=255)              tCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATgg                                                                                                                                     <  1:905219/150‑1 (MQ=255)                  gCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATg                                                                                                                                <  2:538671/151‑1 (MQ=255)                   cTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGCATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATg                                                                                                                                >  1:1573368/1‑150 (MQ=255)                   cTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATg                                                                                                                                >  1:1010615/1‑150 (MQ=255)                   cTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATg                                                                                                                                >  1:1180118/1‑150 (MQ=255)                          cGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGgcag                                                                                                                         <  2:34450/150‑1 (MQ=255)                           gACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGgcagca                                                                                                                       <  1:849446/151‑1 (MQ=255)                            acacGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGgcagcag                                                                                                                      <  1:1312899/151‑1 (MQ=255)                             cacGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGgcagcagc                                                                                                                     >  2:162877/1‑151 (MQ=255)                                   tttGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGcaacaa                                                                                                                <  1:1465480/150‑1 (MQ=255)                                    ttGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCaacaaa                                                                                                               <  1:1204172/150‑1 (MQ=255)                                      gTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCaacaaaca                                                                                                             <  2:1215243/150‑1 (MQ=255)                                        aCCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCa                                                                                                          <  2:859178/151‑1 (MQ=255)                                            atatCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGt                                                                                                      <  1:298859/151‑1 (MQ=255)                                             tatCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTc                                                                                                     >  2:465065/1‑151 (MQ=255)                                               tCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCg                                                                                                    >  1:154492/1‑150 (MQ=255)                                                ccTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAg                                                                                                  <  2:345881/151‑1 (MQ=255)                                                ccTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAAAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGa                                                                                                   <  1:1282792/150‑1 (MQ=255)                                                ccTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAAAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAg                                                                                                  <  2:493291/151‑1 (MQ=255)                                                 cTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGAGTCGAg                                                                                                  >  1:1604800/1‑150 (MQ=255)                                                  ttGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGt                                                                                                 >  1:1196846/1‑150 (MQ=255)                                                   tGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGt                                                                                                 <  2:1426294/149‑1 (MQ=255)                                                    gcgcCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAg                                                                                               <  2:1069243/150‑1 (MQ=255)                                                     cgcCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGAc                                                                                             >  1:947101/1‑151 (MQ=255)                                                         tGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTg                                                                                         >  1:1367491/1‑151 (MQ=255)                                                         tGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTg                                                                                         >  1:1685201/1‑151 (MQ=255)                                                          gCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACttt                                                                                          <  2:1292704/149‑1 (MQ=255)                                                           cGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGt                                                                                       <  2:333354/151‑1 (MQ=255)                                                            gACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTc                                                                                      >  1:84175/1‑151 (MQ=255)                                                             aCGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCa                                                                                     <  1:778268/151‑1 (MQ=255)                                                               ggCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACg                                                                                   <  2:793882/151‑1 (MQ=255)                                                                   aTCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCa                                                                               <  1:75442/151‑1 (MQ=255)                                                                     cAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCAt                                                                              <  2:13832/150‑1 (MQ=255)                                                                     cAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATc                                                                             <  2:1376508/151‑1 (MQ=255)                                                                     cAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATc                                                                             >  2:1234958/1‑151 (MQ=255)                                                                     cAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATc                                                                             >  2:464477/1‑151 (MQ=255)                                                                     cAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATc                                                                             >  1:623645/1‑151 (MQ=255)                                                                     cAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATc                                                                             >  1:1100181/1‑151 (MQ=255)                                                                        cccGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATcac                                                                           >  1:335676/1‑150 (MQ=255)                                                                         ccGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCacac                                                                         <  2:1331861/151‑1 (MQ=255)                                                                         ccGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCacac                                                                         <  1:807134/151‑1 (MQ=255)                                                                           gAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCg                                                                       <  1:1407427/151‑1 (MQ=255)                                                                              cacaACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGat                                                                     >  1:991742/1‑150 (MQ=255)                                                                                  acaaTTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATAtttt                                                                <  2:193270/151‑1 (MQ=255)                                                                                   caaTTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTc                                                               >  2:804102/1‑151 (MQ=255)                                                                                    aaTTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTc                                                               >  2:626917/1‑150 (MQ=255)                                                                                     aTTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTc                                                               >  1:1326198/1‑149 (MQ=255)                                                                                       tGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGc                                                            >  2:1120720/1‑150 (MQ=255)                                                                                         aTGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCAt                                                          >  1:947622/1‑150 (MQ=255)                                                                                            caccacCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAggt                                                      <  2:1045464/151‑1 (MQ=255)                                                                                              ccaccaATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAggtgg                                                    <  2:892978/151‑1 (MQ=255)                                                                                                 ccaATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCa                                                 <  2:165006/151‑1 (MQ=255)                                                                                                  caATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCAc                                                <  1:119577/151‑1 (MQ=255)                                                                                                    aTAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCAcgc                                              <  1:923444/151‑1 (MQ=255)                                                                                                     tAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCAcgc                                              >  2:1021157/1‑150 (MQ=255)                                                                                                      aaCCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCAcgccg                                            <  1:586187/151‑1 (MQ=255)                                                                                                      aaCCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCAcgcc                                             <  1:1249946/150‑1 (MQ=255)                                                                                                          gcgGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTAc                                        >  2:819267/1‑151 (MQ=255)                                                                                                           cgGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACTCCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTAcc                                       <  1:250957/151‑1 (MQ=255)                                                                                                           cgGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTAcc                                       <  2:335702/151‑1 (MQ=255)                                                                                                              ccGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTAccgcc                                    >  1:1369831/1‑151 (MQ=255)                                                                                                               cGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTAccgccg                                   >  1:657180/1‑151 (MQ=255)                                                                                                                ggAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTAccgccg                                   <  2:84179/150‑1 (MQ=255)                                                                                                                 gAATGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGAc                                 >  2:1486747/1‑151 (MQ=255)                                                                                                                  aaTGGTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGAc                                 >  2:326421/1‑150 (MQ=255)                                                                                                                     ggTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCATCGTCGAGTAGACTTGGGTCCCGCCCATCACTCCAATATTTCCAACATAGTTGCCCTCCCCGCTACCGCCGACaga                              <  1:395199/150‑1 (MQ=255)                                                                                                                     ggTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCTAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACagag                             <  1:587815/151‑1 (MQ=255)                                                                                                                     ggTATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACagag                             <  2:355885/151‑1 (MQ=255)                                                                                                                       tATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGcc                           >  2:932024/1‑151 (MQ=255)                                                                                                                       tATGAATCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGcc                           >  1:172307/1‑151 (MQ=255)                                                                                                                           aaTCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGaa                       <  2:750759/151‑1 (MQ=255)                                                                                                                           aaTCAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGa                        >  1:406139/1‑150 (MQ=255)                                                                                                                              cAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAg                     <  2:1222110/150‑1 (MQ=255)                                                                                                                              cAGCGCACCAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGc                    >  2:1695794/1‑151 (MQ=255)                                                                                                                                      cAAATTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGc            >  1:1536940/1‑151 (MQ=255)                                                                                                                                       aaaTTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCACCAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGc            <  1:419275/150‑1 (MQ=255)                                                                                                                                       aaaTTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGcc           >  2:1360444/1‑151 (MQ=255)                                                                                                                                       aaaTTTCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGAAGCAACAAACCCCAGCGTCGAGTAGACTTTGGTCACGGCCATCCCACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGt               <  2:1068126/147‑1 (MQ=255)                                                                                                                                         aTTCCGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACCCCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAg         <  2:640113/151‑1 (MQ=255)                                                                                                                                             cGGTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCa     <  2:154501/151‑1 (MQ=255)                                                                                                                                              ggTGAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCAt    >  1:760963/1‑151 (MQ=255)                                                                                                                                                 gAAGAGCCCAACAGCATGGCGATGCCGGCACAAGCAACAAACACCAGCGTCGAGTAGCCTTTGGCCACGGCCATCCCACCGATATTTTCAGCAGAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCg  <  2:1573465/150‑1 (MQ=255)                                                                                                                                                 gAAGAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCg  <  2:261969/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   ATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAGCGCACCAAATTTCGGTGAAAAGCCCAACAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCG  >  NC_000913/1068683‑1068976

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