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.6% 3978.9 / 8.1 892intergenic (‑2200/‑501)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base . (449/433);  total (452/440)
Fisher's exact test for biased strand distribution p-value = 5.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 936 total aligned reads displayed.

TTTGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_hsa_pgi/4235285‑4235577
                                                                                                                                                    |                                                                                                                                                
tttGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATa                                                                                                                                                  <  1:595961/149‑1 (MQ=255)
tttGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATa                                                                                                                                                  >  1:188643/1‑149 (MQ=255)
tttGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATa                                                                                                                                                  >  2:746088/1‑149 (MQ=255)
tttGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAAGCCAATa                                                                                                                                                  >  2:338724/1‑149 (MQ=255)
  tGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATaat                                                                                                                                                <  2:932105/149‑2 (MQ=255)
   gACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCAGCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATaatc                                                                                                                                               <  2:365030/149‑3 (MQ=255)
      cAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑Aatcatc                                                                                                                                           <  1:741214/149‑1 (MQ=255)
           gggAGTGGAGCTGGGAAAGCGGCTGGCTAATAAAATAGAGCCTGCGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGa                                                                                                                                      <  2:909265/149‑1 (MQ=255)
            ggAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGat                                                                                                                                     >  2:827054/1‑149 (MQ=255)
             gAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGata                                                                                                                                    >  1:112278/1‑149 (MQ=255)
             gAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGata                                                                                                                                    >  1:569258/1‑149 (MQ=255)
                    gCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCCTGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGc                                                                                                                             <  1:242242/149‑1 (MQ=255)
                       gggAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCgg                                                                                                                          >  2:268434/1‑149 (MQ=255)
                        ggAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGACGCCAGATTCCAAT‑AATCATCGTCGATATGTAGGCCGGa                                                                                                                         <  2:962684/149‑1 (MQ=255)
                         gAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGCCCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCAAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGAt                                                                                                                        >  1:765102/1‑149 (MQ=255)
                         gAAAGCAGCCGGCTAAGAAAATAGAGCCTGAGCTTTATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGAt                                                                                                                        <  2:563080/149‑1 (MQ=255)
                             gcagggggCTAGGAACATAGGGGCGGGGCTTGATGGCACGCCTCAAGTTACCTCTCACGGAGCCTCTACCAGCGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAgg                                                                                                                    <  2:111363/143‑1 (MQ=255)
                               agcTGGCTAAGAAAATAGAGCCTGATCTTGATGGCAGTGCCCAAGTGACCTCTCACGACGCTTCTACCAATTGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                                                                  <  2:476095/149‑1 (MQ=255)
                               agcTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                                                                  <  2:284896/149‑1 (MQ=255)
                                gcTGGCTAAGACGATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTTACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGt                                                                                                                 <  1:710968/149‑1 (MQ=255)
                                gcTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGt                                                                                                                 >  2:810315/1‑149 (MQ=255)
                                   ggCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                                                                              <  2:237313/149‑1 (MQ=255)
                                    gCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                                                             >  1:958222/1‑149 (MQ=255)
                                     cTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcg                                                                                                            >  2:457355/1‑149 (MQ=255)
                                      tAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATaatcat                                                                                                                                             >  1:629086/1‑112 (MQ=255)
                                       aaGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGTCCATAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                          <  1:696647/149‑1 (MQ=255)
                                         gAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGatat                                                                                                                                   >  1:33973/1‑122 (MQ=255)
                                              tAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCCCGCCGCATCCg                                                                                                   >  1:53495/1‑149 (MQ=255)
                                               agagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                  >  1:418012/1‑149 (MQ=255)
                                               agagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCCCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                  <  2:591938/149‑1 (MQ=255)
                                                   ccTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                              >  2:21832/1‑149 (MQ=255)
                                                    cTGAGCTGGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                                                                             <  1:154461/149‑1 (MQ=255)
                                                             aTGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                    <  2:41628/149‑1 (MQ=255)
                                                              tGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                   <  1:618976/149‑1 (MQ=255)
                                                              tGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                   <  1:859567/149‑1 (MQ=255)
                                                               ggCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                                  <  2:548068/149‑1 (MQ=255)
                                                                gCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGc                                                                                 >  2:48787/1‑149 (MQ=255)
                                                                  aGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGAt                                                                                                                        >  1:277567/1‑108 (MQ=255)
                                                                      cTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgcg                                                                           >  2:946524/1‑149 (MQ=255)
                                                                      cTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAACCGCTGCCTGATGCGAcgcg                                                                           >  1:30076/1‑149 (MQ=255)
                                                                       tCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGc                                                                                                                             >  2:681109/1‑98 (MQ=255)
                                                                       tCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGc                                                                                                                             <  1:681109/98‑1 (MQ=255)
                                                                        cAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                         <  1:172664/149‑1 (MQ=255)
                                                                         aaGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCTTGATGCGACGCGGTc                                                                        <  1:492599/149‑1 (MQ=255)
                                                                         aaGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                        >  2:731653/1‑149 (MQ=255)
                                                                         aaGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                        <  2:644014/149‑1 (MQ=255)
                                                                         aaGTGACCTCTCACGACGCTCCTACCAATGTGCTCCTCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                        <  2:611459/149‑1 (MQ=255)
                                                                           gTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                      <  1:980024/149‑1 (MQ=255)
                                                                                ctctCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                    >  1:792720/1‑130 (MQ=255)
                                                                                ctctCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGCCGCGGTCGCGTCtt                                                                 >  2:309149/1‑149 (MQ=255)
                                                                                ctccCACGACGCTTCTCCCAATGGGCTCATCAACTTCATCACGCAGCAGCGCGAGGCCAGCGTCCAAT‑AATCCTCGTCGTTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGCCGATGCCTGATGCGCCGCGGCCCCGCCtt                                                                 >  2:524182/1‑149 (MQ=25)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  1:899970/1‑149 (MQ=255)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                >  1:947366/1‑149 (MQ=255)
                                                                                   tcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                              >  1:81354/1‑149 (MQ=255)
                                                                                   tcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                              >  1:401827/1‑133 (MQ=255)
                                                                                   tcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                              >  2:78799/1‑149 (MQ=255)
                                                                                    cACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                             <  1:996639/149‑1 (MQ=255)
                                                                                     acgacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                           >  2:14372/1‑102 (MQ=255)
                                                                                     acgacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                            >  1:779887/1‑149 (MQ=255)
                                                                                      cgacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg                                                           >  1:650952/1‑149 (MQ=255)
                                                                                       gacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGGCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                          >  1:462922/1‑101 (MQ=255)
                                                                                          gCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTa                                                       >  2:316035/1‑149 (MQ=255)
                                                                                          gCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAGGcatc                                                       <  2:712921/149‑4 (MQ=38)
                                                                                          gCTTCTACCAAGGGGCGCACCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTa                                                       <  1:12382/149‑1 (MQ=255)
                                                                                            ttCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCAtcgtc                                                                                                                                        >  2:846848/1‑66 (MQ=255)
                                                                                             tCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                    <  1:507046/149‑1 (MQ=255)
                                                                                             tCTACAAATGGGCTCATCAATTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                    <  1:476670/149‑1 (MQ=255)
                                                                                              cTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAGGcatcggtt                                                   <  2:43543/149‑8 (MQ=37)
                                                                                               tACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                  <  1:912303/149‑1 (MQ=255)
                                                                                               tACCAATGGGCCCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                  <  1:206801/149‑1 (MQ=255)
                                                                                                aCCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                 <  2:83488/133‑1 (MQ=255)
                                                                                                aCCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                 <  1:100478/149‑1 (MQ=255)
                                                                                                  cAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                               >  1:373174/1‑149 (MQ=255)
                                                                                                  cAATGGGCGCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                               <  2:502545/149‑1 (MQ=255)
                                                                                                        gCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATaa                                                                                                                      >  1:881841/1‑72 (MQ=255)
                                                                                                        gCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                         <  2:873728/149‑1 (MQ=255)
                                                                                                          tcatcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                   <  1:662026/73‑1 (MQ=255)
                                                                                                          tcatcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                   >  2:662026/1‑73 (MQ=255)
                                                                                                          tcatcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                              >  1:315575/1‑110 (MQ=255)
                                                                                                            atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                          >  1:824745/1‑128 (MQ=255)
                                                                                                               aaCTTCATCAAGCAGCGGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTAc                                  <  2:5770/149‑1 (MQ=255)
                                                                                                               aaCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                               >  1:286378/1‑120 (MQ=255)
                                                                                                               aaCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTAc                                  >  2:773393/1‑149 (MQ=255)
                                                                                                                aCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACa                                 >  1:721586/1‑149 (MQ=255)
                                                                                                                 cTTCATCAAGCCGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATACGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                <  2:629611/149‑1 (MQ=255)
                                                                                                                  ttCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATACCGATATGTACATc                               >  2:982005/1‑149 (MQ=255)
                                                                                                                      tcaAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAt                           >  1:644012/1‑149 (MQ=255)
                                                                                                                        aaGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                >  2:117318/1‑142 (MQ=255)
                                                                                                                        aaGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                         >  1:450443/1‑149 (MQ=255)
                                                                                                                         agcagcagcGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGCTACGTACGGCGTATTCg                        >  2:156667/1‑149 (MQ=25)
                                                                                                                             gcagcGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTa                                                       >  2:692380/1‑114 (MQ=255)
                                                                                                                             gcagcGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                    <  1:582947/149‑1 (MQ=255)
                                                                                                                                  gcgAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt               >  2:170677/1‑149 (MQ=255)
                                                                                                                                    gAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATACGGCAATTAATAc             >  1:841720/1‑149 (MQ=255)
                                                                                                                                         cAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                        >  1:928614/1‑117 (MQ=255)
                                                                                                                                          agagTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTCCATCGTCTTCGGCAATTAATACATAGCa       >  1:258370/1‑149 (MQ=255)
                                                                                                                                              tCCAAT‑AATCATGGCAGCTATGTCCGCGGGCTACGGCGTTCACGCCGCATCCGCCACCCGATGCCTGTTGCGACGCGGTCGCGTCTTATCAGGCCTACAAGTCGATGGCGATACGTGCATCGTATTAGGCAATTCATCCGTAGCaagct   >  2:194852/1‑145 (MQ=16)
                                                                                                                                              tCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat   <  2:393269/149‑1 (MQ=255)
                                                                                                                                              tCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGATCgat   >  2:254914/1‑149 (MQ=255)
                                                                                                                                              tCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTACCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat   >  2:892267/1‑149 (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 < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑