Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_pgi	526,792	0	A	C	22.1%	122.6 / 24.7	68	S1178R (AGC→CGC)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base C (5/10);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGG  >  NC_000913_3_hsa_pgi/526649‑526938                                                                                                                                                |                                                                                                                                                   ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                               <  2:126535/149‑1 (MQ=255)   gcgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGa                                                                                                                                             <  1:373745/149‑1 (MQ=255)     gcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGagt                                                                                                                                           <  2:256194/149‑3 (MQ=255)       aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                         >  2:1225203/1‑149 (MQ=255)       aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAatat                                                                                                                                         <  2:906918/149‑1 (MQ=255)         tGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                       <  1:474931/149‑1 (MQ=255)         tGCGGGCTGACGGTGGAGCAACAGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatga                                                                                                                                       <  1:514758/149‑1 (MQ=255)          gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                      <  1:1165257/149‑1 (MQ=255)              gCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGCCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGATTATGATGAAt                                                                                                                                  >  2:34677/1‑149 (MQ=255)                 gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                               <  1:570329/149‑1 (MQ=255)                 gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGGGGAGCGCGGAATATGATGAATggg                                                                                                                               <  2:468160/149‑1 (MQ=255)                  aCGGTGGAGCAACTGGTCAGCCAGGGGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                              <  1:102047/149‑1 (MQ=11)                     gTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCaa                                                                                                                           >  1:119986/1‑149 (MQ=255)                      tGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                          >  2:362877/1‑149 (MQ=255)                        gAGCAACTGGCCAGACAGGTGGAGCCGGAATACCCACCGGCGCGAAAAGCTCCTCTGTATCCCTTCGACCACCGGGGACTGCCGCGGGCGCTTATCAGCGAAGACGGCACAACGGCGTGGAGCGCGTAAGAGGAAGACTGGGGCACCCa                                                                                                                        >  1:687534/1‑149 (MQ=17)                          ccAACTGGCCAGACAGGGGGAGCCGGAATACACTCCGGCGCGAAAAGCTCATCTTTATCCCTGCGACCACCGGGGACGGCCGCTGGCGCTTATCACCGAAGACGGCAATACGGCGGGGCGCGGGGAGTATGATGAATGGGGCAACCAGc                                                                                                                      <  2:633006/148‑1 (MQ=11)                           cAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAACCTAATCTTTAACACTGAGACCACCGGGGACTGCCGCTGGCGCTTTTCAGCGAGGACGGCAATTCGTCGTGGCGCGGGGAGTATGAAGAATGGGGCAACCAGcg                                                                                                                     >  1:547254/1‑148 (MQ=11)                              cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGAAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTaa                                                                                                                  <  1:1126188/149‑1 (MQ=255)                               tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCTCCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATAAGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAAt                                                                                                                 <  1:915839/149‑1 (MQ=255)                               tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAAt                                                                                                                 <  1:188240/149‑1 (MQ=255)                                  ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgag                                                                                                              >  1:882053/1‑149 (MQ=255)                                    aGACAGGTGGGGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                            >  1:1113331/1‑149 (MQ=255)                                    aGACAGGTGGAGCCTGAATACCCACCGGCGCGAAACGCTCATCTTTATCACTGCGACCACCGTGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATCTGATGAATGGGGCAACCAGCTTAATgagca                                                                                                            <  1:965902/149‑3 (MQ=17)                                    aGACAGGGGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                            <  1:965904/149‑1 (MQ=255)                                          gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGCCCCCCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                      >  2:83334/1‑149 (MQ=11)                                              aGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcat                                                                                                  <  2:1033459/149‑1 (MQ=255)                                               gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatc                                                                                                 <  1:346662/149‑1 (MQ=255)                                                 cGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcac                                                                                               <  1:125676/149‑2 (MQ=255)                                                 cGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcat                                                                                               >  1:583413/1‑149 (MQ=255)                                                  ggAATACACACCGGCGCGATTAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCTCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTATTGAGGAGAACCCGCATCAtg                                                                                              <  2:447407/149‑1 (MQ=255)                                                   gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGCCGGCAATACGGCGTGGAGCGCGGAATATGTTGAATGGGGCAACCAGCGTAATGAGGAGAACCCGCATCAtgt                                                                                             >  1:111800/1‑149 (MQ=255)                                                   gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                             <  2:484549/149‑1 (MQ=255)                                                   gAATAAACACCGGCGCGAAAAGCTCATCTTTAGCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacct                                                                                             <  1:826858/149‑4 (MQ=1)                                                      tACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACACGcatcacctgca                                                                                          >  2:428297/1‑143 (MQ=1)                                                      tACACACCGGAGCGAAAAGCTCATCTTTATCACTGAGACCACCGTGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGAGGCGCGGGGAGGATGATGAATGGGGCAACCAGCGTAAGTAGGATAACCCGcatcacctgca                                                                                          >  1:639914/1‑143 (MQ=1)                                                       acacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAt                                                                                         <  2:600848/149‑1 (MQ=255)                                                         acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACAGGGGACTGCCGCAGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGACAACCCGCATCATGTGTATcc                                                                                       >  1:456945/1‑148 (MQ=255)                                                             cGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCATCGAATACGGCAATACGACGTCGAGCGATGATTAAGATGAATGTGGCAAACAGCTTATTGATGAGAACCCGCATCATGTGTAtaaaacg                                                                                   >  1:805052/1‑143 (MQ=255)                                                             cGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTAATCATCGAAGACGGCAATACGGCGTGGATCGCGGAAAATGTTGAATGGGGCAACCATCGTAATGAGGAGAACCCGCAACAAGTGTATCAGCCg                                                                                   >  1:271002/1‑149 (MQ=17)                                                              ggCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGt                                                                                  >  1:26376/1‑149 (MQ=255)                                                               gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCCGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTa                                                                                 <  1:1213746/149‑1 (MQ=255)                                                                    aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGTGAACCCGCATCATGTGTATCAGCCGTACCGTc                                                                            <  2:567746/149‑1 (MQ=255)                                                                      aaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTg                                                                          <  1:626049/149‑1 (MQ=255)                                                                          tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCTCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                      <  1:996954/149‑1 (MQ=255)                                                                          tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTGAATGAGGAGACCCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                      >  1:775408/1‑149 (MQ=255)                                                                             tCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgc                                                                   <  1:403236/149‑1 (MQ=255)                                                                                        gCGACCACCGGGTACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCGTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTAGCGGCTGCCAGGGCAGCAGCatgat                                                        <  1:75934/149‑1 (MQ=255)                                                                                           accaccGGGGACTGCCGCTGGCGGTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCATCCAGCTTTATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                     <  2:1132898/149‑1 (MQ=255)                                                                                           accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGGGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                     >  1:619002/1‑149 (MQ=255)                                                                                           accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                     >  1:198995/1‑149 (MQ=255)                                                                                                   ggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAgg                                             <  2:28793/149‑1 (MQ=255)                                                                                                    gACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                            <  2:95570/149‑1 (MQ=255)                                                                                                      cTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCt                                          <  2:1113331/149‑1 (MQ=255)                                                                                                       tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCTAGGCCAGCAGCATGATGAGGATTCAGGGCTg                                         >  1:1229645/1‑149 (MQ=255)                                                                                                         ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                       >  2:717426/1‑149 (MQ=255)                                                                                                              ggCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACtata                                  >  1:1230509/1‑149 (MQ=255)                                                                                                                      tCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCAc                          >  2:34440/1‑149 (MQ=255)                                                                                                                       cAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCAcc                         >  1:458552/1‑149 (MQ=255)                                                                                                                       cAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCAcc                         >  2:631367/1‑149 (MQ=255)                                                                                                                       cAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAGGATGAGGTATCAGGGCTGTACTATAACCGTCAcc                         >  1:494073/1‑149 (MQ=255)                                                                                                                        aGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCg                        >  2:1157226/1‑149 (MQ=255)                                                                                                                         gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCACCCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCTGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACAgg                       >  2:1075117/1‑149 (MQ=255)                                                                                                                         gCGAAGACGGCAATACGGCCTGGAGCCCGGAATATGTTGAATGGGGCACCCAGCTTAATGAGGAGACCCCGCATCATGTGTTTCAGCCCTACCGTCTGCCAGCGCAGCAGCATGATGAGGAATCAGGTCTGTACTAAAACCGGCACAgg                       >  2:546325/1‑149 (MQ=17)                                                                                                                           gAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGta                     <  1:362877/149‑1 (MQ=255)                                                                                                                                 ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAAGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCCGGGCTGTACTATAACCGTCACCGGTACTACGa               >  1:503539/1‑149 (MQ=255)                                                                                                                                   cAATACTGCGTGGATCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATc             <  2:497878/149‑1 (MQ=255)                                                                                                                                    aaTACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATcc            >  1:731494/1‑149 (MQ=255)                                                                                                                                              ggAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGAACTATAACCGTCACCGGTACTACGATCCGTTGCAgggg  <  1:458835/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGG  >  NC_000913_3_hsa_pgi/526649‑526938

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