Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 1,399,449 G→T S26R (AGC→AGA fnr ← oxygen‑sensing anaerobic growth regulon transcriptional regulator FNR; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,399,4490GT97.8% 326.1 / ‑0.8 92S26R (AGC→AGAfnroxygen‑sensing anaerobic growth regulon transcriptional regulator FNR; autorepressor
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (41/49);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

TACCGGAGCGGATGGCATAAAGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGGCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCAAGTAAGTCTAAG  >  NC_000913_3_pae_pgi/1399304‑1399594
                                                                                                                                                 |                                                                                                                                                 
tACCGGAGCGGATGGCATAAAGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGa                                                                                                                                               <  2:406950/150‑1 (MQ=255)
tACCGGAGCGGATGGCATAAAGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGa                                                                                                                                               >  4:390227/1‑150 (MQ=255)
    ggAGCGGATGGCATAAAGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGAtgct                                                                                                                                           <  4:457398/150‑1 (MQ=255)
      aGCGGATGGCATAAAGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGAtgctgc                                                                                                                                         <  3:623757/150‑1 (MQ=255)
          gATGGCATAAAGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATc                                                                                                                                     >  1:99294/1‑150 (MQ=255)
            tGGCATAAAGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCt                                                                                                                                   >  2:579981/1‑150 (MQ=255)
              gCATAAAGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTgg                                                                                                                                 <  1:610339/150‑1 (MQ=255)
               cATAAAGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTATGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGc                                                                                                                                <  3:70678/150‑1 (MQ=255)
                aTAAAGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCa                                                                                                                               >  1:159284/1‑150 (MQ=255)
                  aaaGAGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGTCTTTCTGAATAGGCTTCTTCCGCTCAATGGTATTATCAAGCTGATCAAGCGCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAAGCCTGGCAAt                                                                                                                             <  4:90904/150‑1 (MQ=255)
                   aaGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATg                                                                                                                            <  1:1032270/150‑1 (MQ=255)
                   aaGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATg                                                                                                                            >  2:164693/1‑150 (MQ=255)
                    aGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATgg                                                                                                                           >  4:464295/1‑150 (MQ=255)
                     gCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGa                                                                                                                          <  1:736363/150‑1 (MQ=255)
                       gATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATa                                                                                                                        <  4:635275/150‑1 (MQ=255)
                        aTTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAg                                                                                                                       <  1:367735/150‑1 (MQ=255)
                         tttAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGc                                                                                                                      >  1:1065339/1‑150 (MQ=255)
                            aaGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGcaca                                                                                                                   >  4:723717/1‑150 (MQ=255)
                               ttCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATggg                                                                                                                          >  1:1074012/1‑139 (MQ=255)
                               ttCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATggg                                                                                                                          <  2:1074012/140‑2 (MQ=255)
                               ttCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAAcc                                                                                                                <  2:321991/150‑1 (MQ=255)
                               ttCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAAcc                                                                                                                >  1:86629/1‑150 (MQ=255)
                               ttCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAAcc                                                                                                                >  1:1034656/1‑150 (MQ=255)
                               ttCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAAcc                                                                                                                <  1:1182042/150‑1 (MQ=255)
                                  atcaCCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGcc                                                                                                             >  3:657712/1‑150 (MQ=255)
                                     aCCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGAtgctg                                                                                                                                          <  2:144782/118‑1 (MQ=255)
                                     aCCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGAtgctg                                                                                                                                          >  1:144782/1‑118 (MQ=255)
                                     aCCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGa                                                                                                          >  1:34708/1‑150 (MQ=255)
                                     aCCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGa                                                                                                          <  2:159284/150‑1 (MQ=255)
                                     aCCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGa                                                                                                          >  1:78887/1‑150 (MQ=255)
                                       cAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACt                                                                                                        >  2:313330/1‑150 (MQ=255)
                                           cTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAAt                                                                                                    <  4:18006/150‑1 (MQ=255)
                                               aaCAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATgcgc                                                                                                >  1:23109/1‑150 (MQ=255)
                                                    cGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGtat                                                                                           <  2:34708/150‑1 (MQ=255)
                                                     gTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGtata                                                                                          >  3:623948/1‑150 (MQ=255)
                                                      tCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATaa                                                                                         <  3:723717/150‑1 (MQ=255)
                                                       cTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAAt                                                                                        <  3:464295/150‑1 (MQ=255)
                                                         ggCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTc                                                                                      >  3:461330/1‑150 (MQ=255)
                                                          gCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCg                                                                                     <  2:23122/150‑1 (MQ=255)
                                                              ttCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCttt                                                                                 >  1:731074/1‑150 (MQ=255)
                                                                cTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTc                                                                               <  2:1167728/150‑1 (MQ=255)
                                                                  gAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCg                                                                             <  2:1065339/150‑1 (MQ=255)
                                                                   aaTAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCg                                                                                     >  4:25503/1‑141 (MQ=255)
                                                                   aaTAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCg                                                                                     <  3:25503/141‑1 (MQ=255)
                                                                   aaTAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCgg                                                                            >  4:510306/1‑150 (MQ=255)
                                                                   aaTAGGCTTCTTCCGCTCAATGATATTATCAAGCCGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCgg                                                                            <  1:862627/150‑1 (MQ=255)
                                                                          ttcttcCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATa                                                                     <  3:390227/150‑1 (MQ=255)
                                                                          ttcttcCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATa                                                                     <  3:79836/150‑1 (MQ=255)
                                                                             ttcCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATgcgc                                                                                                <  4:452190/120‑1 (MQ=255)
                                                                             ttcCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATgcgc                                                                                                >  3:452190/1‑120 (MQ=255)
                                                                                 gCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGc                                                              <  4:623948/150‑1 (MQ=255)
                                                                                  cTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACaa                                                                                                                  >  1:444554/1‑97 (MQ=255)
                                                                                  cTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACaa                                                                                                                  <  2:444554/97‑1 (MQ=255)
                                                                                     aaTGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAAt                                                                                                    <  1:426276/108‑1 (MQ=255)
                                                                                     aaTGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAAt                                                                                                    >  2:426276/1‑108 (MQ=255)
                                                                                     aaTGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTc                                                            >  4:368350/1‑148 (MQ=255)
                                                                                     aaTGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTc                                                            <  3:368350/148‑1 (MQ=255)
                                                                                      aTGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAg                                                         >  1:824634/1‑150 (MQ=255)
                                                                                        gATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGcc                                                       <  1:304748/150‑1 (MQ=255)
                                                                                           attatCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTa                                                    <  2:99294/150‑1 (MQ=255)
                                                                                            ttatCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTaa                                                   <  3:401276/150‑1 (MQ=255)
                                                                                              atCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAAtt                                                 <  1:86999/150‑1 (MQ=255)
                                                                                                cAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGa                                               >  3:865/1‑150 (MQ=255)
                                                                                                 aaGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGAAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGat                                              >  4:469653/1‑150 (MQ=255)
                                                                                                  aGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGAtgctgct                                                                                                                                        <  3:92684/59‑2 (MQ=255)
                                                                                                  aGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGAtgctgct                                                                                                                                        >  4:92684/1‑58 (MQ=255)
                                                                                                   gCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGatat                                            <  1:1066671/150‑1 (MQ=255)
                                                                                                     tGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATAttt                                          <  3:350126/150‑1 (MQ=255)
                                                                                                     tGATCAAGCTCATGTTCGTTGAGTGTGAACCGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATAttt                                          <  1:765751/150‑1 (MQ=255)
                                                                                                            gCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCaa                                                          <  2:422289/127‑1 (MQ=255)
                                                                                                            gCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCaa                                                          >  1:422289/1‑127 (MQ=255)
                                                                                                              tCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAAtttt                                 >  1:467693/1‑150 (MQ=255)
                                                                                                              tCATGTTCGTTGAGTGTGAACGGGATGCAAAGATGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTACTTGATATTTGTCAAtttt                                 >  4:290746/1‑150 (MQ=255)
                                                                                                               cATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTa                                <  2:731074/150‑1 (MQ=255)
                                                                                                                aTGTTCGTTGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTaa                               <  2:23109/150‑1 (MQ=255)
                                                                                                                aTGTTCGTTGAGTGTGAACGGGATGCAAAGCTGGCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCg                                                                             <  1:200523/104‑1 (MQ=255)
                                                                                                                aTGTTCGTTGAGTGTGAACGGGATGCAAAGCTGGCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCg                                                                             >  2:200523/1‑104 (MQ=255)
                                                                                                                       ttGAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCttt                        >  1:669916/1‑150 (MQ=255)
                                                                                                                         gAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCttttt                      <  1:1133860/150‑1 (MQ=255)
                                                                                                                         gAGTGTGAACGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCttttt                      >  3:553473/1‑150 (MQ=255)
                                                                                                                                 aCGGGATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCa              >  2:1085744/1‑150 (MQ=255)
                                                                                                                                    ggATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATa                                                                     <  2:683364/92‑1 (MQ=255)
                                                                                                                                    ggATGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATa                                                                     >  1:683364/1‑92 (MQ=255)
                                                                                                                                      aTGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCaagtaa         <  1:313330/150‑1 (MQ=255)
                                                                                                                                      aTGCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCaagtaa         <  2:467693/150‑1 (MQ=255)
                                                                                                                                        gCAAAGCTGTCTGATGGTGTATTCCTGGCAATGGATAGCACAACCGCCAGACTGAATGGGCAGTATAATTCGCTTTTCCGGGATGATAGGTCTGCTCAAGCGGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCaagtaact       <  4:644727/150‑3 (MQ=255)
                                                                                                                                        gCAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCaagtaagt       >  3:277749/1‑150 (MQ=255)
                                                                                                                                         cAAAGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCAAGTAAGTc      <  1:1085744/150‑1 (MQ=255)
                                                                                                                                          aaaGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGAc                                                                                                         <  2:106881/50‑1 (MQ=255)
                                                                                                                                          aaaGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGAc                                                                                                         >  1:106881/1‑50 (MQ=255)
                                                                                                                                          aaaGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCAAGTAAGTCt     >  2:219307/1‑150 (MQ=255)
                                                                                                                                           aaGCTGTCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCAAGTAAGTCTa    >  3:283420/1‑150 (MQ=255)
                                                                                                                                             gctgtctgATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCAAGTAAGTCTAAg  <  2:1141239/150‑1 (MQ=255)
                                                                                                                                             gctgtcggATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCAAGTAAGTCTAAg  <  2:824634/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TACCGGAGCGGATGGCATAAAGCGATTTAAGTTCATCACCAGCCTTAAACAGCGTCTGGCCTTTCTGAATAGGCTTCTTCCGCTCAATGATATTATCAAGCTGATCAAGCTCATGTTCGTTGAGTGTGAACGGGATGCAAAGCTGGCTGATGCTGCAATCCTGGCAATGGATAGCACAACCGCCAGACTGAATGCGCCGTATAATTCGCTTTTCCGGGATCATAGGTCTGCTCAAGCCGTAATTGATATTTGTCAATTTTAACATCTTTTTAGGGAGCAAGTAAGTCTAAG  >  NC_000913_3_pae_pgi/1399304‑1399594

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