Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,354,297 G→T intergenic (‑247/‑428) yhcC ← / → gltB putative Fe‑S oxidoreductase, Radical SAM superfamily protein/glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,354,2970GT99.0% 354.3 / ‑4.3 99intergenic (‑247/‑428)yhcC/gltBputative Fe‑S oxidoreductase, Radical SAM superfamily protein/glutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (51/47);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAA  >  NC_000913_3_pae_pgi/3354152‑3354434
                                                                                                                                                 |                                                                                                                                         
cATAATCAATCAAGATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTgg                                                                                                                                       >  3:144942/1‑150 (MQ=255)
cATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTgg                                                                                                                                       >  4:540712/1‑150 (MQ=255)
cATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGAATTTTAGGCATAAACAGTCGCATGAAACGTCATTACCAATGAAGGCAGTATAAAATTCTgg                                                                                                                                       >  3:144933/1‑150 (MQ=255)
  tAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGtt                                                                                                                                     <  1:748185/150‑1 (MQ=255)
    atcaatcaAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGtttt                                                                                                                                   <  2:596990/150‑1 (MQ=255)
    atcaatcaAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGtttt                                                                                                                                   <  4:291414/150‑1 (MQ=255)
         tcaaAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTgtcgt                                                                                                                              >  1:1113258/1‑150 (MQ=255)
         tcaaAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTgtcgt                                                                                                                              <  2:1004067/150‑1 (MQ=255)
          caaAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTgtcgtc                                                                                                                             <  3:374571/150‑1 (MQ=255)
                 aCCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCaa                                                                                                                      <  3:48363/150‑1 (MQ=255)
                   cGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAgg                                                                                                                    >  4:678407/1‑150 (MQ=255)
                   cGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAgg                                                                                                                    >  3:209853/1‑150 (MQ=255)
                     aaaTTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCa                                                                                                                  >  3:281550/1‑150 (MQ=255)
                     aaaTTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCa                                                                                                                  >  1:324354/1‑150 (MQ=255)
                      aaTTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAg                                                                                                                 >  4:479754/1‑150 (MQ=255)
                      aaTTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAg                                                                                                                 <  2:274017/150‑1 (MQ=255)
                       aTTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAgg                                                                                                                <  3:700816/150‑1 (MQ=255)
                        tttCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGa                                                                                                               <  1:878335/150‑1 (MQ=255)
                           cATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATaa                                                                                                            <  1:817660/150‑1 (MQ=255)
                             tGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAgg                                                                                                          <  4:261887/150‑1 (MQ=255)
                             tGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAgg                                                                                                          <  4:281550/150‑1 (MQ=255)
                                aTAATCATATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGtt                                                                                                       <  4:296271/150‑1 (MQ=255)
                                 tAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTa                                                                                                      <  1:1119990/150‑1 (MQ=255)
                                      acaTAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTgtcgt                                                                                                                              >  1:769820/1‑121 (MQ=255)
                                      acaTAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTgtcgt                                                                                                                              <  2:769820/121‑1 (MQ=255)
                                      acaTAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAacaca                                                                                                 >  4:414142/1‑150 (MQ=255)
                                       caTAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAacacac                                                                                                >  4:649326/1‑150 (MQ=255)
                                        aTAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACAcc                                                                                               >  1:1061935/1‑150 (MQ=255)
                                           aaTCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCttt                                                                                            <  3:678407/150‑1 (MQ=255)
                                             tCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACTTTTAt                                                                                          >  2:154296/1‑150 (MQ=255)
                                              cACTTTTGATTATCGTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACCCCTTTATg                                                                                         <  2:651890/150‑1 (MQ=255)
                                               aCTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGa                                                                                        >  3:93598/1‑150 (MQ=255)
                                               aCTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGa                                                                                        <  2:1061935/150‑1 (MQ=255)
                                                 ttttGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACa                                                                                      >  3:613375/1‑150 (MQ=255)
                                                 ttttGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACa                                                                                      >  3:613535/1‑150 (MQ=255)
                                                      cTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAg                                                                                 >  3:184661/1‑150 (MQ=255)
                                                        tATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGa                                                                               >  4:64224/1‑150 (MQ=255)
                                                        tATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTCTATGAAAGTCAGGa                                                                               >  4:64197/1‑150 (MQ=255)
                                                           cTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAAtt                                                                            <  1:55429/150‑1 (MQ=255)
                                                             tgtgTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCTTGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGa                                                                          <  2:593448/150‑1 (MQ=255)
                                                             tgtgTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGa                                                                          >  3:25720/1‑150 (MQ=255)
                                                              gtgtCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGAc                                                                         <  2:915772/150‑1 (MQ=255)
                                                               tgtCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAACCAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACt                                                                        <  2:177206/150‑1 (MQ=255)
                                                                 tCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGt                                                                      <  1:817067/150‑1 (MQ=255)
                                                                 tCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGt                                                                      >  3:152320/1‑150 (MQ=255)
                                                                   aGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCACTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGCATTGACTGttt                                                                    >  4:469584/1‑150 (MQ=255)
                                                                    gATTTTTTTATTTCCTGATGGAATCTAGGCAAAAACAGTAGCATGATACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAAGCAGGATAAGGGGTAACACACCTGTATTACAGGCAGGAATTGACTGTTtc                                                                   >  2:674216/1‑150 (MQ=255)
                                                                    gATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTtc                                                                   >  1:426570/1‑150 (MQ=255)
                                                                     aTTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTTTTTTTTTTTTGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTtct                                                                  >  4:151651/1‑150 (MQ=255)
                                                                       ttttttATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTctct                                                                >  1:783450/1‑150 (MQ=255)
                                                                          tttATCTCCTGATGGATTTTAGGCAAAAACAGGAGCATGAATCGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAAc                                                             >  2:646072/1‑150 (MQ=255)
                                                                            tATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGa                                                           <  4:664377/150‑1 (MQ=255)
                                                                             aTCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGAc                                                          <  4:536890/150‑1 (MQ=255)
                                                                             aTCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGAc                                                          <  4:536893/150‑1 (MQ=255)
                                                                                tcCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTc                                                       <  3:414142/150‑1 (MQ=255)
                                                                                 ccTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTcc                                                      <  1:541540/150‑1 (MQ=255)
                                                                                 ccTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTcc                                                      >  2:541540/1‑150 (MQ=255)
                                                                                     aTGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCttt                                                  >  4:465969/1‑150 (MQ=255)
                                                                                     aTGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCttt                                                  >  2:1104316/1‑150 (MQ=255)
                                                                                     aTGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCttt                                                  >  1:513000/1‑150 (MQ=255)
                                                                                     aTGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCttt                                                  >  4:234818/1‑150 (MQ=255)
                                                                                       ggATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTa                                                >  2:795050/1‑150 (MQ=255)
                                                                                          ttttAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGcc                                             <  3:573874/150‑1 (MQ=255)
                                                                                           tttAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTaa                                                                                                     <  2:761347/93‑1 (MQ=255)
                                                                                           tttAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTaa                                                                                                     >  1:761347/1‑93 (MQ=255)
                                                                                              aGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTaa                                         >  4:236391/1‑150 (MQ=255)
                                                                                               ggCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTaaa                                        >  1:246982/1‑150 (MQ=255)
                                                                                               ggCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTaaa                                        <  1:62138/150‑1 (MQ=255)
                                                                                               ggCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTaaa                                        <  1:472369/150‑1 (MQ=255)
                                                                                                  aaaaaCAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGAt                                     >  3:694655/1‑150 (MQ=255)
                                                                                                   aaaaCAGTAGCATGAAGCGTCATTACCAATGAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTTAAGGGAGGATAAGGGGTAACACACCTTTATGACAGTTAGGAATTGACTGTTTCTCTAACGTCTTCCCTTTTAGCCTTAAAGATa                                    >  2:613670/1‑150 (MQ=255)
                                                                                                         gtatgaTGAAACGTCATTAACAGTTAAGCCAGTATAAAATGCTGATTTTGTCGGCAGTAAAAGGCTGGATAAAGGTTAACGCGCCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTTTAGATCAAATcc                              <  2:987869/145‑1 (MQ=255)
                                                                                                         gTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATcc                              <  1:1012775/150‑1 (MQ=255)
                                                                                                              aTGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTgt                                                                                                                                 >  1:43877/1‑46 (MQ=255)
                                                                                                              aTGAAACGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTgt                                                                                                                                 <  2:43877/46‑1 (MQ=255)
                                                                                                                   aCGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAAttt                    <  2:783450/150‑1 (MQ=255)
                                                                                                                   aCGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAAttt                    >  3:66398/1‑150 (MQ=255)
                                                                                                                   aCGTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAAttt                    <  3:457845/150‑1 (MQ=255)
                                                                                                                     gTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCa                  <  3:465969/150‑1 (MQ=255)
                                                                                                                     gTCATTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCa                  <  2:1162152/150‑1 (MQ=255)
                                                                                                                        aTTACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTc               >  2:392845/1‑150 (MQ=255)
                                                                                                                         ttACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTaa                                                                                                     >  3:461877/1‑63 (MQ=255)
                                                                                                                         ttACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTaa                                                                                                     <  4:461877/63‑1 (MQ=255)
                                                                                                                          tACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTaaa                                        <  1:388279/123‑1 (MQ=255)
                                                                                                                          tACCAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTaaa                                        >  2:388279/1‑123 (MQ=255)
                                                                                                                            ccAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCAttt           >  3:130834/1‑150 (MQ=255)
                                                                                                                             cAATTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTta          >  3:242685/1‑150 (MQ=255)
                                                                                                                              aaTTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGa                                                                                        >  4:37730/1‑71 (MQ=255)
                                                                                                                              aaTTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGa                                                                                        <  3:37730/71‑1 (MQ=255)
                                                                                                                               aTTAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTtaat        <  3:479754/150‑1 (MQ=255)
                                                                                                                                ttAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGt                >  3:702397/1‑141 (MQ=255)
                                                                                                                                ttAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGt                <  4:702397/141‑1 (MQ=255)
                                                                                                                                ttAAGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTtaata       >  4:452175/1‑150 (MQ=255)
                                                                                                                                  aaGGCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATaaa     <  1:795050/150‑1 (MQ=255)
                                                                                                                                    ggCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGa   >  1:275422/1‑150 (MQ=255)
                                                                                                                                     gCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGaa  <  4:613375/150‑1 (MQ=255)
                                                                                                                                     gCAGTATAAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGaa  <  4:613535/150‑1 (MQ=255)
                                                                                                                                           taAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGACAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTtct                                                                  <  1:958702/80‑1 (MQ=255)
                                                                                                                                           taAAATTCTGGTTTTGTCGTCAGTTCAAGGCAGGACAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTtct                                                                  >  2:958702/1‑80 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
CATAATCAATCAAAATTACCGAAATTTCATGCATAATCACATAAATCACTTTTGCTTATCTTGTGTCAGATTTTTTTATCTCCTGATGGATTTTAGGCAAAAACAGTAGCATGAAACGTCATTACCAATTAAGGCAGTATAAAATGCTGGTTTTGTCGTCAGTTCAAGGCAGGATAAGGGTTAACACACCTTTATGACAGTCAGGAATTGACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAA  >  NC_000913_3_pae_pgi/3354152‑3354434

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