Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_pgi 4,235,433 Δ1 bp 100% intergenic (‑2200/‑501) 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.93.0% 4047.4 / 264.7 1083intergenic (‑2200/‑501)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base . (482/525);  total (520/563)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 1136 total aligned reads displayed.

TTTGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_hsa_pgi/4235285‑4235577
                                                                                                                                                    |                                                                                                                                                
tttGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATTAAGCAGCAGAGCGAGGCCAGAAGCCAATa                                                                                                                                                  >  1:49805/1‑149 (MQ=255)
   gACCAGTTGGGAGGGGAGCTGGGAAAGCAGCAGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACGTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATaatc                                                                                                                                               <  2:237972/149‑3 (MQ=255)
   gACCAGTGGGGAGTGGAGCTTGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGGGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCGCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATactc                                                                                                                                               <  2:396585/149‑4 (MQ=255)
      cAGGGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAatcatc                                                                                                                                           <  1:633623/149‑1 (MQ=255)
         tGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCAtcgtc                                                                                                                                        <  1:922702/149‑1 (MQ=255)
           gggAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGa                                                                                                                                      <  1:173281/149‑1 (MQ=255)
             gAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGata                                                                                                                                    >  1:214588/1‑149 (MQ=255)
                   aGCTGGGAAAGCAGCTGGTTAAGAAAATAGAGCCTGAGCTTGATGGCCGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGg                                                                                                                              >  1:805670‑M2/1‑148 (MQ=255)
                   aGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGg                                                                                                                              <  2:76359‑M2/149‑2 (MQ=255)
                     gtgGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAACCAGCAGCGCGAGGCCAGAGTCCAATAACCATCGTCGATATGTAGgcc                                                                                                                            <  2:517459‑M2/148‑4 (MQ=255)
                      tGGAAAAGCAGCTGGTTAAAAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGgccg                                                                                                                           >  2:740404‑M2/1‑145 (MQ=255)
                        ggAAAGCAGCTGGCTAAGAAAATAGAGCAGGAGCGGGGGGGGGGTGGTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGgccgga                                                                                                                         >  1:437260‑M2/1‑143 (MQ=255)
                         gAAAGCAGCTGGCGAAGAAAATAGAGCCTGAGTTTGATGGCCGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAAAATCGGCTATATGTAGgccggat                                                                                                                        >  2:319488‑M2/1‑142 (MQ=255)
                          aaaGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCCACTTCCTCATGCAGCAGCGCGAGTCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATa                                                                                                                       <  1:534954/149‑1 (MQ=255)
                           aaGGAGCTGGCTAAGAAAATTGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATaa                                                                                                                      <  2:912084/149‑1 (MQ=255)
                           aaGCAGCTGGCTAAGAAAATAGGGCATGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATaa                                                                                                                      <  1:429464/149‑1 (MQ=255)
                           aaGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATaa                                                                                                                      <  2:539949/149‑1 (MQ=255)
                             gcagcTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAgg                                                                                                                    <  1:150139/149‑1 (MQ=255)
                               agcTCGCTAAGAAAATAGAGCCTGAGTTTGATTGTATTGCTTAAGTTACCTCTTACGTCGCTTTTACCAATGGGCTCAGCAACTTCATCAAGCAGCAGCGCGAGGCCATGGTCCAATAATCCTCGTCTATATGTAGGCCGGATAAGGCg                                                                                                                  <  2:34432/149‑1 (MQ=255)
                                    gCTAAGAAAATAGTGCCTGAGCTTGATGTCATTGCTCAAGTTACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                                                             <  1:381430/149‑1 (MQ=255)
                                      tAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAATTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                           <  1:549103/149‑1 (MQ=255)
                                      tAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCCTCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCTATATGTAAGCCGGATAAGGCGTTCAcgc                                                                                                           >  1:400523/1‑149 (MQ=255)
                                      tAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGTATAAGGCGTTCAcgc                                                                                                           >  2:44499/1‑149 (MQ=255)
                                      tAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                           >  1:103288/1‑149 (MQ=255)
                                       aaGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                          <  1:969123/149‑1 (MQ=255)
                                         gAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGATCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTCGCCCGGAAAAGGCGTTCAcgcccc                                                                                                        >  2:60231/1‑147 (MQ=255)
                                            aaTAGAGCCTGAGCCTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                     >  2:198461/1‑149 (MQ=255)
                                              tatagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                   <  1:661103/146‑1 (MQ=255)
                                              tAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACCCCGCATccc                                                                                                   >  2:181492/1‑148 (MQ=255)
                                                gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGAGATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                 <  1:646912/149‑1 (MQ=255)
                                                  gCCGGAGGTTGTTGGCAGGGCCCAGTGGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                               <  2:83311/149‑1 (MQ=38)
                                                   ccTGAGCTTGATGGCAGTGCTCAAGTTACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                              <  1:311924/149‑1 (MQ=255)
                                                   ccTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATaa                                                                                                                      >  2:803988/1‑125 (MQ=255)
                                                   ccTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                              <  1:521462/149‑1 (MQ=255)
                                                    cTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGTGCTTATCAATTTCATCAAGCAGCAGCTCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                                                                             <  1:27585/149‑1 (MQ=39)
                                                     tGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTGCATCAAGCAGCGGCGCCAGGCCAGAGTCCAATAATCATCGTCTATATGTAGGCCGGATAAGGCCTTTACGCCCCATCCGGCAACCg                                                                                            >  2:1625/1‑149 (MQ=255)
                                                      gAGCTTTATGTCAGTTCTTAAGTGACCTCTCACGACGCTTCTACCAATGGGATCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGa                                                                                           <  1:416899/149‑1 (MQ=38)
                                                         cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                        >  2:118057/1‑149 (MQ=255)
                                                         cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCAACTACTGAACCGCGATACGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                        >  1:784580/1‑149 (MQ=38)
                                                          ttGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGcc                                                                                       >  2:254433/1‑149 (MQ=255)
                                                           tGAGGGCAGTGCTCAAGTGACCCCACACGACGCTTCTACCTTTGGGCTCATCAACTTCATCAAGCAGGAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                      <  2:790325/149‑1 (MQ=38)
                                                             aTGGCAGTTCTCAAGTGCCCTCTCACGACGCTTATACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                    <  1:93865/149‑1 (MQ=255)
                                                              tGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                   >  1:677198/1‑149 (MQ=255)
                                                               ggCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGAACCGCGAGGCCAGAGGCGGATAATCGTCCTCCTTACGGAACCCGGATAAGGCATTCACCCCGCATCCGGCAACCGATACCCGAtc                                                                                  >  2:964995/1‑148 (MQ=255)
                                                                gCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGc                                                                                 <  1:729108/149‑1 (MQ=255)
                                                                gCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCAGCAACCGATTCCTGATGc                                                                                 >  1:791866/1‑149 (MQ=255)
                                                                 cAGTGCTTAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCCCACCTCACCGGCGTTCTGCTTCTGCCATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCg                                                                                <  1:740659/149‑1 (MQ=255)
                                                                  aGTGCTCAAGTGACCTCTCACGACGTTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGa                                                                               <  1:71657/149‑1 (MQ=255)
                                                                    tGCTCAAGTGACCTCTCACGACGCGTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                                                                             >  1:901405/1‑149 (MQ=255)
                                                                         aaGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCATAAAGCGAcgcgagc                                                                        >  2:356741/1‑146 (MQ=39)
                                                                             gACCACTCACGACGCTACAACCATAGAGCCCATCAATTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                                                                    <  1:766923/149‑1 (MQ=25)
                                                                              aCCTCTCAGGGCGCTTCTACCTATGGGCTTATCAACTTCTTCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATAGGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                                   <  1:198066/149‑1 (MQ=38)
                                                                              aCCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                                   >  2:682449/1‑149 (MQ=255)
                                                                               ccTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                                                                  >  1:872871/1‑149 (MQ=255)
                                                                                ctctCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                 >  2:776466/1‑149 (MQ=255)
                                                                                ctctCACGACGCTTCTACCAATGGGCTCATCAACTTCACCAAGCAGCAACGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGATCACGCCGCATCCGGCAACCGATGCCTGACTCTCCGCGGACGCATCtc                                                                 >  2:880728/1‑148 (MQ=38)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  1:480793/1‑149 (MQ=255)
                                                                                  ctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGgccg                                                                                                                           <  2:21633‑M2/89‑5 (MQ=255)
                                                                                    cactgcgCTGCTACCAATTGGCTTATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                             <  2:94451/144‑1 (MQ=37)
                                                                                    cACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                             <  1:1000162/149‑1 (MQ=255)
                                                                                     acgacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGACAAGCCTTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                            <  1:697181/149‑1 (MQ=255)
                                                                                     acgacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCACGGGAAGCAAGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTcgcatcat                                                            <  1:986541/149‑8 (MQ=21)
                                                                                       gacgCTTCTACCAATGGGCTCATCAACCTCATCAACCAGCAGCGCGAGCCCAGAGTCACACAATCATCGTCGATATGTAGGCCGGATAAGGCTTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAAAAGGc                                                          <  1:457628/149‑1 (MQ=38)
                                                                                           cTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGgcc                                                                                                                            >  1:99738‑M2/1‑76 (MQ=255)
                                                                                           cTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAAGCGGTCGCGTCTTATCAGGCCtcc                                                      >  1:451539/1‑147 (MQ=255)
                                                                                            ttCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                     <  1:757305/149‑1 (MQ=255)
                                                                                              cTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTCCAgg                                                   >  2:669937/1‑149 (MQ=255)
                                                                                                  cAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATg                                                                                         >  2:570265/1‑107 (MQ=255)
                                                                                                  cAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                               >  2:137684/1‑149 (MQ=255)
                                                                                                  cAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAAAGGTCGa                                               >  1:257673/1‑149 (MQ=255)
                                                                                                     tGGGCTAATCTACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATTGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                            <  1:333803/149‑1 (MQ=255)
                                                                                                      gggCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGTCAACCGATGACTGATACGACGCGGACGCGTCTTATCAGGACTACCGGTCGATGcc                                           >  1:584443/1‑149 (MQ=37)
                                                                                                          tcatcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                       >  2:703161/1‑149 (MQ=255)
                                                                                                           catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                            >  2:71324/1‑111 (MQ=255)
                                                                                                           catcaACTGCATTACGCAGCTGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                      <  2:280385/149‑1 (MQ=25)
                                                                                                            atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                                                                              >  2:980534/1‑76 (MQ=255)
                                                                                                            atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                     >  2:591479/1‑149 (MQ=255)
                                                                                                             tcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTTATTCGACGCGGTCGCGTCTTATAAGGCCTACAGATCGATGCCGCTATGt                                    >  2:629991/1‑149 (MQ=37)
                                                                                                             tcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                    <  2:857205/101‑1 (MQ=255)
                                                                                                                aCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                             >  2:553668/1‑121 (MQ=255)
                                                                                                                aCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGAGGTCGCGTCTTATCAAGCCTACAGGTCGATACCGATATGTACa                                 >  1:23745/1‑149 (MQ=255)
                                                                                                                 cTTCCTCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCTGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                <  1:139953/149‑1 (MQ=255)
                                                                                                                 cTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                >  2:294797/1‑149 (MQ=255)
                                                                                                                 cTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                <  2:374969/149‑1 (MQ=255)
                                                                                                                  ttCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATAAGTACATc                               >  1:119741/1‑149 (MQ=255)
                                                                                                                  ttCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAACAGGCCTACAGGTCGATGCCGATATGTACATc                               >  1:363358/1‑149 (MQ=255)
                                                                                                                  ttCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                               <  1:28545/149‑1 (MQ=255)
                                                                                                                        aaGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGGCTTATCAGGCCTACAGGACGATGCCGATATGTACATCGTATTc                         >  1:614014/1‑149 (MQ=255)
                                                                                                                        aaGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATACGACACGGTCGCATCTTATCAGGCCTACAGGTCGATGACGATATGTACATCGTATTc                         >  1:622189/1‑149 (MQ=38)
                                                                                                                        aaGCAGCAGCGCGACGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                         <  2:650256/149‑1 (MQ=255)
                                                                                                                          gcagcagcGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGGCTACAGGTCGATGCCGATATGTACATCGTATTCgg                       <  2:423307/149‑1 (MQ=255)
                                                                                                                          gcagcagcGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                       <  1:950854/149‑1 (MQ=255)
                                                                                                                            agcagcGCGAGGCAAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCTTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                     <  2:496572/149‑1 (MQ=255)
                                                                                                                                 cgcgcgGCCAGAGTCCAATAATCATCGTCAAGAGGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCGGATGCGACGCGGTCGCGTCGTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa                <  1:824876/149‑1 (MQ=37)
                                                                                                                                  gcgAGGTCAGAGTCCAATAACAATAGTCGATATGTAGGCCGGGTCAGGCGTTCGCGCCGCAGCCTGCGACCGTGGGCGGATGCGACGGGGTCGCTTCTTATGATGCCTGCAGGGCGATGCCGATATGTACATCGTATTCGGCAATTAAt               <  2:89222/149‑1 (MQ=255)
                                                                                                                                  gcgAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt               <  1:561763/149‑1 (MQ=255)
                                                                                                                                      ggCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt           <  1:85606/149‑1 (MQ=255)
                                                                                                                                       gacAGAGTCCAATAATCATCGTCGATAGTGAGGCCGGATATTGTGTTCTGTTCGCATGCGGAAACCGATGCCTGATGCGACGCTGTCGCGTCTGATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa          <  2:461440/147‑1 (MQ=25)
                                                                                                                                       accAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCGTGATGCGACGCGGCCGCTTCTTATCAGGCCTACAGGTCGTTGCCGATATGTACATCGTATACGGCAATTAATACATa          >  2:927158/2‑149 (MQ=37)
                                                                                                                                        ccAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAg         <  1:603626/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TTTGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_hsa_pgi/4235285‑4235577

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