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.0% 481.4 / ‑2.4 101intergenic (‑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 . (47/53);  minor base G (1/0);  total (48/53)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 107 total aligned reads displayed.

GGGTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGA  >  NC_000913_3_pae_pgi/4234517‑4234811
                                                                                                                                                     |                                                                                                                                                 
gggTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGa                                                                                                                                                   <  4:236296/150‑1 (MQ=255)
gggTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGa                                                                                                                                                   <  4:151529/150‑1 (MQ=255)
gggTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGa                                                                                                                                                   >  3:368972/1‑150 (MQ=255)
  gTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatc                                                                                                                                                 >  2:377832/1‑148 (MQ=255)
   tCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCAACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatca                                                                                                                                                >  1:935245/1‑147 (MQ=255)
    cTATATCCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGatcat                                                                                                                                               <  1:254248/150‑5 (MQ=255)
      atataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcg                                                                                                                                            >  3:557749/1‑150 (MQ=255)
      atataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcg                                                                                                                                            >  1:280252/1‑150 (MQ=255)
      atataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcg                                                                                                                                            >  1:988800/1‑150 (MQ=255)
       tataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcgt                                                                                                                                           <  2:607714/150‑1 (MQ=255)
        ataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcgtc                                                                                                                                          <  1:839949/150‑1 (MQ=255)
        ataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcgtc                                                                                                                                          <  1:648356/150‑1 (MQ=255)
          aCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGa                                                                                                                                        >  1:420109/1‑150 (MQ=255)
          aCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGACATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGa                                                                                                                                        >  4:619198/1‑150 (MQ=255)
           cGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCCTGCAGTGATCATCGTCGat                                                                                                                                       >  2:272708/1‑150 (MQ=255)
             cAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGatat                                                                                                                                     <  3:743148/150‑1 (MQ=255)
               gTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGt                                                                                                                                   <  1:403560/150‑1 (MQ=255)
               gTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGt                                                                                                                                   <  2:1178519/150‑1 (MQ=255)
                     aTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCg                                                                                                                             <  4:574440/150‑1 (MQ=255)
                       gCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGa                                                                                                                           <  1:55575/150‑1 (MQ=255)
                         aCCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATa                                                                                                                         <  3:88682/150‑1 (MQ=255)
                          cctcAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATaa                                                                                                                        >  2:350542/4‑150 (MQ=255)
                          ccccAAGGGCGTTTTATCCTTTTTACGAGATGTTGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATaa                                                                                                                        <  3:621475/150‑1 (MQ=255)
                           cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                       >  3:502372/1‑150 (MQ=255)
                           cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                       <  2:757009/150‑1 (MQ=255)
                            ccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCTGATAAgg                                                                                                                      <  4:648510/150‑1 (MQ=255)
                             cAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                     >  1:695056/1‑150 (MQ=255)
                                gggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                  <  2:1130366/150‑1 (MQ=255)
                                    gTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcg                                                                                                              <  2:330408/150‑1 (MQ=255)
                                       ttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                           >  4:430500/1‑150 (MQ=255)
                                         aTCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCCTGCAGCGATCATCGTCGATATGTAGGCCGCATAAGGAGTTCACACCGCa                                                                                                         >  2:128437/1‑150 (MQ=255)
                                          tCCTTTTTACGTGTTTTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                        >  4:336902/1‑150 (MQ=255)
                                           ccTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                       >  1:841307/1‑150 (MQ=255)
                                             tttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGGGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                     >  1:699201/1‑150 (MQ=255)
                                             attactcGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                     <  3:252112/144‑1 (MQ=255)
                                              ttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                    <  1:840814/150‑1 (MQ=255)
                                                  aCGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                                >  1:325194/1‑150 (MQ=255)
                                                    gagaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                            >  2:473436/1‑136 (MQ=255)
                                                    gagaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                            <  1:473436/136‑1 (MQ=255)
                                                    gagaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCg                                                                                              >  3:573910/1‑150 (MQ=255)
                                                    gagaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAACACAAGCACTGAAGAGCCGCCGCCAGCGCTCATCCCCGCTACGCAGCCCGGCCCAGGCGTCCACCCCGCATCCCGCCAccc                                                                                              >  2:901887/1‑149 (MQ=255)
                                                    gagaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACACCGAGAATAGATGCATTGAAGAGGCGCTGGCAGGGATCATCGTCGATATGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCAACCg                                                                                              >  2:286288/1‑150 (MQ=39)
                                                     agaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGa                                                                                             <  4:459850/150‑1 (MQ=255)
                                                     agaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCACGAGAATAGATACATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGa                                                                                             >  4:582786/1‑150 (MQ=255)
                                                       aTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATg                                                                                           <  2:318555/150‑1 (MQ=255)
                                                        tgtgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                          <  1:14076/150‑1 (MQ=255)
                                                         gtgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGcc                                                                                         >  4:329845/1‑150 (MQ=255)
                                                          tgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                        <  1:236212/150‑1 (MQ=255)
                                                          tgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                        <  2:1011306/150‑1 (MQ=255)
                                                                   ccTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGCTGCGAcg                                                                               >  4:463411/1‑150 (MQ=255)
                                                                   ccTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                                                                               <  2:809710/150‑1 (MQ=255)
                                                                    ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                              >  2:85518/1‑150 (MQ=255)
                                                                         tAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                          >  4:336972/1‑149 (MQ=255)
                                                                         tAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                          <  3:336972/149‑1 (MQ=255)
                                                                             aCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                                     <  3:329845/150‑1 (MQ=255)
                                                                               tGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                                     >  2:386566/1‑148 (MQ=255)
                                                                               tGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                   <  3:123188/150‑1 (MQ=255)
                                                                                ggCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                  <  2:1046528/150‑1 (MQ=255)
                                                                                  cAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTc                                                                                                                 <  1:1148370/101‑1 (MQ=255)
                                                                                  cAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTc                                                                                                                 >  2:1148370/1‑101 (MQ=255)
                                                                                   agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                               >  2:560236/1‑150 (MQ=255)
                                                                                     agagTAGTGGCTTAAATCCCCTGGCCGCCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAgg                                                             <  4:466492/150‑1 (MQ=37)
                                                                                     agaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg                                                             <  4:771394/150‑1 (MQ=255)
                                                                                            agGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTGTTATCAGGCCTACAg                                                      >  3:390132/1‑150 (MQ=255)
                                                                                                  aaaTCCCCTGGCCACCCCGAGAATAGATGCATTGATGAGGCGATTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGCAGCGTTTTATCAGGCCTCCAGCTCGAt                                                >  1:767722/1‑150 (MQ=37)
                                                                                                     tCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc                                             >  1:343547/1‑150 (MQ=255)
                                                                                                     tCCCCTGGCCACCACGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc                                             <  3:582786/150‑1 (MQ=255)
                                                                                                      ccccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCg                                            >  1:1152820/1‑150 (MQ=255)
                                                                                                       cccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                            <  3:370388/85‑1 (MQ=255)
                                                                                                       cccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                            >  4:370388/1‑85 (MQ=255)
                                                                                                        ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                    >  3:632155/1‑92 (MQ=255)
                                                                                                        ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                    <  4:632155/92‑1 (MQ=255)
                                                                                                        ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGAAAAGGCGTCTACGCCGCCTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGCCTTATCAGGCCTACAGGTCGAAGCCGat                                          >  4:551149/1‑150 (MQ=255)
                                                                                                          tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                        <  1:757003/150‑1 (MQ=255)
                                                                                                          tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                        >  2:357556/1‑150 (MQ=255)
                                                                                                            gCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                      <  1:873826/150‑1 (MQ=255)
                                                                                                            gCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                      <  1:79258/150‑1 (MQ=255)
                                                                                                            gCCACCCCGAGAATAGATGCATTGAAGAGGCACTTGCAGTGATCATCGTCGATATGTAGGCCGGATACGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                      >  2:700951/1‑150 (MQ=255)
                                                                                                                 cccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                                 <  2:846608/150‑1 (MQ=255)
                                                                                                                  ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc                                             >  2:1106933/1‑137 (MQ=255)
                                                                                                                  ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc                                             <  1:1106933/137‑1 (MQ=255)
                                                                                                                  ccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCg                                <  1:350542/150‑1 (MQ=255)
                                                                                                                   cGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGt                               >  1:1189506/1‑150 (MQ=255)
                                                                                                                       aaTAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                           <  1:85518/150‑1 (MQ=255)
                                                                                                                        aTAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                          <  1:753043/150‑1 (MQ=255)
                                                                                                                         tAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                         <  1:534801/150‑1 (MQ=255)
                                                                                                                            aTGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                      >  2:975927/1‑150 (MQ=255)
                                                                                                                             tGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAt                     >  2:368237/1‑150 (MQ=255)
                                                                                                                                 ttGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt                 <  4:368972/150‑1 (MQ=255)
                                                                                                                                 ttGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt                 >  1:240138/1‑150 (MQ=255)
                                                                                                                                  tGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa                <  4:502372/150‑1 (MQ=255)
                                                                                                                                  tGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa                <  4:573910/150‑1 (MQ=255)
                                                                                                                                    aaGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACa              >  1:1114364/1‑150 (MQ=255)
                                                                                                                                       agGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAg           <  2:988800/150‑1 (MQ=255)
                                                                                                                                         gcgcTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCa         >  3:4356/1‑150 (MQ=255)
                                                                                                                                         gcgcTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCa         >  1:1121256/1‑150 (MQ=255)
                                                                                                                                             ttGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat     <  1:1180066/150‑1 (MQ=255)
                                                                                                                                             ttGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat     >  4:688620/1‑150 (MQ=255)
                                                                                                                                             ttGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAACCGATGAATGATCCGACGCGGTCGCGTTTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat     <  1:548762/150‑1 (MQ=255)
                                                                                                                                                cAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattga  >  1:677670/1‑150 (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: