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.97.2% 303.0 / 0.3 72intergenic (‑1433/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base . (35/34);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAT  >  NC_000913_3_pae_pgi/4234518‑4234808
                                                                                                                                                    |                                                                                                                                              
ggTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGa                                                                                                                                                <  1:1397264/149‑1 (MQ=255)
 gTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGat                                                                                                                                               <  1:279625/149‑2 (MQ=255)
  tCTAGATACGTAGTAGAGATGCTGCGTAAGGGCGTGTTCTCCTGTTTACGAGATGTGGGCATCGCCTCTGTAAAGCTGGCAGAGAGTAGAGGCGTAAGTCGCTTGGCCACCCCGTGAACAGGTGCATGGGAGAGGCGCTTGCAGTGatc                                                                                                                                              <  2:519605/149‑3 (MQ=255)
         aCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGACCTGAAGAGGCGCTTGCAGTG‑ATCAtcgtcg                                                                                                                                      <  1:1048794/149‑1 (MQ=255)
          cGCAGGAGAGATGGACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGa                                                                                                                                     <  2:513473/149‑1 (MQ=255)
            cAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCCGAGATTTGAGGCGTAACTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGata                                                                                                                                   <  1:553135/149‑1 (MQ=255)
              gTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCTCGAATAGATGCATTGAAGAGGCGCTTTCAGTG‑ATCATCGTCGattta                                                                                                                                 >  2:528885/1‑146 (MQ=255)
                 gagaTGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAg                                                                                                                              <  1:704707/149‑1 (MQ=255)
                     tGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCg                                                                                                                          <  2:419422/149‑1 (MQ=255)
                     tGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGAAAATACCCTTTCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCg                                                                                                                          <  1:627353/149‑1 (MQ=255)
                      gCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCgg                                                                                                                         >  1:405965/1‑149 (MQ=255)
                       cACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGa                                                                                                                        >  2:1159438/1‑149 (MQ=255)
                        aCCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAA‑CATCGCCGCTATTTAGGCCgtag                                                                                                                       >  2:323715/1‑146 (MQ=255)
                         ccccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATa                                                                                                                      >  1:1336729/1‑149 (MQ=255)
                              aGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                                                                 >  1:807955/1‑149 (MQ=255)
                              aGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                                                                 <  1:1124196/149‑1 (MQ=255)
                                    ttttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAG‑GGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcg                                                                                                           >  1:567869/1‑149 (MQ=255)
                                           cTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcgt                                                                                                                                        >  1:111569/1‑113 (MQ=255)
                                           cTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCAtcgt                                                                                                                                        <  2:111569/113‑1 (MQ=255)
                                            tttttACGAGATGTGGGCATCGCCTCTGTAAAACTTGCAGAGATTAGAGGCGTAAATCCCCTTGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTCGGCCGGATTAGGCGTTTACGCCGCTTcc                                                                                                   >  2:803729/1‑149 (MQ=255)
                                            tttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATCTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                   >  1:1125639/1‑149 (MQ=255)
                                             ttttACGAGATGTGGGCATCTCCTCTGTAAAATTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCACGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                  <  1:733080/149‑1 (MQ=39)
                                                  tgagaTGGGGGCATCGCCTCTGTAAAACTGGCAGAGATTCGTGGGGTAATTTCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                             <  1:970613/148‑1 (MQ=38)
                                                   gagaTGTGGGCATCGCCTCGGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGAGGCATTGAAGAGGCGCTTGCAGTG‑ATCAACGTCGATATGTAGACCGGAGAAGGCGTTCACCCCGCATCCGGCAAcc                                                                                            >  2:970613/1‑149 (MQ=255)
                                                        gtgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                       >  2:572313/1‑149 (MQ=255)
                                                           ggCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTg                                                                                    >  1:19548/1‑149 (MQ=255)
                                                            gCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                   >  2:957344/1‑149 (MQ=255)
                                                            gCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                   >  1:101521/1‑149 (MQ=255)
                                                             cATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                  <  1:572313/149‑1 (MQ=255)
                                                                  tctctGTAGACCTGGCATAGATTCGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATTCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGTTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                             <  2:1181836/148‑1 (MQ=37)
                                                                  ccTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                             <  1:998282/149‑1 (MQ=255)
                                                                  ccTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                             >  1:1146547/1‑149 (MQ=255)
                                                                  ccTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                             >  1:1287526/1‑149 (MQ=255)
                                                                    attgaaaaaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAAATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                           <  2:486465/144‑1 (MQ=255)
                                                                       gTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCTCATCCGGCAACCGAGGCCTGATGCGCCGCGGt                                                                        >  1:454423/1‑149 (MQ=255)
                                                                       gTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                        >  2:1041882/1‑149 (MQ=255)
                                                                          aaaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                     >  1:1182613/1‑149 (MQ=255)
                                                                          aaaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAAAAGATGCAAAGCAGAAGCGCATACAGT‑AATCATCGTCGATATGTATGCCGGCTAAGGCGGTGACGCGGCAGACGGCGACCGATACCTGCGGGGAcgctatgga                                                                     >  2:974411/1‑143 (MQ=255)
                                                                           aaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATAGACGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAAGCGTTCACACCGCAGCCTGCCACCGATGCCTGATGCTACGGGGTcgcc                                                                    >  2:1152086/1‑148 (MQ=255)
                                                                               ggCAGAAATTAGAGGCGTAAATCCCCTGGGCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                <  1:413900/149‑1 (MQ=255)
                                                                                 cAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                              <  2:19548/149‑1 (MQ=255)
                                                                                   gagaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTGGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                            <  1:1170999/149‑1 (MQ=255)
                                                                                   gagaGTAGAGGCGTGAACCCCCGTGCCACCGCGAGAAGAGAGGCATTGAGGAGGCGTTTGCAGTG‑ATCATCGTCGATGTGTAGGCCGGATAAGGCGTTCACGCCCCATCTTGGAACCGATGCCTGTTGCGACGCGGTCGCGTCTTATCa                                                            <  2:242191/149‑1 (MQ=21)
                                                                                    ataTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGGGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGTGACGCGGTCGCGTCTTATCAg                                                           <  2:640807/147‑1 (MQ=255)
                                                                                    agaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                           <  1:1041882/149‑1 (MQ=255)
                                                                                      aTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAgg                                                          <  2:228990/148‑1 (MQ=38)
                                                                                      aTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAgg                                                          >  1:228990/1‑148 (MQ=38)
                                                                                       ttAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCGTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCGGATGCGACGCGGTCGCGTCTTATCAGGcc                                                        <  1:31042/149‑1 (MQ=255)
                                                                                          ggggCGTAAATCTCTTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                     <  1:71781/147‑1 (MQ=255)
                                                                                          gagGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                     >  1:68387/1‑149 (MQ=255)
                                                                                             gCGTAAATCCCCTGGCCACCCCGAGCATAGATGCATTGACGAGGAGCTTGCAGTG‑ATCATCGTCGATATATAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGCTGCATGATGCGAAGCGGTCGCGTCTTATCAGGCATACCgg                                                  >  1:562807/1‑149 (MQ=37)
                                                                                             gCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                  <  2:1160880/85‑1 (MQ=255)
                                                                                             gCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                                  >  1:1160880/1‑85 (MQ=255)
                                                                                               ggaaaTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                <  2:142209/147‑1 (MQ=255)
                                                                                               gTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                >  1:1145291/1‑149 (MQ=255)
                                                                                                 aaaTCCCCTGGCCACCGCGAGAATTGTTGCATTGAAGAGTCGCTAGGTGTG‑ATCATCGTCGATATGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                              <  1:547620/149‑1 (MQ=25)
                                                                                                  aaTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGCTGCGACGCGGTCGCGTCTTATCAGGCCTAGAGGTCTAt                                             >  2:28033/1‑149 (MQ=255)
                                                                                                    tCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                           <  2:68387/149‑1 (MQ=255)
                                                                                                     ccccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc                                          >  1:345497/1‑149 (MQ=255)
                                                                                                      cccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGAAATGAAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGCCTTATCAGGCCTACAGGTCGATGCCg                                         >  1:478521/1‑149 (MQ=255)
                                                                                                              aCCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCATTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCCGGCCTACAGGTCGCTGCCGGTATGTAc                                 >  2:328667/1‑149 (MQ=255)
                                                                                                                 cagagaATAGATGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                              <  1:165378/147‑1 (MQ=255)
                                                                                                                         cgttGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                      <  1:133846/146‑1 (MQ=255)
                                                                                                                           aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTTCCGTGTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTGTTCGGCa                    >  1:681061/1‑149 (MQ=38)
                                                                                                                           aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                    >  2:1272294/1‑149 (MQ=255)
                                                                                                                           aTGCATTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGGATTCGGCa                    >  1:237474/1‑149 (MQ=255)
                                                                                                                           aTGCAGTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                    <  2:993538/149‑1 (MQ=255)
                                                                                                                               aTTGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGGCGATGCCGATATGTACATCGTATTCGGCAATTa                >  1:444352/1‑149 (MQ=255)
                                                                                                                                ttGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                       >  1:292059/1‑125 (MQ=255)
                                                                                                                                ttGAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                       <  2:292059/125‑1 (MQ=255)
                                                                                                                                  gAAGAGGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa             <  1:1159438/149‑1 (MQ=255)
                                                                                                                                    agagGCGCTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCCGTCGCGTATTATCAGGCCTACAGGTCGACGCCGATATGTACATCGTATTCGGCAATtactacc           >  1:1253088/1‑148 (MQ=25)
                                                                                                                                       ggggcTTGCAGTG‑AGCAGCGTCGGTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAg        <  2:689332/146‑1 (MQ=25)
                                                                                                                                        gcgcTTGCAGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGc       >  2:288539/1‑149 (MQ=255)
                                                                                                                                             tGCAGTG‑ATCAGCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat  <  2:108729/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GGTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAT  >  NC_000913_3_pae_pgi/4234518‑4234808

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