Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,389,934	A→G	Y143H (TAC→CAC)	psd ←	phosphatidylserine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,389,934	0	A	G	100.0%	272.3 / NA	75	Y143H (TAC→CAC)	psd	phosphatidylserine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/38);  total (37/38)

GCAAATCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTAGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATCTTCTTCGATTTT  >  NC_000913/4389800‑4390069                                                                                                                                       |                                                                                                                                        gCAAATCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCa                                                                                                                                     >  2:1470953/1‑139 (MQ=255)      tCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCgg                                                                                                                                <  2:161911/139‑1 (MQ=255)      tCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCgg                                                                                                                                >  2:743411/1‑139 (MQ=255)       cACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCggg                                                                                                                               >  1:587851/1‑139 (MQ=255)       cACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCggg                                                                                                                               >  2:1337463/1‑139 (MQ=255)         cgcgTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGga                                                                                                                             <  2:1703376/139‑1 (MQ=255)             ttcgttACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAgg                                                                                                                         >  1:1813069/1‑139 (MQ=255)                gttACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTaa                                                                                                                      <  1:1470953/139‑1 (MQ=255)                        cAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTcac                                                                                                              >  1:1266113/1‑139 (MQ=255)                                ttCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTAc                                                                                                      <  2:1657611/139‑1 (MQ=255)                                     cACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTg                                                                                                 <  2:587851/139‑1 (MQ=255)                                          tCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGaa                                                                                            <  2:1758738/139‑1 (MQ=255)                                               gCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGt                                                                                       <  2:776988/139‑1 (MQ=255)                                                ccGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                      <  1:1337463/139‑1 (MQ=255)                                                ccGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                      <  1:1620635/139‑1 (MQ=255)                                                   tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGt                                                                                                                        >  1:265020/1‑102 (MQ=255)                                                   tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGt                                                                                                                        <  2:265020/102‑1 (MQ=255)                                                   tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                      >  1:1096731/1‑136 (MQ=255)                                                   tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                      <  2:1096731/136‑1 (MQ=255)                                                   tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGc                                                                                   >  1:1806141/1‑139 (MQ=255)                                                   tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGc                                                                                   >  2:468707/1‑139 (MQ=255)                                                   tGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGc                                                                                   <  2:1266113/139‑1 (MQ=255)                                                          gTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGa                                                                            <  1:743411/139‑1 (MQ=255)                                                              aCGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGt                                                                        <  2:1813069/139‑1 (MQ=255)                                                                 gagaAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACaa                                                                                                            <  1:536753/100‑1 (MQ=255)                                                                 gagaAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACaa                                                                                                            >  2:536753/1‑100 (MQ=255)                                                                  agaagaGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTtgcc                                                                    >  2:863308/1‑139 (MQ=255)                                                                    aagaGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTtg                                                                      >  2:372363/1‑135 (MQ=255)                                                                    aagaGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTtg                                                                      <  1:372363/135‑1 (MQ=255)                                                                     agagATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTtgcctgc                                                                 >  2:391514/1‑139 (MQ=255)                                                                     agagATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTtgcctgc                                                                 >  1:1338139/1‑139 (MQ=255)                                                                         aTCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTAc                                                                                                      >  2:1825006/1‑98 (MQ=255)                                                                         aTCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTAc                                                                                                      <  1:1825006/98‑1 (MQ=255)                                                                            gCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAgg                                                          >  2:723106/1‑139 (MQ=255)                                                                             cccGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAggg                                                         <  1:468707/139‑1 (MQ=255)                                                                                  cACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCg                                                    <  1:723106/139‑1 (MQ=255)                                                                                   aCGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCGATCAGATAGTtgc                                                                     <  1:1233609/121‑1 (MQ=255)                                                                                   aCGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTtgc                                                                     >  2:1233609/1‑121 (MQ=255)                                                                                            atcTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGt                                          >  1:606741/1‑139 (MQ=255)                                                                                             tctcACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGAt                                                                           >  1:972173/1‑105 (MQ=255)                                                                                             tctcACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGAt                                                                           <  2:972173/105‑1 (MQ=255)                                                                                              ctcACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCcagcag                                                           <  1:998782/120‑1 (MQ=255)                                                                                              ctcACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCcagcag                                                           >  2:998782/1‑120 (MQ=255)                                                                                              ctcACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTtg                                        >  2:1000432/1‑139 (MQ=255)                                                                                              ctcACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTtg                                        >  1:1396341/1‑139 (MQ=255)                                                                                               tcACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTtgt                                       <  2:1338139/139‑1 (MQ=255)                                                                                                 aCGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTgg                                     >  1:1824700/1‑139 (MQ=255)                                                                                                    cAGAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCt                                  <  1:1000432/139‑1 (MQ=255)                                                                                                      gAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCg                                                    >  1:142915/1‑119 (MQ=255)                                                                                                      gAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCg                                                    <  2:142915/119‑1 (MQ=255)                                                                                                      gAATACCGTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGACTGTAGTTGTGGCCttt                                <  1:863308/139‑1 (MQ=255)                                                                                                             gTTGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGc                         >  2:1196330/1‑139 (MQ=255)                                                                                                              ttGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTAc                                                                                                      <  1:1649002/61‑1 (MQ=255)                                                                                                              ttGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTAc                                                                                                      >  2:1649002/1‑61 (MQ=255)                                                                                                              ttGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTtg                                                                      <  2:816527/93‑1 (MQ=255)                                                                                                              ttGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTtg                                                                      >  1:816527/1‑93 (MQ=255)                                                                                                              ttGCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCa                        <  1:1196330/139‑1 (MQ=255)                                                                                                                gCACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAgg                      <  1:391514/139‑1 (MQ=255)                                                                                                                 cACGGCATGTGTACGCGGTGGTGGTCACGCGTGGAGAGGTAAGTGGTCACAAACGTAc                                                                                                      >  1:791553/1‑58 (MQ=255)                                                                                                                 cACGGCATGTGTACGCGGTGGTGGTCACGCGTGGAGAGGTAAGTGGTCACAAACGTAc                                                                                                      <  2:791553/58‑1 (MQ=255)                                                                                                                 cACGGCATGTGTACGCGGTGGTGGTCACGCGGGga                                                                                                                             >  1:1291645/1‑35 (MQ=38)                                                                                                                 cACGGCATGTGTACGCGGTGGTGGTCACGCGGGga                                                                                                                             <  2:1291645/35‑1 (MQ=38)                                                                                                                 cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTtgt                                       <  2:743450/121‑1 (MQ=255)                                                                                                                 cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTtgt                                       >  1:743450/1‑121 (MQ=255)                                                                                                                 cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCtt                                 <  2:1169501/127‑1 (MQ=255)                                                                                                                 cACGGCATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGAGCTTCGAGGCTGTAGTTGTGGCCtt                                 >  1:1169501/1‑127 (MQ=255)                                                                                                                      cATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGa                                                   <  2:377409/104‑1 (MQ=255)                                                                                                                      cATGTGTACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGa                                                   >  1:377409/1‑104 (MQ=255)                                                                                                                         gtttACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCAt                                                                                <  1:753565/69‑1 (MQ=255)                                                                                                                         gtgtACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCAt                                                                                >  2:753565/1‑72 (MQ=255)                                                                                                                           gtACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTg                                                                                                 >  1:1611415/1‑53 (MQ=255)                                                                                                                           gtACGCGGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTg                                                                                                 <  2:1611415/53‑1 (MQ=255)                                                                                                                                cgGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGc                             >  1:826009/1‑116 (MQ=255)                                                                                                                                cgGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGc                             <  2:826009/116‑1 (MQ=255)                                                                                                                                cgGTGGTGGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATCTTCTTCGa      >  1:292062/1‑139 (MQ=255)                                                                                                                                  gtggtgGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                      <  2:59196/57‑1 (MQ=255)                                                                                                                                  gtggtgGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                      >  1:59196/1‑57 (MQ=255)                                                                                                                                    ggtggtCACGCGGGGAGAGGTAAGTGGTCACAAACGTAc                                                                                                      <  2:1017726/39‑1 (MQ=255)                                                                                                                                    ggtggtCACGCGGGGAGAGGTAAGTGGTCACAAACGTAc                                                                                                      >  1:1017726/1‑39 (MQ=255)                                                                                                                                    ggtggtCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                      <  1:121127/55‑1 (MQ=255)                                                                                                                                    ggtggtCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTc                                                                                      >  2:121127/1‑55 (MQ=255)                                                                                                                                    ggtggtCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCTATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATCTTCTTCGAtttt  >  2:1695904/1‑139 (MQ=255)                                                                                                                                    ggtggtCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGa                                                   <  2:1149163/90‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        GCAAATCACGCGTTCGTTACGGGCAAACAGATTCGGCACGTTCTGAGCCGTGAGATGGTTAACGGAGAAGAGATCGCCCGGCACGTAGATCATCTCACGCAGAATACCGTTGCACGGCATGTGTACGCGGTGGTAGTCACGCGGGGAGAGGTAAGTGGTCACAAACGTACCGTTGCGGAACAGGTCCGCCATCAGATAGTTGCCTGCCAGCAGGGCTTCGAGGCTGTAGTTGTGGCCTTTGGCTTGCAGGATTTTATCTTCTTCGATTTT  >  NC_000913/4389800‑4390069

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