Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_pgi 3,491,704 G→T 7.2% intergenic (‑84/+21) yhfG ← / ← ppiA putative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,491,7040GT 7.2% 174.9 / 3.7 70intergenic (‑84/+21)yhfG/ppiAputative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base G (27/37);  new base T (3/2);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

GAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAA  >  NC_000913_3_pae_pgi/3491557‑3491852
                                                                                                                                                   |                                                                                                                                                    
gAAAATTACGATTACGCTGAAGCTCCCAGAGACCGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGa                                                                                                                                                     <  1:525814/149‑1 (MQ=255)
      tacgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGccc                                                                                                                                               <  1:882486/149‑1 (MQ=255)
      tacgattacgCGGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGccc                                                                                                                                               <  2:896523/149‑1 (MQ=255)
                      cTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGAAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAtca                                                                                                                               >  2:348548/1‑147 (MQ=255)
                       tCCCAGAGACGGCTCTTTTGCTTATCGGTAAGTTTCGTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGGAATCAttcc                                                                                                                              >  2:344628/1‑147 (MQ=255)
                         ccAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACgg                                                                                                                            >  1:515596/1‑149 (MQ=255)
                                      ttttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCg                                                                                                               >  2:430119/1‑149 (MQ=255)
                                      ttttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCg                                                                                                               >  2:623685/1‑149 (MQ=255)
                                      tttgtcatATCGGGGTGTTTCTTCATTTTATGTTTCCTGTTACTTGTCCATTTGTCACAAGTATATTTGTCTAATCCGGGTTACGTTGGGATAGCAGAAGCGCGTGTAAGATTGCCCGTGCGAGAAATCATTACGGCAGGACTTTAGCg                                                                                                               <  2:28828/142‑1 (MQ=255)
                                         tGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGaa                                                                                                            <  1:404147/149‑1 (MQ=255)
                                         tGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGAAAGATTGCCCGCGCGAGAAATCATTACCGCAGGACTTTAGCGGaa                                                                                                            >  1:931298/1‑149 (MQ=255)
                                          gtttATCGGTTAGTTTCTTCACTTGATGCCTCCTTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGaaa                                                                                                           <  1:286723/147‑1 (MQ=255)
                                             tATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGCGAGCAGAAGCGCGGGCAATCTTGCCcg                                                                                                                                              >  2:28509/1‑111 (MQ=255)
                                             tATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcg                                                                                                                                              <  1:28509/111‑1 (MQ=255)
                                              atgggTAAGTTTCTTCACTTGATGCCTCCGTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGAt                                                                                                       <  2:100759/146‑1 (MQ=255)
                                                  gTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACt                                                                                                   >  1:800663/1‑149 (MQ=255)
                                                   tGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTa                                                                                                  >  2:394495/1‑149 (MQ=255)
                                                    gAGTTTCTTCACTTGATGCCTCCCTGAACGTTTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGAAGAAACGCGGTCAAGGTTGAAAGCGCAAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTAc                                                                                                 <  1:740721/149‑1 (MQ=255)
                                                     aGTTTCTTCACTTGATGCCTCCCTGAACGTGTCGATTTGACACAATTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTAcc                                                                                                <  1:940209/149‑1 (MQ=255)
                                                           ttcACTTGATGCCTCCCTGAACGTGTCCATTTTCCACAAGTATAAGCGGCAAATCCGTGTTACACCTGTAAATCATAATCGCAGGCAAGATTGCCCTCGCTAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGtttt                                                                                          <  1:287037/149‑1 (MQ=255)
                                                             cACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAGGCGGCAAATCGGGGTTAAGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCGGCGCGAGGAATCATTAAGGCAGGCCTTTAGCGGAAAGGATAACTACGGGTTTTGa                                                                                        >  1:296067/1‑149 (MQ=255)
                                                              aCTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGGAGGGATAACTACCGGATTTGAc                                                                                       >  2:792896/1‑149 (MQ=255)
                                                                ttGATGCCTCCCTGAACGTGTCCATTTGCCGCAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACgg                                                                                     >  1:756294/1‑149 (MQ=255)
                                                                  gATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGca                                                                                   <  1:514949/149‑1 (MQ=255)
                                                                   aTGCCTCCCGGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGTAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGcac                                                                                  <  1:832719/149‑1 (MQ=255)
                                                                     gCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAAGATAACTACCGGTTTTGACGGCACAt                                                                                >  1:294380/1‑149 (MQ=255)
                                                                       cTCCCTGAAGGTGTCCATTTGCCACAAGTAGAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTc                                                                              <  2:831789/149‑1 (MQ=255)
                                                                          tctgactgtgtACACTTGCCACAAGTATAGGCGTCAACTCCGGGTTACGCCGGGAGAGTAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGTCTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTgg                                                                           <  2:41112/142‑1 (MQ=255)
                                                                          ccTGAACGTGTTCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTa                                                                                                                               <  1:896217/97‑1 (MQ=255)
                                                                          ccTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTa                                                                                                                               >  2:896217/1‑97 (MQ=255)
                                                                            tGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTa                                                                         >  1:326945/1‑149 (MQ=255)
                                                                            tGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGACTGCCCGCGCGAGAAATCATTACGGCAGGACTATAGCAGAAAGGATAACTACAGGCTTTGACGGCAAATTCTGGTa                                                                         >  1:327830/1‑149 (MQ=255)
                                                                              aaCGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTCGCGGAAAGGATAACTACCGGTTTCGACAGCAGACTCTGGTGCg                                                                       >  2:551370/1‑149 (MQ=255)
                                                                                 gtgtCCATTTGCCACAAGCATAAGCGGCAATTCTGGGTTACGCTGGGTGAGCATAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAGTACGGCAGGACTTTAGCGGAAAGGAGAACTACCGGTTTTGACGGCACATTCTGGTACGGAc                                                                    <  1:5450/149‑1 (MQ=255)
                                                                                   gtCCATTTGCCACAAGTATAAGAGGATCATCCGGGTTACGCCGGTAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTCGTACGGACCa                                                                  >  2:868724/1‑149 (MQ=255)
                                                                                    tCCATTTACCACAAGTATAAGCGGAAAATGGGGGTTTCGGCGGTTGATGAAGAGCGCTGTCAATATTGCTGGGGCGAGGAATCATGACGGCAGGACGTTAGCGTAAAGGATAACTACCGGTTTTGACGGCACATGCTGGTACGGACaaa                                                                 >  2:593114/1‑149 (MQ=255)
                                                                                        tttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTc                                                             >  1:246115/1‑149 (MQ=255)
                                                                                         ttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAGGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCAAATTCTTATACGGACCAAAGtaa                                                            >  1:326118/1‑147 (MQ=255)
                                                                                         ttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCa                                                            <  2:61104/149‑1 (MQ=255)
                                                                                         ttGCCACAAGGATAAGCGGAAAACCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCa                                                            <  2:551398/149‑1 (MQ=255)
                                                                                          tGCGACAATTATAAGCGTGAAATCGAGGTTGCGCTGGGAGAGCAGAGGCGCGGGCAAGAGTGAACGCGAGAGAAATCAGTACGGGAGGTCGTTAGGGGAAAGGATAGCTATCGGTTTTGACGGCACATTCTGGTACGGGCCAACGTCAt                                                           <  2:237732/149‑1 (MQ=255)
                                                                                            ccACTAGTATGAGCTGCAATTCCTGGTTACGCCGGGTGAGCAGAAGGGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGa                                                         <  1:215652/149‑1 (MQ=255)
                                                                                             cacaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAg                                                        >  1:639256/1‑149 (MQ=255)
                                                                                             cacaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAg                                                        <  1:430027/149‑1 (MQ=255)
                                                                                              acaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCAAATTCTGGTAAGGACCAACGTCATGAGt                                                       >  1:489239/1‑149 (MQ=255)
                                                                                               caAGTATAAGCGGCAAATTCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACATTGCCCGCGCGAGAAAGCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTc                                                      <  1:925348/149‑1 (MQ=255)
                                                                                                aaGTATAAGCGGCAAATCCGGGTTACGCCGGGGGAGCAGAAGCGCGGGCAAGAGTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCg                                                     <  2:914354/149‑1 (MQ=255)
                                                                                                           gCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGGAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGAACCTGGGaa                                          >  2:867055/1‑149 (MQ=255)
                                                                                                               aTCCGCGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGACCGCGCGAGAAATCATTACGGCAGGCCTTTCGCGGAAAGGATAACTACAGCCTCTAACCGCACATTCTGGTAAGCAACAACGTCAGGAATCGGAACGTGGGacacat                                      >  2:231320/1‑144 (MQ=255)
                                                                                                                tCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCtt                                     >  1:889470/1‑149 (MQ=255)
                                                                                                                tCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTCAAGAGAAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGCCGGCCCCCGGGAAATCtt                                     >  2:782232/1‑149 (MQ=255)
                                                                                                                  cGGGTTACGCCGGGAGGGCAGAAGCGCGGGCAAGATTGCCCGGGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTAt                                   <  2:12863/149‑1 (MQ=255)
                                                                                                                  cGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTAt                                   <  1:623685/149‑1 (MQ=255)
                                                                                                                     gTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                                <  1:344334/149‑1 (MQ=255)
                                                                                                                             gggAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACGGGGGGTATCTTATCGGCAAGGTcc                        <  2:974734/149‑1 (MQ=255)
                                                                                                                              ggAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATGACTCGCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGGCCCTGGGAAGTCCTGTCGGCGACGTcga                       >  2:2883/1‑147 (MQ=255)
                                                                                                                              ggAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCa                       <  2:852232/149‑1 (MQ=255)
                                                                                                                               gagaGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGACGGCACCTCGGAAATCTTATCGGCCACGTCCAt                      >  2:170820/1‑149 (MQ=255)
                                                                                                                               gagaGCAGAAGCGCGGGCAAGAGTGCCCGCGCGAGAAATCATTACGGCAGGACTCTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCAt                      <  1:609126/149‑1 (MQ=255)
                                                                                                                                aaagCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATg                     <  1:344628/147‑1 (MQ=255)
                                                                                                                                 gagCAGAAGCGCGGGCAAGATTGCCCTCGCGAGAAATCATTACGGCAGGCCTTTAGTGGAAAGGATAACGACCGGTTTTGACGGCACATTCTTGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGc                    <  2:57571/149‑1 (MQ=255)
                                                                                                                                   gCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCt                  <  1:430119/149‑1 (MQ=255)
                                                                                                                                    cAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCtt                 <  1:565588/149‑1 (MQ=255)
                                                                                                                                      gAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTc               <  1:574999/149‑1 (MQ=255)
                                                                                                                                       tcgcgcGGGCAATCTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGca                                                                                   >  1:790695/3‑80 (MQ=38)
                                                                                                                                       tcgcgcGGGCAATCTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGca                                                                                   <  2:790695/78‑1 (MQ=38)
                                                                                                                                          cgcgGGCAAGATTGGTCGCGCGAGAAATCAGTACGGCAGGACTGTGGCGGAAAGGATAACTACCGGTTTTGGCGGCTGATTCTGGTACGCACCAACGTCAGGAGTCGGCACCTGGGAAATCTTATCGGTAACGTCCATGCCTTTcacca           <  2:406990/149‑1 (MQ=255)
                                                                                                                                           gcgGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcaccac          <  2:630274/149‑1 (MQ=255)
                                                                                                                                               gCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa      <  2:246115/149‑1 (MQ=255)
                                                                                                                                                  aGATTGCCCGCGCGAGAAATTATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCa   <  1:509264/149‑1 (MQ=255)
                                                                                                                                                   gATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTATCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTACATGCCTTTAACCACTTTACaaa  >  1:234341/1‑149 (MQ=255)
                                                                                                                                                   gATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAAACTTATCGGCAACGTCCATGCCTTTCACCACTTTACCaa  >  1:123056/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAAAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAA  >  NC_000913_3_pae_pgi/3491557‑3491852

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