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.98.4% 282.4 / ‑3.4 63intergenic (‑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 . (31/31);  minor base C (0/1);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCTATATACGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAT  >  NC_000913_3_pae_pgi/4234518‑4234808
                                                                                                                                                    |                                                                                                                                              
ggactgtataCGCACCAGAGACCACCCCCACGGGCGTACCAGCCTCTTAACGAGACGTGGGCATCGCCTCTGTAAAAATGGCAGCGATTATAGGCATAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGTCGCTTGCAGTGa                                                                                                                                                <  2:477975/143‑1 (MQ=255)
      tataCGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcg                                                                                                                                         <  1:818578/149‑1 (MQ=255)
          cGCAGTAGAGATGAACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGa                                                                                                                                     <  2:435393/149‑1 (MQ=255)
                         ccccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCCGTGATCATCGCCGATATTTAGGCCGGATa                                                                                                                      >  2:463455/1‑149 (MQ=255)
                          cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGcc                                                                                                                           >  2:746630/1‑143 (MQ=255)
                          cccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGcc                                                                                                                           <  1:746630/143‑1 (MQ=255)
                           ccAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                    <  2:217693/149‑1 (MQ=255)
                                ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGat                                                                                                                                    >  2:889696/1‑128 (MQ=255)
                                ggCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGat                                                                                                                                    <  1:889696/128‑1 (MQ=255)
                                     tttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGAt                                                                                                                       >  2:490223/1‑136 (MQ=255)
                                     tttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGAt                                                                                                                       <  1:490223/136‑1 (MQ=255)
                                      ttATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                         <  1:965220/149‑1 (MQ=255)
                                        aTCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATCTGTAAGCCGGATAAGGCGGTCAcgccgc                                                                                                       >  2:969388/1‑149 (MQ=255)
                                            tttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGTATAAGGCGTTCACGCCGCATcc                                                                                                   >  2:512734/1‑149 (MQ=255)
                                            tttttACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAAGCCGCATcc                                                                                                   >  2:963123/1‑149 (MQ=255)
                                             ttttACGAGATGTTGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCAGCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAGCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                  <  1:45128/149‑1 (MQ=255)
                                                 aCGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                              <  1:969388/149‑1 (MQ=255)
                                                  cGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                             >  1:740917/1‑149 (MQ=255)
                                                   gagaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                                                                            >  1:806388/1‑149 (MQ=255)
                                                   gagaTGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCCTGCAGTGATCATCGTCGATATGTAGGCCGGCTAAGACGTTCAAGACGAATCCGACaaga                                                                                            >  2:560244/1‑147 (MQ=255)
                                                        gtgGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                       >  1:908381/1‑149 (MQ=255)
                                                            gCATCGCCTCTGTAAAACTGGGAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                   <  1:787870/149‑1 (MQ=255)
                                                            gCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                   >  2:255477/1‑149 (MQ=255)
                                                             cATCGCCTCTGTAAAATTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                  <  1:315197/149‑1 (MQ=255)
                                                             cATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                  <  1:255477/149‑1 (MQ=255)
                                                              aTCGCCTATGGAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                                 <  2:361137/149‑1 (MQ=255)
                                                                 ggctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGa                                                                              <  2:545390/147‑1 (MQ=255)
                                                                   ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGa                                                                                                                                     >  1:624506/1‑92 (MQ=255)
                                                                   ctctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGa                                                                                                                                     <  2:624506/92‑1 (MQ=255)
                                                                    tctGTAAGACTGGCAGAGATTAGGGGCGTAACCCTACTCGCCGCTGTTTGATTTGTGGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                           <  2:690815/149‑1 (MQ=21)
                                                                    tctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGACCGATGCCTGATGCGAcgc                                                                           <  1:963123/149‑1 (MQ=255)
                                                                    tctGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                           <  1:365176/149‑1 (MQ=255)
                                                                       gTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                        >  1:775143/1‑149 (MQ=255)
                                                                        tAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATCGAAAAGGCGCTAGCAGTGACCATCGTCGATATGTAGGCCGGAAAAAGCGTTCAAGCCGCATCCGGCAACCTATGCCTGATGCGAcgacgtc                                                                       >  2:581240/1‑144 (MQ=37)
                                                                                  agagATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                             >  2:622643/1‑149 (MQ=255)
                                                                                      aTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGCCGATGCCTGATGCGACGCGGTCGGGTGTTATCAGGc                                                         >  1:126660/1‑149 (MQ=255)
                                                                                      aTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                         <  2:775143/149‑1 (MQ=255)
                                                                                      aTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                         <  2:740917/149‑1 (MQ=255)
                                                                                        tAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt                                                       >  1:759078/1‑149 (MQ=255)
                                                                                          gagGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGGCTTATCAGGCTTAc                                                     >  1:884876/1‑149 (MQ=255)
                                                                                               gTAAAGTCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTc                                                <  2:495233/149‑1 (MQ=255)
                                                                                                  aaTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcgtcg                                                                                                                                      <  2:854931/60‑1 (MQ=255)
                                                                                                  aaTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCAtcgtcg                                                                                                                                      >  1:854931/1‑60 (MQ=255)
                                                                                                   aTCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCCTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATg                                            <  1:855533/149‑1 (MQ=255)
                                                                                                    tCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                           >  1:296744/1‑149 (MQ=255)
                                                                                                       ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                        >  1:826353/1‑149 (MQ=255)
                                                                                                       ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCACCGCGGTCGCCTCTTATCAGGCCCACAGGTCGCTGCCGa                                        >  2:680594/1‑149 (MQ=38)
                                                                                                         tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                             >  1:149133/1‑110 (MQ=255)
                                                                                                         tGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                             <  2:149133/110‑1 (MQ=255)
                                                                                                          ggCCACGCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                     <  2:411729/149‑1 (MQ=255)
                                                                                                             cACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGACAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTa                                  >  1:862570/1‑149 (MQ=255)
                                                                                                                cccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                               >  1:967968/1‑149 (MQ=255)
                                                                                                                cccGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                               <  2:826353/149‑1 (MQ=255)
                                                                                                                  cGAGAATAGATGCATTGAAGAGGCGCTTCCAGTGATCATCGTAGATATGTAGGCCGGCTAAGGCGTGCACGCGACATCTGGCAACCGATGCCGGCTTCGACGCAGTCGCGTCTTAAGAAGCCTACATGTCGATAGCGATAGGGACAtag                             >  2:148257/1‑147 (MQ=255)
                                                                                                                   gagaATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGt                            >  2:189633/1‑149 (MQ=255)
                                                                                                                       atcgATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGCCGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                        <  1:962031/146‑1 (MQ=255)
                                                                                                                                ttGAAGAGGCGCTTGCGGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa               <  1:503173/149‑1 (MQ=255)
                                                                                                                                ttGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa               >  1:935865/1‑149 (MQ=255)
                                                                                                                                 ggAAGAGGCGCTTGCACTGCACATCGTCGAGATGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt              <  2:969512/148‑1 (MQ=25)
                                                                                                                                  gAAGCGGCGTTTGCAGTGATCATCGTCGATATGTCGGCCGGCTAAGGCGTTCACGCCGCATCCGGCAACCGATGTCTGTTGCGACGCGGTCGCGTCTTATCAGGCCTACAGCTCGATGCCGATATGTACATCGTATTCGGCAATTACTa             <  2:296744/149‑1 (MQ=25)
                                                                                                                                            ttGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACga   >  1:213409/1‑149 (MQ=255)
                                                                                                                                             tGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgat  <  2:155432/149‑1 (MQ=255)
                                                                                                                                               cAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                           >  2:508498/1‑90 (MQ=255)
                                                                                                                                               cAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                           <  1:508498/90‑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 < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: