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.4% 726.1 / ‑1.1 154intergenic (‑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 . (83/69);  minor base C (1/0);  total (85/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 158 total aligned reads displayed.

GGGTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGA  >  NC_000913_3_pae_pgi/4234517‑4234811
                                                                                                                                                     |                                                                                                                                                 
gggTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGa                                                                                                                                                   <  2:1010570/150‑1 (MQ=255)
 ggTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGat                                                                                                                                                  <  2:329989/150‑2 (MQ=255)
      atataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcg                                                                                                                                            >  3:438216/1‑150 (MQ=255)
      atataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcg                                                                                                                                            >  1:117594/1‑150 (MQ=255)
        ataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcgtc                                                                                                                                          <  1:732115/150‑1 (MQ=255)
            gCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGata                                                                                                                                      <  2:251129/150‑1 (MQ=255)
                tcgagaTCCACCCCAAGGGCGTTTTATCCTTTTTGCAAGATGTAGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGCCCACCCCAAGAATAGAAGCATTGAACAGGCGCTTGCAGTG‑ATCATCATCGAAATgga                                                                                                                                  >  2:1673932/3‑148 (MQ=255)
                 agagATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAg                                                                                                                                 <  3:678598/150‑1 (MQ=255)
                  gagaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAgg                                                                                                                                >  3:875674/1‑150 (MQ=255)
                  gagaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGACCTGTAgg                                                                                                                                >  4:900103/1‑150 (MQ=255)
                  gagaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGACCTGTAgg                                                                                                                                >  4:754080/1‑150 (MQ=255)
                     aTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCg                                                                                                                             <  4:482936/150‑1 (MQ=255)
                       gCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGGTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGa                                                                                                                           <  4:1005872/150‑1 (MQ=255)
                        cACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGAt                                                                                                                          >  2:1133695/1‑150 (MQ=255)
                          ccccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATaa                                                                                                                        <  1:516817/150‑1 (MQ=255)
                           cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                       <  3:351550/150‑1 (MQ=255)
                           cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                       >  2:1259692/1‑150 (MQ=255)
                           cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                       <  1:913137/150‑1 (MQ=255)
                           cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                       >  4:342878/1‑150 (MQ=255)
                            ccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAgg                                                                                                                      <  4:879894/150‑1 (MQ=255)
                            ccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCGTGCAGTG‑ATCATCGTCCATATGAAGGCCGGATAAgg                                                                                                                      >  2:794429/1‑150 (MQ=255)
                                     ttttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                             >  1:826818/1‑150 (MQ=255)
                                     ttttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                             >  1:827214/1‑150 (MQ=255)
                                        tATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                          <  1:559176/150‑1 (MQ=255)
                                          tCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                        >  1:1639010/1‑150 (MQ=255)
                                            cTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                      >  2:1378438/1‑150 (MQ=255)
                                             tttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATACGTAGGCCGGATAAAGCGTTCACGCCGCATCCg                                                                                                     >  2:1583677/1‑150 (MQ=255)
                                                ttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCAGATAAAGCGTTCACGCCGAATCAGGCa                                                                                                  >  2:1345409/1‑150 (MQ=255)
                                                 tACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                                 <  1:147311/150‑1 (MQ=255)
                                                  aCGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                                <  2:1056920/150‑1 (MQ=255)
                                                  aCGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCAGTGAAGAGGCGCTTGCAGTG‑AGCATCGTCGATATGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAAc                                                                                                >  2:127780/1‑150 (MQ=255)
                                                   cGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                                                                               >  4:129451/1‑150 (MQ=255)
                                                     agaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGa                                                                                             >  1:1240626/1‑150 (MQ=255)
                                                      gaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                            <  3:743291/150‑1 (MQ=255)
                                                      gaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                            <  4:415232/150‑1 (MQ=255)
                                                          tgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                        >  1:100884/1‑150 (MQ=255)
                                                          tgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                        <  4:50931/150‑1 (MQ=255)
                                                          tgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGAGGCACTGAAGAGGCGCTGGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCACCGATGCCt                                                                                        >  2:1573059/1‑150 (MQ=39)
                                                            ggCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                      <  2:27800/150‑1 (MQ=255)
                                                              cATCGCCTCTGTCCAACTGGCAGCGATTCGAGGCGTCAGCCTCCTGTTCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCATTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                                    <  1:1596472/150‑1 (MQ=255)
                                                               aTCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGc                                                                                   >  1:1047857/1‑150 (MQ=255)
                                                                  gCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                                >  1:170935/1‑150 (MQ=255)
                                                                     tctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                                <  2:672871/131‑1 (MQ=255)
                                                                     tctGTAAAACTGGCAGAGATTAGAGGCGCAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATCTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGACTGATCCGAcgag                                                                             >  2:700197/1‑148 (MQ=39)
                                                                       tGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGata                                                                                                                                      <  4:151041/91‑1 (MQ=255)
                                                                       tGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGata                                                                                                                                      >  3:151041/1‑91 (MQ=255)
                                                                        gTAAAACTGGCAGAGATTAGAGGCGTAACTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                          <  1:1307787/150‑1 (MQ=255)
                                                                        gTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcgt                                                                                                                                           <  1:1170448/85‑1 (MQ=255)
                                                                        gTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                          >  3:131567/1‑150 (MQ=255)
                                                                           aaaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑GTCATCGGCGGTAGGTAGGCCGGATGAGGCGATCACGCCGCATCCGCCAACGGATGCCTGCTGCGCCGCGGGcgcg                                                                       >  3:778697/1‑150 (MQ=255)
                                                                           aaaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgcg                                                                       >  3:983277/1‑150 (MQ=255)
                                                                           aaaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgcg                                                                       >  1:165509/1‑150 (MQ=255)
                                                                               tGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGACGCTCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                   <  3:808323/150‑1 (MQ=255)
                                                                               tGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCAGTGAAGAGGCGCTTGCACTG‑ACCATCGCCGATACGTCGGCCGGACAAGGCGTTAACGCCGCAGCCGGCCACCGATGCCTGATGCGACGCCGTCGCGTCtt                                                                   >  2:1131023/1‑150 (MQ=255)
                                                                                 gCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                                 >  1:1115231/1‑150 (MQ=255)
                                                                                  cAGGGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                                <  1:1462657/150‑1 (MQ=255)
                                                                                   agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                               >  3:111376/1‑150 (MQ=255)
                                                                                       aTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc                                                           <  1:1259692/150‑1 (MQ=255)
                                                                                         tAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcatc                                                         <  2:1639010/150‑4 (MQ=255)
                                                                                            agGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                      >  4:596300/1‑150 (MQ=255)
                                                                                            agGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                      <  1:706947/150‑1 (MQ=255)
                                                                                              gCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                    >  2:557679/1‑150 (MQ=255)
                                                                                              gCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                    <  4:482774/150‑1 (MQ=255)
                                                                                              gCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                    >  2:682317/1‑150 (MQ=255)
                                                                                                  aaaTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGAt                                                <  3:129451/150‑1 (MQ=255)
                                                                                                     tCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGata                                                                                                                                      <  2:752813/61‑1 (MQ=255)
                                                                                                     tCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc                                             >  2:654560/1‑150 (MQ=255)
                                                                                                       cccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                           >  2:326988/1‑150 (MQ=255)
                                                                                                       cccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGCCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                           >  2:1475428/1‑150 (MQ=255)
                                                                                                        ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATCTGTAGGCGGGATCAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGCCGCGGTCGCGTCTTATCAGGCATAAAGGTCGCTGCCGat                                          >  4:355680/1‑150 (MQ=37)
                                                                                                        ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                          <  2:170935/150‑1 (MQ=255)
                                                                                                        ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                          >  3:596210/1‑150 (MQ=255)
                                                                                                          tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGa                                                                                                                                        >  1:626671/1‑54 (MQ=255)
                                                                                                          tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGCTGCCGatat                                        <  1:1072445/150‑1 (MQ=255)
                                                                                                            gCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                      >  3:1021433/1‑150 (MQ=255)
                                                                                                             ccACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTa                                     >  3:402894/1‑150 (MQ=255)
                                                                                                               aCCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACa                                   >  1:1145445/1‑150 (MQ=255)
                                                                                                               aCCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACa                                   >  4:798471/1‑150 (MQ=255)
                                                                                                                ccccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                                  <  4:875674/150‑1 (MQ=255)
                                                                                                                 cccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATCTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                                 <  3:465746/150‑1 (MQ=255)
                                                                                                                 cccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                                 >  2:695385/1‑150 (MQ=255)
                                                                                                                   cGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                                                                    >  2:1480935/1‑113 (MQ=255)
                                                                                                                         tAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                         <  1:1281610/150‑1 (MQ=255)
                                                                                                                          aGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGc                        <  2:100884/150‑1 (MQ=255)
                                                                                                                            aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                      >  2:513639/1‑150 (MQ=255)
                                                                                                                            aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                      >  1:202440/1‑150 (MQ=255)
                                                                                                                              gCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                    <  2:1175034/150‑1 (MQ=255)
                                                                                                                               cATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                   >  4:794336/1‑150 (MQ=255)
                                                                                                                               cATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa                   >  1:229435/1‑150 (MQ=255)
                                                                                                                                aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCCGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa                  >  3:108606/1‑150 (MQ=255)
                                                                                                                                    aaGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACa              >  3:580243/1‑150 (MQ=255)
                                                                                                                                     agagGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATatcc             >  3:275603/1‑147 (MQ=255)
                                                                                                                                      gagGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                         <  3:49046/137‑1 (MQ=255)
                                                                                                                                      gagGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa            <  1:326988/150‑1 (MQ=255)
                                                                                                                                           gcTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACg       >  3:323595/1‑150 (MQ=255)
                                                                                                                                             ttGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat     <  2:1240626/150‑1 (MQ=255)
                                                                                                                                              tGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgatt    <  1:965547/150‑1 (MQ=255)
                                                                                                                                               gCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGATCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATACCGATATGTACATCGCATTCGGCAATCAATACATAGCACgattg   >  2:1557478/1‑150 (MQ=25)
                                                                                                                                                cAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattga  <  1:700197/150‑1 (MQ=255)
                                                                                                                                                cAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattga  <  2:1210945/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGGTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGA  >  NC_000913_3_pae_pgi/4234517‑4234811

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: ‑