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.

AAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGG  >  NC_000913_3_hsa_pgi/4235324‑4235616
                                                                                                                                               |                                                                                                                                                                  
aaGAAACTAGAGCCTGAGCTGGATGTCAGTGCTCAAGGGACCTATCACGACGCTTCTACCAATGGGCTCATCAAGTTCATCAAGCAGCAGCGCGAGGCCAGAGTACAAT‑AATAATCGTCGATATGTAGGCCGGATAAGGCGTTCAggccg                                                                                                                                                             <  2:1134894/150‑1 (MQ=255)
 aGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCAGGGCGGGTCTACCACGGGGCTCATCATGTTCAGCAAGCAGCCGCGCGAGGCCAGAGT‑CGATGTATGTTCGGCGAGTTGTAGGCCGGCTAAGGCGGTCCcgccgc                                                                                                                                                            >  2:1314229/1‑150 (MQ=255)
    aaaTAGCGCCTGAGCGTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                         <  2:228687/150‑6 (MQ=255)
    aaaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                         <  1:1136482/150‑6 (MQ=255)
    aaaTAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAAcgccttatc                                                                                                                                                         <  1:1479693/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)
       tAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                      <  1:43492/150‑1 (MQ=255)
        agagCCTGCGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCACCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGCTAAGGCGTGCACGCCGCATCCGGc                                                                                                                                                     >  1:1443493/1‑150 (MQ=255)
        agagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                     <  1:1203050/150‑1 (MQ=255)
        agagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                                     <  2:1148505/150‑1 (MQ=255)
         gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                                    <  1:1133119/150‑2 (MQ=255)
         gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                                    <  1:287527/150‑2 (MQ=255)
         gagCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                                    <  2:1396932/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)
          acccTGAGCTGGATGGCAGTGCGCATGTTACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGAAATGTAGGCCGGATGCGTCGTGAACGCCTTATCCGGcct                                                                                                                                                   <  1:363655/148‑3 (MQ=255)
           gCCTGAGCTTGATGGCAGTGATCAAGTGACCTCTAACGACGCTTCTACCAATGGGCTAATCAACTTCAACAAGCACCAGCGCGAGGCCAGAGACCAAT‑AATAATCGTCGATAATTACGACGGATAAGTCCTTAACGCCGCAACCGGC‑AA‑‑‑‑‑‑c                                                                                                                                           >  2:748051/1‑150 (MQ=255)
            ccTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac                                                                                                                                                 <  1:1311225/150‑5 (MQ=255)
            ccTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAAAGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac                                                                                                                                                 <  1:1013755/150‑5 (MQ=255)
             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)
              tGAGCTTGATGGCCGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCCTCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGa                                                                                                                                  <  1:571245/150‑1 (MQ=255)
              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:424624/1‑150 (MQ=255)
                  cTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGcc                                                                                                                              <  1:170939/150‑1 (MQ=255)
                        ggCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACATCTTAAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑TATCATAGTCGATATGTAGGACGGATAAGGCCTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑T‑‑‑‑‑‑CGCTGcacgcagc                                                                                                                        >  1:801819/1‑143 (MQ=37)
                        ggCAGTGCGCAAGTGACCTCTCACGACGCTGCTACCACTGGGCTCATCAACTTCCTCCAGCAGCAGCGCGAGGCCCGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGGTCACGCCGCAACCGGC‑AA‑‑‑‑‑‑T‑‑‑‑‑‑CGGTGCACGATGc                                                                                                                        <  1:824929/150‑1 (MQ=25)
                           aGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑T‑‑‑‑‑‑CGGTGCACGATGcccg                                                                                                                     <  2:473170/150‑4 (MQ=37)
                            gTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGAcg                                                                                                                    >  1:385833/1‑150 (MQ=255)
                               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)
                                       aCCTCACACGACGCTTCTACCAATCGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCt                                                                                                         <  2:653562/150‑1 (MQ=255)
                                        ccTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCtt                                                                                                        <  1:923632/150‑1 (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)
                                                  cgCTTCTACCAATGGGCGCATCAACTTCATCCAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTa                                                                                              <  1:1448280/150‑1 (MQ=255)
                                                    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)
                                                             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                                                                            <  1:1214310/150‑1 (MQ=255)
                                                                    catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                            >  2:402327/1‑150 (MQ=255)
                                                                    catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                            <  2:1340296/150‑1 (MQ=255)
                                                                    catcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATg                                                                            >  2:185954/1‑150 (MQ=255)
                                                                     atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                                                           <  2:422730/150‑1 (MQ=255)
                                                                     atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                                                           >  1:430881/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:666599/1‑150 (MQ=255)
                                                                                 aaGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                                                               >  1:330751/1‑150 (MQ=255)
                                                                                   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: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)
                                                                                                            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)
                                                                                                                                    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‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTAGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa        >  1:1116685/1‑150 (MQ=255)
                                                                                                                                              ttCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCg   <  2:663867/150‑1 (MQ=255)
                                                                                                                                               tCACGCCGCATCCGGC‑AA‑‑‑‑‑‑C‑‑‑‑‑‑CGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCgg  <  1:185954/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                                  
AAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC‑AA‑‑‑‑‑‑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: ‑