Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_pgi	1,507,116	A→G	E112E (GAA→GAG)	ydcP →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_pgi	1,507,116	0	A	G	98.8%	242.2 / NA	84	E112E (GAA→GAG)	ydcP	putative peptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (48/35);  minor base T (0/1);  total (48/36)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATTTTCGTCACGCTTAACACCATTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAACTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCAGGCTACGGACGCGACCA  >  NC_000913_3_hsa_pgi/1506970‑1507258                                                                                                                                                   |                                                                                                                                               aaaaTTTTCGTCACGCTTAACACCATTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGGCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCt                                                                                                                                              <  2:1366530/149‑1 (MQ=255)             aCGCTTAACACCATTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTTACCTCTACAAGACCGGTGTCGATGGGCTGATTGTTCAGTAGATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTAc                                                                                                                                  >  2:1135933/1‑149 (MQ=255)                 ttAACACCATTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAg                                                                                                                              >  1:1087151/1‑149 (MQ=255)                  tAACACCATTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGt                                                                                                                             <  1:8398/149‑1 (MQ=255)                   aaCACCATTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTg                                                                                                                            >  1:1085839/1‑149 (MQ=255)                     cacCATTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCg                                                                                                                          >  2:1173795/1‑149 (MQ=255)                     cacCATTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCg                                                                                                                          >  2:1002207/1‑149 (MQ=255)                        cATTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACa                                                                                                                       >  1:633634/1‑149 (MQ=255)                        cATTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATAACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCTGTACGCAGTGCGACa                                                                                                                       >  1:1032180/1‑149 (MQ=255)                         aTTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACAt                                                                                                                      >  2:708034/1‑149 (MQ=255)                         aTTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACAt                                                                                                                      <  2:1488825/149‑1 (MQ=255)                           tttGCATGATGAGGAGCTGGAACCCGCGCACCGGCTGATTACTGACCTCTACCAGACCGGGGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTc                                                                                                                    <  1:738707/149‑1 (MQ=255)                              gcctgatgatgaGCTGGCACCCGCGCAACGCCTGATTACTGACCTCTACCAGCCCGGGGGCGATGCGCTGCTTGTTCAGGATCTGGGGCTGCTGGTTCTTGATATTCCGCCGATTGAGCTGCAGGCCAGTTCGCAGTGCGACATTCGTa                                                                                                                 <  2:1083026/146‑1 (MQ=255)                               cATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGAACTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTCCGCAGTGAGACATTCGTAc                                                                                                                >  2:420813/1‑149 (MQ=255)                               cATGATGATGAGCTGGAACCCGCGCAACGGCTGAGTACTGACCTCTACCAGACCGGTGTCGATGCGCTGAGTGTTCAGGATATGGGGATTCTGGAACTTGATATTGCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACAATCGtcc                                                                                                                >  2:1011396/1‑147 (MQ=255)                                  gatgatgaGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGt                                                                                                             <  2:1357597/149‑1 (MQ=255)                                  gatgatgaGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGAGCTGATTGTTCAGGATATCGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGt                                                                                                             >  2:790818/1‑149 (MQ=255)                                      atgaGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTAAAGTTGaa                                                                                                         >  2:1440317/1‑149 (MQ=255)                                         ttcTGCAACCCGGGCAACGGCTCATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGaaaaa                                                                                                      <  2:1331974/147‑1 (MQ=255)                                          gCTGTAACCCGCGCAAAGGCTGATTACTGAACTCTACCCGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAg                                                                                                     <  2:428195/149‑1 (MQ=255)                                              gAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCAc                                                                                                                                          >  2:17011/1‑108 (MQ=255)                                              gAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCAc                                                                                                                                          <  1:17011/108‑1 (MQ=255)                                                aCCCGCGCAACGGCTGATTACTGACCTCTACCAGCCCGGTGTGGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCAAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGt                                                                                               >  1:952056/1‑149 (MQ=255)                                                  ccGCGCAACGGCGGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTc                                                                                             <  2:266414/149‑1 (MQ=255)                                                   cgcgCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTcc                                                                                            >  1:1161254/1‑149 (MQ=255)                                                   cgcgCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTAcc                                                                                            >  1:626272/1‑149 (MQ=255)                                                    gcgcAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTGCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCct                                                                                           <  2:1277390/149‑1 (MQ=255)                                                      gcAACGGCTGATTACTGACCTCTACCAGACCGGTGGCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCctct                                                                                         <  1:1417398/149‑1 (MQ=255)                                                        aaCGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGTATATGGTGATTCTGGAACTTGATATTCCGCCTATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCctctct                                                                                       <  1:1407211/149‑1 (MQ=255)                                                        aaCGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCctctct                                                                                       <  2:1161254/149‑1 (MQ=255)                                                         aCGGCTGATTACTGACCTCTACCAGACCGGTTTCGATGCGCTGATTGTTCAGGTTATGATGATTGTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACTCAGAGCGACATTCGTACAGGTGAAAAAGCGAAGTTCCTCTCTg                                                                                      <  1:1109148/149‑1 (MQ=255)                                                              tGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGAATCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGtt                                                                                 <  2:550624/149‑1 (MQ=255)                                                               gATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTg                                                                                >  1:909531/1‑149 (MQ=255)                                                                aTTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTgg                                                                               <  2:517143/149‑1 (MQ=255)                                                                    cTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTc                                                                           >  1:174538/1‑149 (MQ=255)                                                                     tGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCa                                                                          >  2:489594/1‑149 (MQ=255)                                                                      gACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCAc                                                                         <  1:1335099/149‑1 (MQ=255)                                                                         ctctACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCa                                                                      >  2:627325/1‑149 (MQ=255)                                                                           ctACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGa                                                                    >  2:749226/1‑149 (MQ=255)                                                                            tACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGCTGTTGGCTTCACGCAGAt                                                                   >  2:1009924/1‑149 (MQ=255)                                                                             aCCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGAtt                                                                  <  1:1274812/149‑1 (MQ=255)                                                                               cAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTACACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATtgt                                                                >  2:1292450/1‑149 (MQ=255)                                                                                aGACCGGTGTCGATGCGCTGCTTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATtgtg                                                               >  2:1201397/1‑149 (MQ=255)                                                                                 gACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGc                                                              >  2:1116428/1‑149 (MQ=255)                                                                                 gACCGGTGGCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGc                                                              <  1:1293055/149‑1 (MQ=255)                                                                                      gtgtCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGgcg                                                         >  2:873880/1‑149 (MQ=255)                                                                                        gtCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGcgcg                                                       >  2:547512/1‑149 (MQ=255)                                                                                        gtCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGcgcg                                                       >  2:398215/1‑149 (MQ=255)                                                                                           gaagcgcTGATTGTTCCGGATATTGGGATTCTGGAACTTGATATTCCGCCGATTGATCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAATTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCgaga                                                    <  2:633227/146‑1 (MQ=255)                                                                                           gATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCgaga                                                    >  1:644866/1‑149 (MQ=255)                                                                                           gATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCgaga                                                    <  2:53138/149‑1 (MQ=255)                                                                                             tGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGc                                                  <  2:1105518/149‑1 (MQ=255)                                                                                             tGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGc                                                  >  2:777454/1‑149 (MQ=255)                                                                                              gcgcTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCt                                                 >  2:29988/1‑149 (MQ=255)                                                                                                gcTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTTCACGCCAGTACGCAGTGCGACAGTCGAACAGTTGAAAAAGCGAAGTTCCTCTCTTATGTTGGCTTCACGCAGATAGATCTGTCGCGAGAGCTGa                                               >  1:1012596/1‑149 (MQ=255)                                                                                                gcTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGa                                               <  1:186208/149‑1 (MQ=255)                                                                                                gcTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTCGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGa                                               <  2:1085839/149‑1 (MQ=255)                                                                                                   gATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGAAGTGCGCCATTCGTACAGTTGAAAAAGCGATGTTCCTCTCTTAAGTTGTCTAAACTCAGATTGTGCTGGCGCGAGAGCTGAATc                                            >  2:1453701/1‑149 (MQ=255)                                                                                                    aTTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCt                                           >  1:421222/1‑149 (MQ=255)                                                                                                    aTTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCt                                           >  1:228596/1‑149 (MQ=255)                                                                                                    aTTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCt                                           >  1:526147/1‑149 (MQ=255)                                                                                                      tggtCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATGGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTg                                         <  1:500005/146‑1 (MQ=255)                                                                                                        ttCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGAt                                       >  2:34207/1‑149 (MQ=255)                                                                                                                   gggATTCTGTATCTTGATATTCAGCCGATGGAGCTACACTCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGATCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCcgcgc                            <  2:270954/149‑1 (MQ=255)                                                                                                                   gggATTCTGGTACGTGATCTACCGCCCATTGAGCTACACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCGGGCGCGAGAGCGGAATCTTGATCAGATCcgcgc                            <  2:130892/149‑1 (MQ=255)                                                                                                                     gATTCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGa                          >  1:1212490/1‑149 (MQ=255)                                                                                                                       ttCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACAATCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGAtt                        <  1:1467598/149‑1 (MQ=255)                                                                                                                        tCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCATCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTc                       <  1:29988/149‑1 (MQ=255)                                                                                                                        tCTGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGATCCTCTCTGATGTTGGCTTCACGCAGAATGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTc                       >  2:611031/1‑149 (MQ=255)                                                                                                                          tGGAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCAc                     >  1:500492/1‑149 (MQ=255)                                                                                                                           ggAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTcccc                    >  1:482231/1‑149 (MQ=255)                                                                                                                            gAACTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCa                   <  1:547512/149‑1 (MQ=255)                                                                                                                             aaCTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCAg                  >  1:1001430/1‑149 (MQ=255)                                                                                                                               cTTGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGGAGATTGTGCTGGCGCGAGAGCTGCATCTTGATCAGATCCGCGCGATTCACCAGGc                >  2:1175947/1‑149 (MQ=255)                                                                                                                                ttGATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCAGGCt               <  1:708034/149‑1 (MQ=255)                                                                                                                                  gATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGAt                                       >  1:817041/1‑123 (MQ=255)                                                                                                                                  gATATTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGAt                                       <  2:817041/123‑1 (MQ=255)                                                                                                                                   atatTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGAGATGCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCAGGCTACg            <  1:1367291/149‑1 (MQ=255)                                                                                                                                    tatTCCGCCGATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCAGGCTACgg           <  1:1173795/149‑1 (MQ=255)                                                                                                                                        ccgccgATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCAGGCTACGGAcgc       <  1:1135933/149‑1 (MQ=255)                                                                                                                                         cgccgATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCAGGCTACGGAcgcg      >  1:724778/1‑149 (MQ=255)                                                                                                                                         cgccgATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCAGGCTACGGAcgcg      >  2:705531/1‑149 (MQ=255)                                                                                                                                             gATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCAGGCTACGGACGCGACCa  >  1:1398219/1‑149 (MQ=255)                                                                                                                                             gATTGAGCTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCAGGCTACGGACGCGACCa  >  2:605219/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                               AAAATTTTCGTCACGCTTAACACCATTTTGCATGATGATGAGCTGGAACCCGCGCAACGGCTGATTACTGACCTCTACCAGACCGGTGTCGATGCGCTGATTGTTCAGGATATGGGGATTCTGGAACTTGATATTCCGCCGATTGAACTGCACGCCAGTACGCAGTGCGACATTCGTACAGTTGAAAAAGCGAAGTTCCTCTCTGATGTTGGCTTCACGCAGATTGTGCTGGCGCGAGAGCTGAATCTTGATCAGATCCGCGCGATTCACCAGGCTACGGACGCGACCA  >  NC_000913_3_hsa_pgi/1506970‑1507258

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