Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,868,8590C.20.3% 189.3 / 51.9 69intergenic (+17/‑96)yeaG/yeaHprotein kinase, endogenous substrate unidentified; autokinase/UPF0229 family protein
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base . (6/8);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TAACGCCAAAACGTCAACCGACGAGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGG  >  NC_000913_3_pae_pgi/1868710‑1869003
                                                                                                                                                     |                                                                                                                                                 
tAACGCCAAAACGTCAACCGACGAGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTac                                                                                                                                                   <  4:245341/150‑1 (MQ=255)
    gCCAAAACGTCAACCGACGAGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACAt                                                                                                                                               >  3:446659/1‑150 (MQ=255)
    gCCAAAACGTCAACCGACGAGCAGAAGAAACACGACGACTTTGTCGACCGTATGACGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACAt                                                                                                                                               <  3:138452/150‑1 (MQ=255)
     ccAAAACGTCAACCGACGAGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATg                                                                                                                                              >  1:964691/1‑150 (MQ=255)
          aCGTCAACCGACGAGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TAcc                                                                                                                                        <  4:446659/150‑1 (MQ=255)
              cAACCGACGAGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGc                                                                                                                                    <  3:192793/150‑1 (MQ=255)
               aaCCGACGAGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGcc                                                                                                                                   <  3:429479/150‑1 (MQ=255)
                aCCGACGAGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGcca                                                                                                                                  <  2:962560/150‑2 (MQ=255)
                 ccGACGAGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCAGGCCCGGTA‑ACATGTGTACCGGGccag                                                                                                                                 <  3:7206/150‑3 (MQ=255)
                  cgacgaGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTAc                                                                                                                                <  2:964691/150‑1 (MQ=255)
                    acgaGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACaa                                                                                                                              >  1:510254/1‑150 (MQ=255)
                    acgaGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACaa                                                                                                                              >  2:952903/1‑150 (MQ=255)
                         cAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACa                                                                                                                         <  1:180282/150‑1 (MQ=255)
                         cAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACa                                                                                                                         >  1:1015180/1‑150 (MQ=255)
                          agaagaAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAg                                                                                                                        >  4:510601/1‑150 (MQ=255)
                              gaaACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGaa                                                                                                                    <  1:364949/150‑1 (MQ=255)
                                 acacGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCAGGCCCGGTA‑ACATGTGTACCGGGccagg                                                                                                                                >  1:607090/1‑132 (MQ=255)
                                 acacGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCAGGCCCGGTA‑ACATGTGTACCGGGccagg                                                                                                                                <  2:607090/135‑4 (MQ=255)
                                    cgacgaCTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTgg                                                                                                              <  2:584668/150‑1 (MQ=255)
                                      acgacTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGc                                                                                                            <  2:94213/150‑1 (MQ=255)
                                      acgacTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGc                                                                                                            <  3:510601/150‑1 (MQ=255)
                                        gacTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCg                                                                                                                                       >  3:106596/1‑121 (MQ=255)
                                        gacTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCg                                                                                                                                       <  4:106596/121‑1 (MQ=255)
                                           tttGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTgag                                                                                                       >  3:158620/1‑150 (MQ=255)
                                               tCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCt                                                                                                          <  2:916741/143‑1 (MQ=255)
                                               tCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCt                                                                                                          >  1:916741/1‑143 (MQ=255)
                                                   ccGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTg                                                                                                               <  4:183970/134‑1 (MQ=255)
                                                   ccGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTg                                                                                                               >  3:183970/1‑134 (MQ=255)
                                                          atgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCGGGCCCGGTA‑ACATGTGTACCGGGccag                                                                                                                                 >  1:860003/1‑107 (MQ=39)
                                                          atgGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCGGGCCCGGTA‑ACATGTGTACCGGGccag                                                                                                                                 <  2:860003/109‑3 (MQ=39)
                                                            ggAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCg                                                                                      >  2:119081/1‑150 (MQ=255)
                                                                         acCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGGACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTg                                                                         >  3:433299/1‑150 (MQ=255)
                                                                           ccGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGc                                                                       >  1:361172/1‑150 (MQ=255)
                                                                                         gTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACACCCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGCAGTACggg                                                         >  2:477811/1‑150 (MQ=255)
                                                                                            tACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACgggggg                                                      <  2:397155/150‑1 (MQ=255)
                                                                                            tACTGTGCGAATGGTATTTGCGCGTACGTAAAGCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACgggggg                                                      <  2:51619/150‑1 (MQ=255)
                                                                                             aCTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTgag                                                                                                       <  3:276758/100‑1 (MQ=255)
                                                                                             aCTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTgag                                                                                                       >  4:276758/1‑100 (MQ=255)
                                                                                             aCTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGc                                                     <  1:952903/150‑1 (MQ=255)
                                                                                              cTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCa                                                    <  1:572722/150‑1 (MQ=255)
                                                                                                gtgCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCata                                                  >  2:305957/1‑150 (MQ=255)
                                                                                                 tgCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCatat                                                 >  4:378926/1‑150 (MQ=255)
                                                                                                         gTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGt                                         >  2:15552/1‑150 (MQ=255)
                                                                                                           aTTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGttt                                       >  2:1014338/1‑150 (MQ=255)
                                                                                                                  cgtacgtaAATCGTCTTAACAACCCTGGCCCGGTACACATTT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGAcc                                <  4:541292/150‑1 (MQ=255)
                                                                                                                    tacgtaAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCgg                              <  1:1014338/150‑1 (MQ=255)
                                                                                                                            tCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGc                                                                                                                                    <  4:473623/40‑1 (MQ=14)
                                                                                                                            tCGTCTTAACAACCCTGGCCCGGTA‑ACATGTGTACCGGGc                                                                                                                                    >  3:473623/1‑40 (MQ=14)
                                                                                                                               tCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACgg                   >  2:670139/1‑150 (MQ=255)
                                                                                                                               tCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGGAGTACGGGGGGGCTATGACCTGGCTTATTGCCCGGCGTCTGACCgg                   >  2:59065/1‑150 (MQ=255)
                                                                                                                                 ttAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGtt                                                                          >  2:744622/1‑93 (MQ=255)
                                                                                                                                 ttAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGtt                                                                          <  1:744622/93‑1 (MQ=255)
                                                                                                                                   aacaacCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAAt                                                                   <  1:351592/98‑1 (MQ=255)
                                                                                                                                   aacaacCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAAt                                                                   >  2:351592/1‑98 (MQ=255)
                                                                                                                                     caacCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGCCATCTGAACGGCaaaaa             >  3:352062/1‑150 (MQ=255)
                                                                                                                                       ttgtaggCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTAtt                                    <  2:735330/120‑1 (MQ=255)
                                                                                                                                       ttgtaggCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTAtt                                    >  1:735330/6‑125 (MQ=255)
                                                                                                                                        cccTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGc                             >  2:408707/1‑131 (MQ=255)
                                                                                                                                        cccTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGc                             <  1:408707/131‑1 (MQ=255)
                                                                                                                                         gtaggCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGt                                         >  4:285278/4‑118 (MQ=255)
                                                                                                                                         ggaggCCCGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGt                                         <  3:285278/115‑1 (MQ=255)
                                                                                                                                         ccTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACaaa         <  2:361172/150‑1 (MQ=255)
                                                                                                                                         ccTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACa           <  3:480573/148‑1 (MQ=255)
                                                                                                                                         ccTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACa           >  4:480573/1‑148 (MQ=255)
                                                                                                                                          cTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACaaaa        <  1:168696/150‑1 (MQ=255)
                                                                                                                                           tGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTAtt                                    <  1:819158/121‑1 (MQ=255)
                                                                                                                                           tGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTAtt                                    >  2:819158/1‑121 (MQ=255)
                                                                                                                                           tGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAg       >  4:214900/1‑150 (MQ=255)
                                                                                                                                              cccGGTA‑ACATGTGTACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCAt    >  1:689855/1‑150 (MQ=255)
                                                                                                                                                cGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATgg  <  3:378926/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TAACGCCAAAACGTCAACCGACGAGCAGAAGAAACACGACGACTTTGTCGACCGTATGATGGAAAAAGGCTACACCCGTAAACAGGTGCGTTTACTGTGCGAATGGTATTTGCGCGTACGTAAATCGTCTTAACAACCCTGGCCCGGTACACATGT‑TACCGGGCCTACAACGACAGCGAACCGTGGGCCTGAGAAGCGGCAACACAGGCGTAGCATACAGTTGGCAAATGTAGTACGGGGGGCATATGACCTGGTTTATTGACCGGCGTCTGAACGGCAAAAACAAAAGCATGG  >  NC_000913_3_pae_pgi/1868710‑1869003

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