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)
        gTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCAtcgt                                                                                                                                        <  1:718648/149‑1 (MQ=255)
           gggAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGa                                                                                                                                     >  2:866170/1‑149 (MQ=255)
                  gAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAg                                                                                                                              <  2:485179/149‑1 (MQ=255)
                  gAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGAGCGAGGCCAGAGTCCAATAATCATCGTCGATgtgtgg                                                                                                                              >  2:374509/1‑147 (MQ=255)
                        ggAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGCGACCAGANTCCAATAATCCTCGTCGCTATGTNGCCCGGa                                                                                                                        >  2:562167/1‑149 (MQ=255)
                              cagcTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                  >  1:409984/1‑149 (MQ=255)
                               agcTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                                                                 <  2:299145/149‑1 (MQ=255)
                                gcTGGCTAAGAAAATAGACCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGt                                                                                                                <  1:350993/149‑1 (MQ=255)
                                 cTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                               <  2:380391/149‑1 (MQ=255)
                                    gCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCGAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                                                            >  1:691973/1‑149 (MQ=255)
                                      tAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                          >  2:834314/1‑149 (MQ=255)
                                         gAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                       >  2:1051327/1‑149 (MQ=255)
                                             aTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                   >  2:691846/1‑149 (MQ=255)
                                                   ccTGACCTTGATGGAACAGCTCAAGTAATCTCTCACGACCCTTCTACCAATGGGATCATCAACTACATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                             <  1:301450/149‑1 (MQ=37)
                                                         ctcgacggCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGAGATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                       <  1:896976/143‑1 (MQ=255)
                                                         cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                       >  2:142203/1‑149 (MQ=255)
                                                            gATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTg                                                                                    >  1:818228/1‑149 (MQ=255)
                                                             aTGGCACTCCTCAAGTGACCTCCCACGACGCTTCTACCAATGGGCTCATCCACTACATCAATTGTGTGCGCGAGGCCAGAGTCCAATAATCCTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                   <  1:549707/149‑1 (MQ=37)
                                                                    tGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                                                                            >  2:683841/1‑149 (MQ=255)
                                                                       tCAAGTGACATCCCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCgg                                                                         <  1:408222/149‑1 (MQ=255)
                                                                         aaGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                       >  2:319292/1‑149 (MQ=255)
                                                                           gTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                     >  2:107734/1‑149 (MQ=255)
                                                                                 tctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                               >  1:948325/1‑149 (MQ=255)
                                                                                   tcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                             >  1:971187/1‑149 (MQ=255)
                                                                                     acgacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                           >  1:729297/1‑149 (MQ=255)
                                                                                     acgacgCTTCTACCAATCGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                           <  1:192331/149‑1 (MQ=255)
                                                                                       gacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                         <  1:20825/149‑1 (MQ=255)
                                                                                            ttCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                    >  1:63012/1‑149 (MQ=255)
                                                                                               tACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCAGTCGCGTCTTATCAGGCCTACAGGt                                                 >  1:255375/1‑149 (MQ=255)
                                                                                               tACCAATGGGCTCATCAACTTAATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                 <  1:869656/149‑1 (MQ=255)
                                                                                                aCCAAGCGGCTCATCACCTTCACCAAGCAGCAGCGCGACGCCAGACTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                <  1:374509/149‑1 (MQ=25)
                                                                                                 ccAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                  <  2:1127946/114‑1 (MQ=255)
                                                                                                 ccAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                  >  1:1127946/1‑114 (MQ=255)
                                                                                                 ccAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCg                                               >  2:227093/1‑149 (MQ=255)
                                                                                                  cAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                              >  1:654865/1‑149 (MQ=255)
                                                                                                  cAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                              <  2:302228/149‑1 (MQ=255)
                                                                                                   aaTGAGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGAt                                             >  2:104873/1‑149 (MQ=255)
                                                                                                    aTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                   >  2:799815/1‑142 (MQ=255)
                                                                                                    aTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                   <  1:799815/142‑1 (MQ=255)
                                                                                                    aTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATg                                            >  2:802555/1‑149 (MQ=255)
                                                                                                       ggCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCg                                         <  1:488839/149‑1 (MQ=255)
                                                                                                       ggCTAATAAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCg                                         <  2:432826/149‑1 (MQ=255)
                                                                                                        gCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                        <  2:818228/149‑1 (MQ=255)
                                                                                                        gCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCAACCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                        <  1:38247/149‑1 (MQ=255)
                                                                                                           catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                     <  2:691973/149‑1 (MQ=255)
                                                                                                           cagtcaCTTCATCAAGTAGCTGCGCGAGGCCAGAGTCCAATAATCATCGTCCATATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                     <  1:562167/144‑1 (MQ=21)
                                                                                                             tcaACTTCATCAAGCAGCAGCGAGTGACCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACTCCACATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                   >  2:865550/1‑149 (MQ=37)
                                                                                                             tcaACTTCAGCAAGCAGCAGTGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATCAGGCGTTCGCGCCACATCCGGCAACCGATGCCGGATGCGACGCGGTCGCGTCTTATCAGGCGTACAGGTCGATGCCGATATGt                                   <  1:928626/149‑1 (MQ=255)
                                                                                                                 cTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                  <  2:325400/66‑1 (MQ=255)
                                                                                                                 cTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                  >  1:325400/1‑66 (MQ=255)
                                                                                                                  ctcagctagcagcagcGCGAGGCCAGAGTCCAACAATCATCGTCGCTATGTAGGCCGGATAAGGCGTTCCCGCCGCACCCGGCATCCGCTGCCTGCTGCGACGCGGTCGTGTCCTATCAGGCCTACAGGTCGATGCCGATATGTACATc                              <  2:271337/142‑1 (MQ=25)
                                                                                                                       caAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAtt                         <  1:866170/149‑1 (MQ=255)
                                                                                                                        aaGCAGCAGGGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                        <  2:1130662/149‑1 (MQ=255)
                                                                                                                        aaGCAGCAGCGCGAGGCCAGAGGCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                        <  2:706298/149‑1 (MQ=255)
                                                                                                                             gaagcGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGCCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                   <  2:319704/147‑1 (MQ=255)
                                                                                                                                gcgcGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                >  1:825302/1‑149 (MQ=255)
                                                                                                                                gcgcGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGACGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGCCGCATGTGATTAGGTTTACAGGTCGATGCCGATAGATACAACATATTCGGCAAtaa                >  1:317973/1‑147 (MQ=25)
                                                                                                                                 cgcgAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa               <  2:918106/149‑1 (MQ=255)
                                                                                                                                  gcgAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt              >  2:397569/1‑149 (MQ=255)
                                                                                                                                   cgAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa             <  1:319292/149‑1 (MQ=255)
                                                                                                                                    gAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAc            >  2:81649/1‑149 (MQ=255)
                                                                                                                                        ccAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAg        <  1:771996/149‑1 (MQ=255)
                                                                                                                                         cAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGc       >  2:864982/1‑149 (MQ=255)
                                                                                                                                              tCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat  <  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: