Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_pgi4,235,4610AC26.3% 133.8 / 70.8 118intergenic (‑2228/‑473)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (44/43);  new base C (5/26);  total (49/69)
Fisher's exact test for biased strand distribution p-value = 1.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-02
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 122 total aligned reads displayed.

AGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTA  >  NC_000913_3_hsa_pgi/4235313‑4235610
                                                                                                                                                    |                                                                                                                                                                  
tgaatcTGTCTAAGAAAATAGTGCTTGAGTTTTATGGCAGTGATCAAGTGACCTCTCACGACGCTTCTACCAATGGGATAATCAACTTAATAAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑ATTCATCGTCGATATGTTGGCCGTATAAgg                                                                                                                                                                  <  2:1365275/145‑1 (MQ=255)
 gctgcTGGCTAAGAAAATAGAGCCTGAGGGTGATGGCAGTGCTAAAGTGACCTCTCACGACTCTTCTACCAATGGGCGCATCATCTTCATCAAGCAGCAGCGCGAGGCCAGAGGCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                                                                 <  1:1239538/147‑1 (MQ=255)
           aaGAAAATATAGACTGAGATTGATGTCAGTGCTCAAGTGACCTCTCAGGACTCTTCTACCAATTGGCTCATCAAATTCATCAATCAGCATCGCGAGTACAGAGTCCAAT‑AATCATCGTAGATATGTAGGCCTGATAAGGCGTTCAcgccg                                                                                                                                                       <  2:1263425/150‑1 (MQ=255)
               aaaTAGCGCCTGAGCGTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                   <  2:228687/150‑6 (MQ=255)
               aaaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                   <  1:1479693/150‑6 (MQ=255)
               aaaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                   <  1:1136482/150‑6 (MQ=255)
               aaaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                   <  2:733604/150‑6 (MQ=255)
               aaaGAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTGCTACCAATGGGCTCATCAACTTCATCACGCAGCCGCGCGAGGCCAGAGGCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                   <  1:120880/150‑6 (MQ=255)
                aaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCACCGCGAGGCCAGAGTCCAAT‑ACTCCTCGTCGATATGTACGCCGGATACGGCGTGCCCGCCGCCTcc                                                                                                                                                  >  2:1246182/1‑150 (MQ=255)
                aaTAGAGCCTGAGCTTGATGGCAGTGCTAAAGTGCCCTCTAACGACGATTCTAACAATGGGCTCATCAAATTCATCAAGCAGCAGCTCGAAGCCAGAGTCCAAT‑AATCATAGTCGTTATGTTAGCCGGATAAGGCGTACACGCAGCATca                                                                                                                                                  >  2:1431896/1‑149 (MQ=255)
                 aTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                 <  1:859863/150‑1 (MQ=255)
                 aTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                 <  1:75476/150‑1 (MQ=255)
                   agagCCTGCGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCACCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGCTAAGGCGTGCACGCCGCATCCGGc                                                                                                                                               >  1:1443493/1‑150 (MQ=255)
                   agagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                               <  1:1203050/150‑1 (MQ=255)
                    gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              <  2:1396932/150‑2 (MQ=255)
                    gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              <  1:1219836/150‑2 (MQ=255)
                    gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              <  1:1133119/150‑2 (MQ=255)
                    gagCCTGAGCTTGATGGCAGTGATCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              <  2:1233671/150‑2 (MQ=255)
                     cgCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                                             <  1:98684/149‑3 (MQ=255)
                     agCCTGAGCTTTATGTCAGTGCTCAAGTGACCGCTCACGAAGCTTATACGAATGGGCACATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATAAGTAGCCCGGATGCGGCGAGAACGCCTTATCCGGcct                                                                                                                                             <  1:1084448/150‑3 (MQ=255)
                     agCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATCGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                                             <  1:1270098/150‑3 (MQ=255)
                       ccTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAAAGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac                                                                                                                                           <  1:1013755/150‑5 (MQ=255)
                        ctcaGCTTGCGGGCAGTGCTCAAGCGACCTCTCACGACGCGTCTCCCCATGGGCTCAGCAACATCATCAACCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGAAATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑T‑‑‑‑‑‑Cg                                                                                                                             <  1:496295/147‑1 (MQ=37)
                        cTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCTCGACGCTTCTACCAATGGGCTCTTCAACGTCATCAAGCAGCAGCGCGAGGCCAGAGTACAA‑AAATCATCGTCGATATGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGC‑CA‑‑‑‑‑‑C‑‑‑‑‑‑Cg                                                                                                                             >  2:1314638/1‑150 (MQ=255)
                        cTGAGCTTGATGCCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATCGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctaca                                                                                                                                          <  1:911023/150‑6 (MQ=255)
                         tGAGGTTGATGGATGGGCTCAATTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGa                                                                                                                            <  2:724936/150‑1 (MQ=39)
                         tGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAA‑AAATCAGCGTCGATATGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGC‑aatctt                                                                                                                                        >  1:1017049/1‑146 (MQ=39)
                          gAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGCGGCCAGAGTCCAAT‑AATCATCGTCGTTATGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑A‑‑‑‑‑‑Cggg                                                                                                                           >  2:1137700/1‑148 (MQ=39)
                           aggttGATGTCTGTGCTAAAGTGACCTCTCACGACGCTTCTACAAATGGGCTCATCAACTTCATCAAGCAGCATCGCGAGGCCAGAGTCGAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGAGTTCACGCCGCATACGGC‑AA‑‑‑‑‑‑T‑‑‑‑‑‑CGGTg                                                                                                                          <  2:1104500/147‑1 (MQ=255)
                             cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGcc                                                                                                                        <  1:170939/150‑1 (MQ=255)
                             cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGcc                                                                                                                        >  1:424624/1‑150 (MQ=255)
                                   ggCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACATCTTAAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑TATCATAGTCGATATGTAGGACGGATAAGGCCTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑T‑‑‑‑‑‑CGCTGcacgcagc                                                                                                                  >  1:801819/1‑143 (MQ=37)
                                      aGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑T‑‑‑‑‑‑CGGTGCACGATGcccg                                                                                                               <  2:473170/150‑4 (MQ=37)
                                       gTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGAcg                                                                                                              >  1:385833/1‑150 (MQ=255)
                                       gTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTACATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCAACCTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCTGC‑AC‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGAcg                                                                                                              >  2:11442/1‑150 (MQ=38)
                                          cTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCTTCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATTCGACGCgg                                                                                                           <  1:442885/150‑1 (MQ=255)
                                          cTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGG‑AA‑‑‑‑‑‑T‑‑‑‑‑‑CGGTGCAGGATGccggatgc                                                                                                           >  2:810975/1‑143 (MQ=37)
                                          cTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATTCCTGATGCGACGCgg                                                                                                           >  1:1437815/1‑150 (MQ=255)
                                            cAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTc                                                                                                         >  1:689792/1‑150 (MQ=255)
                                             aaGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTcg                                                                                                        <  2:1323464/150‑1 (MQ=255)
                                               gTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AAGCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTcgcg                                                                                                      <  1:865702/150‑1 (MQ=255)
                                                 gACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGAAGCGCGAGGCCAGAGTCCAAT‑AATCATCCTCGATATGTAGGCCGGATACGGCGTTCACGCCGCATCCGGC‑TA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACTCGGTCGCGTc                                                                                                    <  2:1145931/150‑1 (MQ=39)
                                                   ccTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCACAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑A‑‑‑‑‑‑CGAAGCCTGATGCGACGCGGTCGCGTCtt                                                                                                  >  1:631055/1‑150 (MQ=255)
                                                    ctctCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                                                 >  2:924104/1‑150 (MQ=255)
                                                      ctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                                                               <  1:318194/150‑1 (MQ=255)
                                                      ctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                                                               <  1:924104/150‑1 (MQ=255)
                                                       tcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                                                              >  1:663867/1‑150 (MQ=255)
                                                        cACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGACTTATCAg                                                                                             >  1:537953/1‑150 (MQ=255)
                                                           gacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc                                                                                          <  1:581704/150‑1 (MQ=255)
                                                            acgATTCTACCACTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCTGAGGCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt                                                                                         <  1:1000443/150‑1 (MQ=25)
                                                               cTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGCGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AATCGGTGCA‑‑‑‑‑CGATGCCTGATGCGACGCTGGCGCGTCTTATCAg                                                                                             >  1:1207634/1‑150 (MQ=17)
                                                               cTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                                                      >  1:164551/1‑150 (MQ=255)
                                                                ttCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                                                     >  1:178822/1‑150 (MQ=255)
                                                                    aCCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCg                                                                                 >  2:667175/1‑150 (MQ=255)
                                                                    aCCAATGGGCTCATCAACTACATCACGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGACTGATAAGGCGTTCACGCCGCATCCGGC‑AATCGGTGCA‑‑‑‑‑CGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTa                                                                                        <  2:962690/150‑1 (MQ=2)
                                                                     ccAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGtgaac                                                                                                                                                           <  2:586943/88‑5 (MQ=255)
                                                                     ccAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGtgaac                                                                                                                                                           >  1:586943/1‑84 (MQ=255)
                                                                        aTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCGAT‑‑‑‑‑‑CGGTGGACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTAcgtc                                                                                    >  2:794702/1‑147 (MQ=9)
                                                                         tGGGCTCATCCACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGGCTTATCAGGCCTACAGGGCGATGcc                                                                            <  1:1300645/150‑1 (MQ=255)
                                                                               catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                      >  2:402327/1‑150 (MQ=255)
                                                                               catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                      >  2:185954/1‑150 (MQ=255)
                                                                               catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                      <  2:1340296/150‑1 (MQ=255)
                                                                                atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                                                     <  2:422730/150‑1 (MQ=255)
                                                                                atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                                                     >  1:430881/1‑150 (MQ=255)
                                                                                    aCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                                                 >  1:873590/1‑150 (MQ=255)
                                                                                     cTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                                                                <  1:1044750/150‑1 (MQ=255)
                                                                                      ttCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATACTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGAAGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCg                                                               <  2:511007/150‑1 (MQ=255)
                                                                                        catcaAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGCTGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTa                                                             >  2:1136373/1‑150 (MQ=255)
                                                                                            aaGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                                                         >  1:330751/1‑150 (MQ=255)
                                                                                             agcagcagcGCGAGGCCAGAGTCCAA‑GAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTAGCAGGCCGACAGGGCGAGGCCGATATGTACCTCGTATTCgg                                                        >  1:992601/1‑150 (MQ=25)
                                                                                              gcagcagcGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGc                                                       <  1:943324/150‑1 (MQ=255)
                                                                                              gcagcagcGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGc                                                       <  1:316916/150‑1 (MQ=255)
                                                                                               ctgcagcGCGAGGCCAGAGTCAAAT‑ACTCATCGTCGATATGTAGGCCGTATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTTCAGGTCGATGCCGATATGTACATCGTATTCGGCa                                                      <  1:1314638/148‑1 (MQ=37)
                                                                                                agcagcGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                                                     >  2:685925/1‑150 (MQ=255)
                                                                                                  cagcGCGAGGCCTGAGTCCAAT‑AATGATCGTCGATATGTAGGCCGGATAAGGCGTTCAGGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGAGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                                                   <  1:748051/150‑1 (MQ=37)
                                                                                                     cgcgAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGACGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGAGTATTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaaa                                                >  1:1126057/1‑149 (MQ=37)
                                                                                                       cgAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAc                                              <  1:298772/150‑1 (MQ=255)
                                                                                                       cgAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAc                                              <  1:567055/150‑1 (MQ=255)
                                                                                                            ccAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGTC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGACAGATGCGACGCCGTCGAGTCATCACAGGCCTACAGGTCGATGCCGCTATGTACCTCGTATTCGGCAATTAATAAATAGc                                         >  1:856236/1‑150 (MQ=37)
                                                                                                                  gccAAT‑AATCATCGTCGATCTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgatt                                   <  1:1246182/149‑1 (MQ=255)
                                                                                                                      aataaTCATCGTCGATATGTAGGCCGGATAAGGCGATCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATACGACGCGGTCGCGTCTTATCATGCCGACAGGTAGATTCCTATATGTACATCGTATTCGGCATATAAAACCAAGCCCgatttaa                                >  1:1373356/4‑147 (MQ=21)
                                                                                                                       ataaTCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattgatt                               <  1:650494/148‑1 (MQ=255)
                                                                                                                          atcatcGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTaaa                            >  1:42401/1‑150 (MQ=255)
                                                                                                                           tcatcGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAAt                           >  2:1293947/1‑150 (MQ=255)
                                                                                                                               cgtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATGCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAcc                       <  1:11442/150‑1 (MQ=255)
                                                                                                                                       tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAAc               >  1:1186035/1‑150 (MQ=255)
                                                                                                                                        gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGCTTGATTAAATAACCTTAATAACa              >  2:553405/1‑150 (MQ=255)
                                                                                                                                         gttgCCGGATGCGGCGTGAACGCCTTATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACaa             >  2:275128/4‑150 (MQ=38)
                                                                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATg           >  2:1208625/1‑150 (MQ=255)
                                                                                                                                            gCCGGATGCGGCGTGAACGCCTTATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGc          >  2:15707/1‑150 (MQ=38)
                                                                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTTATAACAATGc          >  2:282005/1‑150 (MQ=255)
                                                                                                                                             ccGGATGCGGCGTGAACGCCTTATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGcc         >  2:348706/1‑150 (MQ=38)
                                                                                                                                              cGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCg        <  2:178822/150‑1 (MQ=255)
                                                                                                                                              cGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCg        >  1:417955/1‑150 (MQ=255)
                                                                                                                                               ggATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGa       <  2:666599/150‑1 (MQ=255)
                                                                                                                                               ggATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGa       <  2:330751/150‑1 (MQ=255)
                                                                                                                                                    aGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTTATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa  <  2:247605/150‑1 (MQ=255)
                                                                                                                                                    aGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa  <  2:385833/150‑1 (MQ=255)
                                                                                                                                                    aGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa  >  1:15688/1‑150 (MQ=255)
                                                                                                                                                    aGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTAGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa  >  1:1116685/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                                  
AGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTA  >  NC_000913_3_hsa_pgi/4235313‑4235610

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑