Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,234,667 Δ1 bp intergenic (‑1434/‑499) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,234,6660A.99.2% 573.8 / ‑3.3 120intergenic (‑1433/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (57/62);  minor base T (1/0);  total (58/62)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 124 total aligned reads displayed.

GGTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_pae_pgi/4234518‑4234809
                                                                                                                                                    |                                                                                                                                               
ggTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGat                                                                                                                                                <  2:526726/150‑2 (MQ=255)
    tatataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑atcatc                                                                                                                                           <  4:427123/150‑1 (MQ=255)
     atataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcg                                                                                                                                          >  1:1185232/1‑150 (MQ=255)
      tataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcgt                                                                                                                                         <  4:20771/150‑1 (MQ=255)
        taCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcgtcg                                                                                                                                       <  1:56507/150‑1 (MQ=255)
         aCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGa                                                                                                                                      <  4:761000/150‑1 (MQ=255)
         aCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGa                                                                                                                                      >  4:316899/1‑150 (MQ=255)
          cGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAACAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGat                                                                                                                                     >  4:389283/1‑150 (MQ=255)
           gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGata                                                                                                                                    >  4:836764/1‑150 (MQ=255)
            cAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGatat                                                                                                                                   >  2:317526/1‑150 (MQ=255)
              gTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGt                                                                                                                                 <  1:1081523/150‑1 (MQ=255)
               tAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTa                                                                                                                                >  1:35183/1‑150 (MQ=255)
                 gagaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAgg                                                                                                                              >  2:1326816/1‑150 (MQ=255)
                   gaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGcc                                                                                                                            <  4:53354/150‑1 (MQ=255)
                    aTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCg                                                                                                                           <  1:1165467/150‑1 (MQ=255)
                      gCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGa                                                                                                                         >  3:513021/1‑150 (MQ=255)
                       cACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGAt                                                                                                                        <  1:317526/150‑1 (MQ=255)
                          cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                     >  1:1107901/1‑150 (MQ=255)
                          cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGAGGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                     <  1:1180582/150‑1 (MQ=255)
                           ccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcg                                                                                                                                          >  1:804346/1‑128 (MQ=255)
                           ccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcg                                                                                                                                          <  2:804346/128‑1 (MQ=255)
                            cAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCGTCGCCTCTGTAAAACTGGCAGAGAGTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGGTGCACTGTAGGCCGGATaaagg                                                                                                                   >  4:346880/1‑147 (MQ=255)
                            cAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                   <  3:395857/150‑1 (MQ=255)
                              aGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGt                                                                                                                 <  2:970603/150‑1 (MQ=255)
                              aGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGt                                                                                                                 >  4:24217/1‑150 (MQ=255)
                                ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTc                                                                                                               >  1:343855/1‑150 (MQ=255)
                                    ttttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                           <  3:512427/150‑1 (MQ=255)
                                    ttttATCCTTTTTAAGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                           <  2:710932/150‑1 (MQ=255)
                                       tATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                        >  1:815190/1‑150 (MQ=255)
                                          ccTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGt                                                                                                                                 >  3:222178/1‑122 (MQ=255)
                                           cTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGATGCGCTTGCAGTG‑ATCCTCGTCGATATGTAGTCCGGATGCGGCGTGAACGCTTCATcc                                                                                                    <  4:172268/150‑1 (MQ=255)
                                           cTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                    <  3:178710/150‑1 (MQ=255)
                                             ttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                  <  2:178633/150‑1 (MQ=255)
                                              tttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                 <  2:388018/150‑1 (MQ=255)
                                                tACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                               >  2:1000378/1‑150 (MQ=255)
                                                tACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                               <  1:834590/150‑1 (MQ=255)
                                                        gtgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCGCCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGat                                                                                                                                                <  1:288219/94‑2 (MQ=255)
                                                        gtgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGcc                                                                                       >  2:650283/1‑150 (MQ=255)
                                                           ggCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                    <  2:1185232/150‑1 (MQ=255)
                                                            gCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                   <  2:1212512/150‑1 (MQ=255)
                                                            gCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                   >  2:643150/1‑150 (MQ=255)
                                                            gCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                   >  4:536173/1‑150 (MQ=255)
                                                             cATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                                  <  2:616107/150‑1 (MQ=255)
                                                               tCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCg                                                                                <  1:643150/150‑1 (MQ=255)
                                                                  ccTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                                                                             >  2:216347/1‑150 (MQ=255)
                                                                   ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                >  4:66736/1‑114 (MQ=255)
                                                                   ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                  <  2:851365/128‑1 (MQ=255)
                                                                   ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                  >  1:851365/1‑128 (MQ=255)
                                                                   ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                            >  1:249515/1‑150 (MQ=255)
                                                                      tGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCAGGCATCGGTTGCCGGATGCggcgtgaa                                                                         <  4:410710/150‑8 (MQ=38)
                                                                         aaaaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                      <  2:841574/150‑1 (MQ=255)
                                                                          aaaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgcg                                                                     >  2:267266/1‑150 (MQ=255)
                                                                           aaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                                                                    <  4:19619/150‑1 (MQ=255)
                                                                            aCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                                   <  1:650283/150‑1 (MQ=255)
                                                                              tGGCAGAGATTAGAGGCGTAAATCCCCTGGCCTCCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                 <  1:1053749/150‑1 (MQ=255)
                                                                              tGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                 >  1:114143/1‑150 (MQ=255)
                                                                                gCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                               >  4:389001/1‑150 (MQ=255)
                                                                                 cAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                              <  1:562666/150‑1 (MQ=255)
                                                                                  agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                             >  1:236538/1‑150 (MQ=255)
                                                                                  agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                             >  1:740147/1‑150 (MQ=255)
                                                                                  agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                             <  3:40949/150‑1 (MQ=255)
                                                                                   gagaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                            <  1:267266/150‑1 (MQ=255)
                                                                                      aTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc                                                         >  3:385988/1‑150 (MQ=255)
                                                                                         agagGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                      <  4:449802/150‑1 (MQ=255)
                                                                                           agGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                    >  1:843977/1‑150 (MQ=255)
                                                                                           agGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                    <  3:729825/150‑1 (MQ=255)
                                                                                             gCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                  >  1:719919/1‑150 (MQ=255)
                                                                                             gCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTACCAGTG‑ATCATCGTCGATATGTAGGCCGCATAAGGCGCTCACGCCGCATCCGGCAACCGATGCCTAATGCGACGCGGTCGCGTCTTATCTAGCCTACAGGt                                                  >  2:860250/1‑150 (MQ=255)
                                                                                             gCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATAGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                  >  4:458231/1‑150 (MQ=255)
                                                                                               gTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCg                                                >  4:832145/1‑150 (MQ=255)
                                                                                               gTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCg                                                <  3:836764/150‑1 (MQ=255)
                                                                                                  aaTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATg                                             >  3:10397/1‑150 (MQ=255)
                                                                                                   aTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                            >  1:1175621/1‑150 (MQ=255)
                                                                                                   aTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                            >  2:610989/1‑150 (MQ=255)
                                                                                                      cccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                         <  1:1000378/150‑1 (MQ=255)
                                                                                                      cccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATACGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                         >  3:756904/1‑150 (MQ=255)
                                                                                                       ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                        <  1:859437/150‑1 (MQ=255)
                                                                                                       ccTGGCCACCCAGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGACGCCTGATACGACGCGGTCGCGTCTTATCCGGCCTACAGGTCGATGCCGat                                        <  4:288716/150‑1 (MQ=37)
                                                                                                        cTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                       >  3:514678/1‑150 (MQ=255)
                                                                                                         tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                      <  4:385988/150‑1 (MQ=255)
                                                                                                              aCCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACa                                 <  4:513021/150‑1 (MQ=255)
                                                                                                                cccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                     <  2:435259/128‑1 (MQ=255)
                                                                                                                cccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                     >  1:435259/1‑128 (MQ=255)
                                                                                                                 ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCg                              <  1:1326816/150‑1 (MQ=255)
                                                                                                                   gagaATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTa                            >  3:804452/1‑150 (MQ=255)
                                                                                                                     gaATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAtt                          <  3:94789/150‑1 (MQ=255)
                                                                                                                      aaTAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGCTTGAGTCCGATATGTACATCGTATTc                         <  1:1320122/150‑1 (MQ=25)
                                                                                                                         aGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGc                      >  4:291071/1‑150 (MQ=255)
                                                                                                                           aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                    >  3:632580/1‑150 (MQ=255)
                                                                                                                             gCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                  <  2:773168/150‑1 (MQ=255)
                                                                                                                             gCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                  <  3:19009/150‑1 (MQ=255)
                                                                                                                              cATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                            >  3:740843/1‑123 (MQ=255)
                                                                                                                              cATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                 <  3:389283/150‑1 (MQ=255)
                                                                                                                               aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                              >  1:901815/1‑72 (MQ=255)
                                                                                                                               aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                              <  2:901815/72‑1 (MQ=255)
                                                                                                                               aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa                <  1:216347/150‑1 (MQ=255)
                                                                                                                                ttGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt               >  1:258099/1‑150 (MQ=255)
                                                                                                                                 tGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATACGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa              <  4:756904/150‑1 (MQ=255)
                                                                                                                                    agagGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt           >  2:73281/1‑150 (MQ=255)
                                                                                                                                             tGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgatt  <  2:1141990/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_pae_pgi/4234518‑4234809

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