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.

AGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTA  >  NC_000913_3_hsa_pgi/4235313‑4235610
                                                                                                                                                    |                                                                                                                                                                  
tgaatcTGTCTAAGAAAATAGTGCTTGAGTTTTATGGCAGTGATCAAGTGACCTCTCACGACGCTTCTACCAATGGGATAATCAACTTAATAAAGCAGCAGCGCGAGGCCAGAGTCCAATATTCATCGTCGATATGTTGGCCGTATAAgg                                                                                                                                                                  <  2:1365275/145‑1 (MQ=255)
 gctgcTGGCTAAGAAAATAGAGCCTGAGGGTGATGGCAGTGCTAAAGTGACCTCTCACGACTCTTCTACCAATGGGCGCATCATCTTCATCAAGCAGCAGCGCGAGGCCAGAGGCCAATAATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                                                                 <  1:1239538/147‑1 (MQ=255)
           aaGAAAATATAGACTGAGATTGATGTCAGTGCTCAAGTGACCTCTCAGGACTCTTCTACCAATTGGCTCATCAAATTCATCAATCAGCATCGCGAGTACAGAGTCCAATAATCATCGTAGATATGTAGGCCTGATAAGGCGTTCAcgccg                                                                                                                                                       <  2:1263425/150‑1 (MQ=255)
               aaaTAGCGCCTGAGCGTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                   <  2:228687/150‑6 (MQ=255)
               aaaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                   <  1:1479693/150‑6 (MQ=255)
               aaaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                   <  1:1136482/150‑6 (MQ=255)
               aaaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                   <  2:733604/150‑6 (MQ=255)
               aaaGAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTGCTACCAATGGGCTCATCAACTTCATCACGCAGCCGCGCGAGGCCAGAGGCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                   <  1:120880/150‑6 (MQ=255)
                aaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCACCGCGAGGCCAGAGTCCAATACTCCTCGTCGATATGTACGCCGGATACGGCGTGCCCGCCGCCTcc                                                                                                                                                  >  2:1246182/1‑150 (MQ=255)
                aaTAGAGCCTGAGCTTGATGGCAGTGCTAAAGTGCCCTCTAACGACGATTCTAACAATGGGCTCATCAAATTCATCAAGCAGCAGCTCGAAGCCAGAGTCCAATAATCATAGTCGTTATGTTAGCCGGATAAGGCGTACACGCAGCATca                                                                                                                                                  >  2:1431896/1‑149 (MQ=255)
                 aTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                 <  1:859863/150‑1 (MQ=255)
                 aTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                 <  1:75476/150‑1 (MQ=255)
                   agagCCTGCGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCACCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGCTAAGGCGTGCACGCCGCATCCGGc                                                                                                                                               >  1:1443493/1‑150 (MQ=255)
                   agagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                               <  1:1203050/150‑1 (MQ=255)
                    gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              <  2:1396932/150‑2 (MQ=255)
                    gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              <  1:1219836/150‑2 (MQ=255)
                    gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              <  1:1133119/150‑2 (MQ=255)
                    gagCCTGAGCTTGATGGCAGTGATCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              <  2:1233671/150‑2 (MQ=255)
                     cgCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                                             <  1:98684/149‑3 (MQ=255)
                     agCCTGAGCTTTATGTCAGTGCTCAAGTGACCGCTCACGAAGCTTATACGAATGGGCACATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATAAGTAGCCCGGATGCGGCGAGAACGCCTTATCCGGcct                                                                                                                                             <  1:1084448/150‑3 (MQ=255)
                     agCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATCGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                                             <  1:1270098/150‑3 (MQ=255)
                       ccTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAAAGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac                                                                                                                                           <  1:1013755/150‑5 (MQ=255)
                        ctcaGCTTGCGGGCAGTGCTCAAGCGACCTCTCACGACGCGTCTCCCCATGGGCTCAGCAACATCATCAACCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGAAATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑T‑‑‑‑‑‑Cg                                                                                                                             <  1:496295/147‑1 (MQ=37)
                        cTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCTCGACGCTTCTACCAATGGGCTCTTCAACGTCATCAAGCAGCAGCGCGAGGCCAGAGTACAAAAATCATCGTCGATATGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCCA‑‑‑‑‑‑C‑‑‑‑‑‑Cg                                                                                                                             >  2:1314638/1‑150 (MQ=255)
                        cTGAGCTTGATGCCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATCGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctaca                                                                                                                                          <  1:911023/150‑6 (MQ=255)
                         tGAGGTTGATGGATGGGCTCAATTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGa                                                                                                                            <  2:724936/150‑1 (MQ=39)
                         tGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAAAATCAGCGTCGATATGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGCaatctt                                                                                                                                        >  1:1017049/1‑146 (MQ=39)
                          gAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGCGGCCAGAGTCCAATAATCATCGTCGTTATGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑A‑‑‑‑‑‑Cggg                                                                                                                           >  2:1137700/1‑148 (MQ=39)
                           aggttGATGTCTGTGCTAAAGTGACCTCTCACGACGCTTCTACAAATGGGCTCATCAACTTCATCAAGCAGCATCGCGAGGCCAGAGTCGAATAATCATCGTCGATATGTAGGCCGGATAAGGAGTTCACGCCGCATACGGCAA‑‑‑‑‑‑T‑‑‑‑‑‑CGGTg                                                                                                                          <  2:1104500/147‑1 (MQ=255)
                             cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGcc                                                                                                                        <  1:170939/150‑1 (MQ=255)
                             cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGcc                                                                                                                        >  1:424624/1‑150 (MQ=255)
                                   ggCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACATCTTAAAGCAGCAGCGCGAGGCCAGAGTCCAATTATCATAGTCGATATGTAGGACGGATAAGGCCTTCACGCCGCATCCGGCAA‑‑‑‑‑‑T‑‑‑‑‑‑CGCTGcacgcagc                                                                                                                  >  1:801819/1‑143 (MQ=37)
                                      aGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑T‑‑‑‑‑‑CGGTGCACGATGcccg                                                                                                               <  2:473170/150‑4 (MQ=37)
                                       gTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGAcg                                                                                                              >  1:385833/1‑150 (MQ=255)
                                       gTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTACATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCAACCTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCTGCAC‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGAcg                                                                                                              >  2:11442/1‑150 (MQ=38)
                                          cTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCTTCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATTCGACGCgg                                                                                                           <  1:442885/150‑1 (MQ=255)
                                          cTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGGAA‑‑‑‑‑‑T‑‑‑‑‑‑CGGTGCAGGATGccggatgc                                                                                                           >  2:810975/1‑143 (MQ=37)
                                          cTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATTCCTGATGCGACGCgg                                                                                                           >  1:1437815/1‑150 (MQ=255)
                                            cAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTc                                                                                                         >  1:689792/1‑150 (MQ=255)
                                             aaGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTcg                                                                                                        <  2:1323464/150‑1 (MQ=255)
                                               gTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAGCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTcgcg                                                                                                      <  1:865702/150‑1 (MQ=255)
                                                 gACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGAAGCGCGAGGCCAGAGTCCAATAATCATCCTCGATATGTAGGCCGGATACGGCGTTCACGCCGCATCCGGCTA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACTCGGTCGCGTc                                                                                                    <  2:1145931/150‑1 (MQ=39)
                                                   ccTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCACAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑A‑‑‑‑‑‑CGAAGCCTGATGCGACGCGGTCGCGTCtt                                                                                                  >  1:631055/1‑150 (MQ=255)
                                                    ctctCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                                                 >  2:924104/1‑150 (MQ=255)
                                                      ctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                                                               <  1:318194/150‑1 (MQ=255)
                                                      ctcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                                                               <  1:924104/150‑1 (MQ=255)
                                                       tcACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                                                              >  1:663867/1‑150 (MQ=255)
                                                        cACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGACTTATCAg                                                                                             >  1:537953/1‑150 (MQ=255)
                                                           gacgCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc                                                                                          <  1:581704/150‑1 (MQ=255)
                                                            acgATTCTACCACTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCTGAGGCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt                                                                                         <  1:1000443/150‑1 (MQ=25)
                                                               cTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGCGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCA‑‑‑‑‑CGATGCCTGATGCGACGCTGGCGCGTCTTATCAg                                                                                             >  1:1207634/1‑150 (MQ=17)
                                                               cTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                                                      >  1:164551/1‑150 (MQ=255)
                                                                ttCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                                                     >  1:178822/1‑150 (MQ=255)
                                                                    aCCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCg                                                                                 >  2:667175/1‑150 (MQ=255)
                                                                    aCCAATGGGCTCATCAACTACATCACGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGACTGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCA‑‑‑‑‑CGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTa                                                                                        <  2:962690/150‑1 (MQ=2)
                                                                     ccAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGtgaac                                                                                                                                                           <  2:586943/88‑5 (MQ=255)
                                                                     ccAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGtgaac                                                                                                                                                           >  1:586943/1‑84 (MQ=255)
                                                                        aTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCGAT‑‑‑‑‑‑CGGTGGACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTAcgtc                                                                                    >  2:794702/1‑147 (MQ=9)
                                                                         tGGGCTCATCCACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGGCTTATCAGGCCTACAGGGCGATGcc                                                                            <  1:1300645/150‑1 (MQ=255)
                                                                               catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                      >  2:402327/1‑150 (MQ=255)
                                                                               catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                      >  2:185954/1‑150 (MQ=255)
                                                                               catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                      <  2:1340296/150‑1 (MQ=255)
                                                                                atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                                                     <  2:422730/150‑1 (MQ=255)
                                                                                atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                                                     >  1:430881/1‑150 (MQ=255)
                                                                                    aCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                                                 >  1:873590/1‑150 (MQ=255)
                                                                                     cTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                                                                <  1:1044750/150‑1 (MQ=255)
                                                                                      ttCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATACTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGAAGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCg                                                               <  2:511007/150‑1 (MQ=255)
                                                                                        catcaAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGCTGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTa                                                             >  2:1136373/1‑150 (MQ=255)
                                                                                            aaGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                                                         >  1:330751/1‑150 (MQ=255)
                                                                                             agcagcagcGCGAGGCCAGAGTCCAAGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTAGCAGGCCGACAGGGCGAGGCCGATATGTACCTCGTATTCgg                                                        >  1:992601/1‑150 (MQ=25)
                                                                                              gcagcagcGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGc                                                       <  1:943324/150‑1 (MQ=255)
                                                                                              gcagcagcGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGc                                                       <  1:316916/150‑1 (MQ=255)
                                                                                               ctgcagcGCGAGGCCAGAGTCAAATACTCATCGTCGATATGTAGGCCGTATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTTCAGGTCGATGCCGATATGTACATCGTATTCGGCa                                                      <  1:1314638/148‑1 (MQ=37)
                                                                                                agcagcGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                                                     >  2:685925/1‑150 (MQ=255)
                                                                                                  cagcGCGAGGCCTGAGTCCAATAATGATCGTCGATATGTAGGCCGGATAAGGCGTTCAGGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGAGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                                                   <  1:748051/150‑1 (MQ=37)
                                                                                                     cgcgAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGAGTATTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaaa                                                >  1:1126057/1‑149 (MQ=37)
                                                                                                       cgAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAc                                              <  1:298772/150‑1 (MQ=255)
                                                                                                       cgAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAc                                              <  1:567055/150‑1 (MQ=255)
                                                                                                            ccAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGTCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGACAGATGCGACGCCGTCGAGTCATCACAGGCCTACAGGTCGATGCCGCTATGTACCTCGTATTCGGCAATTAATAAATAGc                                         >  1:856236/1‑150 (MQ=37)
                                                                                                                  gccAATAATCATCGTCGATCTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgatt                                   <  1:1246182/149‑1 (MQ=255)
                                                                                                                      aataaTCATCGTCGATATGTAGGCCGGATAAGGCGATCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATACGACGCGGTCGCGTCTTATCATGCCGACAGGTAGATTCCTATATGTACATCGTATTCGGCATATAAAACCAAGCCCgatttaa                                >  1:1373356/4‑147 (MQ=21)
                                                                                                                       ataaTCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattgatt                               <  1:650494/148‑1 (MQ=255)
                                                                                                                          atcatcGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTaaa                            >  1:42401/1‑150 (MQ=255)
                                                                                                                           tcatcGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAAt                           >  2:1293947/1‑150 (MQ=255)
                                                                                                                               cgtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATGCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAcc                       <  1:11442/150‑1 (MQ=255)
                                                                                                                                       tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAAc               >  1:1186035/1‑150 (MQ=255)
                                                                                                                                        gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGCTTGATTAAATAACCTTAATAACa              >  2:553405/1‑150 (MQ=255)
                                                                                                                                         gttgCCGGATGCGGCGTGAACGCCTTATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACaa             >  2:275128/4‑150 (MQ=38)
                                                                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATg           >  2:1208625/1‑150 (MQ=255)
                                                                                                                                            gCCGGATGCGGCGTGAACGCCTTATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGc          >  2:15707/1‑150 (MQ=38)
                                                                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTTATAACAATGc          >  2:282005/1‑150 (MQ=255)
                                                                                                                                             ccGGATGCGGCGTGAACGCCTTATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGcc         >  2:348706/1‑150 (MQ=38)
                                                                                                                                              cGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCg        <  2:178822/150‑1 (MQ=255)
                                                                                                                                              cGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCg        >  1:417955/1‑150 (MQ=255)
                                                                                                                                               ggATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGa       <  2:666599/150‑1 (MQ=255)
                                                                                                                                               ggATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGa       <  2:330751/150‑1 (MQ=255)
                                                                                                                                                    aGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTTATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa  <  2:247605/150‑1 (MQ=255)
                                                                                                                                                    aGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa  <  2:385833/150‑1 (MQ=255)
                                                                                                                                                    aGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa  >  1:15688/1‑150 (MQ=255)
                                                                                                                                                    aGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTAGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa  >  1:1116685/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                                  
AGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑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: