Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_pgi4,235,4670TG25.6% 159.1 / 65.2 118intergenic (‑2234/‑467)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base T (43/44);  new base G (6/24);  total (50/68)
Fisher's exact test for biased strand distribution p-value = 5.31e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 121 total aligned reads displayed.

AAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGG  >  NC_000913_3_hsa_pgi/4235324‑4235616
                                                                                                                                               |                                                                                                                                                                  
aaGAAACTAGAGCCTGAGCTGGATGTCAGTGCTCAAGGGACCTATCACGACGCTTCTACCAATGGGCTCATCAAGTTCATCAAGCAGCAGCGCGAGGCCAGAGTACAATAATAATCGTCGATATGTAGGCCGGATAAGGCGTTCAggccg                                                                                                                                                             <  2:1134894/150‑1 (MQ=255)
 aGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCAGGGCGGGTCTACCACGGGGCTCATCATGTTCAGCAAGCAGCCGCGCGAGGCCAGAGTCGATGTATGTTCGGCGAGTTGTAGGCCGGCTAAGGCGGTCCcgccgc                                                                                                                                                            >  2:1314229/1‑150 (MQ=255)
    aaaTAGCGCCTGAGCGTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                         <  2:228687/150‑6 (MQ=255)
    aaaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                         <  1:1136482/150‑6 (MQ=255)
    aaaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                         <  1:1479693/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)
       tAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                      <  1:43492/150‑1 (MQ=255)
        agagCCTGCGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCACCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGCTAAGGCGTGCACGCCGCATCCGGc                                                                                                                                                     >  1:1443493/1‑150 (MQ=255)
        agagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                     <  1:1203050/150‑1 (MQ=255)
        agagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                     <  2:1148505/150‑1 (MQ=255)
         gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                                    <  1:1133119/150‑2 (MQ=255)
         gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                                    <  1:287527/150‑2 (MQ=255)
         gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                                    <  2:1396932/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)
          acccTGAGCTGGATGGCAGTGCGCATGTTACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGAAATGTAGGCCGGATGCGTCGTGAACGCCTTATCCGGcct                                                                                                                                                   <  1:363655/148‑3 (MQ=255)
           gCCTGAGCTTGATGGCAGTGATCAAGTGACCTCTAACGACGCTTCTACCAATGGGCTAATCAACTTCAACAAGCACCAGCGCGAGGCCAGAGACCAATAATAATCGTCGATAATTACGACGGATAAGTCCTTAACGCCGCAACCGGCAA‑‑‑‑‑‑c                                                                                                                                           >  2:748051/1‑150 (MQ=255)
            ccTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac                                                                                                                                                 <  1:1311225/150‑5 (MQ=255)
            ccTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAAAGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac                                                                                                                                                 <  1:1013755/150‑5 (MQ=255)
             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)
              tGAGCTTGATGGCCGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCCTCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGa                                                                                                                                  <  1:571245/150‑1 (MQ=255)
              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:424624/1‑150 (MQ=255)
                  cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGcc                                                                                                                              <  1:170939/150‑1 (MQ=255)
                        ggCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACATCTTAAAGCAGCAGCGCGAGGCCAGAGTCCAATTATCATAGTCGATATGTAGGACGGATAAGGCCTTCACGCCGCATCCGGCAA‑‑‑‑‑‑T‑‑‑‑‑‑CGCTGcacgcagc                                                                                                                        >  1:801819/1‑143 (MQ=37)
                        ggCAGTGCGCAAGTGACCTCTCACGACGCTGCTACCACTGGGCTCATCAACTTCCTCCAGCAGCAGCGCGAGGCCCGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGGTCACGCCGCAACCGGCAA‑‑‑‑‑‑T‑‑‑‑‑‑CGGTGCACGATGc                                                                                                                        <  1:824929/150‑1 (MQ=25)
                           aGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑T‑‑‑‑‑‑CGGTGCACGATGcccg                                                                                                                     <  2:473170/150‑4 (MQ=37)
                            gTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGAcg                                                                                                                    >  1:385833/1‑150 (MQ=255)
                               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)
                                       aCCTCACACGACGCTTCTACCAATCGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCt                                                                                                         <  2:653562/150‑1 (MQ=255)
                                        ccTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCtt                                                                                                        <  1:923632/150‑1 (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)
                                                  cgCTTCTACCAATGGGCGCATCAACTTCATCCAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTa                                                                                              <  1:1448280/150‑1 (MQ=255)
                                                    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)
                                                             aTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCGAT‑‑‑‑‑‑CGGTGGACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTAcgtc                                                                                          >  2:794702/1‑147 (MQ=9)
                                                              tGGGCTCATCCACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGGCTTATCAGGCCTACAGGGCGATGcc                                                                                  <  1:1300645/150‑1 (MQ=255)
                                                                    catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                            <  1:1214310/150‑1 (MQ=255)
                                                                    catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                            >  2:402327/1‑150 (MQ=255)
                                                                    catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                            <  2:1340296/150‑1 (MQ=255)
                                                                    catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                            >  2:185954/1‑150 (MQ=255)
                                                                     atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                                                           <  2:422730/150‑1 (MQ=255)
                                                                     atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                                                           >  1:430881/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:666599/1‑150 (MQ=255)
                                                                                 aaGCAGCAGCGCGAGGCCAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                                                               >  1:330751/1‑150 (MQ=255)
                                                                                   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:567055/150‑1 (MQ=255)
                                                                                                 ccAGAGTCCAATAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGTCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGACAGATGCGACGCCGTCGAGTCATCACAGGCCTACAGGTCGATGCCGCTATGTACCTCGTATTCGGCAATTAATAAATAGc                                               >  1:856236/1‑150 (MQ=37)
                                                                                                       gccAATAATCATCGTCGATCTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgatt                                         <  1:1246182/149‑1 (MQ=255)
                                                                                                            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)
                                                                                                                                    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‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTAGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa        >  1:1116685/1‑150 (MQ=255)
                                                                                                                                              ttCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCg   <  2:663867/150‑1 (MQ=255)
                                                                                                                                               tCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCgg  <  1:185954/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                                  
AAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGG  >  NC_000913_3_hsa_pgi/4235324‑4235616

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: