Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,661,849	A→G	V32A (GTC→GCC)	unknown ←	[4Fe‑4S] cluster carrier protein NfuA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,661,849	0	A	G	100.0%	173.1 / NA	62	V32A (GTC→GCC)	unknown	[4Fe‑4S] cluster carrier protein NfuA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/42);  total (20/42)

CGGCGCATCGTCCGCGACTTTTCTCATTTTGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGACTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCCACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAGGGAGCGTTATGTACGGGCAAGCCCCCACACCTGAATATGTAC  >  1/1661620‑1662095                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       cggcgCATCGTCCGCGACTTTTCTCATTTTGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATg                                                                                                                                                                                                                                     <  1:460191/249‑1 (MQ=255) cggcgCATCGTCCGCGACTTTTCTCATTTTGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAg                                                                                                                                                                                                                                   <  1:893059/251‑1 (MQ=255) cggcgCATCGTCCGCGACTTTTCTCATTTTGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAg                                                                                                                                                                                                                                   <  1:85301/251‑1 (MQ=255)           tCCGCGACTTTTCTCATTTTGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCttcttc                                                                                                                                                                                                                         <  1:97053/251‑1 (MQ=255)           tCCGCGACTTTTCTCATTTTGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCttcttc                                                                                                                                                                                                                         <  1:558436/251‑1 (MQ=255)            ccGCGACTTTTCTCATTTTGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTcc                                                                                                                                                                                                                        <  1:213257/251‑1 (MQ=255)                 aCTTTTCTCATTTTGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTg                                                                                                                                                                                                                   <  1:633110/251‑1 (MQ=255)                    tttCTCATTTTGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATc                                                                                                                                                                                                                <  1:463958/251‑1 (MQ=255)                           ttttGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATaaa                                                                                                                                                                                                         >  1:74725/1‑251 (MQ=255)                           ttttGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATaaa                                                                                                                                                                                                         >  1:704122/1‑251 (MQ=255)                            tttGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAAc                                                                                                                                                                                                        <  1:145596/251‑1 (MQ=255)                            tttGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAAc                                                                                                                                                                                                        <  1:561098/251‑1 (MQ=255)                                          cTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAatat                                                                                                                                                                                                                                      <  1:225126/207‑1 (MQ=255)                                          cTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAAc                                                                                                                                                                                                        >  1:97025/1‑237 (MQ=255)                                          cTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATgcg                                                                                                                                                                                            >  1:263914/1‑249 (MQ=255)                                          cTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATgcg                                                                                                                                                                                            >  1:363194/1‑249 (MQ=255)                                              tAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAAc                                                                                                                                                                                                        >  1:362414/1‑233 (MQ=255)                                               aaGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTg                                                                                                                                                                                                                   <  1:729778/221‑1 (MQ=255)                                               aaGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTgcg                                                                                                                                                                                      >  1:899883/1‑250 (MQ=255)                                                            cgcCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAAt                                                                                                                                                                         <  1:859166/250‑1 (MQ=255)                                                                     cttcGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTCGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATttgt                                                                                                                                                               <  1:436342/251‑1 (MQ=255)                                                                       tcGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATaa                                                                                                                                                                                                          >  1:710597/1‑206 (MQ=255)                                                                       tcGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGaa                                                                                                                                                                           <  1:208551/237‑1 (MQ=255)                                                                          gTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGaaa                                                                                                                                                                          >  1:643042/1‑235 (MQ=255)                                                                              cGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATg                                                                                                                                                                                                                                     <  1:242596/172‑1 (MQ=255)                                                                                tAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAActc                                                                                                                                                     <  1:113480/250‑1 (MQ=255)                                                                                  aTCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACtctc                                                                                                                                                   >  1:420740/1‑250 (MQ=255)                                                                                     aaTTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCAGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTcctac                                                                                                                                               <  1:630359/251‑4 (MQ=255)                                                                                      aTTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCtaca                                                                                                                                              >  1:504657/1‑251 (MQ=255)                                                                                       tttCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACaa                                                                                                                                             <  1:241811/251‑1 (MQ=255)                                                                                                        aaCGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAAtt                                                                                                                            <  1:858716/251‑1 (MQ=255)                                                                                                               ggCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTgcgc                                                                                                                                                                                     >  1:876136/1‑187 (MQ=255)                                                                                                                gCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCgg                                                                                                                      <  1:705569/249‑1 (MQ=255)                                                                                                                     tttCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGaaa                                                                                                                                                                          <  1:353950/192‑1 (MQ=255)                                                                                                                        cGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACg                                                                                                            >  1:5407/1‑251 (MQ=255)                                                                                                                          tCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGcc                                                                                                          <  1:134867/251‑1 (MQ=255)                                                                                                                                  ggCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTcctac                                                                                                                                               <  1:13429/206‑4 (MQ=255)                                                                                                                                      gagaAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTcctac                                                                                                                                               <  1:413847/202‑4 (MQ=255)                                                                                                                                      gagaAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGtttt                                                                                              >  1:790803/1‑251 (MQ=255)                                                                                                                                         aaaCTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGtttttt                                                                                            <  1:715497/250‑1 (MQ=255)                                                                                                                                              gTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGAt                                                                                                               <  1:167884/226‑1 (MQ=255)                                                                                                                                              gTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTa                                                                                       >  1:556632/1‑250 (MQ=255)                                                                                                                                                 gTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAAtt                                                                                                                                                                  >  1:351937/1‑172 (MQ=255)                                                                                                                                                                gtcgtcGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTAc                                                                               <  1:276236/240‑1 (MQ=255)                                                                                                                                                                   gtcgGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATCCGGCGTGCCGGGATAGGCTACGAAAATACGAATATGAGGGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACGCTcctac                                                                                                                                               >  1:237137/1‑170 (MQ=255)                                                                                                                                                                        cTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAata                                                                                                                                                                                                                                       <  1:406722/80‑1 (MQ=255)                                                                                                                                                                                       gTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCcttttc                                                   <  1:144035/245‑1 (MQ=255)                                                                                                                                                                                          ggCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAgg                                           <  1:748637/250‑1 (MQ=255)                                                                                                                                                                                           gCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTccta                                                                                                                                                <  1:375576/148‑3 (MQ=255)                                                                                                                                                                                           gCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATa                                                                                                              >  1:583049/1‑182 (MQ=255)                                                                                                                                                                                                    aaCGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCATGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAGGGAGCGTTATGt                                <  1:179338/251‑1 (MQ=255)                                                                                                                                                                                                    aaCGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAGGGAGCGTTATGt                                <  1:240774/251‑1 (MQ=255)                                                                                                                                                                                                      cGCCACATTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAGGGAGCGTTATGTAc                              <  1:556972/251‑1 (MQ=255)                                                                                                                                                                                                          acaTTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATAc                                                                                                             >  1:397423/1‑168 (MQ=255)                                                                                                                                                                                                          acaTTCGGCATTCGGCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGc                                                                                                           <  1:168526/170‑1 (MQ=255)                                                                                                                                                                                                           cattcggcattcggCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAGGGAGCGTTATGTACGGGCa                         <  1:308384/251‑1 (MQ=255)                                                                                                                                                                                                           cattcggcattcggCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAGGGAGCGTTATGTACGGGCa                         <  1:196848/251‑1 (MQ=255)                                                                                                                                                                                                           cattcggcattcggCGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAGGGAGCGTTATGTACGGGCa                         <  1:442788/251‑1 (MQ=255)                                                                                                                                                                                                             ttcggcattcggcGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAGGGAGCGTTATGTACGGGCAAg                       >  1:842373/1‑251 (MQ=255)                                                                                                                                                                                                                        gcGTGCCGGGATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAGGGAGCGTTATGTACGGGCAAGCCCCCACAcc             <  1:460267/250‑1 (MQ=255)                                                                                                                                                                                                                                 gATTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAggg                                          <  1:234697/212‑1 (MQ=255)                                                                                                                                                                                                                                  aTTGGCTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCTACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAGGGAGCGTTATGTACGGGCAAGCCCCCACACCTGAATATGTAc  <  1:793545/251‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       CGGCGCATCGTCCGCGACTTTTCTCATTTTGGCATTCGGCGCTTTTAAGGTGAGCTGCGCGCCCAGTTCTTCGGTGACGTAATCAATTTCCGCATCTTCCAGAAACGGCAGGCTGATTTCGTCTACAAAGGCGGAGAAACTGTCGTATTTGATTTCCATGTCGTCGGCTTCAACGGAATTCGGTGGGCAGTAAGAAACGCCACATTCGGCATTCGGCGTGCCGGGATTGACTACGAAAATACGAATATGAGTGCCTTCTTCCTGTTGATCCAATAAACGGCGAAAATGCGCTTGCGCCGCTTCTGAAATAGTAATTTGTTGCATAACTCTCCCACAATAAAATCAGATTAAATTTGTCGGATATGATACGCCTCACGCCGTTTTTTTTCTACCTTTTACTATAATAAAAAAGCGGATTAGCCTTTTCTTTTAAGGGAGCGTTATGTACGGGCAAGCCCCCACACCTGAATATGTAC  >  1/1661620‑1662095

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