Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_pgi 3,724,829 T→C 6.7% T164A (ACC→GCC)  glyQ ← glycine tRNA synthetase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,724,8290TC 6.7% 161.8 / 3.3 63T164A (ACC→GCC) glyQglycine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (31/25);  new base C (2/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

TATTCGAAGTTGTAAGTGGACTGCTCCACTTCGTTCTGATGGAACACGTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACG  >  NC_000913_3_pae_pgi/3724681‑3724976
                                                                                                                                                    |                                                                                                                                                   
tATTCGAACGTGTAAGGGGACCGCTCCCTTTCGTTCTGATGGAACACGTCGCCGTAGGTGGTTTTACCCAGCGTGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGt                                                                                                                                                     <  1:414949/149‑1 (MQ=255)
        gTTGTAAGTGGACTGCTCCACTTCGTTCTGATGGAACACGTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGc                                                                                                                                             <  1:790331/149‑1 (MQ=255)
                  gACTGCTCCCCTTCGTTCGGATGGAACACGTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCgg                                                                                                                                   >  2:122042/1‑149 (MQ=255)
                    cTGCTCCACTTCGTTCTGATGGAACACGTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCGGAATGTACTTGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGtt                                                                                                                                 <  1:645804/149‑1 (MQ=255)
                       cTCCACTTCGTTCTGCTGAAACACGTCGCCGTAGGTGGATTTCCCCAGCGGGCCGTACCTCCACACGAGGGCGAAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGCAACCGGTTTac                                                                                                                              >  2:966442/1‑149 (MQ=255)
                       cTCCACTTCGTTCTGATGGAACACGTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCCGACCGTAGGTGATCTCGCCGGTAACCGGTTTac                                                                                                                              >  2:729656/1‑149 (MQ=255)
                            cTTCGTTCTGATGGAACACGTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTcc                                                                                                                         <  1:648745/149‑1 (MQ=255)
                             ttCGTTCTGATGGAACACATCGCCGTAGGTGGTTCCGTCCTGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAACGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTCCACTCCa                                                                                                                        >  2:925667/1‑149 (MQ=255)
                                gggcgggtGGAACACGCCGCGGTAGGTGGTTGTACCTAGGGGGCCGTCGTTCCAGCCCAGGTCGTAAACGCTGTCTCCGCCCTGCATGTACATGGGCAGACGTTCGAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGac                                                                                                                     <  2:42887/142‑1 (MQ=255)
                                     gATGGAACACGTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCCGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCAGTTTACACTCAAGActacca                                                                                                                >  1:471744/1‑149 (MQ=255)
                                          aacgtgtcgtcGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGGGATCTCGCCTGCACCCGGTTTACACGCCCGACCACCAACctg                                                                                                           >  2:877207/6‑149 (MQ=255)
                                          aaCACGTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACctg                                                                                                           <  1:337576/149‑1 (MQ=255)
                                           acacGTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACctgc                                                                                                          >  1:614045/1‑149 (MQ=255)
                                             gcGTCGCCGCAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCCGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACctgctg                                                                                                        <  1:11991/148‑1 (MQ=255)
                                             acGTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACctgctg                                                                                                        <  2:854511/149‑1 (MQ=255)
                                               gTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGa                                                                                                      >  2:637977/1‑149 (MQ=255)
                                               gTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGa                                                                                                      >  1:575046/1‑149 (MQ=255)
                                                tCGGCGTAGGTGGTTTTACCCGGTGTGCCGTCTATCCAGGCAAGGTCGTAAACGCGGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGaa                                                                                                     <  2:891701/149‑1 (MQ=255)
                                                  gCCGTAGGTGGTTTTACCCATCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGCGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGaagt                                                                                                   >  2:396466/1‑149 (MQ=255)
                                                  gCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGCGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGaagt                                                                                                   >  2:392514/1‑149 (MQ=255)
                                                      tAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCTTGAATGTACATGGAAAGAAACTCCAGACCGAAGAGAAACCCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGaagtaagt                                                                                               >  1:291088/1‑149 (MQ=255)
                                                         gtggtTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGaa                                                                                            >  2:833516/1‑149 (MQ=255)
                                                         gcggtGTTACCCAGCGGGCCGTCGCGCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGaa                                                                                            <  2:153468/147‑1 (MQ=255)
                                                          tggtTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAAACCACCAACCTGCTGGAAGTAAGTGAAc                                                                                           >  2:307614/1‑149 (MQ=255)
                                                             ttttACCCAGCGGGCCGTCGCTCCAGACCAGGTCGGAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGc                                                                                        <  1:374958/149‑1 (MQ=255)
                                                             ttttACCCAGCGGGCCGTCGCTCAAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCGACCTGCTGGTAGTAAGTGAACTGc                                                                                        <  1:870649/149‑1 (MQ=255)
                                                                 aCCCAGCGGGCCGTAGCTCCCGCCCAGTTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAATCGTTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCa                                                                                    <  2:5754/149‑1 (MQ=255)
                                                                    cAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTAACtt                                                                                 >  1:495398/1‑149 (MQ=255)
                                                                     aGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGCACCTGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACATTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTc                                                                                >  1:948441/1‑149 (MQ=255)
                                                                       ccggCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTCCGTCATTTCCa                                                                              <  2:569903/147‑1 (MQ=255)
                                                                        gggCCGTCGCTTCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCAt                                                                             <  1:390118/149‑1 (MQ=255)
                                                                                 cTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTTCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAg                                                                    <  2:354039/149‑1 (MQ=255)
                                                                                  tCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGACCAGGCGGACCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCACCCTGCTGGAAGTAAGTGAACTGCGCCACTTCCATGCCGTTCAGc                                                                   >  2:508377/1‑149 (MQ=255)
                                                                                   ccAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGCGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCACCCGGCTGGACGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGcc                                                                  <  1:396466/149‑1 (MQ=255)
                                                                                   ccAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGCGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGcc                                                                  <  1:392514/149‑1 (MQ=255)
                                                                                    cAGACCAGGTCGTAAACGCTGTCTACGCCCTGAACGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGATTTACACTCCAGACCACCAACCTGCTGGAAGAAAGTGAACTGCCTCACCTCCCCGCCCTTCAcccc                                                                 >  2:985529/1‑148 (MQ=255)
                                                                                         cAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCCTGCCGTTCAGCCCCCCtt                                                            >  2:353179/1‑149 (MQ=255)
                                                                                             tCGTAAACGCTGTCTACGCCCTGAATGTACAGGGCCAGACGTTCCAGACCGTAGGTGAACTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTccca                                                        <  1:725269/149‑1 (MQ=255)
                                                                                              ggTAAATGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTGCCATGCCGTTCAGCCACACTTcccag                                                       <  1:157469/148‑1 (MQ=255)
                                                                                                aaaaCGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGCTTACACTCCAGACTACCAACCTGCTGGAAGTAAGTGAACCGCGTCACTTCCATGCCGTTCTGCCACACTTcccagcc                                                     >  2:928907/2‑149 (MQ=255)
                                                                                                  aaCGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTcccagccca                                                   <  1:637977/149‑1 (MQ=255)
                                                                                                  aaCGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTcccagccca                                                   >  2:587049/1‑149 (MQ=255)
                                                                                                    cGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACctgc                                                                                                          >  2:256555/1‑92 (MQ=255)
                                                                                                    cGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACctgc                                                                                                          <  1:256555/92‑1 (MQ=255)
                                                                                                              gcgcTGAATGTACAGAGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCCGGAAGTACGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCAc                                       <  1:575274/146‑1 (MQ=255)
                                                                                                                ccTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCAccc                                     >  1:431371/1‑149 (MQ=255)
                                                                                                                ccTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTACATGCCGTTCAGCCACACTTCCCAGCCCAGTCACCAGGCAccc                                     >  2:735945/1‑149 (MQ=255)
                                                                                                                ccTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTCCGTCACTTCCATGCCGTTCATCCACACTTCCCAACCCAATCACAAGGCAccc                                     >  1:473991/1‑149 (MQ=255)
                                                                                                                  tGAATGTACATGGCCAGACGTTCCAGACGGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAGCTGCGTCACTTCCATGCCGTTCAACAACACTACCCAGCGCACTACCCAGGCACCCaa                                   >  1:445304/1‑148 (MQ=255)
                                                                                                                        tACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCgg                             <  2:291088/149‑1 (MQ=255)
                                                                                                                        tACATGGCCAGACGTTCCAGACCGGAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGTAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCgg                             <  1:470895/149‑1 (MQ=255)
                                                                                                                        tACATGGCAAGACTTGCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCgg                             <  1:125237/149‑1 (MQ=255)
                                                                                                                         aCATGGCCAGAAGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCggg                            >  1:615296/1‑149 (MQ=255)
                                                                                                                           aTGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCGGGtt                          >  1:409376/1‑149 (MQ=255)
                                                                                                                                cAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCGGGTTTTccc                     >  1:382730/1‑149 (MQ=255)
                                                                                                                                       tCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCGGGTTTTCCCAGTTATc              >  1:968266/1‑149 (MQ=255)
                                                                                                                                        ccAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCGGGTTTTCCCAGTTATCt             >  1:807482/1‑149 (MQ=255)
                                                                                                                                         cAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCAAACTTCCCAGCCCATTCCCCAGGCACCCAGCGGCGGGTTTTCCCAGTTATCtc            >  2:434600/1‑148 (MQ=255)
                                                                                                                                         cAGACCGTAGGATAATTCGCCGGTAACCGGTTTACACTCCAGACCACCCACCTGCTTGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCGGGTTTTCCCAGTTATCtt            <  1:131112/149‑1 (MQ=255)
                                                                                                                                          aGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCTCAGTCCCCATGCACCCAGCGTCGGGTTTTCCCATTTATCTTc           >  1:24242/1‑149 (MQ=255)
                                                                                                                                          aGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAACGACGGGTTTTCCCAGATATCTTc           >  1:90701/1‑149 (MQ=255)
                                                                                                                                           gACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCGGGTTTTCCCAGTTATCTTcc          <  2:431371/149‑1 (MQ=255)
                                                                                                                                             ccGTAGGTGATCTCGCCGGTAACCGGTTTACTCTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCAc        <  1:762203/149‑1 (MQ=255)
                                                                                                                                                   gTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACg  >  2:522759/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TATTCGAAGTTGTAAGTGGACTGCTCCACTTCGTTCTGATGGAACACGTCGCCGTAGGTGGTTTTACCCAGCGGGCCGTCGCTCCAGACCAGGTCGTAAACGCTGTCTACGCCCTGAATGTACATGGCCAGACGTTCCAGACCGTAGGTGATCTCGCCGGTAACCGGTTTACACTCCAGACCACCAACCTGCTGGAAGTAAGTGAACTGCGTCACTTCCATGCCGTTCAGCCACACTTCCCAGCCCAGTCCCCAGGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACG  >  NC_000913_3_pae_pgi/3724681‑3724976

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