Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,026,802	A→T	6.7%	S22C (AGT→TGT)	ocd →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,026,802	0	A	T	6.7%	280.5 / 8.6	89	S22C (AGT→TGT)	ocd	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base T (2/4);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

TCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCGCTATCA  >  NC_002947/5026660‑5026949                                                                                                                                                |                                                                                                                                                    tCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGa                                                                                                                                              <  1:850627/150‑1 (MQ=255) tCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGa                                                                                                                                              <  2:782637/150‑1 (MQ=255)  cGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTAAGTGCGTAACGCTCGATACTGTCTGTATCGag                                                                                                                                             <  1:1770418/150‑2 (MQ=255)   gACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCCTGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGagc                                                                                                                                            <  2:466334/150‑3 (MQ=255)   gACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCCTGAACTGTGAGTGCGTAACGCTCGATACAGACAGTATCGATg                                                                                                                                            >  1:1833705/1‑150 (MQ=255)    aCATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                              <  1:442247/147‑1 (MQ=255)    aCATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                              >  2:442247/1‑147 (MQ=255)    aCATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGt                                                                                                                                           >  2:1216503/1‑150 (MQ=255)      aTGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAg                                                                                                                                         <  1:2349034/150‑1 (MQ=255)        ggCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAACACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTg                                                                                                                                       >  2:1073045/1‑150 (MQ=255)         gCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTgg                                                                                                                                      <  1:711712/150‑1 (MQ=255)             aaCACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGaaaa                                                                                                                                  <  1:504335/150‑1 (MQ=255)                acCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                              >  2:90394/1‑135 (MQ=255)                acCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                              <  1:90394/135‑1 (MQ=255)                  cAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCAt                                                                                                                             >  2:617286/1‑150 (MQ=255)                    aCGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATAAGTGAGTGCGTAACGCTCGATACAGACAGTATCGCTGTACTGGTTAACTCAttt                                                                                                                           <  2:703315/150‑1 (MQ=255)                       catcatCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCa                                                                                                                        >  1:1664237/1‑150 (MQ=255)                            tcCTGGAAGCTGTGCGCGACCGCCGCCCAACACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCt                                                                                                                   <  1:1073045/150‑1 (MQ=255)                             ccTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTg                                                                                                                  <  1:236167/150‑1 (MQ=255)                               tGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGc                                                                                                                <  1:1216503/150‑1 (MQ=255)                                   aGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCaga                                                                                                            <  1:2145521/150‑1 (MQ=255)                                    gCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCagag                                                                                                           <  2:2058041/150‑1 (MQ=255)                                    gCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCagag                                                                                                           <  1:313952/150‑1 (MQ=255)                                      tgtgCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTAt                                                                                                                                                 >  2:1606887/1‑110 (MQ=255)                                      tgtgCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTAt                                                                                                                                                 <  1:1606887/110‑1 (MQ=255)                                      tgtgCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTCAGCGAAGT‑CCTCACGCTAACCCAACACGACATTCTTCCGTTCTTAACGGTGGACACGGACAGTATCGTGGTAGTGTAAAACGCATTTGCATCGCTGGCCAGAGGc                                                                                                         <  2:274907/150‑1 (MQ=255)                                             gACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTgg                                                                                                  >  1:2075877/1‑150 (MQ=255)                                               ccgccgCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTc                                                                                                <  2:801935/150‑1 (MQ=255)                                               cagccgccCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGAAAAGTGGTc                                                                                                >  2:2337981/3‑150 (MQ=255)                                                 gccgccCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGt                                                                                                    >  2:113598/1‑144 (MQ=255)                                                 gccgccCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGt                                                                                                    <  1:113598/144‑1 (MQ=255)                                                    gccCCGCACGTGCTTAAAGAGGTGCCCTGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGcc                                                                                           >  2:2220216/1‑150 (MQ=255)                                                      ccAGCACGTGCTTAAAGATGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAAAGCGAGATCCGTGAGTGCGCAACGCGCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGt                                                                                         >  1:1499198/1‑150 (MQ=255)                                                         gCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAGCGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCaa                                                                                      >  2:596051/1‑150 (MQ=255)                                                             gTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAAtcct                                                                                  >  2:112338/1‑150 (MQ=255)                                                             gTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAAtcct                                                                                  >  2:1013120/1‑150 (MQ=255)                                                               gCTTAAAGAGGTGCCATGTAAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAAAGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAAtactcc                                                                                >  2:495151/1‑150 (MQ=255)                                                               gCTTAAAGAGGTGCAATGTGAGCGAAGT‑CCTTCTGCTTAACCACCGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAAtcctcc                                                                                <  1:951607/150‑1 (MQ=255)                                                                cTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCg                                                                               <  1:617286/150‑1 (MQ=255)                                                                 ttAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATATTGACAGTATCGATGTAGTGGGAAACCCATTTGCATCGCTGGCCAGACGCACAGTGGTCATGCCGTCAATCCTCCgg                                                                              >  2:1389823/1‑150 (MQ=255)                                                                 ttAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCgg                                                                              <  1:112338/150‑1 (MQ=255)                                                                  tAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAAACAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGt                                                                             >  1:440850/1‑150 (MQ=255)                                                                    aaGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTa                                                                           >  1:2067544/1‑150 (MQ=255)                                                                    aaGAGGTGCCATGTGA‑CCAAGTGTTTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTa                                                                           <  1:1328582/150‑1 (MQ=255)                                                                            ccATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCaa                                                                                      <  2:1974836/131‑1 (MQ=255)                                                                            ccATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCaa                                                                                      >  1:1974836/1‑131 (MQ=255)                                                                               tgtgAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGa                                                                >  1:524864/1‑150 (MQ=255)                                                                               tgtgAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGa                                                                <  2:1664237/150‑1 (MQ=255)                                                                               tgtgAGCGAAGT‑CCTTCTGCTTAAACAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGa                                                                <  2:440850/150‑1 (MQ=255)                                                                                gtgAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCaga                                                                                                            >  2:1122437/1‑105 (MQ=255)                                                                                gtgAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCaga                                                                                                            <  1:1122437/105‑1 (MQ=255)                                                                                 tgAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGc                                                              >  2:315221/1‑150 (MQ=255)                                                                                    gCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCg                                                                               <  2:1909792/130‑1 (MQ=255)                                                                                    gCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCg                                                                               >  1:1909792/1‑130 (MQ=255)                                                                                     cGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACaa                                                          <  2:1499198/150‑1 (MQ=255)                                                                                         gT‑CCTTCTTCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGt                                                                                         <  1:1047241/115‑1 (MQ=255)                                                                                         gT‑CCTTCTTCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGt                                                                                         >  2:1047241/1‑115 (MQ=255)                                                                                          t‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGtt                                                                            <  2:835236/127‑1 (MQ=255)                                                                                          t‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGtt                                                                            >  1:835236/1‑127 (MQ=255)                                                                                          t‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGa                                                                         >  2:734339/1‑130 (MQ=255)                                                                                          t‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGa                                                                         <  1:734339/130‑1 (MQ=255)                                                                                                cTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTg                                                <  2:524864/150‑1 (MQ=255)                                                                                                 tGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTACGAGCACAATGGCGAGGTgg                                               >  2:1097213/1‑150 (MQ=255)                                                                                                  gCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGa                                              <  2:1277340/150‑1 (MQ=255)                                                                                                   cTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTa                                                                           >  1:206520/1‑120 (MQ=255)                                                                                                   cTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTa                                                                           <  2:206520/120‑1 (MQ=255)                                                                                                   cTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGAc                                             <  1:596051/150‑1 (MQ=255)                                                                                                   cTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGAc                                             <  2:2067544/150‑1 (MQ=255)                                                                                                        ccAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAg                                                                                                     >  1:109339/1‑89 (MQ=255)                                                                                                        ccAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAg                                                                                                     <  2:109339/89‑1 (MQ=255)                                                                                                         cAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCaaa                                       >  2:1349386/1‑150 (MQ=255)                                                                                                         cAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCaaa                                       >  2:816877/1‑150 (MQ=255)                                                                                                          aaCGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCaaaa                                      >  1:1579918/1‑150 (MQ=255)                                                                                                               gAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCt                                 <  1:1013120/150‑1 (MQ=255)                                                                                                                      gTGAGTGCGTAACGGTCAAGACAGACAGTATCGATGTAGTGGAAAACGCGTTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGATTACCTGc                          >  1:695122/1‑150 (MQ=255)                                                                                                                      gTGAGTGCATAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGc                          <  1:816877/150‑1 (MQ=255)                                                                                                                         aGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGt                                                                                         <  1:1486977/84‑1 (MQ=255)                                                                                                                         aGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGt                                                                                         >  2:1486977/1‑84 (MQ=255)                                                                                                                              gTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCtt                  >  2:1818040/1‑150 (MQ=255)                                                                                                                              gTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCtt                  >  2:1608634/1‑150 (MQ=255)                                                                                                                              gTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCAAACTTGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCtt                  >  1:634691/1‑150 (MQ=255)                                                                                                                                aaCGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTga                <  1:1349386/150‑1 (MQ=255)                                                                                                                                     tCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCaa                                                                                      >  1:1840567/1‑75 (MQ=255)                                                                                                                                     tCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCaa                                                                                      <  2:1840567/75‑1 (MQ=255)                                                                                                                                     tCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGCTCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAAtt           >  1:658982/1‑150 (MQ=255)                                                                                                                                       gctACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGCCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAACAAGGCGTACCTGCCAGACCACGGGAATTTc         >  1:506680/3‑150 (MQ=255)                                                                                                                                       gATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTc         <  2:1579918/150‑1 (MQ=255)                                                                                                                                        aTACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCg        <  2:2075877/150‑1 (MQ=255)                                                                                                                                        aTACACACAGTACCGATGTAGTGGAAAACGCATTTCCATCGCTGGCCAGAGGCAAAGTGCTCATGTCGCCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGCACGTCAAAACCGCTTACATGCCAGACCTAGAGACTTTCg        <  2:506680/150‑1 (MQ=255)                                                                                                                                              acagTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCGCTATCa  >  1:1194060/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                    TCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGT‑CCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCGCTATCA  >  NC_002947/5026660‑5026949

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