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.98.4% 290.0 / ‑3.4 64intergenic (‑2200/‑501)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (31/32);  minor base C (0/1);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAT  >  NC_000913_3_hsa_pgi/4235285‑4235576
                                                                                                                                                    |                                                                                                                                               
tttGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATa                                                                                                                                                 >  1:380391/1‑149 (MQ=255)
   gACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGCCGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATaatc                                                                                                                                              <  1:996430/149‑3 (MQ=255)
     ccAGGGGGGAGTGGATGTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATaatcat                                                                                                                                            <  2:925519/149‑5 (MQ=255)
        gTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCAtcgt                                                                                                                                        <  1:718648/149‑1 (MQ=255)
           gggAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGa                                                                                                                                     >  2:866170/1‑149 (MQ=255)
                  gAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAg                                                                                                                              <  2:485179/149‑1 (MQ=255)
                  gAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGAGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATgtgtgg                                                                                                                              >  2:374509/1‑147 (MQ=255)
                        ggAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGCGACCAGANTCCAAT‑AATCCTCGTCGCTATGTNGCCCGGa                                                                                                                        >  2:562167/1‑149 (MQ=255)
                              cagcTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                  >  1:409984/1‑149 (MQ=255)
                               agcTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                                                                 <  2:299145/149‑1 (MQ=255)
                                gcTGGCTAAGAAAATAGACCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGt                                                                                                                <  1:350993/149‑1 (MQ=255)
                                 cTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                               <  2:380391/149‑1 (MQ=255)
                                    gCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCGAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                                                            >  1:691973/1‑149 (MQ=255)
                                      tAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                          >  2:834314/1‑149 (MQ=255)
                                         gAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                       >  2:1051327/1‑149 (MQ=255)
                                             aTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                   >  2:691846/1‑149 (MQ=255)
                                                   ccTGACCTTGATGGAACAGCTCAAGTAATCTCTCACGACCCTTCTACCAATGGGATCATCAACTACATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                             <  1:301450/149‑1 (MQ=37)
                                                         ctcgacggCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGAGATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                       <  1:896976/143‑1 (MQ=255)
                                                         cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                       >  2:142203/1‑149 (MQ=255)
                                                            gATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTg                                                                                    >  1:818228/1‑149 (MQ=255)
                                                             aTGGCACTCCTCAAGTGACCTCCCACGACGCTTCTACCAATGGGCTCATCCACTACATCAATTGTGTGCGCGAGGCCAGAGTCCAAT‑AATCCTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                   <  1:549707/149‑1 (MQ=37)
                                                                    tGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                                                                            >  2:683841/1‑149 (MQ=255)
                                                                       tCAAGTGACATCCCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCgg                                                                         <  1:408222/149‑1 (MQ=255)
                                                                         aaGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                       >  2:319292/1‑149 (MQ=255)
                                                                           gTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                     >  2:107734/1‑149 (MQ=255)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                               >  1:948325/1‑149 (MQ=255)
                                                                                   tcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                             >  1:971187/1‑149 (MQ=255)
                                                                                     acgacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                           >  1:729297/1‑149 (MQ=255)
                                                                                     acgacgCTTCTACCAATCGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                           <  1:192331/149‑1 (MQ=255)
                                                                                       gacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                         <  1:20825/149‑1 (MQ=255)
                                                                                            ttCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                    >  1:63012/1‑149 (MQ=255)
                                                                                               tACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCAGTCGCGTCTTATCAGGCCTACAGGt                                                 >  1:255375/1‑149 (MQ=255)
                                                                                               tACCAATGGGCTCATCAACTTAATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                 <  1:869656/149‑1 (MQ=255)
                                                                                                aCCAAGCGGCTCATCACCTTCACCAAGCAGCAGCGCGACGCCAGACTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                <  1:374509/149‑1 (MQ=25)
                                                                                                 ccAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                  <  2:1127946/114‑1 (MQ=255)
                                                                                                 ccAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                  >  1:1127946/1‑114 (MQ=255)
                                                                                                 ccAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCg                                               >  2:227093/1‑149 (MQ=255)
                                                                                                  cAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                              >  1:654865/1‑149 (MQ=255)
                                                                                                  cAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                              <  2:302228/149‑1 (MQ=255)
                                                                                                   aaTGAGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGAt                                             >  2:104873/1‑149 (MQ=255)
                                                                                                    aTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                   >  2:799815/1‑142 (MQ=255)
                                                                                                    aTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                   <  1:799815/142‑1 (MQ=255)
                                                                                                    aTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATg                                            >  2:802555/1‑149 (MQ=255)
                                                                                                       ggCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCg                                         <  1:488839/149‑1 (MQ=255)
                                                                                                       ggCTAATAAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCg                                         <  2:432826/149‑1 (MQ=255)
                                                                                                        gCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                        <  2:818228/149‑1 (MQ=255)
                                                                                                        gCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCAACCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                        <  1:38247/149‑1 (MQ=255)
                                                                                                           catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                     <  2:691973/149‑1 (MQ=255)
                                                                                                           cagtcaCTTCATCAAGTAGCTGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCCATATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                     <  1:562167/144‑1 (MQ=21)
                                                                                                             tcaACTTCATCAAGCAGCAGCGAGTGACCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACTCCACATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                   >  2:865550/1‑149 (MQ=37)
                                                                                                             tcaACTTCAGCAAGCAGCAGTGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATCAGGCGTTCGCGCCACATCCGGCAACCGATGCCGGATGCGACGCGGTCGCGTCTTATCAGGCGTACAGGTCGATGCCGATATGt                                   <  1:928626/149‑1 (MQ=255)
                                                                                                                 cTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                  <  2:325400/66‑1 (MQ=255)
                                                                                                                 cTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                  >  1:325400/1‑66 (MQ=255)
                                                                                                                  ctcagctagcagcagcGCGAGGCCAGAGTCCAA‑CAATCATCGTCGCTATGTAGGCCGGATAAGGCGTTCCCGCCGCACCCGGCATCCGCTGCCTGCTGCGACGCGGTCGTGTCCTATCAGGCCTACAGGTCGATGCCGATATGTACATc                              <  2:271337/142‑1 (MQ=25)
                                                                                                                       caAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAtt                         <  1:866170/149‑1 (MQ=255)
                                                                                                                        aaGCAGCAGGGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                        <  2:1130662/149‑1 (MQ=255)
                                                                                                                        aaGCAGCAGCGCGAGGCCAGAGGCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                        <  2:706298/149‑1 (MQ=255)
                                                                                                                             gaagcGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGCCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                   <  2:319704/147‑1 (MQ=255)
                                                                                                                                gcgcGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                >  1:825302/1‑149 (MQ=255)
                                                                                                                                gcgcGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGACGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGCCGCATGTGATTAGGTTTACAGGTCGATGCCGATAGATACAACATATTCGGCAAtaa                >  1:317973/1‑147 (MQ=25)
                                                                                                                                 cgcgAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa               <  2:918106/149‑1 (MQ=255)
                                                                                                                                  gcgAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt              >  2:397569/1‑149 (MQ=255)
                                                                                                                                   cgAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa             <  1:319292/149‑1 (MQ=255)
                                                                                                                                    gAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAc            >  2:81649/1‑149 (MQ=255)
                                                                                                                                        ccAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAg        <  1:771996/149‑1 (MQ=255)
                                                                                                                                         cAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGc       >  2:864982/1‑149 (MQ=255)
                                                                                                                                              tCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat  <  1:802555/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TTTGACCAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAT  >  NC_000913_3_hsa_pgi/4235285‑4235576

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