Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	37,311	T→A	12.5%	T172S (ACT→TCT) ‡	caiC ←	carnitine‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	37,311	0	T	A	12.5%	236.3 / 21.8	88	T172S (ACT→TCT) ‡	caiC	carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (39/38);  new base A (8/3);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

CATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGAT  >  NC_000913/37162‑37457                                                                                                                                                       |                                                                                                                                                   cATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGAGCAGGTTGTAATGGGTAATCACCACAACTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGACG‑TATCGTCAGt                                                                                                                                                    >  2:174568/1‑150 (MQ=255)   ttACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTCCCGCCG‑TATCGTCAGTCg                                                                                                                                                  <  2:435155/150‑1 (MQ=255)     aCCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGAt                                                                                                                                                <  1:1474607/150‑1 (MQ=255)            ggTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAgcggcg                                                                                                                                         >  1:359834/1‑150 (MQ=255)              tAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATA‑cggcgga                                                                                                                                       <  2:456615/150‑2 (MQ=255)               aGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTg                                                                                                                                      <  1:73676/150‑1 (MQ=255)               aGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTg                                                                                                                                      <  1:1599791/150‑1 (MQ=255)               aGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTg                                                                                                                                      <  2:1788547/150‑1 (MQ=255)                gACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGc                                                                                                                                     <  1:1475228/150‑1 (MQ=255)                 aCGTCATCGTCACGCAGTGCACACTGCTAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATA‑cggcggaaat                                                                                                                                    <  1:1264376/150‑5 (MQ=255)                    tCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAg                                                                                                                                 >  2:1864072/1‑150 (MQ=255)                         gTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATA‑cggcggaaat                                                                                                                                    >  1:1092044/1‑138 (MQ=255)                         gTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATA‑cggcggaaat                                                                                                                                    <  2:1092044/142‑5 (MQ=255)                         gTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACaa                                                                                                                            >  1:581863/1‑150 (MQ=255)                           cACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAgg                                                                                                                          <  1:590714/150‑1 (MQ=255)                             cGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTg                                                                                                                        <  2:359834/150‑1 (MQ=255)                                aGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCTTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCa                                                                                                                     <  2:1310117/150‑1 (MQ=255)                                 gTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAg                                                                                                                    <  2:920533/150‑1 (MQ=255)                                   gCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgt                                                                                                                  >  2:1394145/1‑150 (MQ=255)                                   gCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgt                                                                                                                  >  2:1204766/1‑150 (MQ=255)                                       acTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgttgtt                                                                                                              >  1:1258313/1‑150 (MQ=255)                                        cTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATTGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGttgttg                                                                                                             <  1:1452635/150‑1 (MQ=255)                                          gCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGAt                                                                                                           >  1:1529274/1‑150 (MQ=255)                                            cAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGAttt                                                                                                         <  2:422210/150‑1 (MQ=255)                                            cAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGAttt                                                                                                         <  1:588059/150‑1 (MQ=255)                                               gCGGTGTAATATACAGCGTATCTCATGTTGTTATTGGTAATCACCACTGCTTTCTGTCGGGTGGTGGTGCCGGAGGTGAAGTGAATTTCCTCCG‑TATATTCAGTCGATCGCGGCGGTGCATATCACAAGGTGGCAGGTTGTTGATTTTTc                                                                                                      <  2:1045909/150‑1 (MQ=255)                                                 ggAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAg                                                                                                    <  2:1399527/150‑1 (MQ=255)                                                  gAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGt                                                                                                   >  1:942529/1‑150 (MQ=255)                                                   aGTAATATCCATCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGtt                                                                                                  <  2:1446684/150‑1 (MQ=255)                                                    gTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTg                                                                                                 <  1:1315475/150‑1 (MQ=255)                                                    gTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTg                                                                                                 <  1:1392324/150‑1 (MQ=255)                                                    gTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTg                                                                                                 >  2:1147164/1‑150 (MQ=255)                                                       atatCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGt                                                                                              <  2:1216623/150‑1 (MQ=255)                                                          tCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTaaa                                                                                           >  1:141654/1‑150 (MQ=255)                                                                      cAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGcc                                                                               <  1:1685748/150‑1 (MQ=255)                                                                           tGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCcatcat                                                                          >  2:490314/1‑150 (MQ=255)                                                                              aaTGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGa                                                                                        >  2:144926/1‑133 (MQ=255)                                                                              aaTGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGa                                                                                        <  1:144926/133‑1 (MQ=255)                                                                                  ggTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCggg                                                                   >  1:1162297/1‑150 (MQ=255)                                                                                        caccacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGc                                                             >  1:260006/1‑150 (MQ=255)                                                                                        caccacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGc                                                             >  1:267295/1‑150 (MQ=255)                                                                                         accacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAg                                                                                                    >  2:665552/1‑110 (MQ=255)                                                                                         accacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAg                                                                                                    <  1:665552/110‑1 (MQ=255)                                                                                            acacCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCcaca                                                         <  2:492212/150‑1 (MQ=255)                                                                                               ccTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCt                                                      >  2:1859379/1‑150 (MQ=255)                                                                                               ccTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCt                                                      >  2:540514/1‑150 (MQ=255)                                                                                                cTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTg                                                     >  2:717992/1‑150 (MQ=255)                                                                                                  ttCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTc                                                   >  1:1181181/1‑150 (MQ=255)                                                                                                  ttCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTc                                                   <  1:1147164/150‑1 (MQ=255)                                                                                                   tcggtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCa                                                  <  1:1864072/150‑1 (MQ=255)                                                                                                    cggtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAg                                                 <  2:1181181/150‑1 (MQ=255)                                                                                                     ggtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAgg                                                <  2:942529/150‑1 (MQ=255)                                                                                                      gtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGc                                               <  2:5780/150‑1 (MQ=255)                                                                                                      gtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGc                                               <  2:141654/150‑1 (MQ=255)                                                                                                         gggAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCaaa                                            >  1:475541/1‑150 (MQ=255)                                                                                                         gggAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCaaa                                            >  2:1857711/1‑150 (MQ=255)                                                                                                               tggtgCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGcc                                      >  1:834812/1‑150 (MQ=255)                                                                                                                ggtgCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCg                                     >  1:415543/1‑150 (MQ=255)                                                                                                                  tgCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCa                                   >  1:222248/1‑150 (MQ=255)                                                                                                                   gCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCaa                                  >  1:570022/1‑150 (MQ=255)                                                                                                                     cGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAAtt                                >  2:1549671/1‑150 (MQ=255)                                                                                                                      ggAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTg                               >  1:560677/1‑150 (MQ=255)                                                                                                                      ggAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTg                               <  2:646478/150‑1 (MQ=255)                                                                                                                         ggTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGt                            <  1:844362/150‑1 (MQ=255)                                                                                                                          gTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTg                           <  2:1529274/150‑1 (MQ=255)                                                                                                                           tGAAGAGAATTTCCGCCG‑TATCGTTAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTgg                          >  2:1449233/1‑150 (MQ=255)                                                                                                                           tGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTgg                          >  2:69306/1‑150 (MQ=255)                                                                                                                              agagAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCAt                       >  1:571982/1‑150 (MQ=255)                                                                                                                              agagAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCAt                       >  1:523665/1‑150 (MQ=255)                                                                                                                               gagaATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGGAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATc                      >  1:1033968/1‑150 (MQ=255)                                                                                                                                agaATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCt                     >  2:1388091/1‑150 (MQ=255)                                                                                                                                 gaATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGCTGGAAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACTCCATCATCAGCGGGAAGTGCCACATCTGTCCGGCAAATGTGCCGCAATTGAGTGGCATCtt                    >  2:299501/1‑150 (MQ=255)                                                                                                                                      accgccgcTATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctgc               >  2:1631377/2‑150 (MQ=255)                                                                                                                                       ccgccgcTATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctgct              >  1:1892132/1‑150 (MQ=255)                                                                                                                                        cgccgctATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctgctg             >  1:1964465/7‑150 (MQ=255)                                                                                                                                        cgccg‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctgctg             <  1:717992/150‑1 (MQ=255)                                                                                                                                         gccgctATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa            >  2:406176/6‑150 (MQ=255)                                                                                                                                         gccgctATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa            >  2:1085774/6‑150 (MQ=255)                                                                                                                                         gccgctATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa            >  2:712572/6‑150 (MQ=255)                                                                                                                                         gccgctATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGCGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa            >  2:1878353/6‑150 (MQ=255)                                                                                                                                         gccg‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa            >  1:173812/1‑150 (MQ=255)                                                                                                                                            g‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATc         <  1:490314/150‑1 (MQ=255)                                                                                                                                               aTCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTg       <  1:1204766/150‑1 (MQ=255)                                                                                                                                               aTCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTg       >  2:1257340/1‑150 (MQ=255)                                                                                                                                                 cGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGtt     <  2:560677/150‑1 (MQ=255)                                                                                                                                                 cGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGtt     <  1:1131303/150‑1 (MQ=255)                                                                                                                                                 cGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATATGtt     >  2:1131303/1‑150 (MQ=255)                                                                                                                                                  gTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg    >  2:1341603/1‑150 (MQ=255)                                                                                                                                                    cAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGAt  <  1:1341603/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                   CATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGAT  >  NC_000913/37162‑37457

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: ‑