Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_pgi 4,235,435 Δ1 bp intergenic (‑2202/‑499) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_pgi4,235,4330A.99.0% 433.4 / ‑0.4 103intergenic (‑2200/‑501)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (66/35);  total (66/37)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

TTTGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_hsa_pgi/4235285‑4235577
                                                                                                                                                    |                                                                                                                                                
   gACCCGTGGGGAGTGGAGCGGGGAAAGCAGCTGGCTAAGAAAATCGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATaatc                                                                                                                                               <  2:152932/149‑3 (MQ=255)
       agggcggAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCAtcg                                                                                                                                          <  2:1495690/144‑1 (MQ=255)
          ggggAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCAtcgtcg                                                                                                                                       >  1:875008/1‑149 (MQ=255)
          ggggAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCAtcgtcg                                                                                                                                       <  2:688227/149‑1 (MQ=255)
           gggAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGa                                                                                                                                      >  2:1232048/1‑149 (MQ=255)
            ggAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGat                                                                                                                                     <  2:587951/149‑1 (MQ=255)
            ggAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGat                                                                                                                                     >  2:708750/1‑149 (MQ=255)
                  gAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTGGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGCGTCCAATAATCATCGACGATATGTAg                                                                                                                               >  1:116791/1‑149 (MQ=255)
                  gAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTGGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCGCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAg                                                                                                                               <  2:1435272/149‑1 (MQ=255)
                   aGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAgg                                                                                                                              >  1:697669/1‑149 (MQ=255)
                    gCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATAATCGTCGATAAGTAGGc                                                                                                                             >  2:1133980/1‑149 (MQ=255)
                           aaGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATaa                                                                                                                      >  1:1108150/1‑149 (MQ=255)
                           aaGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATaa                                                                                                                      <  2:1188972/149‑1 (MQ=255)
                            agcagcTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                     <  2:570614/149‑1 (MQ=255)
                               agcTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCTATATGTAGGCCGGATAAGGCg                                                                                                                  >  2:488662/1‑149 (MQ=255)
                               agcTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                                                                  >  1:966571/1‑149 (MQ=255)
                                 cTGTCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                <  2:1383482/149‑1 (MQ=255)
                                       aaGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGCGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGCATAAGGCGTTCAcgcc                                                                                                          >  1:405198/1‑149 (MQ=255)
                                       aaGAAAATAGAGCCTGAGCTTGAGGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCACAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                          <  1:1367462/149‑1 (MQ=255)
                                        aGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                         <  1:262742/149‑1 (MQ=255)
                                         gAAAATAGCGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGAAGGCCGGATAAGGCGTTCAcgccgc                                                                                                        <  1:90034/149‑1 (MQ=255)
                                         gAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                        >  2:215601/1‑149 (MQ=255)
                                           aaaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                      >  2:841419/1‑149 (MQ=255)
                                            aaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                     >  2:649576/1‑149 (MQ=255)
                                               agagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                  >  1:630754/1‑149 (MQ=255)
                                                gagCCTGAGCTTTATTGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGTATGCGGCGTGAACGCCTTATCCGGc                                                                                                 <  1:913207/149‑1 (MQ=255)
                                                gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                 >  2:740472/1‑149 (MQ=255)
                                                     tGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCCGATAAGGCGTTCACGCCGCATCCGGCAAccc                                                                                            >  2:203586/1‑148 (MQ=255)
                                                       aGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                          >  1:50332/1‑149 (MQ=255)
                                                        gCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATg                                                                                         <  1:266959/149‑1 (MQ=255)
                                                         cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                  >  2:1319085/1‑139 (MQ=255)
                                                         cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                  <  1:1319085/139‑1 (MQ=255)
                                                          ttGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGcc                                                                                       >  2:142841/1‑149 (MQ=255)
                                                          ttGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGcc                                                                                       >  2:591311/1‑149 (MQ=255)
                                                          ttGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAACCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGcc                                                                                       <  2:115827/149‑1 (MQ=255)
                                                            gATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTg                                                                                     >  1:872303/1‑149 (MQ=255)
                                                             aTGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCTATAATCATCGTCGAAATGTAGGCCGGATAAGGCGTTCACGCCGCATCCCGCAACCGCTGCCTGa                                                                                    >  2:1260723/1‑149 (MQ=39)
                                                              tGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                   >  1:942601/1‑149 (MQ=255)
                                                                gCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGc                                                                                 <  1:637309/149‑1 (MQ=255)
                                                                      cTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgcg                                                                           <  2:875008/149‑1 (MQ=255)
                                                                      cTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGAGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgcg                                                                           <  1:1304449/149‑1 (MQ=255)
                                                                       tCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCCTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCgg                                                                          >  2:718497/1‑149 (MQ=255)
                                                                       tCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTACAACAACCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCCGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCAGATGCGAcgcgc                                                                          >  2:1480548/1‑148 (MQ=38)
                                                                           gTGCCCTCTCACGCCGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGCAGTTCACTCCGCATCCCGCAACCGATGCCAGATGCGACGCGCTcgc                                                                      >  1:378414/1‑149 (MQ=255)
                                                                           gTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                      <  1:364441/149‑1 (MQ=255)
                                                                             gACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGAGGTCGCGt                                                                    >  1:121592/1‑149 (MQ=255)
                                                                                ctctCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCCGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                 >  1:376634/1‑149 (MQ=255)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCACGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCCCAGCCGGCAGCCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  2:1383869/1‑149 (MQ=38)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACCCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  1:579712/1‑149 (MQ=255)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtga                                                                >  1:916345/1‑147 (MQ=255)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  2:1498748/1‑149 (MQ=255)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  2:1320838/1‑149 (MQ=255)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  2:1279335/1‑149 (MQ=255)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAACGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  1:309623/1‑149 (MQ=255)
                                                                                     ccgacgCTTCTACCACTGGGCTCATCAAATTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAg                                                            <  2:1335824/148‑1 (MQ=37)
                                                                                     acgacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                            >  1:760056/1‑149 (MQ=255)
                                                                                       gacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                          <  2:942601/149‑1 (MQ=255)
                                                                                        acgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc                                                         >  1:889850/1‑149 (MQ=255)
                                                                                           cTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                      <  1:196044/149‑1 (MQ=255)
                                                                                            ttCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                     >  2:398667/1‑149 (MQ=255)
                                                                                              cTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAgg                                                   >  1:357582/1‑149 (MQ=255)
                                                                                              cTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAgg                                                   >  2:1374126/1‑149 (MQ=255)
                                                                                               tACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                  >  1:1449058/1‑149 (MQ=255)
                                                                                               aaCCAATCGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                  <  2:1564523/148‑1 (MQ=255)
                                                                                                aCCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCGAAAGGTc                                                 >  2:213637/1‑149 (MQ=255)
                                                                                                  cAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                               >  2:718461/1‑149 (MQ=255)
                                                                                                  cAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                               >  1:1267651/1‑149 (MQ=255)
                                                                                                  cAATGGGCTCATAAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCAGCATCCGCCAACCCATGCCTGATGCGACGCTGTCGCGTCTAATCAGGCCTACAGGACGa                                               >  1:379137/1‑149 (MQ=25)
                                                                                                   aaTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGAt                                              <  2:630754/149‑1 (MQ=255)
                                                                                                    aTGGGCTCATCAACTTCATCAAGCAGCTGCGCGAGGCCAGATTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTGTTTTCAGGCCTACAGGTCGATg                                             <  2:1090364/149‑1 (MQ=37)
                                                                                                     tGGGCTCATCAACTTCCTCTAGCAGCAGCGCGTGGCCAGAGTCCATTAATCCTCGGCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGGCGATGc                                            <  2:757103/149‑1 (MQ=21)
                                                                                                     tGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                            <  1:1232048/149‑1 (MQ=255)
                                                                                                     tGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                            <  2:1125982/149‑1 (MQ=255)
                                                                                                         cTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                        <  2:405198/149‑1 (MQ=255)
                                                                                                          tcatcaCCTTCATCTAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGGCGATGCCGata                                       <  2:164569/149‑1 (MQ=25)
                                                                                                            atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCGGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                     >  2:1125074/1‑149 (MQ=255)
                                                                                                            atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACCCGATGCCTGATGCGACGCGGGCGCGTCTTATCATGCCGACAGGTCGATGCCGagatt                                     >  1:112548/1‑148 (MQ=25)
                                                                                                               aaCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCGTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGtcc                                  >  1:1059372/1‑147 (MQ=255)
                                                                                                                 cTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                >  1:781322/1‑149 (MQ=255)
                                                                                                                  ttCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATTAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCGTAGCCGGCCTACAGGTCGATGCCGATATGTACATc                               >  2:1428774/1‑149 (MQ=38)
                                                                                                                  ttCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                               <  2:1221124/149‑1 (MQ=255)
                                                                                                                  ttCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCAGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                               >  2:286302/1‑149 (MQ=255)
                                                                                                                     atcaAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTa                            >  1:106286/1‑149 (MQ=255)
                                                                                                                       caAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAtt                          >  1:436544/1‑149 (MQ=255)
                                                                                                                        aaGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                         >  1:358611/1‑149 (MQ=255)
                                                                                                                         agcagcagcGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                        <  1:649576/149‑1 (MQ=255)
                                                                                                                          gcagcagcGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                       >  1:1244414/1‑149 (MQ=255)
                                                                                                                           cagcagcGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCTACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGc                      <  1:213637/149‑1 (MQ=255)
                                                                                                                            agcagcGCGAGGCCAGAGTCCAATAATCAGCGTGGATATGTAGGCCGGAGAAGGCGTTCACGCCGCATCCGGCACCCGGTGCCTGATGCGACGCGGTCGCGGCGTATCAGGCGTACAGGTCGATGCCGATATGTAGATCGTATTGGGCa                     >  2:866293/1‑149 (MQ=25)
                                                                                                                             gcagcGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                    >  1:35260/1‑149 (MQ=255)
                                                                                                                                 cgcgAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa                <  2:872303/149‑1 (MQ=255)
                                                                                                                                  gcgAGGCCAGAGTCGAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCAGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCGGATTTGTACATCGTATTCGGCAATTAAt               <  2:12214/149‑1 (MQ=25)
                                                                                                                                  gcgAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt               <  2:239907/149‑1 (MQ=255)
                                                                                                                                   cgAGGCCAGAGTCCTATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa              <  1:1138753/149‑1 (MQ=255)
                                                                                                                                    gAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTATTATCAGGCCTACAGGTCGATGCCGATATGTAAATCGTATTCGGCAATTAAtta             >  1:522759/1‑147 (MQ=255)
                                                                                                                                       gCCAGAGGCCAATACTCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa          <  1:708750/149‑1 (MQ=255)
                                                                                                                                         cAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGc        >  1:84181/1‑149 (MQ=255)
                                                                                                                                         cAGAGTCCAATAATAATCGTCGATATGTAGGACGGATAAGGCGTTCACGCCGCATACGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGACTAAAGGTCGATGACTATATGTACATCGTATTCTGCAATTAATACTTAGc        >  1:428102/1‑149 (MQ=25)
                                                                                                                                            agTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAGCGTATTCGGCAATTAATACATAGCACg     >  1:1173578/1‑149 (MQ=255)
                                                                                                                                               ccAATAATCATCGTCGCTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTGTTATCAGGCCTACAGGTCGATGCCGATAGGTACATCGTATTCGGCAATTAATACATAGCAGgatt  >  2:959281/1‑149 (MQ=25)
                                                                                                                                                    |                                                                                                                                                
TTTGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_hsa_pgi/4235285‑4235577

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: